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

%下面是显示结果
第1组解：
α=0
β=-0.25
γ=-1
δ=0.1875
ζ=4
第2组解：
α=-0.37903
β=-4.2708
γ=0.11342
δ=-4.2708
ζ=-1.8751
第3组解：
α=-0.31049-0.69335i
β=0.10415-0.061534i
γ=-0.43171+0.60389i
δ=0.10415-0.061534i
ζ=1.9376+0.71569i
第4组解：
α=-0.31049+0.69335i
β=0.10415+0.061534i
γ=-0.43171-0.60389i
δ=0.10415+0.061534i
ζ=1.9376-0.71569i