矩阵积分

Question

我想找到一个矩阵的积分

T=[1 1 1]';
R=[1 1];
I=eye(3);
A=diag([1 1 1]');
B=[1 1 1]';
F=[1 1 1];
syms x
G=F*((exp(i*x)*I-A-B*F)\B);
f = T*G*(R*R')*G'*T';
int(f,[0 2.*pi])

我期望输出3 x 3的实数矩阵，但是我得到3 x 3的整数和未知x矩阵。

Answer 1

尝试以下操作：

T=[1 1 1]';
R=[1 1];
I=eye(3);
A=diag([1 1 1]');
B=[1 1 1]';
F=[1 1 1];
G=@(x) F*((exp(1i*x)*I-A-B*F)\B);
f = @(x) T*G(x)*(R*R')*G(x)'*T';
integral(f,0,2*pi,'ArrayValued',true)

输出

7.5398    7.5398    7.5398
7.5398    7.5398    7.5398
7.5398    7.5398    7.5398