dmctuningwboralatiqi.m

to calculate a tuning parameter of dmc and mpc controller

% DMC TUNING for Wood and Berry Distillation Column
display('DMC TUNING for Wood and Berry Distillation Column')
%kp1=12.8;thau1=16.7;theta1=1;
%kp2=-18.9;thau2=21;theta2=3;
%kp3=6.6;thau3=10.9;theta3=7;
%kp4=-19.4;thau4=14.4;theta4=3;
kp=[12.8 -18.9;6.6 -19.4];
thau=[16.7 21;10.9 14.4];
theta=[1 3;7 3];
for r=1:2
    for s=1:2
        T(r,s)=max(0.1*thau(r,s),0.5*theta(r,s));
    end
end
T=[T(1,:) T(2,:)];
T=min(T)
for r=1:2
    for s=1:2
        P(r,s)=(5*thau(r,s)/T)+(theta(r,s)/T)+1;
    end
end
P=[P(1,:) P(2,:)];
P=max(P)
for r=1:2
    for s=1:2
        M(r,s)=(thau(r,s)/T)+(theta(r,s)/T)+1;
    end
end
M=[M(1,:) M(2,:)];
M=max(M)
gamma=1;
for r=1:2
    for s=1:2
        Lambda(r,s)=gamma*kp(r,s).^2*((P-((theta(r,s)/T)+1)-1.5*(thau(r,s)/T)+2-((M-1)/2)));
    end
end
Lambda=[Lambda(1,:) Lambda(2,:)];
Lambda=sum(Lambda);
Lambda=(M/10)*Lambda

% DMC TUNING for OR 3x3
display('DMC TUNING for OR 3x3')
%kp1=0.66;thau1=6.7;theta1=2.6;
%kp2=-0.61;thau2=8.64;theta2=3.5;
%kp3=-0.0049;thau3=9.06;theta3=1;
%kp4=1.11;thau4=3.25;theta4=6.5;
%kp5=-2.3;thau5=5;theta5=3;
%kp6=-0.01;thau6=7.09;theta6=1.2;
%kp7=-34.68;thau7=8.15;theta7=9.2;
%kp8=46.2;thau8=10.9;theta8=9.4;
%kp9=0.87;thau9=9.12;theta9=0.29;
kp=[0.66 -0.61 -0.0049;1.11 -2.3 -0.01;-34.68 46.2 0.87];
thau=[6.7 8.64 9.06;3.25 5 7.09;8.15 10.9 9.12];
theta=[2.6 3.5 1;6.5 3 1.2;9.2 9.4 0.29];
for r=1:3
    for s=1:3
        T(r,s)=max(0.1*thau(r,s),0.5*theta(r,s));
    end
end
T=[T(1,:) T(2,:) T(3,:)];
T=min(T)
for r=1:3
    for s=1:3
        P(r,s)=(5*thau(r,s)/T)+(theta(r,s)/T)+1;
    end
end
P=[P(1,:) P(2,:) P(3,:)];
P=max(P)
for r=1:3
    for s=1:3
        M(r,s)=(thau(r,s)/T)+(theta(r,s)/T)+1;
    end
end
M=[M(1,:) M(2,:) M(3,:)];
M=max(M)
gamma=1;
for r=1:3
    for s=1:3
        Lambda(r,s)=gamma*kp(r,s).^2*((P-((theta(r,s)/T)+1)-1.5*(thau(r,s)/T)+2-((M-1)/2)));
    end
end
Lambda=[Lambda(1,:) Lambda(2,:) Lambda(3,:)];
Lambda=sum(Lambda);
Lambda=(M/10)*Lambda

% DMC TUNING for ALATIQI 4x4
display('DMC TUNING for ALATIQI 4x4')
%kp1=2.22;thau1=49.15;theta1=17.74;
%kp2=-2.94;thau2=35.56;theta2=5.27;
%kp3=0.017;thau3=30.04;theta3=7.54;
%kp4=-0.64;thau4=48.02;theta4=34.26;
%kp5=-2.33;thau5=57.08;theta5=22.53;
%kp6=3.46;thau6=32;theta6=1.01;
%kp7=-0.51;thau7=52.64;theta7=23.35;
%kp8=1.68;thau8=41.999;theta8=6.19;
%kp9=-1.06;thau9=22.87;theta9=22.43;
%kp10=3.511;thau10=20.05;theta10=18.856;
%kp11=4.41;thau11=16.2;theta11=1.01;
%kp12=-5.38;thau12=17;theta12=0.5;
%kp13=-5.73;thau13=48.63;theta13=10.79;
%kp14=4.32;thau14=25;theta14=0.01;
%kp15=-1.25;thau15=43.95;theta15=11.39;
%kp16=4.737;thau16=46.65;theta16=6.38;
kp=[2.22 -2.94 0.017 -0.64;-2.33 3.46 -0.51 1.68;-1.06 3.511 4.41 -5.38;-5.73 4.32 -1.25 4.737];
thau=[49.15 35.56 30.04 48.02;57.08 32 52.64 41.999;22.87 20.05 16.2 17;48.63 25 43.95 46.65];
theta=[17.74 5.27 7.54 34.26;22.53 1.01 23.35 6.19;22.43 18.856 1.01 0.5;10.79 0.01 11.39 6.38];
for r=1:4
    for s=1:4
        T(r,s)=max(0.1*thau(r,s),0.5*theta(r,s));
    end
end
T=[T(1,:) T(2,:) T(3,:) T(4,:)];
T=min(T)
for r=1:4
    for s=1:4
        P(r,s)=(5*thau(r,s)/T)+(theta(r,s)/T)+1;
    end
end
P=[P(1,:) P(2,:) P(3,:) P(4,:)];
P=max(P)
for r=1:4
    for s=1:4
        M(r,s)=(thau(r,s)/T)+(theta(r,s)/T)+1;
    end
end
M=[M(1,:) M(2,:) M(3,:) M(4,:)];
M=max(M)
gamma=1;
for r=1:4
    for s=1:4
        Lambda(r,s)=gamma*kp(r,s).^2*((P-((theta(r,s)/T)+1)-1.5*(thau(r,s)/T)+2-((M-1)/2)));
    end
end
Lambda=[Lambda(1,:) Lambda(2,:) Lambda(3,:) Lambda(4,:)];
Lambda=sum(Lambda);
Lambda=(M/10)*Lambda