24小时热门版块排行榜    

查看: 1649  |  回复: 5

sdcdzt

新虫 (初入文坛)

[求助] 求解二阶非线性常微分方程!!!解决之后另有重谢!!! 已有2人参与

dz(2)和dz(4)是z1,z2的二阶导,z(2),z(4)是z1,z2的一阶导。方程里最高阶项是耦合的。



dz(2) + (8165768846995289*z(1))/43980465111040 + 260*z(2) + (729274196212461*z(1)^2*z(3)^2)/657666867200 - (810461751221981*z(1)^2*z(3)^3)/770703360000 + (1056424101821091*z(1)^3*z(3)^2)/2055208960000 + (8722212816045573*dz(4)*z(3))/3777893186295716170956800 - (6306482750015563*z(1)*z(3))/858993459200 + (7086797913037029*z(3)*z(4))/11805916207174113034240 - (756465654738930894864411090880955416576*v^2*((5583899589019513*z(1)^2)/28823037615171174400 - (1535799397223231*z(1)*z(3))/2882303761517117440 + (1705094118577331*z(3)^2)/1801439850948198400))/1726659758420692373277280879638240234375 + (6240407844064763*dz(4)*z(3)^2)/483570327845851669882470400 - (3584792003555446609710148086013952*v^2*((5616914817747799*z(1))/219902325555200 - (4981953625221941*z(3))/219902325555200))/1726659758420692373277280879638240234375 + (4800904038295969*z(1)*z(3)^2)/140928614400 - (4926513293399861*z(1)^2*z(3))/493250150400 - (7448269298026523*z(1)*z(3)^3)/3946001203200 - (1754174731218613*z(1)^3*z(3))/8418135900160 + (1349947142680613*z(1)*z(3)^4)/1233125376000 - (4778439099191479*z(1)^4*z(3))/49325015040000 + (253516568665131*z(3)^2*z(4))/75557863725914323419136 - (756465654738930894864411090880955416576*v^2*((5584448069578993*z(1)^2)/115292150460684697600 - (6186436149201411*z(1)*z(3))/28823037615171174400 + (5375473397785631*z(3)^2)/11529215046068469760))/1726659758420692373277280879638240234375 + (41486480788177196822439296062044766208*v^2)/115110650561379491551818725309216015625 + (1375585107332589*z(1)^2)/2748779069440 + (5931535843587817*z(1)^3)/6313601925120 + (3741979822281889*z(1)^4)/101017630801920 + (4577882713395029*z(3)^2)/214748364800 + (6818002503693121*z(1)^5)/526133493760000 - (262863617127891*z(3)^3)/46976204800 + (3438474728478809*z(3)^4)/1973000601600 - (811723297025401*z(3)^5)/112742891520000=0
dz(4) + (297872865246189*z(1))/5629499534213120 + 260*z(4) - (563948109211209*z(1)^2*z(3)^2)/858993459200 + (1509940229868547*z(1)^2*z(3)^3)/2055208960000 - (899018850540299*z(1)^3*z(3)^2)/3758096384000 + (3584792003555446609710148086013952*v^2*((311265982272027*z(1))/13743895347200 - (5439582840008959*z(3))/27487790694400))/1726659758420692373277280879638240234375 - (756465654738930894864411090880955416576*v^2*((7021818610175179*z(1)*z(3))/7205759403792793600 - (5967775428393763*z(1)^2)/57646075230342348800 + (681586795427753*z(3)^2)/1441151880758558720))/1726659758420692373277280879638240234375 + (3298744320994199*dz(4)*z(3))/236118324143482260684800 + (4647990245028827*z(1)*z(3))/214748364800 + (5360459521615573*z(3)*z(4))/1475739525896764129280 + (4720249406991477*dz(4)*z(3)^2)/60446290980731458735308800 - (2922476536051257*z(1)*z(3)^2)/164416716800 + (774060836814941*z(1)^2*z(3))/35232153600 + (19783176809747*z(1)*z(3)^3)/11744051200 + (1753869418212413*z(1)^3*z(3))/5261334937600 - (6448419081252833*z(1)*z(3)^4)/31568009625600 + (2956279016688527*z(1)^4*z(3))/32883343360000 + (767040528636115*z(3)^2*z(4))/37778931862957161709568 + (31759674132034291766388478409227370496*v^2)/191851084268965819253031208848693359375 - (4870310409024751*z(1)^2)/43980465111040 - (4158104683784381*z(1)^3)/6313601925120 - (7911407438985527*z(1)^4)/1010176308019200 + (4164490152882159*z(3)^2)/274877906944 - (2403471885628813*z(1)^5)/1578400481280000 + (1454845076733109*z(3)^3)/24662507520 + (890364628444807*z(3)^4)/1683627180032 + (1162785006804349*z(3)^5)/2466250752000 - (756465654738930894864411090880955416576*v^2*((1612342847780083*z(1)*z(3))/115292150460684697600 + (3202078393359277*z(1)^2)/57646075230342348800 + (3408158844499937*z(3)^2)/28823037615171174400))/1726659758420692373277280879638240234375=0
回复此楼

» 猜你喜欢

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

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

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
方程很复杂,不过解二阶非线性方程的方法还是很成熟的:
比如  y''=f
令  y1=y, y2=y'
原方程化为常微分方程组:
y1'=y2
y2'=f

加上初值或者边值条件,MATLAB的ode45或者ode15s ode23s即可解。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2014-04-15 21:35:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
2楼: Originally posted by 月只蓝 at 2014-04-15 21:35:43
方程很复杂,不过解二阶非线性方程的方法还是很成熟的:
比如  y''=f
令  y1=y, y2=y'
原方程化为常微分方程组:
y1'=y2
y2'=f

加上初值或者边值条件,MATLAB的ode45或者ode15s ode23s即可解。

补充一下
MATLAB的ode45或者ode15s ode23s适合于初值问题哈
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
3楼2014-04-15 21:39:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sdcdzt

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by 月只蓝 at 2014-04-15 21:39:23
补充一下
MATLAB的ode45或者ode15s ode23s适合于初值问题哈...

这些都用过。都不适用。因为我这个方程里连最高阶项都是二次的。ode45根本解不了。
4楼2014-04-17 19:22:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sdcdzt

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by 月只蓝 at 2014-04-15 21:39:23
补充一下
MATLAB的ode45或者ode15s ode23s适合于初值问题哈...

说错了,最高阶项是非线性的
5楼2014-04-17 19:24:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onesupeng

金虫 (职业作家)

【答案】应助回帖

自己编程用迭代法去求解看看

另外,两个未知量一个方程?
长期招收博士生,参见http://fsl-unsw.com
6楼2014-05-01 13:23:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sdcdzt 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳又变了 +11 wuchongjun 2026-08-20 16/800 2026-08-20 19:26 by wuchongjun
[基金申请] 科研孤儿太难了 +9 我4大白菜 2026-08-20 9/450 2026-08-20 18:23 by 2549091707
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 22/1100 2026-08-20 17:03 by yaoyewhu2008
[基金申请] 放榜前的不淡定 40+3 snowwithsea 2026-08-19 12/600 2026-08-20 15:56 by lfy8008
[基金申请] 言之凿凿,明天官网会通知放榜? +14 医学老男孩 2026-08-20 16/800 2026-08-20 15:47 by LNP@mRNA
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见