matlab实现幅频响应
若h(t)={1-t/Ts 0<t<Ts;1+t/Ts -Ts<t<0;
0 其他
怎么求出它的频域响应啊
用matlab写出了h(t)的程序
dt=0.001;
t=-2:dt:2;
Ts=1;
ht=zeros(size(t));
for i=1:length(t)
if (t(i)>0)&(t(i)<Ts)
ht(i)=1-t(i)/Ts;
elseif (t(i)>-Ts)&(t(i)<=0)
ht(i)=1+t(i)/Ts;
else ht(i)=0;
end
end
那怎么画出它的幅频特性图啊？
用abs(fft(ht))吗？
画出来的图形怎么跟书本不一样？
f的区间怎么弄 本帖最后由 wdhd 于 2016-9-12 13:08 编辑
原帖由 liumilan 于 2008-6-14 20:31 发表
那怎么画出它的幅频特性图啊？
用abs(fft(ht))吗？
画出来的图形怎么跟书本不一样？
f的区间怎么弄
把程序稍修改，看是否满足楼主要求：
dt=0.001; fs=1/dt;
t=-2:dt:2;
Ts=1;
ht=zeros(size(t));
for i=1:length(t)
if (t(i)>0)&(t(i)<Ts)
ht(i)=1-t(i)/Ts;
elseif (t(i)>-Ts)&(t(i)<=0)
ht(i)=1+t(i)/Ts;
else ht(i)=0;
end
end
subplot 211; plot(t,ht); grid;
H=fft(ht);
N=length(ht);
HH=H/max(abs(H));
N2=fix(N/2);
ff=(-N2:N2)*fs/N;
n2=1:N/2+1;
subplot 212; plot(ff,20*log10(abs(fftshift(HH)))); grid;
axis([-10 10 -100 0]);
[ 本帖最后由 songzy41 于 2008-6-15 16:27 编辑 ]
页:
[1]