syszh_fc2.m

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

%用于画最好点响应曲线
function y=SYSZH_FC2(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)
hold on
plot(t,y1(1,:))
text(4.4,0.67,'\leftarrow 极好点曲线','FontSize',10,'color','b')
xlabel('时间（秒）')
ylabel('输出')
title('单纯形法寻优（极坏点和极优点曲线比较）')