u1_4.m

来自「matlab 实现的离散信号分析」· M 代码 · 共 12 行

%卷积计算
n=1:50;
hb=zeros(1,50);
hb(1)=1;hb(2)=2.5;hb(3)=2.5;hb(4)=1;
close all;subplot(3,1,1);stem(hb);title('系统hb[n])');
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);
subplot(3,1,2);stem(x);title('输入信号x[n]');
y=conv(x,hb);
subplot(3,1,3);stem(y);title('输出信号y[n]');