syszl_fc1.m

来自「MatLab环境下编程实现简单交互式单纯形法进行寻优并显示寻优结果。」· M 代码 · 共 33 行

%用于画最坏点响应曲线
function y=SYSZL_FC1(K1,T1)
num=[10*K1*T1 0];
den=[2*T1 2+T1 1 0];
[a b c d]=tf2ss(num,den);%传递函数转换成状态空间表示形式

%四阶龙格库塔法计算程序
hs = 0.003;
num = 10000;
y1=ones(1,num);
x=[0;0;0];
y1(1,1)=0;
t=zeros(1,num);
for i=2:num
    u=1-y1(1,i-1);
    k1=a*x+b*u;
    k2=a*(x+hs*k1/2)+b*u;
    k3=a*(x+hs*k2/2)+b*u;
    k4=a*(x+hs*k3)+b*u;
    x1=x+hs/6*(k1+2*k2+2*k3+k4);
    y=c*x+d*u;
    x=x1;
    y1(1,i)=y;
    t(1,i)=hs*i;
end
figure(1);

plot(t,y1(1,:),'r')

grid on

text(1.8,1.7,'\leftarrow 极坏点曲线','FontSize',10,'color','r')