MATLAB如何求传函,FeedBack system matlab整理

一 求传函

G(s) = (2S^2+3S+1)/(S^3+3.2S^2+2.4S)

den = [1 3.2 2.4 0] %分母

sys = tf(num,den) %tf求传函函数

利用conv( )求多项式系数

(S^2+2S+3)(S^2+0.1S+0.2)

num = conv([1 2 3],[1 0.1 0.2])

二 求取频率特性

函数:freqs(num,den,w)

num:传递函数分子系数

den:传递函数分母系数

w:构造的图表空间

G(s) = (2S^2+3S+1)/(S^3+3.2S^2+2.4S)

num = [2 3 1] %分子

den = [1 3.2 2.4 0] %分母

sys = tf(num,den) %tf求传函函数

w = logspace(-1,1)

freqs(num,den,w)

于是幅频和相频特性曲线就画出来了

三 特殊信号输入下的响应曲线

Matlab中已经封装好了一些特殊的信号模型

例如

impulse

单位脉冲响应

step

单位阶跃相应

G(s) = (2S^2+3S+1)/(S^3+3.2S^2+2.4S)下单位阶跃信号输入的系统响应曲线

step(sys)

即： 信号关键字(sys)

四 根轨迹

G(s) = K(s+2)/(s(s+1)(s+3))

num = [1 2];

den = conv([1 0],conv([1 1],[1 3]));

sys =tf(num,den);

rlocus(sys) %根轨迹

title('根轨迹图')

axis equal

rlocus()

五 奈奎斯特曲线

函数：

nyquist()

;

nyquist(sys); %sys为计算出的传函

axis equal

grid

六 波特图

函数:

bode()

;

bode(sys);%sys为传函

七 对数相图

函数:

nichols()

;

G(s) = (0.1S^3+0.3S^2+0.11S+1.7)/(2.01S^3+1.9S^2+3.12S+1)

num = [0.1 0.3 0.11 1.7];

den = [2.01 1.9 3.12 1];

sys = tf(num,den);

nichols(sys)

Simulink 控制仿真模块

在标准模块中，应用的一些主要模块位于 Continuous. Discontinuous. Discrete 模块中

1 Continuous模块

2 Discontinuous模块

提供模拟非线性环节的模块

Backlash 间隙环节

Dead Zone 死区环节

Saturation 饱和环节

3 Discrete模块

模块与Continuous类似，但是都是离散的。其中还有 Unit Delay 单位延迟，和Zero-Order Hold零阶保持器