实验二离散时间系统的时域和频域分析(学生).docVIP

实验二 离散时间系统的时域和频域分析 相关MATLAB函数 1. 求解零状态响应：filter()函数 y=filter(b,a,x) 由向量b和a组成的系统对输入x进行滤波 2. 求解单位序列响应：impz()函数 （1）h=impz(b,a) 计算单位序列响应的序列值，取样点个数由MATLAB自动选取 （2）h=impz(b,a,n) 计算指定范围内（0: n-1）的单位序列响应的序列值 （3）impz(b,a) 绘制单位序列响应的时域波形 3. 求解频率响应：freqz()函数 （1）[h,w]=freqz(b,a,n) 可得到n点频率响应，这n个点均匀地分布在上半单位圆（即），并将这n点频率记录在w中，相应的频率响应记录在h中。n最好能取2的幂次方，如果缺省，则n=512。 （2）[h,w]=freqz(b,a,n,whole) 在之间均匀选取n个点计算频率响应。 （3）[h,w]=freqz(b,a,n,Fs) Fs为采样频率（以Hz为单位），在0~Fs/2频率范围内选取n个频率点，计算相应的频率响应。 （4）[h,w]=freqz(b,a,n,whole,Fs) 在0~Fs之间均匀选取n个点计算频率响应。 （5）freqz(b,a) 可以直接得到系统的幅频和相频特性曲线。其中幅频特性以分贝的形式给出，频率特性曲线的横轴采用的是归一化频率，即Fs/2=1。 4. 系统函数的零极点图的绘制：zplane(b,a) 实验内容： 1. 已知某系统的系统函数为，，要求：（1）从理论上求解系统的单位冲激响应和零状态响应，并根据求解结果用MATLAB绘制其时域波形；（2）试分别用MATLAB的impz()函数和filter()函数绘制系统的单位冲激响应和零状态响应。 N=20; n=0:N-1; h=-7/4*(-0.6).^n+11/4*(0.2).^n; yzs=-11/16*(0.2).^n-21/32*(-0.6).^n+75/32; subplot(221) stem(n,h,.);title(理论计算的单位冲激响应); subplot(222) stem(n,yzs,.);title(理论计算的零状态响应); b=[1 2]; a=[1 0.4 -0.12]; x=ones(1,N); yzs1=filter(b,a,x); subplot(224) stem(n,yzs1,.);title(系统的零状态响应); subplot(223) y=impz(b,a,20); stem(n,y,.);title(系统的单位冲击响应); 已知某系统的系统函数为 （1）绘制其零极点图 （2）用freqz()函数绘出该系统的幅频特性曲线和相频特性曲线，并说明该系统的作用。 b=[1]; a=[1 -0.5]; subplot(221) zplane(b,a);title(系统的零极点分布); n=(0:50)*pi/50; [h,w]=freqz(b,a,n); subplot(223) plot(w/pi,abs(h));title(幅频特性曲线); subplot(224) plot(w/pi,angle(h));title(相频特性曲线); 已知某系统的差分方程为，其中，，，设输入信号为，绘出该系统的幅频特性曲线和相频特性曲线，以及系统的输入信号和响应，并说明该系统的性能。 b=[1 -1.1 0.6]; a=[0.6 -1.1 1]; n=(0:50)*pi/50; [h,w]=freqz(b,a,n); subplot(221) plot(w/pi,abs(h));title(幅频特性曲线); subplot(222) plot(w/pi,angle(h));title(相频特性曲线); n=0:20; x=sin(0.1*pi*n); subplot(223) plot(n,x);title(输入信号); y=impz(b,a,21); subplot(224) stem(n,y,.);title(冲击响应);