动态矩阵控制算法(DMC)

% DMC.m 动态矩阵控制（DMC）

num=1;

den=[1 1 1];

G=tf(num,den); %连续系统

Ts=0.2; %采样时间Ts

G=c2d(G,Ts); %被控对象离散化[num,den,]=tfdata(G,'v');

N=60; %建模时域N

[a]=step(G,1*Ts:Ts:N*Ts); %计算模型向量a

M=2; %控制时域

P=15; %优化时域

for j=1:M

for i=1:P-j+1

A(i+j-1,j)=a(i,1);

end

end %动态矩阵A

Q=1*eye(P); %误差权矩阵Q

R=1*eye(M); %控制权矩阵R

C=[1,zeros(1,M-1)]; %取首元素向量C 1*M

E=[1,zeros(1,N-1)]; %取首元素向量E 1*N

d=C*(A'*Q*A+R)^(-1)*A'*Q; %控制向量d=[d1 d2 ...dp] h=1*ones(1,N); %校正向量h(N维列向量) I=[eye(P,P),zeros(P,N-P)]; %Yp0=I*YNo

S=[[zeros(N-1,1) eye(N-1)];[zeros(1,N-1),1]]; %N*N移位阵S sim('DMCsimulink') %运行siumlink文件

subplot(2,1,1); %图形显示plot(y,'LineWidth',2);

hold on;

plot(w,':r','LineWidth',2);

xlabel('\fontsize{15}k');

ylabel('\fontsize{15}y,w');

legend('输出值','设定值')

grid on;

subplot(2,1,2);

plot(u,'g','LineWidth',2);

xlabel('\fontsize{15}k');

ylabel('\fontsize{15}u');

grid on;