r5r10.m

来自「dsp学习中matlab试验源代码」· M 代码 · 共 33 行

%先运行R5.m程序生成一个宽度为5的门函数

%Ha的长度为十
n=0:9; 
x=sign(sign(10-n)+1); 
close all; 
subplot(3,2,1); 
stem(n,x);
title('Ha(t)或Rn(t)信号序列');
k=-25:25; 
X=x*(exp(-j*pi/25)).^(n'*k);
magX=abs(X);                     
subplot(3,2,2);
stem(magX);
title('Ha(t)或Rn(t)信号的幅度谱');axis([0,50,0,20]);



subplot(3,2,3); 
stem(m,y);
title('Xc(t)或Rn(t)信号序列');
subplot(3,2,4);
stem(magY);
title('Xc(t)或Rn(t)信号的幅度谱');axis([0,50,0,20]);

%Xc和Ha进行卷积
z=conv(x,y); 
subplot(3,2,5);stem(z);title('输出信号 y[n]');
p=0:13;
Z=z*(exp(-j*pi/12.5)).^(p'*k);
magZ=abs(Z);       %绘制x(n)的幅度谱 
subplot(3,2,6);stem(magZ);title('输出信号 y[n]的幅度谱');axis([0,51,0,150]);