用MATLAB画零极点图
电源 零极点 matlab
电源零极点 matlab介绍在电力系统中，电源是一个重要的组成部分，它为电力设备提供所需的电能。
电源的性能直接影响到电力系统的稳定性和可靠性。
为了更好地了解电源的特性和行为，我们可以使用Matlab来进行分析和模拟。
什么是零极点在电源的分析中，零极点是非常重要的概念。
零点是电源传递函数中使得输出为零的输入值，而极点则是使得传递函数无穷大的输入值。
了解电源的零极点分布可以帮助我们更好地理解电源的频率响应和稳定性。
使用Matlab分析电源的零极点Matlab是一个功能强大的数学计算软件，它提供了丰富的工具和函数来进行电源的分析。
以下是使用Matlab分析电源零极点的一般步骤：1. 定义电源传递函数首先，我们需要定义电源的传递函数。
传递函数是描述输入和输出之间关系的数学表达式。
在Matlab中，可以使用tf函数来定义传递函数。
例如，如果我们有一个传递函数为s/(s+1)的电源，可以使用以下代码进行定义：num = [1];den = [1, 1];sys = tf(num, den);2. 绘制零极点图接下来，我们可以使用zpk函数来获取电源的零极点信息，并使用zplane函数来绘制零极点图。
零极点图可以帮助我们直观地了解电源的频率响应和稳定性。
以下是绘制零极点图的代码示例：[z, p, k] = zpkdata(sys);zplane(z, p);3. 分析零极点的位置和数量通过观察零极点图，我们可以分析电源的零极点的位置和数量。
零极点的位置可以告诉我们电源的频率响应特性，而零极点的数量则与电源的稳定性相关。
通常情况下，零极点越远离原点，电源的频率响应越宽，但稳定性可能会受到影响。
4. 评估电源的稳定性除了通过零极点图来评估电源的稳定性外，我们还可以使用Matlab提供的稳定性分析工具来进行评估。
Matlab提供了一些函数，如isstable和margin，可以帮助我们快速评估电源的稳定性。
以下是使用isstable函数来评估电源稳定性的代码示例：isStable = isstable(sys);5. 优化电源设计根据对电源零极点的分析和评估结果，我们可以进一步优化电源的设计。
自动控制原理MATLAB实验报告
实验一 典型环节的MATLAB 仿真一、实验目的1．熟悉MATLAB 桌面和命令窗口，初步了解SIMULINK 功能模块的使用方法。
2．通过观察典型环节在单位阶跃信号作用下的动态特性，加深对各典型环节响应曲线的理解。
3．定性了解各参数变化对典型环节动态特性的影响。
二、实验原理1．比例环节的传递函数为 K R K R R R Z ZsG 200,1002)(211212==-=-=-=其对应的模拟电路及SIMULINK 图形如图1所示。
2．惯性环节的传递函数为uf C K R K R s C R R R Z Z s G 1,200,10012.021)(121121212===+-=+-=-=其对应的模拟电路及SIMULINK 图形如图2所示。
图1 比例环节的模拟电路及SIMULINK 图形3．积分环节(I)的传递函数为uf C K R s s CR Z Z s G 1,1001.011)(111112==-=-=-=其对应的模拟电路及SIMULINK 图形如图3所示。
4．微分环节(D)的传递函数为uf C K R s s C R Z Z s G 10,100)(111112==-=-=-= uf C C 01.012=<<其对应的模拟电路及SIMULINK 图形如图4所示。
5．比例+微分环节（PD ）的传递函数为)11.0()1()(111212+-=+-=-=s s C R R R Z Z s G uf C C uf C K R R 01.010,10012121=<<===其对应的模拟电路及SIMULINK 图形如图5所示。
图3 积分环节的模拟电路及及SIMULINK 图形图4 微分环节的模拟电路及及SIMULINK 图形6．比例+积分环节（PI ）的传递函数为 )11(1)(11212sR s C R Z Z s G +-=+-=-= uf C K R R 10,100121=== 其对应的模拟电路及SIMULINK 图形如图6所示。
matlab zplane函数用法
matlab zplane函数用法
`zplane`是MATLAB中用于绘制数字滤波器零极点图的函数。
以下是关于`zplane`函数的一些基本用法：
```matlab
%语法:
zplane(b,a)
%示例:
b=[1,-1];%分子多项式系数
a=[1,-0.5];%分母多项式系数
zplane(b,a)
title('零极点图')%添加图的标题
xlabel('实部')%添加x轴标签
ylabel('虚部')%添加y轴标签
```
这里，`b`是数字滤波器的分子多项式系数，`a`是数字滤波器的分母多项式系数。
`zplane`函数会绘制这些系数所对应的零极点图。
-`b`和`a`的系数是多项式的系数，例如，`[1,-1]`对应于多项式`1-z^(-1)`。
-在零极点图中，零点用"o"表示，极点用"x"表示。
-通过在`zplane`函数后面添加标题和轴标签，可以使图更具可读性。
此外，你还可以使用其他参数和选项来调整图的外观。
可以通过输入`help zplane`来查看更多详细的用法和选项。
用matlab绘制差分方程Z变换-反变换-zplane-residuez-tf2zp-zp2tf-tf2sos-sos2tf-幅相频谱等等
《数字信号处理》（一） 实验目的使用ztrans,iztrans 函数分别求出离散时间信号的Z 变换和Z 反变换的结果,并用pretty 函数进行结果美化。
编写函数时养成良好的注释习惯,有利于对函数的理解。
复习MATLAB 的基本应用，如：help,可以帮助查询相关的函数的使用方法,巩固理论知识中的离散时间信号的传递函数与二次项式之间的转换，以及使用zplane 函数画出相关系统的零极点分布图,根据零极点的分布情况估计系统的滤波特性。
（二） 程序的运行与截图实验项目一Z 变换（1）求)(])31()21[()(n u n x nn += Z 变换clear all;close all;clc;syms nf=0.5^n+(1/3)^n; %定义离散信号F=ztrans(f) %z 变换pretty(F); 运算结果F（2）4)(n n x = Z 变换clear all ;close all ;clc;syms nf=n^4; %定义离散信号F=ztrans(f) %Z 变换 pretty(F)运算结果（3）)sin()(b an n x += Z 变换clear all;close all;clc;syms a b nf = sin(a*n+b) %定义离散信号F=ztrans(f) %Z 变换pretty(F)运算结果实验项目二Z 反变换（1）2)2(2)(-=z z z X Z 反变换 clear all;close all;clc;syms k zFz=2*z/(z-2)^2; %定义Z 反变换表达式fk=iztrans(Fz,k) %Z 反变换pretty(fk);运算结果（2）12)1()(2++-=z z z z z X Z 反变换 clear all;close all;clc;syms k zFz=z*(z-1)/(z^2+2*z+1); %定义Z 反变换表达式fk=iztrans(Fz,k) %Z 反变换pretty(fk);运算结果f（3） 211cos 211)(---+-+=zz z z X ω Z 反变换 clear all;close all;clc;syms k z wFz=(1+z^(-1))/(1-2*z^-1*cos(w)+z^-2); %定义Z 反变换表达式fk=iztrans(Fz,k) %Z 反变换pretty(fk);运算结果实验项目三各种模型之间的变换2)2)(1(10)(--=z z z z H =4851023-+-z z z z (1)clear all;close all;clc;b=[0 0 10 0];%分子的系数数组a=[1 -5 8 -4]; %分母的系数数组zplane(b,a)% 使用zplane 函数绘制如下系统的零极点分布图 运算结果(2)clear all ;close all ;clc;b=[0 0 10 0]; %分子的系数数组a=[1 -5 8 -4]; %分母的系数数组[r,p,c]=residuez(b,a) %使用matlab 中的residuez 函数，将)(z H 分解成为多个简单有理分式之和运算结果r =-15.00005.000010.0000p =2.00002.00001.0000c =(3)clear all;close all;clc;b=[0 0 10 0]; %分子的系数数组a=[1 -5 8 -4]; %分母的系数数组[z,p,k]=tf2zp(b,a) %使用tf2zp求出系统函数的零、极点和增益运算结果z =p =2.00002.00001.0000k =10(4)clear all;close all;clc;z=[1;-3];%零点,列向量p=[2; -4];%极点,列向量k=5; %增益[b,a] = zp2tf(z,p,k) %根据求出的零、极点和增益，然后自学使用zp2tf还原出分子和分母的系数运算结果(5)clear all;close all;clc;b=[0 0 10 0]; %分子的系数数组a=[1 -5 8 -4]; %分母的系数数组[sos,g]=tf2sos(b,a) %使用tf2sos将系统函数分解成一系列二阶子系统的级联形式运算结果sos =0 1.0000 0 1.0000 -2.0000 00 1.0000 0 1.0000 -3.0000 2.0000g =10(6)clear all;close all;clc;sos=[0 1.0000 0 1.0000 -2.0000 0;0 1.0000 0 1.0000 -3.0000 2.0000];g=10;%增益[b,a]=sos2tf(sos,g) %根据求出的一系列二阶子系统，使用sos2tf还原出H分子和分母的系数)(z运算结果b =0 0 10 0a =1 -5 8 -4(7)clear all;close all;clc;b=[0 0 10 0]; %分子的系数数组a=[1 -5 8 -4]; %分母的系数数组n=(0:500)*pi/500; %在pi范围内取501个采样点[h,w]=freqz(b,a,n);%求系统的频率响应subplot(2,1,1),plot(n/pi,abs(h));grid %作系统的幅度频响图axis([0,1,1.1*min(abs(h)),1.1*max(abs(h))]);ylabel(‘幅度’);subplot(2,1,2),plot(n/pi,angle(h));grid %作系统的相位频响图axis([0,1,1.1*min(angle(h)),1.1*max(angle(h))]);ylabel(‘相位’);xlabel(‘以pi为单位的频率’);运行结果2221)(232+++++=z z z z z z H (1)clear all ;close all ;clc;b=[0 1 1 1]; %分子的系数数组a=[1 2 2 2]; %分母的系数数组zplane(b,a)% 使用zplane 函数绘制如下系统的零极点分布图运行结果(2)clear all ;close all ;clc;b=[0 1 1 1]; %分子的系数数组a=[1 2 2 2]; %分母的系数数组[r,p,k]=residuez(b,a) %使用matlab 中的residuez 函数，将)(z H 分解成为多个简单有理分式之和运行结果r =-0.4006-0.0497 - 0.1609i-0.0497 + 0.1609ip =-1.5437-0.2282 + 1.1151i-0.2282 - 1.1151ik =0.5000(3)clear all;close all;clc;b=[0 1 1 1]; %分子的系数数组a=[1 2 2 2]; %分母的系数数组[z,p,k]=tf2zp(b,a) %使用tf2zp求出系统函数的零、极点和增益运行结果z =-0.5000 + 0.8660i-0.5000 - 0.8660ip =-1.5437-0.2282 + 1.1151i-0.2282 - 1.1151ik =1(4)clear all;close all;clc;z=[-0.5000 + 0.8660i-0.5000 - 0.8660i];p=[-1.5437-0.2282 + 1.1151i-0.2282 - 1.1151i ];k=1;[b,a]=zp2tf(z,p,k) %根据求出的零、极点和增益，使用zp2tf 还原出)(z H 分子和分母的系数运行结果b =0 1.0000 1.0000 1.0000a =1.00002.0001 2.0001 1.9999(5)clear all ;close all ;clc;b=[0 1 1 1]; %分子的系数数组a=[1 2 2 2]; %分母的系数数组[sos,g]=tf2sos(b,a) %使用tf2sos 将系统函数分解成一系列二阶子系统的级联形式运行结果sos =0 1.0000 0 1.0000 1.5437 01.0000 1.0000 1.0000 1.0000 0.4563 1.2956g =1(6)clear all ;close all ;clc;sos=[ 0 1.0000 0 1.0000 1.5437 0;1.0000 1.0000 1.0000 1.0000 0.4563 1.2956];g=1;[b,a]=sos2tf(sos,g) %根据求出的一系列二阶子系统，自学使用sos2tf 还原出)(z H 分子和分母的系数运行结果b =0 1 1 1a =1.00002.0000 2.0000 2.0000(7)clear all;close all;clc;b=[0 1 1 1]; %分子的系数数组a=[1 2 2 2]; %分母的系数数组n=(0:500)*pi/500; %在pi范围内取501个采样点[h,w]=freqz(b,a,n);%求系统的频率响应subplot(2,1,1),plot(n/pi,abs(h));grid %作系统的幅度频响图axis([0,1,1.1*min(abs(h)),1.1*max(abs(h))]);ylabel('幅度');subplot(2,1,2),plot(n/pi,angle(h));grid %作系统的相位频响图axis([0,1,1.1*min(angle(h)),1.1*max(angle(h))]);ylabel('相位');xlabel('以pi为单位的频率');运行结果实验项目四根据零极点分布图估计系统的滤波特性。
数字信号处理-原理与实践(方勇)习题答案(1-2)
第一章1-1 有一个连续信号)2cos()(ψπ+=ft t x a ，式中Hz f 20=，2πψ=，（１） 求出)(t x a 的周期；（２） 用采样间隔s T 02.0=对)(t x a 进行采样，写出采样信号)(ˆt xa 的表达式； （３） 画出对应)(ˆt xa 的时域离散信号（序列）)(n x 的波形，并求出)(n x 的周期。
解：（1）)(t x a 的周期是s fT a 05.01==（2）∑∞-∞=-+=n a nT t fnT t x)()2cos()(ˆδψπ∑∞-∞=-+=n nT t nT )()40cos(δψπ（3）)(n x 的数字频率为πω8.0=，252=ωπ周期5=N 。
)28.0cos()(ππ+=n n x ，画出其波形如题1-1图所示。
题1-1图 1-2 设)sin()(t t x a π=，()()sin()a s s x n x nT nT π==，其中s T 为采样周期。
（1）)(t x a 信号的模拟频率Ω为多少？ （2）Ω和ω的关系是什么？（3）当s T s 5.0=时，)(n x 的数字频率ω为多少？ 解：（1）)(t x a 的模拟频率s rad /π=Ω。
（2）Ω和ω的关系是：s T ⋅Ω=ω。
（3）当s T s 5.0=时，rad πω5.0=。
1-3 判断下面的序列是否是周期的，若是周期的，确定其周期。
（1）)873cos()(ππ-=n A n x ，A 为常数；（2）)81()(π-=n j e n x 。
解： （1）πω73=，3142=ωπ，这是有理数，因此是周期序列，周期是14=T ； （2）81=ω，πωπ162=，这是无理数，因此是非周期序列。
1-4 研究一个线性时不变系统，其单位脉冲响应为指数序列)()(n u a n h n =，10<<a 。
对于矩阵输入序列，1,01()0N n N R n ≤≤-⎧=⎨⎩,其他 求出输出序列，并用MA TLAB 计算，比较其结果。
根轨迹分析的MATLAB实现
函数命令调用格式：
[k,poles]=rlocfind(sys) [k,poles]=rlocfind(sys,p)
函数命令使用说明：
rlocfind（）函数命令可计算出与根轨迹上极点（位 置为poles）相对应的根轨迹增益（k）。rlocfind（）函 数既适用于连续时间系统，也适用于离散时间系统。
pzmap(a,b,c,d)函数可以在复平面内绘制用状态空间 模型描述系统的零极点图，对于MIMO系统，可绘制所有 输入到输出间的传递零点。
pzmap(sys)函数可以在复平面里绘制以传递函数模 型sys表示开环系统的零极点。传递函数模型sys即 G(s)=num(s)/den(s)。
pzmap(p,z)函数可在复平面里绘制零极点图，其中 行矢量p为极点位置，列矢量z为零点位置。这个函数命令 用于直接绘制给定的零极点图。
[k,poles]=rlocfind(sys,p)函数可对指定根计算对应 的增益与根矢量p。
[k,poles]=rlocfind(sys)函数输入参量sys可以是由函 数tf（）、zpk（）、ss（）中任何一个建立的LTI对象模 型。函数命令执行后，可在根轨迹图形窗口中显示十字形 光标，当用户选择根轨迹上一点时，其相应的增益由k记 录，与增益相关的所有极点记录在poles中。
2．求系统根轨迹的函数rlocus()
函数命令调用格式：
[r,k]=rlocus(a,b,c,d) [r,k]=rlocus(sys) [r,k]=rlocus(a,b,c,d,k) [r,k]=rlocus(num,den,k)
函数命令使用说明：
rlocus()函数命令用来绘制SISO系统的根轨迹图。 rlocus()函数既可适用于连续时间系统，也适用于离散时 间系统。
实验三 利用MATLAB进行时域分析
实验三
一、实验目的
(1)学会使用MATLAB编程绘制控制系统的单位阶跃响应曲线；
(2)研究二阶控制系统中，、n对系统动态特性和时域指标的影响；
(3)掌握准确读取动态特性指标的方法；
lab2='zunibi = 0.25';text(0.3,1.5,lab2),
>> lab3='zunibi = 0.5';text(0.3,1.2,lab3),
lab4='zunibi = 0.75';text(0,3,1.05,lab4),
lab5='zunibi = 1';text(0,35,0.9,lab5),
den=[1,4,25];
step(num,den)
grid%绘制网格线。
title(¹Unit-Step Response of G(s)=25/(s^2+4s+25) ¹) %图像标题
图3-2MATLAB绘制的响应曲线
还可以用下面的语句来得出阶跃响应曲线
>> G=tf([0,0,25],[1,4,25]);
den = [1 2*sigma*10 100];
damp(den)
sys = tf(num,den);
i=i+1;
step(sys,2)
hold on
end
Eigenvalue Damping Freq. (rad/s)
0.00e+000 + 1.00e+001i 0.00e+000 1.00e+001
MATLAB应用-传递函数与画图
用MATLAB求系统的零点、极点 及特征多项式
• 部分分式展开
n−1 +⋅⋅⋅+bn−1s +bn M(s) num b0sn +bs 1 = = 考虑传递函数：N(s) den a sn +asn−1 +⋅⋅⋅+a s +a 0 1 n−1 n
带有选项的曲线绘制命令的调用格式： plot(x1,y1,s1,x2,y2,s2,…)
用MATLAB绘制二维图形
• 子图的命令
MATLAB允许将一个图形窗口按矩阵形式分成多个子 窗口，分别显示多个图形，需要用到subplot( )函数。 该函数把一个图形窗口分割成m*n个子绘图区域。调 用格式：subplot (m, n, k) 例如：subplot(4,3,6)表示将窗口分割成4*3个部 分，在第6 个部分上绘制图像。 通过参考数 k 可以调用各个绘图区域，子绘图区域按 注意：MATLAB最多允许9*9的分割。 行从左到右从上到下编号。
10s + 5 G2 ( s ) = s
H (s) =
1 0.01s + 1
用MATLAB建立传递函数模型
• 多项式模型TF和零极点模型ZPK的转换 • TF对象转换成ZPK对象 调用格式为：G1=zpk(G) 例
6.8s 2 + 61.2 s + 95.2 G (s) = 4 s + 7.5s 3 + 22 s 2 + 19.5s
用MATLAB建立传递函数模型
• 多项式模型—TF对象（单入—单出系统） 线性时不变（LTI）系统的传递函数模型：
用MATLAB画零极点图.
>> num=[1 2 1]; >> den=[1 5 3 8 9]; >> t=0.2; >> G1=tf(num,den) %没有采样时间 t，则显示为传递函数,自变量为 s
Transfer function: s^2 + 2 s + 1
----------------------------s^4 + 5 s^3 + 3 s^2 + 8 s + 9
• s^4 + 3 s^3 + 8 s^2 + 4 s + 2
•
• >> G1=zpk(G);
• >> z=G1.z;
• >> p=G1.p;
• >> Z=z{:};
• >> P=p{:};
• >> k=G1.k;
• >> pzmap(G);
• >> pzmap(G1);
• >> grid on
将传递函数化为零极点增益模型 并绘制零极点图
>> num2=[1 7 12]; >> den2=[1 3 2]; >> G2=tf(num2,den2)；
>> G=series(G1,G2) 或 >> G=G1*G2 %串联
Transfer function: 2 s^4 + 20 s^3 + 71 s^2 + 107 s + 60
---------------------------------------s^5 + 7 s^4 + 19 s^3 + 25 s^2 + 16 s + 4
Matlab绘制零极点图、判断稳定
4
Xidian University, ICIE. All Rights Reserved
Matlab绘制零极点图、判断稳定
系统的零极图，如图所示，可见系统有4个极点，
2个零点，其中在虚轴上有一对共轭极点，故该系统是
不稳定的。
1.5 0.89
不稳定系统
0.81
0.7 0.56 0.38 0.2
1 0.95
Imaginary Axis
0.5 0.988
03
2.5
2Hale Waihona Puke 1.510.5
-0.5 0.988
-1 0.95
0.89
0.81
0.7 0.56 0.38 0.2
-1.5
-3
-2.5
-2
-1.5
-1
-0.5
0
0.5
Real Axis
5
Xidian University, ICIE. All Rights Reserved
Matlab绘制零极点图、判断稳定
K1.19 Matlab绘制零极点图、判断稳定
例：利用MATLAB画出系统
H
(s)
s4
s2 3s3
4s 3 4s2 6s
4
的零极点图，并判断系统的稳定性。
解：b=[1 4 3]; %分子系数，按降幂顺序排列 a=[1 3 4 6 4]; %分母系数，按降幂顺序排列 sys=tf(b,a) pzmap(sys);sgrid; azp=roots(a); %求出极点azp，在左半平面即为稳定
title('不稳定系统'); elseif wd==1
title('稳定系统'); end
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
• >> num=[1 2 1]; • >> den=[1 5 3 8 9]; • >> [z,p,k]=tf2zp(num,den)%化传递函数为零极点增益形式
z= -1 -1
p= -4.6347 0.3174 + 1.3569i 0.3174 - 1.3569i -1.0000
k= 1
化零极点增益形式为传递函数
将传递函数化为零极点增益模型 并绘制零极点图
• • • • •
Transfer function: 3 s^2 + 2 s + 8 ----------------------------s^4 + 3 s^3 + 8 s^2 + 4 s + 2 将上式化为零极点增益模型，并绘制零极 点图。
• >> num=[3 2 8]; • >> den=[1 3 8 4 2]; • >> G=tf(num,den)
>> Z=z{:} Z= -0.3333 + 1.5986i -0.3333 - 1.5986i 显示极点 显示零点
>> P=p{:}
P= -1.2496 + 2.2082i -1.2496 - 2.2082i -0.2504 + 0.4980i -0.2504 - 0.4980i
>> pzmap(G1);%绘制零极点图 >> grid on %打开绘图网络
ζwn
M e p
/ 1 2
化传递函数为零极点形式
Transfer function: z^2 + 2 z + 1 ----------------------------z^4 + 5 z^3 + 3 z^2 + 8 z + 9 ，Sampling time: 0.2
>> num=[1 2 1]; >> den=[1 5 3 8 9]; >> t=0.2; >> G2=tf(num,den,t)
鼠标指向某零点或极点会显示该方 框。 极点、阻尼比ζ、超调量、固有频率 wn.
2 w n 1
Pole-Zero Map 2.5 2 1.5 1
wn
Imaginary Axis
0.5 0 -0.5 -1 -1.5 -2 -2.5 -1.4
ζ
由斜边和实部求出。
超调量由ζ求出。
-1.2 -1 -0.8 -0.6 -0.4 -0.2 0 Real Axis
Transfer function: 2 s^4 + 20 s^3 + 71 s^2 + 107 s + 60 ---------------------------------------s^5 + 7 s^4 + 19 s^3 + 25 s^2 + 16 s + 4 >> Gp=parallel(G1,G2) 或>> Gp=G1+G2 %并联 Transfer function: s^5 + 13 s^4 + 57 s^3 + 112 s^2 + 101 s + 34 -------------------------------------------s^5 + 7 s^4 + 19 s^3 + 25 s^2 + 16 s + 4
Transfer function: 3 s^2 + 2 s + 8 ----------------------------s^4 + 3 s^3 + 8 s^2 + 4 s + 2
求传递函数。命令后无 “;”则屏幕立即显示，否 则不显示。
• • • •
>> G1=zpk(G);%化为零极点增益形式 >> z=G1.z;%将G1零点存入z %下面语句也可获得零极点。 >> p=G1.p; %将G1极点点存入p >> [z,p,k]=tf2zp(num,den) >> k=G1.k;%将G1增益存入k
将传递函数化为零极点增益模型 并绘制零极点图
• • • • • • • • • • • • • • • • • >> num=[3 2 8]; >> den=[1 3 8 4 2]; >> G=tf(num,den) Transfer function: 3 s^2 + 2 s + 8 ----------------------------s^4 + 3 s^3 + 8 s^2 + 4 s + 2 >> G1=zpk(G); >> z=G1.z; >> p=G1.p; >> Z=z{:}; >> P=p{:}; >> k=G1.k; >> pzmap(G); >> pzmap(G1); >> grid on
• s^2 + 7 s + 12 • G2= -------------• s^2 + 3 s + 2 • 分别求G1、G2串、并联后的传递函数。
>> num1=[2 6 5]; >> den1=[1 4 5 2]; >> G1=tf(num1,den1)；
>> num2=[1 7 12]; >> den2=[1 3 2]; >> G2=tf(num2,den2)； >> G=series(G1,G2) 或 >> G=G1*G2 %串联
%没有采样时间 t，则显示为传递函数,自变量为 s
Transfer function: s^2 + 2 s + 1 ----------------------------s^4 + 5 s^3 + 3 s^2 + 8 s + 9
>> G11=zpk(G1) %求连续系统G1零极点
Zero/pole/gain: (s+1)^2 ---------------------------------------(s+4.635) (s+1) (s^2 - 0.6347s + 1.942)
Transfer function: z^2 + 2 z + 1 ----------------------------z^4 + 5 z^3 + 3 z^2 + 8 z + 9 Sampling time: 0.2
பைடு நூலகம்
%有采样时间 t，则显示为脉冲传递函数，自变量为 z
>> G22=zpk(G2) %求离散系统G2零极点
>> [num1,den1]=zp2tf(z,p,k) %化零极点增益形式为传递函数
num1 = 0 0 1 2 1
den1 = 1.0000 5.0000 3.0000 8.0000 9.0000
系统的串、并联的传递函数
• 2 s^2 + 6 s + 5 • G1= --------------------• s^3 + 4 s^2 + 5 s + 2
Zero/pole/gain: (z+1)^2 ---------------------------------------(z+4.635) (z+1) (z^2 - 0.6347z + 1.942) Sampling time: 0.2
>> num=[1 2 1]; >> den=[1 5 3 8 9]; >> t=0.2; >> G1=tf(num,den)