exprience1_3.m

来自「matlab exprience,M文件,实验代码」· M 代码 · 共 16 行

clear all;
b=[1 2 1];
a=[1 -0.5 0.25];
zplane(b,a);
N=100;
pulse=[1 zeros(1,N-1)];
hn=impz(b,a,N);
figure
stem(hn)
title('pulse response hn')
n=0:1:200;
xn=5+3*cos(0.2*pi*n)+4*sin(0.6*pi*n);
yn=filter(b,a,xn);
figure
stem(n,yn)
title('xn response')