MATLAB频域分析
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
求该模拟滤波器的频率响应。 MATLAB 源程 序如下。 B=1;A=[1 2.6131 3.4142 2.6131 1]; W=0:0.1:pi; freqs(B,A,W) 注：如果单一用freqs则系统自动画图， 可能将没有价值的信息省略。
若要获得幅频特性和相频特性，则分别取其模(abs)和角度 （angle） clear 1.5 clc 1 B=1 A=[1 2.6131 3.4142 2.6131 1] 0.5 W=0:0.1:pi 0 0 h=freqs(B,A,W) 4 H=abs(h) 2 s=angle(h) 0 subplot(2,1,1) -2 plot(W,H) -4 0 subplot(2,1,2) plot(W,s)
4.6 线性时不变系统的频率响应
1．求模拟滤波器Ha(s)的频率响应函数freqs( ) 格式：H＝freqs(B,A,W) 功能：计算由向量W(rad/s)指定的频率点上模拟滤器系统函数Ha(s) 的频率响应Ha(jΩ )，结果存于H向量中。
[例4-31] 已知某模拟滤波器可表示为：
y(4) (t ) 2.6131y(3) (t ) 3.4141y(2) (t ) 2.6131y(1) (t ) y(t ) f (t )
0.5
1
1.5百度文库
2
2.5
3
3.5
0.5
1
1.5
2
2.5
3
3.5
2．求数字滤波器H(z)的频率响应函数freqz( ) 格式：H=freqz(B,A,W) 功能：计算由向量W(rad)指定的数字频率点上(通常指在H(z)的频率 响应H(ejw )。 [例4-32] 已知某滤波器可表示为
y(n 8) y(n) x(n)
求该滤波器的频率响应。MATLAB源程序为： B=[1 0 0 0 0 0 0 0 –1]; A=1; freqz(B,A) 该程序运行所绘出的幅频与 相频性曲线如图4.31所示。
图4.31滤波器幅度和相位曲线
若要获得幅频特性和相频特性，则分别取其模(abs)和角度 （angle） clear 1.5 clc 1 B=1 A=[1 2.6131 3.4142 2.6131 1] 0.5 W=0:0.1:pi 0 0 h=freqs(B,A,W) 4 H=abs(h) 2 s=angle(h) 0 subplot(2,1,1) -2 plot(W,H) -4 0 subplot(2,1,2) plot(W,s)
4.6 线性时不变系统的频率响应
1．求模拟滤波器Ha(s)的频率响应函数freqs( ) 格式：H＝freqs(B,A,W) 功能：计算由向量W(rad/s)指定的频率点上模拟滤器系统函数Ha(s) 的频率响应Ha(jΩ )，结果存于H向量中。
[例4-31] 已知某模拟滤波器可表示为：
y(4) (t ) 2.6131y(3) (t ) 3.4141y(2) (t ) 2.6131y(1) (t ) y(t ) f (t )
0.5
1
1.5百度文库
2
2.5
3
3.5
0.5
1
1.5
2
2.5
3
3.5
2．求数字滤波器H(z)的频率响应函数freqz( ) 格式：H=freqz(B,A,W) 功能：计算由向量W(rad)指定的数字频率点上(通常指在H(z)的频率 响应H(ejw )。 [例4-32] 已知某滤波器可表示为
y(n 8) y(n) x(n)
求该滤波器的频率响应。MATLAB源程序为： B=[1 0 0 0 0 0 0 0 –1]; A=1; freqz(B,A) 该程序运行所绘出的幅频与 相频性曲线如图4.31所示。
图4.31滤波器幅度和相位曲线