User:Mtxu

From Wikimization

(Difference between revisions)
Jump to: navigation, search
Line 1: Line 1:
Hello,
Hello,
 +
 +
[[Beginning with CVX|My question in Control theory: Click here]].
I'm Maria Calle, I have an answer abot CVX and how to work with. I hope here is the place for answer about.
I'm Maria Calle, I have an answer abot CVX and how to work with. I hope here is the place for answer about.
Line 8: Line 10:
So I'm begining whit CVX program and now I'm lost. Could you help me? I don't know how to continue.
So I'm begining whit CVX program and now I'm lost. Could you help me? I don't know how to continue.
-
 
-
I have a system, A, B and C. It have to move poles to one region, for it, I make a transform to the original system. (A_alfa, B_alfa, C_alfa).
 
-
I have to know the gain matrix, with is:
 
-
<math>K=R*W^(-1)</math>
+
[http://www.convexoptimization.com/wikimization/index.php?title=Special:Userlogin&returnto=Special:Userlogout Login] required to [http://www.convexoptimization.com/wikimization/index.php/Special:Emailuser/Mtxu E-mail this user].<br>
-
 
+
If you do not already have an account, [http://www.convexoptimization.com/wikimization/index.php?title=Special:Userlogin&returnto=Special:Userlogout create one here].
-
It's knew that:
+
-
 
+
-
<math>v(x)=x'*P*x=x'*W^(-1)*x</math>.
+
-
 
+
-
<math>x'*(A'*P+P*A)*x-2*x'P*B*K*x<0</math>.
+
-
 
+
-
<math>u(t)=-K*x(t)</math>.
+
-
 
+
-
 
+
-
<math>Hij(W,R)=W*A/alfa_i '+A/alfa_i*W-B/alfa_i*R-R'*B/alfa_j '<0</math>
+
-
 
+
-
and
+
-
 
+
-
<math>W=eye(2,2)*w, w \in{} R^(2x2)</math>
+
-
 
+
-
<math>R=eye(2,2)*w, w \in{} R^(4x2)</math>.
+
-
+
-
The ecuation must: <math>min f(W,R)=min(p1)+min(p2)</math>
+
-
 
+
-
subject to:
+
-
 
+
-
<math>w_nn \leq{} p1</math>
+
-
 
+
-
<math>-p2 \leq{} r_qn\leq{} p2</math>
+
-
+
-
with
+
-
 
+
-
<math>n=1,....4, q=1...2.</math>
+
-
 
+
-
<math>W=w*eye(4)</math>
+
-
 
+
-
<math>R=r*eye(2,4)</math>
+
-
+
-
and:
+
-
 
+
-
<math>W \geq{} \epsilon1*eye(4)</math>
+
-
 
+
-
<math>Hij(W,R) \leq{} -\epsilon2*eye(4)</math>
+
-
+
-
Now the article says that aplaing the convex programm, it posible to solve the ecuation.
+
-
+
-
They present the next algoritm:
+
-
1) Initialization
+
-
 
+
-
<math>l=0, W1=W1=eye(4), R1=R2=zeros(2,4)</math>
+
-
 
+
-
2)Calculate
+
-
 
+
-
<math>\lambda(W)=\lambda min(W)</math>
+
-
 
+
-
<math>\lambda_Hij=\lambda max(Hij(Wi,Ri))</math>
+
-
+
-
3)if
+
-
 
+
-
<math>\lambda(W) \geq \epsilon1</math>
+
-
 
+
-
<math>lambda_H(ij) \leq -\epsilon2</math>
+
-
 
+
-
STOP
+
-
+
-
Else
+
-
 
+
-
<math>l=l+1</math>
+
-
 
+
-
calculate eigenvalues(v_w y v_Hij).
+
-
 
+
-
calculate the constraint linear: C1(W,R)
+
-
 
+
-
<math>v_w'*W1*v_w \geq \epsilon1</math> if <math>(\epsilon1-\lambda_w)>(\lambda_Hij+\epsilon2)</math>
+
-
 
+
-
<math>v_Hij'*Hij*v_Hij \leq -\epsilon2</math> if <math>(\epsilon1-\lambda_w)\leq (\lambda_Hij+\epsilon2)</math>
+
-
+
-
4)Solve
+
-
 
+
-
<math>min(p1+p2)</math>
+
-
 
+
-
under
+
-
 
+
-
<math>\epsilon1 \leq w_nn \leq p1</math>
+
-
 
+
-
<math>-p2 \leq r_qn\leq p2</math>
+
-
 
+
-
<math>C_k(W,R)</math>
+
-
 
+
-
<math>k=1,...z</math>
+
-
 
+
-
Now I have no idea of how to continue, or how to program the algoritm.
+
-
+
-
They said the method for solving it's mathematical convex program.
+
-
+
-
Thanks a lot for all, Greetings.
+

Revision as of 13:51, 2 February 2009

Hello,

My question in Control theory: Click here.

I'm Maria Calle, I have an answer abot CVX and how to work with. I hope here is the place for answer about.

I'm one student from Spain. I'm studing engienering in automatic control. I'm in my last year and I'm making a proyect of investigacion. It's about LQR and how to sintonice it by pole assignment.

I have some problems with the initialitation of the gain matrix. I'm reading one article for know how to make it. In the article they said that it's necessary to solve by convex programm. Until last week I didn't knew anything about convex program.

So I'm begining whit CVX program and now I'm lost. Could you help me? I don't know how to continue.


Login required to E-mail this user.
If you do not already have an account, create one here.

Personal tools