高分求救！！！急啊~~~MATLAB频域校正，关于子函数的调用，大家都来帮帮忙
书上给出的程序是这样的
子函数程序代码：
function Gc=cqjz_frequency(G,kc,yPm)
G=tf(G)
[mag,pha,w]=bode(G*kc);Mag=20*log10(mag)
[Gm,Pm.Wcg,Wcp]=margin(G*kc)
phi=(yPm-getfield(Pm,'Wcg'))*pi/180;alpha=(1+sin(phi))/(1-sin(phi))
Mn=-10*log10(alpha);Wcgn=spline(Mag,w,Mn)
T=1/Wcgn/sqrt(alpha);Tz=alpha*T
Gc=tf([Tz,1],[T,1])
主程序代码：
num=2;den=conv([1,0],conv([0.3,1],[0.1,1]));G=tf(num,den)
kc=3;yPm=45+12
Gc=cqjz_frequency(G,kc,yPm)
G=G*kc;GGc=G*Gc
Gy_close=feedback(G,1)
Gx_close=feedback(GGc,1)
figure(1);step(Gx_close,'b');hold on
step(Gy_close,'r');grid
gtext('校正前的');gtext('校正后的')
figure(2);bode(G,'b');hold on
bode(GGc,'r');grid
gtext('校正前的');gtext('校正后的');gtext('校正前的');gtext('校正后的');
figure(3);nyquist(Gx_close,'b');hold on
nyquist(Gy_close,'r');grid
gtext('校正前的');gtext('校正后的')
然后、然后该怎么用啊？把程序输进去都得不出结果的。子程序到底该怎么用？哪位高手给详细解答一下啊，就针对这个题吧，谢谢了，急需啊！
好像没有问题啊。这是画三个图，闭环的阶跃曲线，波德图，奈奎斯特图。
吧函数保存，把主程序直接拷到命令下，回车即可。
看这图，
