题目

解方程 以矩阵元素作为变量
程序如下：
clear;clc;
a=sym('a',[3,1]);
a=sym(a,'real');
b=[1 2 3;4 5 6;7 8 9;];
c=a'*b*a;
d=sym('d',[3,1]);
d=[diff(c,a(1,1)); diff(c,a(2,1)); diff(c,a(3,1))];
s=solve(d(1,1),d(2,1),d(3,1),'a(1,1)','a(2,1)','a(3,1)');
e=[s.a(1,1); s.a(2,1); s.a(3,1)];
最后s.a(1,1)是输出不出来的

提问时间：2021-12-20

答案

clear;clc;for i=1:3a(i,1)=sym(['a' num2str(i)],'real');endb=[1 2 3;4 5 6;7 8 9;];c=a'*b*a;for i=1:3d(i,1)=sym(['d' num2str(i)]);endd=[diff(c,a(1,1)); diff(c,a(2,1)); diff(c,a(3,1))];s=solve(d(1,1),d(2...