matlab 离散时间信号,离散时间信号和离散时间系统Matlab实现B.ppt

本章相关MATLAB函数 要求a(1)＝1，若不为1，则程序会自动将其归一化为1。 1、filter.m filter可用来求一个离散系统的输出。调用格式： y=filter(b,a,x); x(t)=sin(2*pi*200*t)+ sin(2*pi*10*t)T=1/1000;n=0:199;x(n)=sin(2*pi*200*n*T)+ sin(2*pi*10*n*T); n = 0:199; %取200个点 T=1/1000; %采样频率1KHz x = sin(2*pi*200*n*T)+ sin(2*pi*10*n*T); b=[1,1]; a=[1,-0.9]; y=filter(b,a,x); subplot(2,1,1); stem(n, x); grid on; title(' x'); subplot(2,1,2); stem(n, y); grid on; title(' y'); 2、impz.m impz可用来求一个离散系统的h(n)。调用格式： h=impz(b,a,N); [h,t]=impz(b,a,N); 注：其中，N是h(n)所需的长度。 n = 0:199; %x(n)取200个点 T=1/1000; %采样频率1KHz x = sin(2*pi*200*n*T)+ sin(2*pi*10*n*T); b=[1,1]; a=[1,-0.9]; h= impz(b,a,200); ny=0:398; y=conv(x,h); subplot(2,1,1); stem(n, x); grid on; title(' x'); subplot(2,1,2); stem(ny, y); grid on; title(' y'); axis([0 200 -20 20]); stem(n, h) 3、freqz.m freqz可用来求一个离散系统的频率响应。调用格式： freqz(b,a,N,’whole’,Fs); [H,w] = freqz(b,a,N,’whole’,Fs); 1) N是频率轴分点数，建议N为2的整次幂。2) w返回频率轴坐标向量供绘图用。3) Fs是采样频率，若Fs=1，频率轴给出归一化频率。4) whole指定计算的频率范围从0～Fs，缺省时从 0～Fs/2。 b=[1,1]; a=[1,-0.9]; [H,w]=freqz(b,a,512,'whole',1000); subplot(2,1,1); plot(w,abs(H)) ; grid on; subplot(2,1,2); plot(w,angle(H)); grid on; 4、zplane.m zplane可用来显示离散系统的零极点图。调用格式： zplane(b,a); b=[1,1]; a=[1,-0.9]; zplane(b,a);