Solving a set of simultaneous equations

Question

MarshallSc 2022-4-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1691500-solving-a-set-of-simultaneous-equations

评论： Alex Sha 2022-4-9

I have the following set of simultaneous equations that I want to solve. I'm just confused as to what solver is the best choice!

Can it be even solved?

a = 0.5;
sqrt(x1^2 + a*x1*x2 + a*x1*x3 + a*x1*x4) = 30;
sqrt(a*x2*x1 + x2^2 + a*x2*x3 + a*x2*x4) = 33;
sqrt(a*x3*x1 + a*x3*x2 + x3^2 + a*x3*x4) = 25;
sqrt(a*x4*x1 + a*x4*x2 + a*x4*x3 + x4^2) = 10;

Thank you!

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Alex Sha 2022-4-9

在 MATLAB Online 中打开

There will be multi-solution if no range limition:

No.	x1	x2	x3	x4
-27.5584616528484	-31.4642304277135	-21.2013312551081	42.4667927458073
27.5584616528484	31.4642304277135	21.2013312551081	-42.4667927458073
21.092393187031	24.5327728282962	15.6460319211707	2.9752314576261
-35.124387658461	55.4211403974195	-28.2016268231437	-8.21718324886969
35.124387658461	-55.4211403974195	28.2016268231437	8.21718324886969
46.8115886629648	-42.6760633957222	-31.4217676113028	18.926670615613
-39.7163903825374	49.5555904803269	-32.6637611270595	17.2196122200305
39.7163903825374	-49.5555904803269	32.6637611270595	-17.2196122200305
-46.8115886629648	42.6760633957222	31.4217676113028	-18.926670615613
31.8811243735484	35.9706867855451	-49.7196043306907	6.44640901641614
-31.8811243735484	-35.9706867855451	49.7196043306907	-6.44640901641614
36.4849054651774	40.6839971426023	-42.3597267431898	-21.95861020684
-36.4849054651774	-40.6839971426023	42.3597267431898	21.95861020684
53.1637277627363	-38.0038803965751	-26.9964048344445	-7.46949936450619
-53.1637277627363	38.0038803965751	26.9964048344445	7.46949936450619

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Torsten 2022-4-8

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1691500-solving-a-set-of-simultaneous-equations#answer_937595

在 MATLAB Online 中打开

a = 0.5;
fun =@(x1,x2,x3,x4) [x1^2 + a*x1*x2 + a*x1*x3 + a*x1*x4 - 30^2;...
                     a*x2*x1 + x2^2 + a*x2*x3 + a*x2*x4 - 33^2;...
                     a*x3*x1 + a*x3*x2 + x3^2 + a*x3*x4 - 25^2;...
                     a*x4*x1 + a*x4*x2 + a*x4*x3 + x4^2 - 10^2];
x = fsolve(@(x)fun(x(1),x(2),x(3),x(4)),[1 1 1 1])
fun(x(1),x(2),x(3),x(4))