yanzheng.m

dsp学习中matlab试验源代码

%
k=-25:25; 
n=1:50;         
A=444.128;       %设置信号有关的参数 
a=50*sqrt(2.0)*pi;      
T=0.001;        %采样率 
w0=50*sqrt(2.0)*pi;     %ω符号在MatLab中不能输入，用w代替 
x=A*exp(-a*n*T).*sin(w0*n*T);   %pi是MATLAB定义的π，信号乘可采用“.*” 
X=x*(exp(-j*pi/12.5)).^(n'*k); 
magX=abs(X);        %绘制x(n)的幅度谱 
subplot(4,1,1);stem(magX);title('输入信号的幅度谱');

hb=zeros(1,50);       %注意：MATLAB中数组下标从1开始 
hb(1)=1; hb(2)=2.5; hb(3)=2.5; hb(4)=1; 
Hb=hb*(exp(-j*pi/12.5)).^(n'*k);
magHb=abs(Hb);       %绘制hb(n)的幅度谱 
subplot(4,1,2);stem(magHb);title('系统响应的幅度谱');

y=conv(x,hb);
n=1:99; k=1:99; 
Y=y*(exp(-j*pi/12.5)).^(n'*k);
magY=abs(Y);        %绘制y(n)的幅度谱 
subplot(4,1,3);stem(magY);title('输出信号的幅度谱');axis([0,50,0,8000]);

%以下将验证的结果显示 
XHb=X.*Hb; 
Subplot(4,1,4);stem(abs(XHb));title('(R)x(n的幅度谱与 hb(n)幅度谱相乘');
axis([0,50,0,8000]);