使用matlab求离散系统的频率响应分析和零、极点分布

程序及相关说明自提
链接：https://pan.baidu.com/s/1rGWwINx4uOmwS3OuvJor0w?pwd=pdhj
提取码：pdhj

编写代码如下：

% 创建一个传递函数
num = [0.0528 0.797 0.1295 0.1295 0.797 0.0528]; % 分子系数
den = [1 -1.8007 2.4947 -1.8821 0.9537 -0.2336]; % 分母系数

% 将传递函数转换为零极点形式
[z, p, K] = tf2zp(num, den);

% 输出结果
disp('零点:');
disp(z);
disp('极点:');
disp(p);
disp('增益:');
disp(K);

%绘制零极点分布图
figure(1)
zplane(num,den)
% figure(2)
% zplane(z,p)

%计算部分分式展开
[r, p, k] = residuez(num,den);

% 输出结果
disp('残差项系数:');
disp(r);
disp('极点:');
disp(p);
disp('常数项:');
disp(k);

% 将传递函数转换为SOS形式
[z, p, K] = tf2zp(num, den); % 获取零点、极点和增益
sos = zp2sos(z, p, K); % 将传递函数转换为SOS形式

% 输出结果
disp('SOS形式的系数:');
disp(sos);

实验结果：

>> Z_transform_related_exercises
零点:
 -14.9370 + 0.0000i
   0.4546 + 0.8907i
   0.4546 - 0.8907i
  -1.0000 + 0.0000i
  -0.0669 + 0.0000i

极点:
   0.2807 + 0.9097i
   0.2807 - 0.9097i
   0.3743 + 0.6206i
   0.3743 - 0.6206i
   0.4907 + 0.0000i

增益:
    0.0528

残差项系数:
  -0.0011 - 0.5172i
  -0.0011 + 0.5172i
  -1.4372 + 0.4195i
  -1.4372 - 0.4195i
   3.1554 + 0.0000i

极点:
   0.2807 + 0.9097i
   0.2807 - 0.9097i
   0.3743 + 0.6206i
   0.3743 - 0.6206i
   0.4907 + 0.0000i

常数项:
   -0.2260

SOS形式的系数:
    0.0528    0.7887         0    1.0000   -0.4907         0
    1.0000    1.0669    0.0669    1.0000   -0.7486    0.5252
1.0000   -0.9092    1.0000    1.0000   -0.5614    0.9064

零极点图：