MATLAB中滤波函数、频率响应函数以及频率响应函数不同表达形式的转换-Toy模板网

这篇具有很好参考价值的文章主要介绍了MATLAB中滤波函数、频率响应函数以及频率响应函数不同表达形式的转换。希望对大家有所帮助。如果存在错误或未考虑完全的地方，请大家不吝赐教，您也可以点击"举报违法"按钮提交疑问。

一、理论公式

频率响应函数的表达式： $MATLAB中滤波函数、频率响应函数以及频率响应函数不同表达形式的转换$

对应的z变换的多项表达式： $MATLAB中滤波函数、频率响应函数以及频率响应函数不同表达形式的转换$

Z变换的零极点表达形式： $MATLAB中滤波函数、频率响应函数以及频率响应函数不同表达形式的转换$

Z变换的二阶因子级联形式： $MATLAB中滤波函数、频率响应函数以及频率响应函数不同表达形式的转换$

二、滤波函数filter

filter函数，仅可以用于零状态响应系统。

y=filter(b,a,x) ；

%b为z变换多项表达式公式中[b0,b1...bM]的矩阵

%a为z变换多项表达式公式中[a0,a1...aN]的矩阵

%x为输入的原始信号

移动平均滤波的matlab程序示例：

%求解零状态差分方程函数--LTI系统
% y=filter(b,a,x);


N=201;
n=rand(1,N)-0.5;%噪声信号
k=0:N-1;%
x=2*k.*(0.9.^k)+2.0*cos(0.02*pi*k)+5.0;%输入信号
y=x+n;%包含噪声的信号
plot(k,n,'g--',k,x,'b--',k,y,'r-');%显示三个函数图像，
xlabel('Time index k');
legend('n[k]','x[k]','y[k]');%添加图例


M=10;
b=ones(M,1)/M;%移动平均滤波，10个点相加，再除10，算出平均值
a=[1];
filter_y=filter(b,a,y);
plot(k,x,'b-',k,y,'g--',k,filter_y,'r-');%显示三个函数的波形
xlabel('Time index k');
legend('x[k]','y[k]','filter_y[k]');%添加图例

得出的波形如下图：

MATLAB中滤波函数、频率响应函数以及频率响应函数不同表达形式的转换

三、频率响应函数

H=freqz（b,a,w）；

%b为z变换多项表达式公式中[b0,b1...bM]的矩阵

%a为z变换多项表达式公式中[a0,a1...aN]的矩阵

%w为输入的角频率

matlab程序示例

%freqz用以分析离散系统的频率响应

b=[1]; %分子矩阵
a1=[1,-0.9]; a2=[1,0.9];%分母矩阵
w=linspace(0,pi,512);%在0-π范围内等间隔分512份
h1=freqz(b,a1,w);%计算频率响应
h2=freqz(b,a2,w);%计算频率响应
% plot(w,abs(h1),w,abs(h2),':');
plot(w/pi,abs(h1),w/pi,abs(h2),':');
legend('a=0.9','a=-0.9');

MATLAB中滤波函数、频率响应函数以及频率响应函数不同表达形式的转换

可以看出在分母矩阵为[1,-0.9]时，系统为低通滤波器；分母矩阵为[1,0.9]时，系统为高通滤波器。

四、频率响应函数不同形式的转换

b=[1,4];%Z变换的分子矩阵
a=[1,0.1,-0.2];%Z变换的分母矩阵
z=zeros(1,2);%2*1的矩阵，值为0
% [z,p,K]=tf2zp(b,a);%多项式转零极点表达式
[z,p,K]=tf2zpk(b,a);%多项式转零极点表达式
sos=zp2sos(z,p,K);%零极点转二阶因子级联形式

多项式表达式为 $MATLAB中滤波函数、频率响应函数以及频率响应函数不同表达形式的转换$

经过matlab计算，得到其零极点表达式为 $MATLAB中滤波函数、频率响应函数以及频率响应函数不同表达形式的转换$

二阶因子级联形式： $MATLAB中滤波函数、频率响应函数以及频率响应函数不同表达形式的转换$ （只有1级）

需要注意的是：tf2zp计算出来的结果中只有一个零点，而 tf2zpk计算出来的结果中有两个零点（多出来一个位置在0处的零点）

tf2zpk是用于离散系统传递函数，基于z变换；而tf2zp是用于连续系统传递函数，基于拉普拉斯变换。依据二者的函数公式，计算该系统函数的零极点出来，确实会存在零点的细微差别。此处不做讨论，注意即可。文章来源地址https://www.toymoban.com/news/detail-468158.html

到了这里，关于MATLAB中滤波函数、频率响应函数以及频率响应函数不同表达形式的转换的文章就介绍完了。如果您还想了解更多内容，请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章，希望大家以后多多支持TOY模板网！