求助:如何用Matlab模拟扫频激励信号下的响应信号?

chris-chen

各位高手们,小弟最近编程遇到点问题:
我看到一篇老外的文献中有模拟扫频信号和响应信号的,响应信号是通过扫频信号和脉冲响应函数h(t)卷积得到,可是用matlab怎么也弄不出来响应信号,还请各位帮帮忙,真诚感谢!
附文献主要公式:





结果如下:

chris-chen

附程序:请各位看看

1. clear all;
   
2. close all;
   
3. clc;
   
4. 
   
5. 
   
6. % -------------------------------------------------------------------------
   
7. %       设置参数
   
8. % -------------------------------------------------------------------------
   
9. fs=256;
   
10. f0=1;
    
11. fm=17;
    
12. T=64;                               %时间长度
    
13. C=1/T *log(fm/f0);
    
14. t=0:1.0/fs:T-1.0/fs;
    
15. x=sin(((2*pi*f0)/C).*(exp(C*t)-1)); %扫频信号
    
16. 
    
17. figure(1);
    
18. plot(t,x);         
    
19. xlabel('Time[sec]');
    
20. ylabel('Amplitude');
    
21. 
    
22. % -------------------------------------------------------------------------
    
23. %脉冲响应函数
    
24. gamma=[0.005,0.05,0.005,0.05,0.005,0.05] %阻尼比
    
25. omega=2*pi.*[4,6,8,10,12,14]             %频率
    
26. h=0;
    
27. for k=1:6
    
28.   h=h+exp((-gamma(k)+i*sqrt(1-gamma(k)^2))*omega(k)*t);  
    
29.    end

复制代码

chris-chen

这儿没有高手帮忙解决吗?

meiyongyuandeze

运行了下没有报错信息啊，不知道什么问题啊、！

chris-chen

程序还没有算完,还有一点,最终得不到响应信号的波形,不知道怎么回事

1. H=fft(h);
   
2. X=fft(x);
   
3. Y=X.*H;
   
4. y=abs(ifft(Y));
   
5. figure;
   
6. plot(t,y);

复制代码

chris-chen

这儿有高手能把这个响应信号模拟出来吗?

cff190

请问这篇文章的题目是什么？

ChaChing

没空细看! 请版主先移信号版块!?

xiazunping

加一句

1. Y=[Y conj(Y(end-1:-1:2))];

复制代码

就可以了，不要说不知道为什么！

白起

期待能够有所帮助