|
|
A1=rand(2,2)+i*rand(2,2)
A2=rand(2,2)+i*rand(2,2)
A1=A1*A1'
A2=A2*A2'
B1=rand(2,2)+i*rand(2,2)
B2=rand(2,2)+i*rand(2,2)
B1=B1*B1'
B2=B2*B2'
t=0.3
a1=t*A1+(1-t)*B1
a2=t*A2+(1-t)*B2
ans1=log2(det(expm(a1)+expm(a2)))
ans2=t*log2(det(expm(A1)+expm(A2)))+(1-t)*log2(det(expm(B1)+expm(B2)))
这是仿真代码 |
|