xbhb.m

dsp学习中matlab试验源代码

n=0:50;              %定义序列的长度是50 
hb=zeros(1,51);      
hb(1)=1; hb(2)=2.5; hb(3)=2.5; hb(4)=1; 
close all; subplot(3,2,1);stem(n,hb);title('系统 hb[n]');axis([0,10,0,4]);
k=-25:25; 
HB=hb*(exp(-j*pi/12.5)).^(n'*k);
magX=abs(HB);       %绘制x(n)的幅度谱 
subplot(3,2,2);stem(magX);title('系统 hb[n]的幅度谱');

m=0:50;              %定义序列的长度是50 
x=zeros(1,51);       %注意：MATLAB中数组下标从1开始 
x(1)=1; 
subplot(3,2,3);stem(m,x);title('输入信号 xb[n]');axis([0,10,0,1]);
k=-25:25; 
X=x*(exp(-j*pi/12.5)).^(n'*k);
magX=abs(X);       %绘制x(n)的幅度谱 
subplot(3,2,4);stem(magX);title('输入信号 xb[n]的幅度谱');




y=conv(x,hb); 
subplot(3,2,5);stem(y);title('输出信号 y[n]');axis([0,10,0,4]);
p=0:100;
Z=y*(exp(-j*pi/12.5)).^(p'*k);
magZ=abs(Z);       %绘制x(n)的幅度谱 
subplot(3,2,6);stem(magZ);title('输出信号 y[n]的幅度谱');