马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
function [sys,x0,str,ts]=PF_1(t,x,u,flag,Ts)
global A,B,C,D,G,H
A=[0,1,0,0;0,0,-3.8354,0;0,0,0,1;0,0,28.0160,0];
B=[0;1.7392;0;3.5735];
C=[1,0,0,0;0,0,1,0];
D=[0;0];
G=exp(A*Ts);
H=quad(@fun,0,Ts);
switch flag,
case 0,
[sys,x0,str,ts]=mdlInitializeSizes;
case 2,
sys=mdlUpdate(t,x,u,Ts);
case 3,
sys=mdlOutputs(t,x,u,Ts);
case 9,
sys=[];
otherwise,
error=(['Unhandled flag=',num2str(flag)]);
end
function [sys,x0,str,ts]=mdlInitializeSizes
sizes=simsizes;
sizes.NumConStates=0;
sizes.NumDiscStates=1;
sizes.NumOutputs=1;
sizes.NumInputs=3;
sizes.DirFeedthrough=1;
sizes.NumSampleTimes=1;
sys=simsizes(sizes);
x0=u(3);
str=[];
ts=[1 0];
function sys=mdlUpdates(t,x,u,Ts)
sys=G*x+H*u(2);
function sys=mdlOutPuts(t,x,u,Ts)
sys=C*G*G*x+C*G*H*u(2)+C*H*u(1);
function y=fun(t);
y=exp(A*Ts-A*t)*B;
以上是我自己编的程序,编辑S-函数模块时出现如下问题
??? Undefined function or variable "B".
Error in ==> PF_1 at 2
global A,B,C,D,G,H
高手帮我指点下怎么改谢谢了~~~~~ |