matlab：G(s)=30(2s+1)/[s(0.1s+1)(0.2s+1)(26s+1)],要求画bode图，并且求出截止频率，相角裕度，幅值裕度

Question

求完整的matlab程序

Answer 1

很简单 你先对此方程计算展开 得出s的开方次数 因为你这个我没计算过 我只举个简单例子
G（s)=(3S^3+S^2+1)/(2S^4+S^3+9S^2+S^1)
num=[3 1 0 1];%分子
den=[2 1 9 1 0];%分母
figure(1)
nyquist(num,den);
figure(2)
bode(num,den);
figure(3)
margin(num,den)
figure(4)
nichols(num,den); ngrid
运行就能出图 这么说应该很简单了吧

追问

谢谢，关键是怎么写剪切频率，幅频特性，相频特性的计算程序，能再帮一下么？

追答

不用计算程序啊 只有有一个计算式就行
G(s)=分子/分母
[ ]内是各个的系数比如[2 1 9 1 0]就是表示2S^4+S^3+9S^2+S^1
只要你把上面那一串G(s)=30(2s+1)/[s(0.1s+1)(0.2s+1)(26s+1)]具体展开列出式子就行
直接就能出图

追问

厄，我知道。这样说吧，bode图上的Gm和Pm两组数据是对应的剪切频率和幅频特性，相频特性吗?是的话为什么和我计算的不一样

追答

我个人觉得吧。。如果学校里完成任务的话能出结果就可以了，因为是算式什么都是对的。。如果楼主是个人兴趣的话。。就只好你钻研一下了。。

Answer 2

num=[60 30];
den=[13/25 391/50 263/10 1 0];
f=tf(num,den)
bode(f)
S= allmargin(f)
S中给出了各种参数