我有一个复杂的矩阵A(NxN)。在Matlab中,eig(A)将给出矩阵的所有复杂特征值。现在我很有兴趣找到每个复特征值的绝对值(r)和参数(\ phi)(每个特征值都有自己的r = abs(Z)和\ phi = arg(Z))。 如何编写以下产品表达式:
\ prod_j(sin(\ phi_j)+ r ^(1/2)_j其中索引j遍历矩阵A的所有特征值。
答案 0 :(得分:0)
要获得r
和phi
,只需使用Matlab函数abs
和angle
,就像这样......
z = eig(rand(5));
r = abs(z)
phi = angle(z)
然后你可以做任何你需要做的事情来得到矢量。
例如,所有sin(phi) + sqrt(r)
和phi
对的数量r
的乘积为:
prod( sin(phi) + sqrt(r) )
(注意,sin
和sqrt
函数的矢量化消除了对任何循环的需要。)