facter3.m

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

function facter3(a,w0,N1)
N=50;
n=1:N;
A=1;
T=1;
xa=A*exp(-a*n*T).*sin(w0*n*T);
subplot(3,3,1);
stem(xa);
title('输入序列');
k=1:N;
w=2*pi/N;
Xa=xa*(exp(-j*w)).^(n'*k);
subplot(3,3,2);
stem(abs(Xa));
title('输入信号序列的幅度谱');
subplot(3,3,3);
stem(angle(Xa));
title ('输入信号序列的相位谱');
ha=R(N1,N);
subplot(3,3,4);
stem(ha);
title('冲激序列');
k=1:N;
n=1:N;
w=2*pi/N;
Ha=ha*(exp(-j*w)).^(n'*k); 
subplot(3,3,5);
stem(abs(Ha));
title('冲激响应序列的幅度谱');
subplot(3,3,6);
stem(angle(Ha)); 
title ('冲激响应序列的相位谱');
Y=Xa.*Ha;
subplot(3,2,5);
stem(abs(Y));
title('输出信号序列的幅度谱');
subplot(3,2,6);
stem(angle(Y)); 
title ('输出信号序列的相位谱');