标签：subplot 频率响应 title jw 对系统 xlabel ylabel 极点
系列文章目录
数字信号处理（DSP：Digital Signal Process）是电子通信领域非常重要的研究方向，博主汇总了数字信号处理（DSP）中常用的经典案例分析，主要基于算法分析、MATLAB程序实现、信号图像显示，对数字信号处理的实际应用进行详细分析。第一篇 常见离散信号的MATLAB产生和图形显示
常见离散信号的MATLAB产生和图形显示
第二篇 零极点分布对系统频率响应的影响
零极点分布对系统频率响应的影响
文章目录
一、研究目的
1.加深对离散系统的频率响应分析和零极点分布概念的理解。
2.学习用零极点分布的几何方法分析研究系统的频率响应。
二、原理分析
如果知道系统的系统函数 ，可以得到它的零极点分布，由零极点分布可以方便地对系统的频率响应进行定性分析。
按照教材（2.6.8）式和（2.6.9）式，系统的幅度特性由零点矢量长度之积除以极点矢量长度之积，当频率 从 变化到 时，观察零点矢量长度和极点矢量长度的变化，重点观察那些矢量长度较短的情况。另外，由分析知道，极点主要影响频率响应的峰值，极点愈靠近单位圆，峰值愈尖锐；零点主要影响频率响应的谷值，零点愈靠近单位圆，谷值愈深，如果零点在单位圆上，那么该频率点的频率特性为零。根据这些规律可以定性画出频率响应的幅度特性。
峰值频率和谷值频率可以近似用响应的极点和零点的相角表示，例如极点 ，峰值频率近似为 ，极点愈靠近单位圆，估计法结果愈准确。
本实验借助计算机分析系统的频率响应，目的是掌握用零极点分布的几何方法分析研究系统的频率响应，实验时需要将 代入系统函数 中，再在 之间，等间隔选择若干点，并计算它的频率响应。
三、实验内容
1.情况一
假设系统用下面差分方程描述：
假设a = 0.7 ，a = 0.8 ，a = 0.9 ，分别在三种情况下分析系统的频率特性，并打印幅度特性曲线。
2.情况二
假设系统用下面差分方程描述：
假设a = 0.7 ，a = 0.8 ，a = 0.9 ，分别在三种情况下分析系统的频率特性，并打印幅度特性曲线。
3.情况三
假设系统用下面差分方程描述：
试分析它的频率特性，打印其幅度特性曲线，并求出峰值频率和谷值频率，并与零极点的相角比较。
四、实验代码（matlab）
1.情况一
%%3（1）
B=[1];
a=0.7; A=[1,-a];
[H,w]=freqz(B,A,512,'whole');
figure; subplot(1,3,1); zplane(B,A);
xlabel ('Re(z)'); ylabel('Im(z)');title('零极点图');
subplot(1,3,2); plot(w,abs(H));
xlabel ('w'); ylabel('|H(jw)|'); title('幅频特性曲线');
subplot(1,3,3); plot(w,angle(H));
xlabel ('w'); ylabel ('φ(jw)'); title('相频特性曲线');
B=[1];
a=0.8; A=[1,-a];
[H,w]=freqz(B,A,512,'whole');
figure(2); subplot(1,3,1); zplane(B,A);
xlabel ('Re(z)'); ylabel('Im(z)');title('零极点图');
subplot(1,3,2); plot(w,abs(H));
xlabel ('w'); ylabel('|H(jw)|'); title('幅频特性曲线');
subplot(1,3,3); plot(w,angle(H));
xlabel ('w'); ylabel ('φ(jw)'); title('相频特性曲线');
B=[1];
a=0.9; A=[1,-a];
[H,w]=freqz(B,A,512,'whole');
figure(3); subplot(1,3,1); zplane(B,A);
xlabel ('Re(z)'); ylabel('Im(z)');title('零极点图');
subplot(1,3,2); plot(w,abs(H));
xlabel ('w'); ylabel('|H(jw)|'); title('幅频特性曲线');
subplot(1,3,3); plot(w,angle(H));
xlabel ('w'); ylabel ('φ(jw)'); title('相频特性曲线');
2.情况二
%%3（2）
a=0.7;
B=[1,a];
A=[1];
[H,w]=freqz(B,A,512,'whole');
figure; subplot(1,3,1); zplane(B,A);
xlabel ('Re(z)'); ylabel('Im(z)');title('零极点图');
subplot(1,3,2); plot(w,abs(H));
xlabel ('w'); ylabel('|H(jw)|'); title('幅频特性曲线');
subplot(1,3,3); plot(w,angle(H));
xlabel ('w'); ylabel ('φ(jw)'); title('相频特性曲线');
a=0.8;
B=[1,a];
A=[1];
[H,w]=freqz(B,A,512,'whole');
figure(2); subplot(1,3,1); zplane(B,A);
xlabel ('Re(z)'); ylabel('Im(z)');title('零极点图');
subplot(1,3,2); plot(w,abs(H));
xlabel ('w'); ylabel('|H(jw)|'); title('幅频特性曲线');
subplot(1,3,3); plot(w,angle(H));
xlabel ('w'); ylabel ('φ(jw)'); title('相频特性曲线');
a=0.9;
B=[1,a];
A=[1];
[H,w]=freqz(B,A,512,'whole');
figure(3); subplot(1,3,1); zplane(B,A);
xlabel ('Re(z)'); ylabel('Im(z)');title('零极点图');
subplot(1,3,2); plot(w,abs(H));
xlabel ('w'); ylabel('|H(jw)|'); title('幅频特性曲线');
subplot(1,3,3); plot(w,angle(H));
xlabel ('w'); ylabel ('φ(jw)'); title('相频特性曲线');
3.情况三
%%3（3）
B=[1,1];
A=[1,-1.273,0.81];
[H,w]=freqz(B,A,512,'whole');
figure; subplot(1,3,1); zplane(B,A);
xlabel ('Re(z)'); ylabel('Im(z)');title('零极点图');
subplot(1,3,2);
y=abs(H);
plot(w,y);
xlabel ('w'); ylabel('|H(jw)|'); title('幅频特性曲线');
subplot(1,3,3); plot(w,angle(H));
xlabel ('w'); ylabel ('φ(jw)'); title('相频特性曲线');
f=find(y==max(y));
w(f)%得到幅度最大时所对应的频率点（即峰值频率）
f=find(y==min(y));
w(f)%得到幅度最小时所对应的频率点（即谷值频率）
五、结果图像
1.情况一
a=0.7
a=0.8
a=0.9
2.情况二
a=0.7
a=0.8
a=0.9
3.情况三
标签：subplot,频率响应,title,jw,对系统,xlabel,ylabel,极点 来源： https://blog.csdn.net/weixin_43048780/article/details/121622235
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。