24小时热门版块排行榜    

查看: 1400  |  回复: 7

malin877

新虫 (小有名气)

[求助] 非线性方程求解问题,求助各位兄弟姐妹们~~

本人非数学专业,但是最近实验过程中遇到了一个非线性方程的求解问题,思考数日,实在弄不懂了,烦请诸位帮忙看下以下哪个解是正确的,或者都不正确?
多谢多谢!!!
回复此楼

» 猜你喜欢

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

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

yongcailiu

金虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
malin877: 金币+2, ★★★很有帮助 2013-04-14 18:05:25
首先,区分一下符号,c^(k)表示第k次迭代的结果,应该和第一行公式中出现的k区分开。其次,如果k是常数,即第一行公式成立的话,那么迭代公式除了第二个,其他三个(其中第三个和第四个应该是N-1次方,而不是N次方)都是对的(同样要区分k)。最后,第三个格式是隐式格式,第四个格式是显式格式,一般而言,前者稳定性好,需要迭代求解。
2楼2013-04-11 10:02:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

malin877

新虫 (小有名气)

引用回帖:
2楼: Originally posted by yongcailiu at 2013-04-11 10:02:26
首先,区分一下符号,c^(k)表示第k次迭代的结果,应该和第一行公式中出现的k区分开。其次,如果k是常数,即第一行公式成立的话,那么迭代公式除了第二个,其他三个(其中第三个和第四个应该是N-1次方,而不是N次方) ...

非常感谢你的回答,是要区分一下符号 两个k 分开 一个是时间 一个是常数。不过我可不可以继续问一下,我用了matlab里的fzero 公式 写了公式里的第一行 但是结果很怪 不知是不是迭代的问题,还有就是如果用第四个 显式格式的话  是不是要有一个特定 for---end 去定义时间 不同于之前的线性?

谢谢!
3楼2013-04-14 18:05:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

malin877

新虫 (小有名气)

引用回帖:
2楼: Originally posted by yongcailiu at 2013-04-11 10:02:26
首先,区分一下符号,c^(k)表示第k次迭代的结果,应该和第一行公式中出现的k区分开。其次,如果k是常数,即第一行公式成立的话,那么迭代公式除了第二个,其他三个(其中第三个和第四个应该是N-1次方,而不是N次方) ...

还有就是如果1),4)的离散是正确的话并且都是用的显式的有限差分法,1)和4)结果能一样吗?
如果就1)和4)都用显式方法来解 到最后 就是解一个一元的不定式,完全可以用fzero作为内置函数 ???

多谢多谢!
4楼2013-04-14 23:55:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongcailiu

金虫 (小有名气)

【答案】应助回帖

★ ★
malin877: 金币+2, ★★★很有帮助 2013-04-18 23:27:29
引用回帖:
4楼: Originally posted by malin877 at 2013-04-14 23:55:41
还有就是如果1),4)的离散是正确的话并且都是用的显式的有限差分法,1)和4)结果能一样吗?
如果就1)和4)都用显式方法来解 到最后 就是解一个一元的不定式,完全可以用fzero作为内置函数 ???

多谢多谢 ...

首先,帖子中只是针对某个格式给出了它的差分形式,单单这些是不能求解的(如3x^2+x), 因此除了需要把这个微分方程的格式详细的给出,即除了左端项以及等价项(如帖子中给出的那样)外,还需要给出右端项才能使用数值方法迭代求解。如果两种数值格式都收敛到精确解(当然要要求解唯一)的话,那么(1)(4)最后的结果应该是一样的,虽然它们的迭代格式略有不同。
5楼2013-04-15 12:19:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

malin877

新虫 (小有名气)

引用回帖:
5楼: Originally posted by yongcailiu at 2013-04-15 12:19:13
首先,帖子中只是针对某个格式给出了它的差分形式,单单这些是不能求解的(如3x^2+x), 因此除了需要把这个微分方程的格式详细的给出,即除了左端项以及等价项(如帖子中给出的那样)外,还需要给出右端项才能使用数 ...

这样子对么?
6楼2013-04-18 23:27:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongcailiu

金虫 (小有名气)

【答案】应助回帖


malin877: 金币+1, ★★★很有帮助 2013-04-22 18:31:34
引用回帖:
6楼: Originally posted by malin877 at 2013-04-18 23:27:17
这样子对么?
...

只说第一个式子,我觉得,右端的第二项两个数相乘的第一个分母应该是r_j(当然如果r从0开始,步长为delta t,那么这个结果从形式上来看没问题,至于收敛性令当别论,数值结果好就接受,不好的话就应该考虑构造其他的格式) 。后面的部分你代入自己推导一下吧,貌似问题不大
7楼2013-04-19 12:44:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

malin877

新虫 (小有名气)

引用回帖:
7楼: Originally posted by yongcailiu at 2013-04-19 12:44:26
只说第一个式子,我觉得,右端的第二项两个数相乘的第一个分母应该是r_j(当然如果r从0开始,步长为delta t,那么这个结果从形式上来看没问题,至于收敛性令当别论,数值结果好就接受,不好的话就应该考虑构造其他的 ...

非常感谢~~~
8楼2013-04-22 18:31:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 malin877 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 25/1250 2026-08-16 10:53 by hmhminy
[基金申请] 时间戳又变了8-15 +12 archvillain 2026-08-15 23/1150 2026-08-15 20:47 by foyo404
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +7 Ldrop2023 2026-08-13 7/350 2026-08-15 15:57 by mzhh2000
[基金申请] 关于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
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 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
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见