声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1458|回复: 3

[其他相关] S函数问题求助

[复制链接]
发表于 2007-11-25 13:56 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
我要编写一个关于参数自整定的S函数
function [sys,x0,str,ts] = para(t,x,u,flag)
global Ke Kec Ku;
Ke=18;Kec=1;Ku=0.5;
switch flag,
case 0,
    [sys,x0,str,ts]=mdlInitializeSizes;
case 3,
    sys=mdlOutputs(t,x,u);
case {1,2,4,9},
    sys=[];
  otherwise
    error(['Unhandled flag = ',num2str(flag)]);
end
function [sys,x0,str,ts]=mdlInitializeSizes
sizes = simsizes;
sizes.NumContStates  = 0;
sizes.NumDiscStates  = 0;
sizes.NumOutputs     = 3;
sizes.NumInputs      = 2;
sizes.DirFeedthrough = 0;
sizes.NumSampleTimes = 1;   % at least one sample time is needed
sys = simsizes(sizes);
function sys=mdlOutputs(t,x,u)
%u(1)和u(2)分别为偏差e和偏差变化率ec
%sys(1),sys(2),sys(3)分别为调整后的Ke,Kec,Ku
if abs(u(1))>0.5|abs(u(2))>0.3
    sys(1)=0.2*Ke;sys(2)=0.3*Kec;sys(3)=Ku+0.5;
elseif abs(u(1))<0.2|abs(u(2))<0.1
    sys(1)=0.5*Ke;sys(2)=0.5*Kec;sys(3)=Ku+0.1;
end
sys=[];
主要就是最后这个输出,我想输出调整后的Ke,Kec,Ku,如何设置sys=[]这里面呢?上面红色的地方的地方要如何修改?
有高手可以提示一下吗?谢谢!
回复
分享到:

使用道具 举报

发表于 2007-11-28 10:02 | 显示全部楼层
看上去好像没什么问题,并不需要什么特别的设置

function [sys,x0,str,ts] = para(t,x,u,flag)
global Ke Kec Ku;
Ke=18;Kec=1;Ku=0.5;
switch flag,
case 0,
    [sys,x0,str,ts]=mdlInitializeSizes;
case 3,
    sys=mdlOutputs(t,x,u);
case {1,2,4,9},
    sys=[];
  otherwise
    error(['Unhandled flag = ',num2str(flag)]);
end
function [sys,x0,str,ts]=mdlInitializeSizes
sizes = simsizes;
sizes.NumContStates  = 0;
sizes.NumDiscStates  = 0;
sizes.NumOutputs     = 3;
sizes.NumInputs      = 2;
sizes.DirFeedthrough = 0;
sizes.NumSampleTimes = 1;   % at least one sample time is needed
sys = simsizes(sizes);
x0= [];
str=[];
ts=[-1 0];  %这个要有的

function sys=mdlOutputs(t,x,u)
%u(1)和u(2)分别为偏差e和偏差变化率ec
%sys(1),sys(2),sys(3)分别为调整后的Ke,Kec,Ku
if abs(u(1))>0.5||abs(u(2))>0.3
    sys(1)=0.2*Ke;sys(2)=0.3*Kec;sys(3)=Ku+0.5;
elseif abs(u(1))<0.2||abs(u(2))<0.1
    sys(1)=0.5*Ke;sys(2)=0.5*Kec;sys(3)=Ku+0.1;
end
%sys=[];

也可以把Ke,Kec,Ku写在function里面
function [sys,x0,str,ts] = para(t,x,u,flag,Ke,Kec,Ku)
%global Ke Kec Ku;
%Ke=18;Kec=1;Ku=0.5;  %这三个变量是输入量
switch flag,
case 0,
    [sys,x0,str,ts]=mdlInitializeSizes;
case 3,
    sys=mdlOutputs(t,x,u,Ke,Kec,Ku);
case {1,2,4,9},
    sys=[];
  otherwise
    error(['Unhandled flag = ',num2str(flag)]);
end
function [sys,x0,str,ts]=mdlInitializeSizes
sizes = simsizes;
sizes.NumContStates  = 0;
sizes.NumDiscStates  = 0;
sizes.NumOutputs     = 3;
sizes.NumInputs      = 2;
sizes.DirFeedthrough = 0;
sizes.NumSampleTimes = 1;   % at least one sample time is needed
sys = simsizes(sizes);
x0= [];
str=[];
ts=[-1 0];
function sys=mdlOutputs(t,x,u,Ke,Kec,Ku)
%u(1)和u(2)分别为偏差e和偏差变化率ec
%sys(1),sys(2),sys(3)分别为调整后的Ke,Kec,Ku
if abs(u(1))>0.5||abs(u(2))>0.3
    sys(1)=0.2*Ke;sys(2)=0.3*Kec;sys(3)=Ku+0.5;
elseif abs(u(1))<0.2||abs(u(2))<0.1
    sys(1)=0.5*Ke;sys(2)=0.5*Kec;sys(3)=Ku+0.1;
end

试一下看看吧
 楼主| 发表于 2007-11-29 15:14 | 显示全部楼层
谢谢!我会试试
发表于 2007-12-1 11:43 | 显示全部楼层

值得看看

正在使用s函数,来学习一下
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-25 13:31 , Processed in 0.055387 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表