[MATLAB混合编程] Matlab关于设计一个控制系统通用分析和校正程序

[jhxjsy]

求助高手，急~~~
当已知开环传递函数可以按不同的分析方法对闭环系统的性能进行分析，具体内容包括：
时域分析：阶跃响应曲线；判稳定性；求稳态误差；求上升时间、调节时间、超调量；
频率响应法分析：作出伯德图；求相角稳定裕量及幅值交接频率与幅值稳定裕量及相角交接频率；
要求：
以上两种不同的分析方法可以在一个程序中任意选择。
子功能可用自编函数实现，由程序文件调用。
设计一个控制系统通用校正程序
要求：
通用程序包括频率响应法校正方法
频率响应法包括超前校正与滞后校正二种补偿器
根据不同的具体情况，不同补偿器可以任意选择
子功能可用自编函数实现，由程序文件调用

我做了一部分，但是不对，给高手参考，谢谢啦：
function [pos,tr,ts2,tp]=stepchar(y,t)
[mp,ind]=max(y);
dimt=length(t);
yss=y(dimt);
pos=100*(mp-yss)/yss;
tp=t(ind);
i=1;j=1;k=1;q=1;
while y(i)<0.1
i=i+1;
end
t1=t(i);
while y(j)<0.9
j=j+1;
end
t2=t(j);
tr=t2-t1;
i=dimt+1;n=0;
while n==0
i=i-1;
if i==1
n=1;
elseif y(i)>=1.02
n=1;
end
end
t1=t(i);i=dimt+1;n=0;
while n==0
i=i-1;
if y(i)<=0.98
n=1;
end
t2=t(i);
if t1>t2
ts2=t1;
else
ts2=t2;
end
end