用MATLAB求Lyapunove方程,matlab方程的求解.doc

matlab方程的求解.doc

MATLAB的LYAPUNOV、SYLVESTER和RICCATI方程的MATLAB求解一、连续LYAPUNOV方程连续LYAPUNOV方程可以表示为LYAPUNOV方程来源与微分方程稳定性理论，其中要求C为对称正定的NN方阵，从而可以证明解X亦为NN对称矩阵，这类方程直接求解比较困难，不过有了MATLAB中LYAP函数，就简单多了。A123456780A123456780C1054567479C1054567479XLYAPA,CX394443888903889388892777802222038890222201111二、LYAPUNOV方程的解析解利用KRONCECKER乘积的表示方法，可以将LYAPUNOV方程写为FUNCTIONXLYAP2A,CLYAPUNOV方程的符号解法NSIZEC,1A0KRONA,EYENKRONEYEN,ACCX0INVA0CXRESHAPEX0,N,N例子A123456780C1054567479XLYAP2SYMA,SYMCX71/18,35/9,7/1835/9,25/9,2/97/18,2/9,1/9三、离散LYAPUNOV方程离散LYAPUNOV方程的一般形式为MATLAB中直接提供了DLYAP函数求解该方程XDLYAPA,Q其实，如果A矩阵非奇异，则等式两边同时右乘得到就可以将其变换成连续的SYLVESTER方程而SYLVESTER方程是广义LYAPUNOV方程，故离散的LYAPUNOV方程还可以使用下面的方法求解BINVA’CQINVA’XLYAPA,B,C下面总结下我们上面的讲到的知识点XLYAPA,C连续LYAPUNOV方程数值解法XLYAP2A,C连续LYAPUNOV方程符号解法XLYAPA,B,C广义LYAPUNOV方程，即SYLVESTER方程XDLYAPA,Q或者XLYAPA,INVA’,QINVA’离散LYAPUNOV方程SYLVESTER方程MATLAB求解SYLVESTER方程的一般形式为该方程又称为广义的LYAPUNOV方程，式中A是NN方阵，B是MM方阵，X和C是NM矩阵。MATLAB控制工具箱提供了直接的求解该方程的LYAP函数A816357492B2345C123456XLYAPA,B,CA816357492B2345C123456X020110201600393015540642808966同理，我们使用KRONECKER乘机的形式将原方程进行如下变化故可以编写SYLVESTER方程的解析解函数FUNCTIONXLYAP3A,B,CSYLVESTER方程的解析解法REWRITEDBYDYNAMICMOREINATIONHTTP//WWWILOVEMATLABCNIFNARGIN2,CBBA ENDNR,NCSIZECA0KRONA,EYENCKRONEYENR,B TRYC1C X0INVA0C1XRESHAPEX0,NC,NRCATCHERROR MATLABSKY提醒您矩阵奇异 END用上面的数据，我们试验下该解析解法的XLYAP3SYMA,B,CX2853/14186,557/14186,9119/1418611441/56744,8817/56744,50879/56744RICCATI方程的MATLAB求解RICCATI方程是一类很著名的二次型矩阵形式，其一般形式为由于含有矩阵X的二次项，所有RICCATI方程求解要LYAPUNOV方程更难，MATLAB控制工具箱提供了ARE函数，可以直接求解该函数A213102012B222152112C544104115XAREA,B,CA213102012B222152112C544104115X098740798304189057740130805775028400073006924如何用MATLAB求解LYAPUNOV指数我是需要分析计算LOGISTIC数据,都是用来说明对初值的敏感以下是LOGISTIC求解的程序,希望得到LYAPUNOV的程序CLCCLEARCLOSEALLLAMBDA35E44X04ONES1,LENGTHLAMBDAN1400前面的迭代点数N2100后面的迭代点数FZEROSN1N2,LENGTHLAMBDAFORI1N1N2XLAMBDAX1XFI,XENDFFN11END,PLOTLAMBDA,F, K , MARKERSIZE ,1XLABEL \LAMBDA YLABEL X 慢慢看吧，很有用