matlab求解线性方程组请问这个怎么解啊?
matlab求解线性方程组请问这个怎么解啊?
matlab求解线性方程组
matlab求解线性方程组请问这个怎么解啊?
% 这显然不是线性方程组,上述解析解目测需要直接化解推导,比如方程（3）-（2）
% 方程（4）/ζ-（1）
% 用matlab求解如下
syms malpha mbeta mgamma mdelta mzeta;
rs=solve('1+mdelta*(4*malpha+4*mgamma-2)+2*(2*malpha-1)*mbeta*(1+4*mgamma*mdelta)',...
'1+2*malpha+(2+4*mbeta+8*malpha*mbeta)*mgamma',...
'1+3*malpha+(3+8*mbeta+20*malpha*mbeta)*mgamma',...
'(1+mdelta*(4*malpha+4*mgamma+2)+2*(2*malpha+1)*mbeta*(1+4*mgamma*mdelta))*mzeta-2',...
'1+(4*mbeta-2)*mgamma-2*malpha*(1+4*mgamma*mdelta)-mzeta',...
'malpha','mbeta','mgamma','mdelta','mzeta');
% 五个方程,五个未知量,貌似是定解问题,求解结果共有四组：
for n=1:4
disp(['第' num2str(n) ']);
disp(['α=' num2str(eval(rs.malpha(n)))]);
disp(['β=' num2str(eval(rs.mbeta(n))) ]);
disp(['γ=' num2str(eval(rs.mgamma(n))) ]);
disp(['δ=' num2str(eval(rs.mdelta(n))) ]);
disp(['ζ=' num2str(eval(rs.mzeta(n))) ]);
end
%下面是显示结果
α=0
β=-0.25
γ=-1
δ=0.1875
ζ=4
α=-0.37903
β=-4.2708
γ=0.11342
δ=-4.2708
ζ=-1.8751
α=-0.31049-0.69335i
β=0.10415-0.061534i
γ=-0.43171+0.60389i
δ=0.10415-0.061534i
ζ=1.9376+0.71569i
α=-0.31049+0.69335i
β=0.10415+0.061534i
γ=-0.43171-0.60389i
δ=0.10415+0.061534i
ζ=1.9376-0.71569i