在Matlab7.0里怎样用riccati方程方法求解H∞控制器
\u6025\u6c42matlab R2007b\u7248\u672c\u4e2driccati\u65b9\u7a0b\u7684\u6c42\u89e3\u547d\u4ee4\u662f\u4ec0\u4e48\uff1f\uff01\uff01\uff01\uff01\uff01\uff01care Solve continuous-time algebraic Riccati equation\u8fde\u7eed\u65f6\u95f4\u578b
dare Solve discrete-time algebraic Riccati equations (DAREs)\u79bb\u6563\u65f6\u95f4\u578b
\u540e\u9762\u4e24\u4e2a\u662f\u5e7f\u4e49\u89e3\uff1a
gcare Generalized solver for continuous-time algebraic Riccati equation
gdare Generalized solver for discrete-time algebraic Riccati equation
\u5177\u4f53\u8bf7\u67e5\u9605 help \u6587\u6863
\u662f\u7684\u3002
\u89e3: \u8bbey=y(x), \u5728R\u4e0a\u53ef\u5fae, \u6709\u91cc\u5361\u8482\u65b9\u7a0b(Riccati Equation):
y' = f1(x) + f2(x)y + f3(x)y^2;
\u53ef\u4ee4v=v(x) = y*f3(x)\u4ee3\u5165, \u5f97
v'(x) = v^2 + R(x)v + S(x)
\u5176\u4e2d\u6709R(x)=f2+[f3'(x)/f3(x)], S(x)=f1(x)*f3(x);
\u518d\u4ee4v(x)=−u'(x)/u(x), u=u(x)\u5728R\u4e0a\u53ef\u5fae\u4e14\u4e0d\u4e3a\u96f6, \u4ee3\u5165\u4e0a\u5f0f, \u5373\u53ef\u4ee5\u5316\u4e3a\u5173\u4e8eu\u7684\u4e8c\u9636\u7ebf\u6027\u5fae\u5206\u65b9\u7a0b:
u''−R(x)u'+S(x)u=0.
函数是system Function的简称，用它来写自己的simulink模块。（够简单吧，^_^， 详细的概念介绍大伙看帮助吧）可以用matlab、C、C++、Fortran、Ada等语言来写， 这儿我只介绍怎样用matlab语言来写吧（主要是它比较简单）
先讲讲为什么要用s函数，我觉得用s函数可以利用matlab的丰富资源，而不仅仅局 限于simulink提供的模块，而用c或c++等语言写的s函数还可以实现对硬件端口的操作， 还可以操作windows API等的
先介绍一下simulink的仿真过程（以便理解s函数），simulink的仿真有两个阶段： 一个为初始化，这个阶段主要是设置一些参数，像系统的输入输出个数、状态初值、 采样时间等；第二个阶段就是运行阶段，这个阶段里要进行计算输出、更新离散状态、 计算连续状态等等，这个阶段需要反复运行，直至结束。
在matlab的workspace里打edit sfuntmpl(这是matlab自己提供的s函数模板)，我们 看它来具体分析s函数的结构。 它的第一行是这样的：function [sys,x0,str,ts]=sfuntmpl(t,x,u,flag)
先讲输入与输出变量的含义：t是采样时间，x是状态变量，u是输入（是做成simulink模 块的输入）,flag是仿真过程中的状态标志（以它来判断当前是初始化还是运行等）；sys输出根据flag的不同而不同（下面将结合flag来讲sys的含义），x0是状态变量的初始值，str是保留参数（mathworks公司还没想好该怎么用它，嘻嘻，一般在初始化中将它置空就可以了,str=[])，ts是一个1×2的向量，ts(1)是采样周期，ts(2)是偏移量。
下面结合sfuntmpl.m中的代码来讲具体的结构：
switch flag, %判断flag，看当前处于哪个状态
case 0,
[sys,x0,str,ts]=mdlInitializeSizes;
flag=0表示处于初始化状态，此时用函数mdlInitializeSizes进行初始化，此函数在 sfuntmpl.m的149行
我们找到他，在初始化状态下，sys是一个结构体，用它来设置模块的一些参数，各个参 数详细说明如下
size = simsizes;%用于设置模块参数的结构体用simsizes来生成
sizes.NumContStates = 0;%模块连续状态变量的个数
sizes.NumDiscStates = 0;%模块离散状态变量的个数
sizes.NumOutputs = 0;%模块输出变量的个数
sizes.NumInputs = 0;%模块输入变量的个数
sizes.DirFeedthrough = 1;%模块是否存在直接贯通（直接贯通我的理解是输入能 %直接控制输出）
sizes.NumSampleTimes = 1;%模块的采样时间个数，至少是一个
sys = simsizes(sizes); %设置完后赋给sys输出
举个例子，考虑如下模型：
dx/dt=fc(t,x,u) 也可以用连续状态方程描述：dx/dt=A*x+B*u
x(k+1)=fd(t,x,u) 也可以用离散状态方程描述：x(k+1)=H*x(k)+G*u(k)
y=fo(t,x,u) 也可以用输出状态方程描述：y=C*x+D*u
设上述模型连续状态变量、离散状态变量、输入变量、输出变量均为1个，我们就只需改上面那一段代码为：
（一般连续状态与离散状态不会一块用，我这儿是为了方便说明）
sizes.NumContStates=1;sizes.NumDiscStates=1;sizes.NumOutputs=1;sizes.NumInpu
ts=1;
其他的可以不变。继续在mdlInitializeSizes函数中往下看：
x0 = []; %状态变量设置为空，表示没有状态变量，以我们上面的假设，可改 %为x0=[0,0](离散和连续的状态变量我们都设它初值为0)
str = []; %这个就不用说了，保留参数嘛，置[]就可以了，反正没什么用，可 %能7.0会给它一些意义
ts = [0 0]; %采样周期设为0表示是连续系统，如果是离散系统在下面的mdlGet %TimeOfNextVarHit函数中具体介绍
嘻嘻，总算讲完了初始化，后面的应该快了
在sfuntmpl的106行继续往下看：
case 1,
sys=mdlDerivatives(t,x,u);
flag=1表示此时要计算连续状态的微分，即上面提到的dx/dt=fc(t,x,u)中的dx/dt,找到 mdlDerivatives函数（在193行）如果设置连续状态变量个数为0，此处只需sys=[]; 就可以了（如sfuntmpl中一样），按我们上述讨论的那个模型，此处改成 sys=fc(t,x(1),u)或sys=A*x(1)+B*u %我们这儿x(1)是连续状态变量，而x(2)是离散的 ，这儿只用到连续的，此时的输出sys就是微分
继续，在sfuntmpl的112行：
case 2,
sys=mdlUpdate(t,x,u);
flag=2表示此时要计算下一个离散状态，即上面提到的x(k+1)=fd(t,x,u)，找到mdlUpd ate函数（在206行）它这儿sys=[];表示没有离散状态，我们这而可以改成 sys=fd(t,x(2),u)或sys=H*x(2)+G*u;%sys即为x(k+1)
看来后面几个一两句话就可了，呵呵，在sfuntmpl的118行
case 3,
sys=mdlOutputs(t,x,u);
flag=3表示此时要计算输出，即y=fo(t,x,u),找到mdlOutputs函数（在218行），如上，如果sys=[]表示没有输出，我们改成sys=fo(t,x,u)或sys=C*x+D*u %sys此时为输出y
好像快完了，嘻嘻，在sfuntmpl的124行
case 4,
sys=mdlGetTimeOfNextVarHit(t,x,u);
flag=4表示此时要计算下一次采样的时间，只在离散采样系统中有用(即上文的mdlInit ializeSizes中提到的ts设置ts(1)不为0)
连续系统中只需在mdlGetTimeOfNextVarHit函数中写上sys=[];这个函数主要用于变步长的设置，具体实现大家可以用edit vsfunc看vsfunc.m这个例子
最后一个，在sfuntmpl的130行
case 9,
sys=mdlTerminate(t,x,u);
flag=9表示此时系统要结束，一般来说写上在mdlTerminate函数中写上sys=[]就可，如 果你在结束时还要设置什么，就在此函数中写
关于sfuntmpl这个s函数的模板讲完了。
s函数还可以带用户参数，下面给个例子，和simulink下的gain模块功能一样，大伙自己 看吧，我睡觉去了，累了
function [sys,x0,str,ts] = sfungain(t,x,u,flag,gain)
switch flag,
case 0,
sizes = simsizes;
sizes.NumContStates = 0;
sizes.NumDiscStates = 0;
sizes.NumOutputs = 1;
sizes.NumInputs = 1;
sizes.DirFeedthrough = 1;
sizes.NumSampleTimes = 1;
sys = simsizes(sizes);
x0=[];
str=[];
ts=[0,0];
case 3,
sys=gain*u;
case {1,2,4,9},
sys = [];
end
做好了s函数后，simulink--user-defined function下拖一个S-Function到你的模型，就可以用了 在simulink----user-defined function还有个s-Function Builder，他可以生成用c语 言写的s函数 在matlab的workspace下打sfundemos，可以看到很多演示s函数的程序
绛旓細浣犺鐨r鎸囩殑鏄釜浠涔堜笢瑗匡紵simulink閲屾病鏈夌壒瀹氬搴斾簬r鐨勪俊鎭傚彟澶7.0鐗堟湰澶佷簡锛鐢鏂扮増鐨勫惂銆
绛旓細more 鍦ㄥ懡浠ょ獥鍙ｄ腑鎺у埗鍒嗛〉杈撳嚭 鍚姩鍜岄鍑篗ATLAB Quit 閫鍑篗ATLAB Startup 寮曠敤MATLAB鏃舵墍鎵ц鐨凪鏂囦欢 Matlabrc 涓诲惎鍔∕鏂囦欢 涓鑸俊鎭 Info MATLAB绯荤粺淇℃伅鍙奙athworks鍏徃淇℃伅 Subscribe 鎴愪负MATLAB鐨璁㈣喘鐢ㄦ埛 hostid MATLAB涓绘湇鍔＄▼搴忕殑璇嗗埆浠ｅ彿 Whatsnew 鍦ㄨ鏄庝功涓湭鍖呭惈鐨勬柊淇℃伅 Ver 鐗堟湰淇℃伅 鎿嶄綔绗﹀拰鐗规畩瀛楃 + ...
绛旓細鈥r鈥濅竴鑸槸绾㈣壊鐨勬剰鎬 涓嶅ソ鎰忔 娌＄湅鍒扳 ' 鈥濓紝杩欎釜浠ｈ〃瀵箁鐭╅樀鐨勮浆缃
绛旓細matlab7.0鏄疪14鐗堟湰銆俶atlab7.0鏄竴娆剧敤浜庣畻娉曞紑鍙戙佹暟鎹彲瑙嗗寲銆佹暟鎹垎鏋愪互鍙婃暟鍊艰绠楃殑楂樼骇鎶鏈绠楄瑷鍜屼氦浜掑紡鐜锛屼富瑕佸寘鎷琈ATLAB鍜孲imulink涓ゅぇ閮ㄥ垎銆
绛旓細鐩稿綋绠鍗曪紝鐢–trl+R杩涜鏃嬭浆锛
绛旓細7.0灞炰簬R2004 R2007鏄7.4鍜7.5 鐐瑰嚮about MATLAB閲岄潰鏈夈傚弬鑰冭祫鏂欙細http://baike.baidu.com/view/10598.htm?fr=ala0_1_1
绛旓細鍗曞嚮鍏冧欢锛屽湪鑿滃崟涓湁涓涓猣ormat 锛閲岄潰鏈変竴涓猺otate block锛屽彲浠ユ棆杞殑銆備篃鍙互ctrl+R瀹炵幇 绁濅綘瀛︿範鎰夊揩锛
绛旓細鎸変互涓6涓楠わ細1.棣栧厛鍦ㄦ闈笂榧犳爣鍙抽敭鐐瑰嚮matlab7.0蹇嵎鏂瑰紡 2.杩涘叆蹇嵎鏂瑰紡灞炴х晫闈紝鐐瑰嚮鍏煎鎬ч夐」鍗 3.鍕鹃夛細浠ュ吋瀹规ā寮忚繍琛岃繖涓▼搴 4.閫夋嫨windows2000杩欎釜閫夋嫨 5.鐐瑰嚮搴旂敤鎸夐挳锛屽畬鎴愯缃 6.閲嶆柊鍙屽嚮matlab7.0蹇嵎鏂瑰紡锛屽嵆鍙垚鍔熻繍琛屻
绛旓細quit 鎺ㄥ嚭Matlab 鐜 quiver 浜岀淮鏂瑰悜绠ご鍥 quiver3 涓夌淮鏂瑰悜绠ご鍥 R r rand 浜х敓鍧囧寑鍒嗗竷闅忔満鏁 randn 浜х敓姝ｆ佸垎甯冮殢鏈烘暟 randperm 闅忔満缃崲鍚戦噺 range 鏍锋湰鏋佸樊 rank 鐭╅樀鐨勭З rats 鏈夌悊杈撳嚭 rcond 鐭╅樀鍊掓潯浠舵暟浼拌 real 澶嶆暟鐨勫疄閮 reallog 鍦ㄥ疄鏁板煙鍐呰绠楄嚜鐒跺鏁 realpow 鍦ㄥ疄鏁板煙鍐呰绠椾箻鏂 realsqrt ...
绛旓細鍦╩atlab鐨剋orkspace閲屾墦edit sfuntmpl(杩欐槸matlab鑷繁鎻愪緵鐨剆鍑芥暟妯℃澘)锛屾垜浠 鐪嬪畠鏉ュ叿浣撳垎鏋恠鍑芥暟鐨勭粨鏋勩 瀹冪殑绗竴琛屾槸杩欐牱鐨勶細function [sys,x0,str,ts]=sfuntmpl(t,x,u,flag)鍏堣杈撳叆涓庤緭鍑哄彉閲忕殑鍚箟锛歵鏄噰鏍锋椂闂达紝x鏄姸鎬佸彉閲忥紝u鏄緭鍏ワ紙鏄仛鎴恠imulink妯 鍧楃殑杈撳叆锛,flag鏄豢鐪熻繃绋嬩腑鐨...
