python牛顿法解非线性方程组_萌新请教牛顿法求解三元非线性方程组
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
请问牛顿迭代法求解三元非线性方程组,不同迭代初值貌似有很多不同结果,如何求解到满足要求的解,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}
python牛顿法解非线性方程组_萌新请教牛顿法求解三元非线性方程组相关推荐
- python数据分析常用的算法_萌新向Python数据分析及数据挖掘 第三章 机器学习常用算法 第二节 线性回归算法 (上)理解篇...
理解 以a b为变量,预测值与真值的差的平方和为结果的函数 参数学习的基本方法:找到最优参数使得预测与真实值差距最小 假设可以找到一条直线 y = ax+b 使得预测值与真值的差的平方和最小 故事 假 ...
- 如何系统地自学python100天_Github上发布了一个Python学习秘笈,从萌新到王者的100天Python学习之旅...
北京千锋互联科技有限公司成都分公司骆昊(jackfrued)在Github上发布了一个Python学习秘笈,从萌新到王者的100天Python学习之旅. 简单的说,Python是一个"优雅& ...
- ctfshow_萌新_萌新隐藏题
https://harvey-blog.com/Safety/525 0x01 萌新认证 加群召唤flag 0x02 萌新_密码 萌新不会 53316C6B5A6A42684D3256695A4456 ...
- python新手入门到放弃_python萌新:从零基础入门到放弃
原标题:python萌新:从零基础入门到放弃 不管是在什么领域,自学者都占绝大多数,你说自学可以吗?可以,没问题的,只需要你具备以下几点最基础的能力: 第一点:天赋.对于python而言其实是非常需要 ...
- python语法报错原因_python萌新最常见5大错误
Python 小萌新最常见 5 大错误的: 1. invalid syntax 报错原因:语法错误 首先这个时候需要检查语法是否有用错,如果报错行没有发现错误,通常在相 邻行会出现问题,这个时候就要检 ...
- python input与返回值-Python 详解基本语法_函数_返回值
Python 详解基本语法 概要: 函数的返回值是函数重要的组成部分.函数的根本在于实现程序的部分功能,所以很多时候我们需要将函数执行后的结果返回给程序再由程序作出进一步的操作.可以说是函数的返回值令 ...
- 组装台式电脑配置清单_萌新攒机必备!多价位台式电脑配置清单!
首先需要声明的是,配置单仅做参考,为了给一些攒机萌新提供可以对比的选择,没有绝对完美的配置,都要根据需求和预算进行取舍.不同人对于品牌也有不同程度的偏好,可能因为之前用过类似的产品得到了糟糕的体验,从 ...
- 声卡loopback有什么用_萌新做音乐那点事 | 外置专业声卡的选择方法与推荐
[本文长期更新,感兴趣可以点收藏] 前言:作为一名独立音乐人,接触了这么多年,自认为了解一点声卡的常识,为正在准备购买声卡[或是音频接口]的萌新音乐人做一篇介绍文. 市面常见的声卡品牌 常见的品牌真的 ...
- excel自动保存_萌新必看!python处理excel实例
电脑上必安装的办公软件三大软件:word.excel.ppt,其中excel是在处理数据效率最高,也是最为繁琐的,因此高效快捷使用excel尤为重要,以下为大家介绍使用python更自动化处理exce ...
最新文章
- 8086汇编-实验9-字符打印
- win7 安装mysql 5.7.9记录
- 论文浅尝 | AAAI2020 - 基于生成对抗的知识图谱零样本关系学习
- 数据结构实验之二叉树的建立与遍历
- jquery+ajax+ashx
- wsl2无法使用systemctl_Win 10 更新,Linux 内核的 WSL 2 开始上线
- 【codevs1228】苹果树,哦
- Facebook斥资5亿美元 建设全风电数据中心
- mysql之使用python操作mysql需要commit的情况
- 错误未找到引用源_你好,C++(77)12.1 用右值引用榨干C++的性能
- 190616每日一句
- SPSS--Friedman检验步骤
- android als传感器,环境光传感器(ALS)背光控制系统解决方案
- openstack虚拟机配置vip
- 给爱车安装一台GPS,贴身“保镖”保护爱车
- html怎么用wps打开pdf文件,怎么使用电脑wps对PDF文件进行修改
- Mac Sublime Text安装可以使xml格式化的插件
- 小曾WRF自学日记(3)渐入佳境 ——WRF实例-数据下载与WPS前处理
- 给 Java 初学者的学习路线建议
- 静态IP和动态IP有什么区别?