三元方程组牛顿 java_萌新请教牛顿法求解三元非线性方程组

最新推荐文章于 2021-02-28 01:58:39 发布

惢安

最新推荐文章于 2021-02-28 01:58:39 发布

阅读量423

点赞数

文章标签：三元方程组牛顿 java

本文链接：https://blog.csdn.net/weixin_35625530/article/details/114574436

版权

该楼层疑似违规已被系统折叠隐藏此楼查看此楼

请问牛顿迭代法求解三元非线性方程组，不同迭代初值貌似有很多不同结果，如何求解到满足要求的解，0

FindRoot[{-6.565159793181527`*^-6 (0.008850684113090817` -

0.7905181223796485` (0.008850684113090817` - t0)) z1 -

1.3980263449472176`*^-10 (1.438702749038092`*^8 +

1.2840021082598146`*^7 t0) z1 +

0.00011254833333333333` z1^3 - (

0.0005367544235159816` (E^(-0.4` z1) + E^(

0.4` z1)) (-3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 - 3.696` E^(0.4` z2) z2))/(

E^(0.8` z1) - E^(0.8` z2)) + (

0.0005367544235159816` (-E^(-0.4` z1) + E^(

0.4` z1)) (3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 + 3.696` E^(0.4` z2) z2))/(

E^(0.8` z1) - E^(0.8` z2)) - t0 == 0,

1.3980263449472176`*^-10 (1.438702749038092`*^8 +

1.2840021082598146`*^7 t0) +

0.000338625` z1^2 - (-6.565159793181527`*^-6 \

(0.008850684113090817` -

0.7905181223796485` (0.008850684113090817` - t0)) +

9.8`*^-7 z1^2 - (

0.00021470176940639266` (-E^(-0.4` z1) + E^(

0.4` z1)) (-3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 - 3.696` E^(0.4` z2) z2))/(

E^(0.8` z1) - E^(0.8` z2)) + (

0.00021470176940639266` (E^(-0.4` z1) + E^(

0.4` z1)) (3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 + 3.696` E^(0.4` z2) z2))/(

E^(0.8` z1) - E^(0.8` z2))) == 0,

0.00031751999999999996` -

6.565159793181527`*^-6 (0.008850684113090817` -

0.7905181223796485` (0.008850684113090817` - t0)) +

9.8`*^-7 z1^2 - 9.8`*^-7 z2^2 - (

6.606735159817352`*^-7 (-E^(-0.4` z1) + E^(

0.4` z1)) (-3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 - 3.696` E^(0.4` z2) z2))/(

E^(0.8` z1) - E^(0.8` z2)) + (

6.606735159817352`*^-7 (-E^(-0.4` z2) + E^(

0.4` z2)) (-3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 - 3.696` E^(0.4` z2) z2))/(

E^(0.8` z1) - E^(0.8` z2)) + (

6.606735159817352`*^-7 (E^(-0.4` z1) + E^(

0.4` z1)) (3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 + 3.696` E^(0.4` z2) z2))/(

E^(0.8` z1) - E^(0.8` z2)) - (

6.606735159817352`*^-7 (E^(-0.4` z2) + E^(

0.4` z2)) (3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 + 3.696` E^(0.4` z2) z2))/(

E^(0.8` z1) - E^(

0.8` z2)) - (-0.000031536` (-0.65638188` +

0.00011320166666666666` z2^3 - (

0.0005367544235159816` (E^(-0.4` z2) + E^(

0.4` z2)) (-3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 - 3.696` E^(0.4` z2) z2))/(

E^(0.8` z1) - E^(0.8` z2)) - (

0.0005334510559360729` (-E^(-0.4` z2) + E^(

0.4` z2)) (3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 + 3.696` E^(0.4` z2) z2))/(

E^(0.8` z1) - E^(0.8` z2)) +

18 (-6.565159793181527`*^-6 (0.008850684113090817` -

0.7905181223796485` (0.008850684113090817` - t0)) +

9.8`*^-7 z1^2 - 9.8`*^-7 z2^2 - (

6.606735159817352`*^-7 (-E^(-0.4` z1) + E^(

0.4` z1)) (-3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 -

3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (

6.606735159817352`*^-7 (-E^(-0.4` z2) + E^(

0.4` z2)) (-3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 -

3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (

6.606735159817352`*^-7 (E^(-0.4` z1) + E^(

0.4` z1)) (3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 +

3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) - (

6.606735159817352`*^-7 (E^(-0.4` z2) + E^(

0.4` z2)) (3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 +

3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2))) +

z2 (-6.565159793181527`*^-6 (0.008850684113090817` -

0.7905181223796485` (0.008850684113090817` - t0)) +

9.8`*^-7 z1^2 - 9.8`*^-7 z2^2 - (

6.606735159817352`*^-7 (-E^(-0.4` z1) + E^(

0.4` z1)) (-3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 -

3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (

6.606735159817352`*^-7 (-E^(-0.4` z2) + E^(

0.4` z2)) (-3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 -

3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (

6.606735159817352`*^-7 (E^(-0.4` z1) + E^(

0.4` z1)) (3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 +

3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) - (

6.606735159817352`*^-7 (E^(-0.4` z2) + E^(

0.4` z2)) (3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 +

3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2))) +

18 (-6.565159793181527`*^-6 (0.008850684113090817` -

0.7905181223796485` (0.008850684113090817` - t0)) +

9.8`*^-7 z1^2 + 0.000338625` z2^2 - (

6.606735159817352`*^-7 (-E^(-0.4` z1) + E^(

0.4` z1)) (-3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 -

3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) - (

0.00021404109589041095` (-E^(-0.4` z2) + E^(

0.4` z2)) (-3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 -

3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (

6.606735159817352`*^-7 (E^(-0.4` z1) + E^(

0.4` z1)) (3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 +

3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (

0.00021404109589041095` (E^(-0.4` z2) + E^(

0.4` z2)) (3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 +

3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2))) -

z2 (-6.565159793181527`*^-6 (0.008850684113090817` -

0.7905181223796485` (0.008850684113090817` - t0)) +

9.8`*^-7 z1^2 + 0.000338625` z2^2 - (

6.606735159817352`*^-7 (-E^(-0.4` z1) + E^(

0.4` z1)) (-3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 -

3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) - (

0.00021404109589041095` (-E^(-0.4` z2) + E^(

0.4` z2)) (-3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 -

3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (

6.606735159817352`*^-7 (E^(-0.4` z1) + E^(

0.4` z1)) (3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 +

3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (

0.00021404109589041095` (E^(-0.4` z2) + E^(

0.4` z2)) (3.696` E^(0.4` z1) z1 +

3.696` E^(0.4` z1 + 0.8` z2) z1 +

3.696` E^(0.8` z1 + 0.4` z2) z2 +

3.696` E^(0.4` z2) z2))/(

E^(0.8` z1) - E^(0.8` z2))))) == 0}, {z1, 3, 6}, {z2, 9,

17}, {t0, 0, 0.00837}]

{z1 -> 10.9876, z2 -> 14.529, t0 -> -0.169798}

惢安

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
三元方程组牛顿 java_萌新请教牛顿法求解三元非线性方程组

该楼层疑似违规已被系统折叠隐藏此楼查看此楼请问牛顿迭代法求解三元非线性方程组，不同迭代初值貌似有很多不同结果，如何求解到满足要求的解，0FindRoot[{-6.565159793181527`*^-6 (0.008850684113090817` -0.7905181223796485` (0.008850684113090817` - t0)) z1 -1.3980263449472176`...
复制链接

扫一扫