MATLAB的答案

解决一个多元方程(方程,4变量)n行数据使用循环

37视图(30天)
maryam西文
maryam西文 2023年5月15日12:22
评论道: maryam西文2023年5月16日在7:08
我写的这段代码解决一个多元方程(4方程,4变量)虽然aprim0, bprim0, xmid, ymid是可用的
为i = 1: n
eqn1 =是的(我,1)- (aprim0(我,1)* xup(我,1))——bprim0(我,1)= = 0;
eqn2 = ydown (1) - (aprim0(我,1)* xdown(我,1))——bprim0(我,1)= = 0;
eqn3 =√(xup(我,1)- xmid(我,1))^ 2 +(是的(我,1)- ymid(我,1))^ 2)- 0.008 = = 0;
eqn4 =√(xdown(我,1)- xmid(我,1))^ 2 + (ydown(我,1)- ymid(我,1))^ 2)- 0.008 = = 0;
命令= [eqn1 eqn2、eqn3 eqn4);
var = [xup(1),我是的(我,1),xdown(我,1),ydown(我,1)];
索尔=解决(“命令”,“var”);
结束
在解决函数,但它没有工作我的错误在哪里?

在回答这个问题。

接受的答案

Torsten
Torsten 2023年5月15日在51
编辑:Torsten 2023年5月15日在7:30
跑:
信谊Aprim0 Bprim0 Xmid Ymid R真实
信谊Xup Xdown是的Ydown
%定义方程组
eqn1 =是的- Aprim0 * Xup Bprim0 = = 0;
eqn2 = Ydown - Aprim0 * Xdown Bprim0 = = 0;
eqn3 = (Xup - Xmid) ^ 2 +(是的- Ymid) ^ 2 - R ^ 2 = = 0;
eqn4 = (Xdown - Xmid) ^ 2 + (Ydown - Ymid) ^ 2 - R ^ 2 = = 0;
%解方程组
索尔=解决([eqn1, eqn2],[是的,Ydown]);
eqn3 =潜艇(eqn3,是的,sol.Yup);
sol_Xup =简化(解决(eqn3 Xup));
eqn4 =潜艇(eqn4 Ydown sol.Ydown);
sol_Xdown =简化(解决(eqn4 Xdown));
sol_Yup =简化(潜艇(sol.Yup, Xup sol_Xup));
sol_Ydown =简化(潜艇(sol.Ydown, Xdown sol_Xdown));
%的插入数值参数来得到数值
%的未知数
n = 1;
aprim0 (1) = 5.0;
bprim0 (1) = 1.5;
xmid (1) = 6;
ymid (1) = -2.5;
r = 0.008;
为i = 1: n
xup(我)=双(潜艇(sol_Xup, [Aprim0 Bprim0 Xmid Ymid R], [Aprim0 (i) Bprim0 (i) Xmid (i) Ymid(我)R]));
xdown(我)=双(潜艇(sol_Xdown, [Aprim0 Bprim0 Xmid Ymid R], [Aprim0 (i) Bprim0 (i) Xmid (i) Ymid(我)R]));
是的(我)=双(潜艇(sol_Yup, [Aprim0 Bprim0 Xmid Ymid R], [Aprim0 (i) Bprim0 (i) Xmid (i) Ymid(我)R]));
ydown(我)=双(潜艇(sol_Ydown, [Aprim0 Bprim0 Xmid Ymid R], [Aprim0 (i) Bprim0 (i) Xmid (i) Ymid(我)R]));
结束
xup
xup =
-0.5385 - 1.3077我-0.5385 + 1.3077
xdown
xdown =
-0.5385 - 1.3077我-0.5385 + 1.3077
是的
是的=
-1.1923 - 6.5385我-1.1923 + 6.5385
ydown
ydown =
-1.1923 - 6.5385我-1.1923 + 6.5385
4评论
显示3年长的评论

登录置评。

更多的答案(0)

类别

找到更多的在矩阵索引在帮助中心和文件交换

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!