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

迷南。 2023-10-05 19:25 11阅读 0赞

## 【 实验目的】 ##

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

## 【 实验内容】 ##

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0_size_16_color_FFFFFF_t_70_pic_center]

## 【 实现方法】 ##

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0_size_16_color_FFFFFF_t_70_pic_center 1]

## 【 实验报告要求】 ##

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

## 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;

## 运行结果一： ##

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0_size_16_color_FFFFFF_t_70_pic_center 2]

## 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;

## 运行结果二： ##

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0_size_16_color_FFFFFF_t_70_pic_center 3]

## 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;

## 运行结果三： ##

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0_size_16_color_FFFFFF_t_70_pic_center 4]

## 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;

## 运行结果四： ##

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0_size_16_color_FFFFFF_t_70_pic_center 5]

## 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 磊');

## 运行结果五： ##

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0_size_16_color_FFFFFF_t_70_pic_center 6]

## 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 磊');

## 运行结果六： ##

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0_size_16_color_FFFFFF_t_70_pic_center 7]

## 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 磊');

## 运行结果七： ##

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0_size_16_color_FFFFFF_t_70_pic_center 8]

## 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 磊');

## 运行结果八： ##

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0_size_16_color_FFFFFF_t_70_pic_center 9]

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

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

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

[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0_size_16_color_FFFFFF_t_70_pic_center]: https://img-blog.csdnimg.cn/20210207212728981.bmp?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0,size_16,color_FFFFFF,t_70#pic_center
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0_size_16_color_FFFFFF_t_70_pic_center 1]: https://img-blog.csdnimg.cn/20210207212817425.bmp?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0,size_16,color_FFFFFF,t_70#pic_center
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0_size_16_color_FFFFFF_t_70_pic_center 2]: https://img-blog.csdnimg.cn/20210207213114952.jpg?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0,size_16,color_FFFFFF,t_70#pic_center
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0_size_16_color_FFFFFF_t_70_pic_center 3]: https://img-blog.csdnimg.cn/20210207213147479.jpg?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0,size_16,color_FFFFFF,t_70#pic_center
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0_size_16_color_FFFFFF_t_70_pic_center 4]: https://img-blog.csdnimg.cn/20210207213216594.jpg?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0,size_16,color_FFFFFF,t_70#pic_center
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0_size_16_color_FFFFFF_t_70_pic_center 5]: https://img-blog.csdnimg.cn/20210207213248599.jpg?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0,size_16,color_FFFFFF,t_70#pic_center
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0_size_16_color_FFFFFF_t_70_pic_center 6]: https://img-blog.csdnimg.cn/20210207213319404.jpg?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0,size_16,color_FFFFFF,t_70#pic_center
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0_size_16_color_FFFFFF_t_70_pic_center 7]: https://img-blog.csdnimg.cn/20210207213350445.jpg?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0,size_16,color_FFFFFF,t_70#pic_center
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0_size_16_color_FFFFFF_t_70_pic_center 8]: https://img-blog.csdnimg.cn/20210207213420307.jpg?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0,size_16,color_FFFFFF,t_70#pic_center
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0_size_16_color_FFFFFF_t_70_pic_center 9]: https://img-blog.csdnimg.cn/2021020721344947.jpg?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQyMDc4OTM0,size_16,color_FFFFFF,t_70#pic_center