信号与系统仿真实验——实验四 离散系统频率响应特性分析及系统函数极零点分布特性

【 实验目的】

1.系统函数极、零点分布特性分析
2.离散系统频率响应特性分析及其特性曲线的绘制

【 实验内容】

【 实现方法】

【 实验报告要求】

记录实验波形并分析波形特征，根据实验要求得出相应的结论；对系统函数极、零点图的应用做总结。

Matlab程序一：

方程式1和2
b=[1,0];  %1式                 
a=[1,-1];
k=0:10;
subplot(2,2,1);
zplane(b,a);%用zplane()函数绘制系统零极点图
title('零极点图Make by 磊');
grid on;
h=impz(b,a,k);%用impz()函数来绘制单位函数响应。
subplot(2,2,2);
stem(h);%离散序列绘图用stem()函数
title('单位函数响应Make by 磊');
grid on;
b=[1,0];  %2式                 
a=[1,-0.8];
k=0:10;
subplot(2,2,3);
zplane(b,a);%用zplane()函数绘制系统零极点图
title('零极点图Make by 磊');
grid on;
h=impz(b,a,k);%用impz()函数来绘制单位函数响应。
subplot(2,2,4);
stem(h);%离散序列绘图用stem()函数
title('单位函数响应Make by 磊');
grid on;

运行结果一：

Matlab程序二：

方程式3和4
b=[1,0];  %3式                 
a=[1,-1.2];
k=0:10;
subplot(2,2,1);
zplane(b,a);%用zplane()函数绘制系统零极点图
title('零极点图Make by 磊');
grid on;
h=impz(b,a,k);%用impz()函数来绘制单位函数响应。
subplot(2,2,2);
stem(h);%离散序列绘图用stem()函数
title('单位函数响应Make by 磊');
grid on;
b=[1,0];  %4式                 
a=[1,-4*sqrt(2)/5,0.64];
k=0:10;
subplot(2,2,3);
zplane(b,a);%用zplane()函数绘制系统零极点图
title('零极点图Make by 磊');
grid on;
h=impz(b,a,k);%用impz()函数来绘制单位函数响应。
subplot(2,2,4);
stem(h);%离散序列绘图用stem()函数
title('单位函数响应Make by 磊');
grid on;

运行结果二：

Matlab程序三：

方程式5和6
b=[1,0];  %5式                 
a=[1,-2*cos(pi/8),1];
k=0:10;
subplot(2,2,1);
zplane(b,a);%用zplane()函数绘制系统零极点图
title('零极点图Make by 磊');
grid on;
h=impz(b,a,k);%用impz()函数来绘制单位函数响应。
subplot(2,2,2);
stem(h);%离散序列绘图用stem()函数
title('单位函数响应Make by 磊');
grid on;
b=[1,0];  %6式                 
a=[1,-6*sqrt(2)/5,1.44];
k=0:10;
subplot(2,2,3);
zplane(b,a);%用zplane()函数绘制系统零极点图
title('零极点图Make by 磊');
grid on;
h=impz(b,a,k);%用impz()函数来绘制单位函数响应。
subplot(2,2,4);
stem(h);%离散序列绘图用stem()函数
title('单位函数响应Make by 磊');
grid on;

运行结果三：

Matlab程序四：

方程式7
b=[1,1];  %7式                 
a=[3,-1,0,0,0,1];
k=0:10;
subplot(1,2,1);
zplane(b,a);%用zplane()函数绘制系统零极点图
title('零极点图Make by 磊');
grid on;
h=impz(b,a,k);%用impz()函数来绘制单位函数响应。
subplot(1,2,2);
stem(h);%离散序列绘图用stem()函数
title('单位函数响应Make by 磊');
grid on;

运行结果四：

Matlab程序五：

方程1
b=[1];                                            
a=[1,-0.5];
subplot(2,2,1);
zplane(b,a);
title('零极点图Make by 磊');
h=freqz(b,a);%用freqz()函数来分析系统的频率响应特性
subplot(2,2,2);
plot(h);
title('频率响应特性曲线Make by 磊');
subplot(2,2,3);
plot(abs(h));
title('幅频特性曲线Make by 磊');
subplot(2,2,4);
plot(angle(h));
title('相频特性曲线Make by 磊');

运行结果五：

Matlab程序六：

方程2
b=[1,0];                                            
a=[1,-0.5];
subplot(2,2,1);
zplane(b,a);
title('零极点图Make by 磊');
h=freqz(b,a);%用freqz()函数来分析系统的频率响应特性
subplot(2,2,2);
plot(h);
title('频率响应特性曲线Make by 磊');
subplot(2,2,3);
plot(abs(h));
title('幅频特性曲线Make by 磊');
subplot(2,2,4);
plot(angle(h));
title('相频特性曲线Make by 磊');

运行结果六：

Matlab程序七：

方程3
b=[1,0];                                            
a=[1,0.5];
subplot(2,2,1);
zplane(b,a);
title('零极点图Make by 磊');
h=freqz(b,a);%用freqz()函数来分析系统的频率响应特性
subplot(2,2,2);
plot(h);
title('频率响应特性曲线Make by 磊');
subplot(2,2,3);
plot(abs(h));
title('幅频特性曲线Make by 磊');
subplot(2,2,4);
plot(angle(h));
title('相频特性曲线Make by 磊');

运行结果七：

Matlab程序八：

方程4
b=[1,0.5];                                            
a=[1,0];
subplot(2,2,1);
zplane(b,a);
title('零极点图Make by 磊');
h=freqz(b,a);%用freqz()函数来分析系统的频率响应特性
subplot(2,2,2);
plot(h);
title('频率响应特性曲线Make by 磊');
subplot(2,2,3);
plot(abs(h));
title('幅频特性曲线Make by 磊');
subplot(2,2,4);
plot(angle(h));
title('相频特性曲线Make by 磊');

运行结果八：

本人能力有限，解释尚不清楚明了，如遇任何问题，大家可留言或私信。后续将程序文件打包上传，供大家学习使用。

本文希望对大家有帮助，当然上文若有不妥之处，欢迎指正。

分享决定高度，学习拉开差距