/
W_cvx_22.m
57 lines (31 loc) · 1.42 KB
/
W_cvx_22.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
function [W_r W_f]=W_cvx_22(Mn,Mr,N0,v_1,v_2 ,R_1,R_2 , H_1_1,H_1_2,H_2_1,H_2_2 ,H_1_r,H_2_r,H_r_1,H_r_2 ,H_1_f,H_2_f,H_f_1,H_f_2,D1,D2 )
cvx_begin
variable W_r(Mn,Mn)
variable W_f(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_2_r_1=H_2_r*W_r*H_r_1;
H_2_r_2=H_2_r*W_r*H_r_2;
H_1_f_1=H_1_f*W_f*H_f_1;
H_1_f_2=H_1_f*W_f*H_f_2;
H_2_f_1=H_2_f*W_f*H_f_1;
H_2_f_2=H_2_f*W_f*H_f_2;
H_p_1_1=H_1_1+H_1_r_1+H_1_f_1 ;
H_p_1_2=H_1_2+H_1_r_2+H_1_f_2 ;
H_p_2_1=H_2_1+H_2_r_1+H_2_f_1 ;
H_p_2_2=H_2_2+H_2_r_2+H_2_f_2 ;
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_1*[R_1']* N0)+trace(R_2*R_2'*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_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(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_f'*H_1_f'*R_1')'*vec(W_f'*H_1_f'*R_1')*N0+...
vec(W_f'*H_2_f'*R_2')'*vec(W_f'*H_2_f'*R_2')*N0);
subject to
vec(W_r)'*vec(W_r) <= 1;
vec(W_f)'*vec(W_f) <= 1;
cvx_end