matlab求解线性方程组请问这个怎么解啊?
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/27 20:25:28
![matlab求解线性方程组请问这个怎么解啊?](/uploads/image/z/7622474-50-4.jpg?t=matlab%E6%B1%82%E8%A7%A3%E7%BA%BF%E6%80%A7%E6%96%B9%E7%A8%8B%E7%BB%84%E8%AF%B7%E9%97%AE%E8%BF%99%E4%B8%AA%E6%80%8E%E4%B9%88%E8%A7%A3%E5%95%8A%3F)
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