matlab中基于传递函数或者状态方程的幅频特性分析
来源：CSDN博客，作者：我不是瓜皮。
随便举一个传递函数为
在matlab中写入代码如下。
代码有两种写法，一种是直接写传递函数G 的表达式，简单粗暴。tf 其实就是transfer function的缩写，即传递函数。
s=tf('s');
G=(s^2 1)*(s^2 s 1)/((s 1)*(s^2 5*s 7));
bode(G);
grid
xlabel('频率f');title(' ');
另一种写代码的方法：
num=conv([1,0,1],[1,1,1]);
den=conv([1,1],[1,5,7]);
bode(tf(num,den));
grid
xlabel('频率f');title(' ');
这里conv命令就是将分子分母的因式乘出来，然后按照s的幂次降序得到num 和den 两个一维矩阵，num=[1,1,2,1,1]，den=[1,6,12,7]。再直观一点说，这时候
此时，tf (num, den) 和第一种中的G 就是等价的。运行得：
有的时候我们不需要bode图中的相角图，只想要频率图，只需要将bode(G) 换成
h = bodeplot(G);
setoptions(h,'FreqUnits','Hz','PhaseVisible','off');
得出来：
对于一般的状态方程
传递函数G=Y/U。
[num,den]=ss2tf(A,B,C,D)
这里求出的num 和den 与上述的一样，只不过此时的传递函数就取决于你所设定的Y 有几个输出量了，如果Y 有三个量，那么，num 和den 将会是三维矩阵，对应着三个传递函数，剩下的工作就和上面的一样了。
至于将传递函数转成状态方程，就是上面的逆过程 [A,B,C,D]=tf2ss(num,den)。不过我目前还没用到，因为我觉得根据自己所需要的状态变量最后得出的状态方程会不一样，用matlab中求出的状态空间只是其中一种。
原文链接：
