24小时热门版块排行榜    

查看: 2052  |  回复: 8

949497669

新虫 (初入文坛)

[求助] 求助如何用matlab求非线性方程组的解 已有1人参与

本人matlab菜鸟一枚,现遇到一个瓶颈问题,就是求解一条曲线与曲面的交点,已知条件和要求如下:
1、已知曲线和曲面的参数方程,各参数的范围已知;
2、两者可能有交点,可能没有交点,需要首先判断有无交点;
3、对于有交点的情况,问题转化为曲线和曲面x,y,z坐标对应相等组成的非线性方程组的求解问题;可能有两个交点,也可能有一个交点;
现在尝试了solve,fsolve,和牛顿迭代法,都不行
solve解非线性方程组的能力较弱,
fsolve需要人为建立m文件,由于解方程组的过程放在了一个循环里,因此方程组求解完之后,需要改变方程再次进行求解,因此使用fsolve也不现实。
牛顿迭代法不容易确定初值,而且只能求出一个接,不适用于有多个交点的情况。
哪位大牛能给我出出主意,能方便的求出交点。谢谢!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

参数方程转换为隐函数方程:F(x, y, z)=0;G(x, y, z)=0
新方程  F(x, y, z) + k*G(x, y, z)=0
解就是交点。
2楼2015-01-04 10:02:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

949497669

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by feixiaolin at 2015-01-04 10:02:33
参数方程转换为隐函数方程:F(x, y, z)=0;G(x, y, z)=0
新方程  F(x, y, z) + k*G(x, y, z)=0
解就是交点。

关键是新方程怎么解
3楼2015-01-04 11:51:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

引用回帖:
3楼: Originally posted by 949497669 at 2015-01-04 11:51:34
关键是新方程怎么解...

solve函数

[ 发自小木虫客户端 ]
4楼2015-01-04 12:41:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

949497669

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by feixiaolin at 2015-01-04 12:41:04
solve函数
...

如何转化为隐函数
5楼2015-01-04 13:03:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

引用回帖:
5楼: Originally posted by 949497669 at 2015-01-04 13:03:01
如何转化为隐函数

参数方程转化。
比如 x=t; y=t^2; 方程就是 y-x^2=0

[ 发自小木虫客户端 ]
6楼2015-01-04 13:45:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wurongjun

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
949497669(feixiaolin代发): 金币+3 2015-01-07 08:38:15
这个问题没有大一统的方法!
可以具体问题具体分析!
建议用牛顿迭代法!如你所说这个与初值有关!你可以根据你的问题给出合适的初值!
另外,还可以用最优化算法,这个比较多!
善恶到头终有报,人间正道是沧桑.
7楼2015-01-04 20:10:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

949497669

新虫 (初入文坛)

自己顶一下
8楼2015-01-14 10:43:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

949497669

新虫 (初入文坛)

我这个问题的最大特点是在程序的循环过程中,每次循环都有一个对应的方程组,进而解出方程组的解,因此如果需要将方程组存在另一个m文件里不太现实,或者有我不知道的处理方法。另外每次循环中方程组解的范围不确定,因此用牛顿迭代的方法也有问题。
9楼2015-01-14 10:52:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 949497669 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +10 Ldrop2023 2026-08-13 10/500 2026-08-17 09:59 by 小豌豆_发芽
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +11 家与远方 2026-08-10 16/800 2026-08-16 10:28 by ray43
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
信息提示
请填处理意见