找回密码
 欢迎注册
查看: 11346|回复: 3

[求助] 非线性方程组 牛顿法中某个方程乘上倍数a后对 雅克比矩阵有多大影响。

[复制链接]
发表于 2016-3-21 15:33:23 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?欢迎注册

×
我本来有“一批”方程组,只是常数项不同。有少量解不出来。昨天我把方程组其中一个乘上差不多 10^10(分母消了公约数)之后,就大部分解不出来了。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2016-3-22 11:06:51 | 显示全部楼层
你用的什么软件?可能是精度的问题。
我的一个项目中的需要求解大量的(可多达几亿个)整数系数 的线性方程组的 整数解,程序是自己编写的,大数计算部分用到GMP大整数库。这些方程组中,绝大部分是无解的。但是随便挑一个无解的线性方程组用·matlab或者matamatica来算,是有解的,因为他们都是基于浮点数的。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2016-3-22 19:31:17 | 显示全部楼层
liangbch 发表于 2016-3-22 11:06
你用的什么软件?可能是精度的问题。
我的一个项目中的需要求解大量的(可多达几亿个)整数系数 的线性方 ...

非线性,带指数和对数,必须浮点。
50位和100位的迭代记录,差不多相等。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2016-3-23 16:07:54 | 显示全部楼层
zeroieme 发表于 2016-3-22 19:31
非线性,带指数和对数,必须浮点。
50位和100位的迭代记录,差不多相等。


要看具体情况。10^10,万一数值算法的计算过程中,还要做一次幂运算,就更大了。数值计算最忌 非同等数量级的两数发生二元运算了
====
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2024-11-21 21:15 , Processed in 0.034163 second(s), 16 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表