xahb.m

dsp学习中matlab试验源代码

n=0:50;              %定义序列的长度是50 
hb=zeros(1,51);      %注意：MATLAB中数组下标从1开始 
hb(1)=1; hb(2)=2.5; hb(3)=2.5; hb(4)=1; 
close all; subplot(3,1,1);stem(n,hb);title('系统 hb[n]');axis([0,50,-4,4]);

m=1:50; T=0.001;      %定义序列的长度是和采样率 
A=444.128; a=50*sqrt(2.0)*pi;   %设置信号有关的参数 
w0=50*sqrt(2.0)*pi; 
x=A*exp(-a*m*T).*sin(w0*m*T);   %pi是MATLAB定义的π，信号乘可采用“.*” 
subplot(3,1,2);stem(x);title('输入信号 xa[n]');axis([0,50,-200,200]);
y=conv(x,hb); 
subplot(3,1,3);stem(y);title('输出信号 y[n]');axis([0,50,-1000,1200]);