1. 求助 matlab中LMI工具箱gevp的详细用法,最好举个程序例子,万分感激
Solving gevp Problem in Matlab
30/10/10 13:51 Filed in: Matlab Code
http://arashfotouhi.com/files/-0.html
In order to solve the following optimization problem:
minimize a
subject to
[A’+PA+C-aP,(BP+CQ)’;
BP+CQ,-I] < 0
&
P=P’>0
in Matlab use the code below:
%Solving a generalized eigenvalue problem in Matlab
%Date: October 30, 2010
%Author: Arash Fotouhi
setlmis([]);
p=lmivar(1,[n 1]);
q=lmivar(2,[m m]);
y=lmivar(1,[n 1]);
lmiterm([1 1 1 0],0);
lmiterm([-1 1 1 p],1,1);
lmiterm([2 1 1 p],1,A,'s');
lmiterm([2 1 1 0],C);
lmiterm([2 2 1 p],B,1);
lmiterm([2 2 1 q],1,C);
lmiterm([2 2 2 0],-eye(m));
lmiterm([-2 1 1 y],1,1);
lmiterm([-2 2 1 0],0);
lmiterm([-2 2 2 0],0);
lmiterm([3 1 1 y],1,1);
lmiterm([-3 1 1 p],1,1);
lmis=getlmis;
[lambda,decvars]=gevp(lmis,1);
Popt = dec2mat(lmis,decvars,p);
Qopt = dec2mat(lmis,decvars,q);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2. 如何使用matlab中的工具箱
1、首先给出对应的拟合数据:>> x=1:100;>> y=2*x;一条直线。
3. matlab中LMI工具箱中lmiterm中怎么表示(PBK-Y)的转置,P、B、K、Y都是适应维数的矩阵。
1.因为令F=K*X1后,若所得不等式中中只有X1,F,则该不等式为LMI;若不等式中出现X1,F,K三个变量,则不等式变为受等式约束的LMI形式,无法直接通过Matlab求解。
2.不等式中出现X1^-1时,可以采取如下处理方法:
a)若X1^-1为矩阵对角元素,则可根据不等式-X1^-1<=X1-2I将原矩阵不等式替换成LMI形式;(此方法简便易行,但保守性增大)
b)利用锥补线性化算法(cone complementary linearisation),可参考以下文献
郭亚锋, 李少远. 网络控制系统的H∞状态反馈控制器设计. 控制理论与应用
G. Zhou, D. Wang, P. Chen, X. Zhao. Guaranteed Cost Control For Networked Control System With Interval Time-varying Delay. //Proceedings of the 21st Chinese Control and Decision Conference
4. 求高手用matlab的LMI工具箱 或者 YALMIP工具包 求下列矩阵不等式中η的最优解(最大值)。
和可惜不会啊
5. matlab的lmi工具箱求解线性矩阵不等式。
第二个lmiterm和第三个多打了一个括号。在Y]),这个)不要。
6. A和B等于什么,用matlab中的lmi工具箱求解
是优化工具性(Optimization Toolbox)的函数,但从2014b这个函数被删除了,官方建议用intlinprog代替。intlinprog是2014a引入的函数,可用于求解混合整数线性规划(MILP)问题。
7. MATLABLMI工具箱使用方法
所谓“工具箱”只是说提供了一组函数,这些函数跟matlab自带的其他函数本质上没有任回何区别,因此答编程方法上也没有什么特别的技巧。
先写m文件,在里面的代码里调用LMI工具箱的函数,然后运行就可以了。附件里是一个示例m文件,可以直接在matlab里运行了。matlab2013b运行通过。
8. matlab中的LMI工具箱
拜托,是Lmiterm,不是Imiterm!
9. matlab中的LMI工具箱如何求解带有克罗内克积的线性矩阵不等式
设函数 φ (x)连续且满足 φ (x)=e^x+ ∫(x,0)(t-x) φ(t)dt,求φ(x) 解: φ (x)=e^x+ ∫[0→x] (t-x) φ(t)dt =e^x+ ∫[0→x] tφ(t)dt-x∫[0→x] φ(t)dt 两边对x求导得: φ'(x)=e^x+ xφ(x)-∫[0→x] φ(t)dt-xφ(x) =e^x-∫[0→x] φ(t)dt (1) 两边再对导: φ''(x)...
10. 如何用matlab的LMI工具箱或yalmip工具箱求解下面矩阵不等式:
P = double (P)即可求出