%目標函數(shù)計算程序
function Q=SYS_FUC(K1,T1)
num=[10*K1*T1 0];
den=[2*T1 2+T1 1 0];
[a b c d]=tf2ss(num,den);%傳遞函數(shù)轉換成狀態(tài)空間表示形式

%四階龍格庫塔法計算程序
h=0.003;
n=10000;
y1=ones(1,n);
x=[0;0;0];
y1(1,1)=0;
for i=2:n
    u=1-y1(1,i-1);
    k1=a*x+b*u;
    k2=a*(x+h*k1/2)+b*u;
    k3=a*(x+h*k2/2)+b*u;
    k4=a*(x+h*k3)+b*u;
    x1=x+h/6*(k1+2*k2+2*k3+k4);
    y=c*x+d*u;
    x=x1;
    y1(1,i)=y;
end
%計算Q的值
Q=0;
for j=1:n
    f=abs(1-y1(1,j));
    %Q=Q+h*f*j*(150-0)/n;
    Q=Q+j*h*h*f;
end