chris-chen 发表于 2011-4-26 15:08

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

本帖最后由 chris-chen 于 2011-4-26 15:20 编辑

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





结果如下:


chris-chen 发表于 2011-4-26 15:34

本帖最后由 牛小贱 于 2015-3-31 09:21 编辑

附程序:请各位看看
clear all;
close all;
clc;


% -------------------------------------------------------------------------
%       设置参数
% -------------------------------------------------------------------------
fs=256;
f0=1;
fm=17;
T=64;                               %时间长度
C=1/T *log(fm/f0);
t=0:1.0/fs:T-1.0/fs;
x=sin(((2*pi*f0)/C).*(exp(C*t)-1)); %扫频信号

figure(1);
plot(t,x);         
xlabel('Time');
ylabel('Amplitude');

% -------------------------------------------------------------------------
%脉冲响应函数
gamma= %阻尼比
omega=2*pi.*             %频率
h=0;
for k=1:6
h=h+exp((-gamma(k)+i*sqrt(1-gamma(k)^2))*omega(k)*t);
   end

chris-chen 发表于 2011-4-26 16:57

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

meiyongyuandeze 发表于 2011-4-26 17:04

本帖最后由 meiyongyuandeze 于 2011-4-26 17:06 编辑

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

chris-chen 发表于 2011-4-26 17:10

本帖最后由 牛小贱 于 2015-3-31 09:21 编辑

程序还没有算完,还有一点,最终得不到响应信号的波形,不知道怎么回事
H=fft(h);
X=fft(x);
Y=X.*H;
y=abs(ifft(Y));
figure;
plot(t,y);

chris-chen 发表于 2011-4-26 20:05

本帖最后由 chris-chen 于 2011-4-26 22:47 编辑

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

cff190 发表于 2011-5-17 16:49

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

ChaChing 发表于 2011-5-18 00:40

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

xiazunping 发表于 2011-7-22 12:59

本帖最后由 牛小贱 于 2015-3-31 09:22 编辑

加一句Y=;就可以了，不要说不知道为什么！

白起 发表于 2015-3-31 07:32

期待能够有所帮助
页: [1]
查看完整版本: 求助:如何用Matlab模拟扫频激励信号下的响应信号?