matlab实现幅频响应

liumilan

若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的区间怎么弄

songzy41

原帖由 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 编辑 ]