DMC最基础的simulink仿真程序

function y=my_fun1213() P=5; M=4; N=8; C=[1 0 0 0 0 0 0 0 ]; H=[0.01 ;0.01; 0.01 ;0.01;0.01; 0.01 ;0.01;0.01 ]; A=[1 0.5 0.75 0.625 0.7 0.65 0.675 0.67]; a=[ A(1) 0 0 0 ; A(2) A(1) 0 0 ; A(3) A(2) A(1) 0 ; A(4) A(3) A(2) A(1); A(5) A(4) A(3) A(2)];%P*M的矩阵 q=[0.01 0 0 0 0; 0 0.01 0 0 0; 0 0 0.01 0 0; 0 0 0 0.01 0; 0 0 0 0 0.01];%对角矩阵P*P 也可以用0.01*eye(P) r=[0.01 0 0 0 ; 0 0.01 0 0 ; 0 0 0.01 0 ; 0 0 0 0.01 ];%对角矩阵M*M 也可以用0.01*eye(M) c= [1; 0; 0; 0 ];% M*1 也可以设置 c1=zeros(1,M-1) c=[1,c1] y=c'*((a'*q*a+r)^-1)*a'*q