/
W_cvx_3.m
64 lines (40 loc) · 2.23 KB
/
W_cvx_3.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
function W_r=W_cvx_3(Mn,Mr,N0,v_1,v_2,v_3,R_1,R_2,R_3, H_1_1,H_1_2,H_2_1,H_2_2,H_1_3,H_2_3,H_3_3,H_r_3,H_3_1,H_3_2, H_3_r,H_1_r,H_2_r,H_r_1,H_r_2,D1,D2,D3)
cvx_begin
variable W_r(Mn,Mn)
H_1_r_1=H_1_r*W_r*H_r_1;
H_1_r_2=H_1_r*W_r*H_r_2;
H_1_r_3=H_1_r*W_r*H_r_3;
H_2_r_1=H_2_r*W_r*H_r_1;
H_2_r_2=H_2_r*W_r*H_r_2;
H_2_r_3=H_2_r*W_r*H_r_3;
H_3_r_1=H_3_r*W_r*H_r_1;
H_3_r_2=H_3_r*W_r*H_r_2;
H_3_r_3=H_3_r*W_r*H_r_3;
H_p_1_1=H_1_1+H_1_r_1 ;
H_p_1_2=H_1_2+H_1_r_2 ;
H_p_1_3=H_1_3+H_1_r_3 ;
H_p_2_1=H_2_1+H_2_r_1 ;
H_p_2_2=H_2_2+H_2_r_2 ;
H_p_2_3=H_2_3+H_2_r_3 ;
H_p_3_1=H_3_1+H_3_r_1 ;
H_p_3_2=H_3_2+H_3_r_2 ;
H_p_3_3=H_3_3+H_3_r_3 ;
minimize( trace( -R_1*H_p_1_1*v_1-v_1'*H_p_1_1'*R_1'+eye(D1,D1))+ ...
+trace( -R_2*H_p_2_2*v_2-v_2'*H_p_2_2'*R_2'+eye(D2,D2))+...
+trace( -R_3*H_p_3_3*v_3-v_3'*H_p_3_3'*R_3'+eye(D3,D3))+...
+trace(R_1*[R_1']* N0)+trace(R_2*R_2'*N0)+trace(R_3*R_3'*N0)+...
vec(v_1'*H_p_1_1'*R_1')'*vec(v_1'*H_p_1_1'*R_1')+...
vec(v_2'*H_p_1_2'*R_1')'*vec(v_2'*H_p_1_2'*R_1')+...
vec(v_3'*H_p_1_3'*R_1')'*vec(v_3'*H_p_1_3'*R_1')+...
vec(v_1'*H_p_2_1'*R_2')'*vec(v_1'*H_p_2_1'*R_2')+...
vec(v_2'*H_p_2_2'*R_2')'*vec(v_2'*H_p_2_2'*R_2')+...
vec(v_3'*H_p_2_3'*R_2')'*vec(v_3'*H_p_2_3'*R_2')+...
vec(v_1'*H_p_3_1'*R_3')'*vec(v_1'*H_p_3_1'*R_3')+...
vec(v_2'*H_p_3_2'*R_3')'*vec(v_2'*H_p_3_2'*R_3')+...
vec(v_3'*H_p_3_3'*R_3')'*vec(v_3'*H_p_3_3'*R_3')+...
vec(W_r'*H_1_r'*R_1')'*vec(W_r'*H_1_r'*R_1')*N0+...
vec(W_r'*H_2_r'*R_2')'*vec(W_r'*H_2_r'*R_2')*N0+...
vec(W_r'*H_3_r'*R_3')'*vec(W_r'*H_3_r'*R_3')*N0);
subject to
vec(W_r)'*vec(W_r) <= 1;
cvx_end