我怎样才能解决这个?

1视图(30天)
CHENG WEI李
CHENG WEI李 2022年10月23日
回答: Torsten 2022年10月24日
clc
清晰的
信谊PL c k
eqn = (PL * exp (- c * exp (- k * 0)) = = 179323;PL * exp (- c * exp (- k * 10)) = = 203302;PL * exp (- c * exp (- k * 20)) = = 226542]
S =解决(eqn (PL; c k))
S =
结构体字段:
PL(0×1符号):
c: [0×1 sym]
凯西:[0×1 sym]
2的评论
沃尔特·罗伯森
沃尔特·罗伯森 2022年10月23日
在实数方程没有解决方案。

登录置评。

答案(2)

穆罕默德·乌斯曼
穆罕默德·乌斯曼 2022年10月23日
%解决方程组起始点(0,0,0)。
% PL = x (1);c = x (2);k = x (3);
%初始猜测是(0,0,0),你可以改变它accordingily
有趣= @root2d;
x0 = (0, 0, 0);
x0, x = fsolve(有趣)
警告:Trust-region-dogleg FSOLVE算法不能处理方阵系统;使用Levenberg-Marquardt算法。
没有找到解决方案。fsolve停止因为当前步骤的相对大小小于步长值的宽容,但函数值的向量不是接近零的值函数的宽容。
x = 1×3
1.0 e + 5 * 2.0192 - -0.0000 0
函数F = root2d (x)
F (1) = x (1) * exp (- x (2) * exp (- x (3) * 0)) - 179323;
F (2) = x (1) * exp (- x (2) * exp (- x (3) * 10)) - 203302;
F (2) = x (1) * exp (- x (2) * exp (- x (3) * 20)) - 226542;
结束
2的评论
亚历克斯·沙
亚历克斯·沙 2022年10月24日
实际上有数值解像下图:
x1: 446505.431672107
x2: 0.912262916225993
x3: 0.0148006249649759

登录置评。


Torsten
Torsten 2022年10月24日
信谊PL c k
eqn1 = PL * exp (- c * exp (- k * 0)) = = 179323;
eqn2 = PL * exp (- c * exp (- k * 10)) = = 203302;
eqn3 = PL * exp (- c * exp (- k * 20)) = = 226542;
sPL =解决([eqn1 eqn2], [c k]);
警告:只有在特定条件下有效的解决方案。在解决方案,包括参数和条件指定“ReturnConditions”作为“真实”的价值。
PLsol =解决(潜艇(eqn3 [c k], [sPL。c sPL.k])、PL);
sck =解决([潜艇(PL, eqn1 PLsol),潜艇(eqn2, PL PLsol)], [c k]);
csol = sck.c;
ksol = sck.k;
vpa (PLsol)
ans =
446505.43167210849515476976968181
vpa (csol)
ans =
0.91226291622599614437018856917074
vpa (ksol)
ans =
0.014800624964975891465693800088887

类别

找到更多的在非线性方程组帮助中心文件交换

标签

世界杯预选赛小组名单社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!世界杯预选赛小组名单

开始狩猎!

Baidu
map