应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 运筹学与控制论 » 求助解复杂非线性方程组的好的方法
251/3返回列表123下一页
查看: 3611  |  回复: 24
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

wangfumin

金虫 (小有名气)

[求助] 求助解复杂非线性方程组的好的方法

本人是做空间坐标定位系统的，是用光学的方法利用CCD采集干涉图，然后分析得出相位分布，根据相位分布反过来求解得到发出相干光点光源的空间坐标。
主要是求解出两个点的空间三维坐标，是6个未知数。请问有没有哪位知道好的方法求解非线性方程组呢？本人试过很多方法，发现L-M算法不错的，但是仍然有些问题，主要是将所得的相位加入噪声之后，就没办法算出空间坐标值了，不知道为什么这么微小的一个变化会引起达不到迭代精度的问题。细节的东西感觉一句两句说不清楚，先请大家给出一个指导思路。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-11-18 23:18:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

feixiaolin

荣誉版主 (文坛精英)

优秀版主
引用回帖:
14楼: Originally posted by cooooldog at 2013-11-20 07:49:27
这种算法跟普通牛顿法相比有什么特点?...

Newton流线法是近两年提出的一种从Newton场的观点出发，沿流线的追踪算法。与传统求解非线性方程组的方法类似，Newton流线法基于迭代方向的收敛性进行求解，但普通的求解方法只在根附近有二次收敛性，对初值比较敏感。而Newton流线法建立了与问题维数和函数的具体性质无关的4个判断参数：减幅指标 ，降速指标 ，最大斜率 ，斜率比 ，这4个参数是可计算量，基于这些可计算的判断参数，沿Newton流可以自适应地追踪或者识别奇点，在计算过程中可以跳过奇点，所以降低了算法对初始值的敏感度；在很大程度上克服了常规方法的局部收敛性，是一种全局收敛的算法；并且可以自适应实时确定合适的步长，使得计算更精确。
赞一下(3人) 回复此楼
15楼2013-11-20 14:58:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
引用回帖:
16楼: Originally posted by cooooldog at 2013-11-21 08:24:34
专业! 非常感谢; 请问这个方法最早是国内还是国外专家提出来的? 有没有放出原始代码的网站? 英文关键词是什么?

Newton法似乎一直有一种对步长做 damping factor修正的方法,以确保全局收敛, 两者又有什么区别?...

据陈文所述，是他们提出来的。
网络可见，该陈曾到中科院数学所做有关专题报告。
赞一下(2人) 回复此楼
17楼2013-11-21 09:28:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangfumin

金虫 (小有名气)
引用回帖:
12楼: Originally posted by cooooldog at 2013-11-20 07:42:44
相位是什么方式求的? 从什么数据计算? 大概是多少维的矩阵?

很好奇, (x_k,y_k,z_k)是CCD的像的齐次坐标吗? 从图上看, F是一个标量(三个平方根),如果c是矩阵,这种减肥是怎么计算的?

没看懂...

定义以CCD平面为xoy平面，这样z_k默认就是0了。
c是和x_k,y_k,z_k同维的列向量。
相位是移相算法求得的，矩阵的维数就是CCD的网格数，640*240，从640*240中取点代入F中进行计算。
赞一下(2人) 回复此楼
18楼2013-11-21 10:17:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
非线性方程组的Newton流线法.pdf
见
http://www.cnki.com.cn/Article/CJFDTotal-JSSX201203003.htm
赞一下(2人) 回复此楼
20楼2013-11-21 10:27:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangfumin

金虫 (小有名气)
引用回帖:
22楼: Originally posted by cooooldog at 2013-11-21 12:04:43
因为c是矩阵, 前面根号下的都是标量;
我是不太明白标量怎么跟矩阵相减;
还是标量跟矩阵的每个元素elementwise的减法?

相位是如何从照片得到?用的不是普通的CCD相机?可以得到相位信息的,难道是激光全息照相...

其实方程组中每个方程的形式是一模一样的，只是选的点不同而已。相位是通过对干涉图分析得到

[ 发自小木虫客户端 ]
赞一下(2人) 回复此楼
23楼2013-11-21 12:20:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
wangfumin: 金币+5 2013-11-19 14:50:35
wangfumin: 金币+5, ★有帮助 2013-11-19 14:51:14
牛顿流线法
陈传淼, 胡宏伶, 雷蕾, 等. 非线性方程组的 Newton 流线法[J]. 计算数学, 2012, 34(3): 235-258.
赞一下(1人) 回复此楼
2楼2013-11-18 23:24:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็
看上去很高端,有没有直接一点的数学模型?
给两个方程组的例子出来
赞一下(1人) 回复此楼
ส็็็็็็็็็็็็็็็็็็็็
4楼2013-11-19 09:54:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็
引用回帖:
5楼: Originally posted by wangfumin at 2013-11-19 13:36:41
说白了就是空间点的三维坐标计算
其实数学模型还是比较简单，只不过描述起来不太方便。

问题.jpg
...

相位是什么方式求的? 从什么数据计算? 大概是多少维的矩阵?

很好奇, (x_k,y_k,z_k)是CCD的像的齐次坐标吗? 从图上看, F是一个标量(三个平方根),如果c是矩阵,这种减肥是怎么计算的?

没看懂
赞一下(1人) 回复此楼
ส็็็็็็็็็็็็็็็็็็็็
12楼2013-11-20 07:42:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็
引用回帖:
7楼: Originally posted by feixiaolin at 2013-11-19 14:57:19
我们用了，效果不错。...

你们用什么语言实现算法的?
赞一下(1人) 回复此楼
ส็็็็็็็็็็็็็็็็็็็็
13楼2013-11-20 07:43:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็
引用回帖:
2楼: Originally posted by feixiaolin at 2013-11-18 23:24:57
牛顿流线法
陈传淼, 胡宏伶, 雷蕾, 等. 非线性方程组的 Newton 流线法. 计算数学, 2012, 34(3): 235-258.

这种算法跟普通牛顿法相比有什么特点?
赞一下(1人) 回复此楼
ส็็็็็็็็็็็็็็็็็็็็
14楼2013-11-20 07:49:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็
引用回帖:
15楼: Originally posted by feixiaolin at 2013-11-20 14:58:14
Newton流线法是近两年提出的一种从Newton场的观点出发，沿流线的追踪算法。与传统求解非线性方程组的方法类似，Newton流线法基于迭代方向的收敛性进行求解，但普通的求解方法只在根附近有二次收敛性，对初值比较敏 ...

专业! 非常感谢; 请问这个方法最早是国内还是国外专家提出来的? 有没有放出原始代码的网站? 英文关键词是什么?

Newton法似乎一直有一种对步长做 damping factor修正的方法,以确保全局收敛, 两者又有什么区别?
赞一下(1人) 回复此楼
ส็็็็็็็็็็็็็็็็็็็็
16楼2013-11-21 08:24:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangfumin

金虫 (小有名气)
引用回帖:
17楼: Originally posted by feixiaolin at 2013-11-21 09:28:05
据陈文所述，是他们提出来的。
网络可见，该陈曾到中科院数学所做有关专题报告。...

这个算法一年前提出来的，感觉相关资料不多，网上也没有相关代码。
有没有相关资料可以提供一下呢？
赞一下(1人) 回复此楼
19楼2013-11-21 10:21:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็
引用回帖:
20楼: Originally posted by feixiaolin at 2013-11-21 10:27:59
非线性方程组的Newton流线法.pdf
见
http://www.cnki.com.cn/Article/CJFDTotal-JSSX201203003.htm

我个人感觉, 因为Numerical recipes 这样的书上,2007第三版,已经提到类似的damping因子改善Newton方法全局收敛性的办法,很可能这种改进的意义并不显著;

否则,数学方面的文章应该是相对容易投到国外的好一些的刊物的;
从第一作者通讯作者的资历看,应该也不是写英文有困难的.

没有经过外国同行认证过的创新性,其作用的显著性和创新的大小,可能都得存疑.
赞一下(1人) 回复此楼
ส็็็็็็็็็็็็็็็็็็็็
21楼2013-11-21 12:01:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็
引用回帖:
18楼: Originally posted by wangfumin at 2013-11-21 10:17:05
定义以CCD平面为xoy平面，这样z_k默认就是0了。
c是和x_k,y_k,z_k同维的列向量。
相位是移相算法求得的，矩阵的维数就是CCD的网格数，640*240，从640*240中取点代入F中进行计算。...

因为c是矩阵, 前面根号下的都是标量;
我是不太明白标量怎么跟矩阵相减;
还是标量跟矩阵的每个元素elementwise的减法?

相位是如何从照片得到?用的不是普通的CCD相机?可以得到相位信息的,难道是激光全息照相
赞一下(1人) 回复此楼
ส็็็็็็็็็็็็็็็็็็็็
22楼2013-11-21 12:04:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

wangfumin

金虫 (小有名气)
引用回帖:
2楼: Originally posted by feixiaolin at 2013-11-18 23:24:57
牛顿流线法
陈传淼, 胡宏伶, 雷蕾, 等. 非线性方程组的 Newton 流线法. 计算数学, 2012, 34(3): 235-258.

谢谢 明天找来看一下先

[ 发自小木虫客户端 ]
赞一下 回复此楼
3楼2013-11-18 23:58:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangfumin

金虫 (小有名气)
引用回帖:
4楼: Originally posted by cooooldog at 2013-11-19 09:54:48
看上去很高端,有没有直接一点的数学模型?
给两个方程组的例子出来

说白了就是空间点的三维坐标计算
其实数学模型还是比较简单，只不过描述起来不太方便。
求助解复杂非线性方程组的好的方法
问题.jpg
赞一下 回复此楼
5楼2013-11-19 13:36:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangfumin

金虫 (小有名气)
引用回帖:
2楼: Originally posted by feixiaolin at 2013-11-18 23:24:57
牛顿流线法
陈传淼, 胡宏伶, 雷蕾, 等. 非线性方程组的 Newton 流线法. 计算数学, 2012, 34(3): 235-258.

这位仁兄，此文献数学知识比较广泛，没怎么看懂。请问你有没有实施过这个方法求解复杂的非线性方程，效果怎样？
赞一下 回复此楼
6楼2013-11-19 14:50:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
引用回帖:
6楼: Originally posted by wangfumin at 2013-11-19 14:50:28
这位仁兄，此文献数学知识比较广泛，没怎么看懂。请问你有没有实施过这个方法求解复杂的非线性方程，效果怎样？...

我们用了，效果不错。
赞一下 回复此楼
7楼2013-11-19 14:57:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangfumin

金虫 (小有名气)
引用回帖:
7楼: Originally posted by feixiaolin at 2013-11-19 14:57:19
我们用了，效果不错。...

你们求的是什么？也是空间坐标吗？求解的是几个未知数呢
赞一下 回复此楼
8楼2013-11-19 20:27:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
引用回帖:
8楼: Originally posted by wangfumin at 2013-11-19 20:27:48
你们求的是什么？也是空间坐标吗？求解的是几个未知数呢...

复现了论文算法，与L-M, BFGS等方法做了比较。
赞一下 回复此楼
9楼2013-11-19 21:24:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangfumin

金虫 (小有名气)
引用回帖:
9楼: Originally posted by feixiaolin at 2013-11-19 21:24:46
复现了论文算法，与L-M, BFGS等方法做了比较。...

发论文没，发了我下下来自己钻研
赞一下 回复此楼
10楼2013-11-19 22:59:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wangfumin 的主题更新
251/3返回列表123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +7 scalable 2026-08-24 7/350 2026-08-24 23:01 by anata1209
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-24 21:59 by BZKMTicpDhFj
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 0XLacIJUOj8D 2026-08-24 5/250 2026-08-24 21:51 by BZKMTicpDhFj
[基金申请] 人气不行了 +10 fansofjerry 2026-08-21 10/500 2026-08-24 21:03 by zhanghaozhu
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 明天应该可查了！？ +4 chengyan1220 2026-08-23 4/200 2026-08-24 08:47 by gloomy6159
[基金申请] 什么时候开奖？ +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
信息提示
关闭
请填处理意见
关闭