facter1.m

有关dsp的,有好多子程序,还是不错的,欢迎使用

function facter1(T)
n=0:50;
x=Xa(444.128,50.*sqrt(2.0).*pi,50.*sqrt(2.0).*pi,T);
subplot(3,1,1);
stem(x);
title('理想采样信号序列');
k=-25:25;
W=(pi./12.5).*k;
X=x*(exp(-j*pi./12.5)).^(n'*k);
magX=abs(X);
subplot(3,1,2);stem(magX);
title('幅频响应');
angX=angle(X);
subplot(3,1,3);stem(angX);
title('相频响应');