mpc建模仿真2

%%%%%%系统模型建立与转换函数
%例8.5传递函数转为mpc状态空间模型
num=[1 3 1];den=[1 2 2 1];[A B C D]=tf2ss(num,den);%传递函数转为状态空间方程
pmod=ss2mod(A,B,C,D)%状态空间方程转为mpc状态空间方程

%例8.6传递函数模型转为mpc传递函数模型
num=[1 1];
den=[1 4 4];
g=poly2tfd(num,den,0,0.5)

%例8.7mpc传递函数转为mpc状态空间模型
num=[1 1];den=[1 3 6];
g=poly2tfd(num,den,0,0);%传递函数转为mpc传递函数
mod1=tfd2mod(0.1,1,g)

%例8.8传递函数转为阶跃响应模型
num=[1 2];
den=[1 3 1];tf1=poly2tfd(num,den,0,0);%传递函数转为mpc传递函数
plant=tfd2step(5,0.1,1,tf1);
plotstep(plant)%系统阶跃响应曲线

%例8.9多变量系统给传递函数转换为阶跃响应模型
clear
clc
tf1=poly2tfd(2,[1,1],0,0);tf2=poly2tfd(3,[1,2],0,0);
tf3=poly2tfd(1,[1,4],0,0);tf4=poly2tfd(1,[1,1],0,0);
plant=tfd2step(5,0.1,2,tf1,tf2,tf3,tf4)%转为阶跃响应模型
plotstep(plant)
%为什么只有两幅图

%例8.10状态空间转为mpc状态空间模型
A=[-2,-1;1,0]; B=[1;0];C=[1 2];D=0
plant=ss2step(A,B,C,D,5,0,0.2,1);
plotstep(plant)

%例8.11mpc状态空间模型改变采样周期
num=[1 1];den=[1 3 6];
g=poly2tfd(num,den,0,0);%将传递函数转为Mpc传递函数
mod1=tfd2mod(0.1,1,g);%采样周期，输出个数，传递函数%mpc传递函数转为mpc状态空间
plant1=mod2step(mod1,5,0.1);
plotstep(plant1)
%mod2=mod2mod(mod1,0.5)%改变系统采样周期
%plant2=mod2step(mod2,5,0.1);%mpc状态空间模型，截止时间，采样时间%计算阶跃响应模型
%plotstep(plant2)