matlab中分数阶线性系统的bode图程序看不懂，求高手指点，越详细越好，多谢啦！

function H=f_bode(a,eta,b,gamma,w,Gc)

%分数笑旦阶线性系统的bode图程序

%a,b,eta,gamma均为整数，w实值矢量，Gc权重，输出H是frd模型，例碰腊扰如：

%a=3b=1eta=3gamma=2freq = logspace(1,2)w = 0.05*(freq)f_bode(a,eta,b,gamma,w)

%a=3b=1eta=3gamma=2freq = linspace(0.01,0.2)w = 0.05*(freq)f_bode(a,eta,b,gamma,w)

H=zeros(size(w))

for i=1:length(w)

P=b*((sqrt(-1)*w(i)).^gamma.')

Q=a*((sqrt(-1)*w(i)).^eta.')

H(i)=P/Q

end

H=frd(H,w)%生成频率响应局野数据模型

if nargin==6,

H=H*Gc

end

if nargout==0,

bode(H)

end

前两天刚好回答了一个相同的问题（编号689391549523114324），把答案转贴过来，供参考。

这个题直接把传递函数表示出来然后调用相应的函数即可。

参考代码：

G1=tf([1 2],[1 31 229 旦码198])

s=tf('s')

G=feedback(G1,1)/s

bode(1000*G)

figure,margin(1000*G)

figure,rlocus(G)

（1）从margin的绘图结果可知幅值裕度为氏旁14.1dB，相角裕度为40.6度。

（2）从根轨迹图可见，闭歼迟橡环系统稳定的条件是K<=4890（近似值）。

---