24小时热门版块排行榜    

查看: 3024  |  回复: 11
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zzstephen

新虫 (初入文坛)

[求助] 【紧急求助】 带有微分平方项的ODE求解

dx/(1+(dx)^2) + x = 0.1sin(2t)   

想用 Matlab 解,自己弄了很久了  dsolve 解不出来, 又不知道怎么转化成ode45 能解的格式,请大家帮帮忙,解析解或者数值解的方法都好
回复此楼

» 猜你喜欢

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

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

yimeng123456

铁虫 (小有名气)

【答案】应助回帖


soliton923(金币+1): 谢谢参与讨论~~~ 2011-11-06 18:01:55
换成微分方程组呗 令dx=y ,再对y求导
3楼2011-11-06 17:37:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

guojunhang

金虫 (正式写手)


soliton923(金币+1): 谢谢参与讨论~~~ 2011-11-06 18:01:45
现在感觉这个问题属于“提法不当”,但是又看不出来。
我认为你这个方程的形式要稍微变换一下,再写出来。
变成一般的微分方程,大家好讨论
比如说
f(t,x,dx/dt,....) = 0
的形式。
2楼2011-11-06 14:40:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meveny

银虫 (初入文坛)

【答案】应助回帖

我试了下,貌似没法直接用MATLAB的ode45求解,因为求不出dx.是否可以先用符号函数求解工具(maple?)解二次方程求出dx,再用matlab的ode45做?
4楼2011-11-06 19:48:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

zzstephen(金币+1): 多谢提供解题思路!后来用matlab 的 ode15i 可以解,对比之后跟你的答案相同。 2011-11-07 00:37:59
有点类似隐式微分方程了,用1stOpt应该可以求解,假设x初值为0.1,t范围=【0,2*pi】:
CODE:
Variable t=[0,2*pi], x=0.1;
Plot x;
ODEFunction x' = (0.1*sin(2*t)-x)*(1+(x')^2);

t        x(t)
0        0.1
0.20943951        0.0850711745007771
0.41887902        0.0801305233227716
0.62831853        0.0813026887426947
0.83775804        0.0846100642268956
1.0471976        0.0864196821522688
1.2566371        0.0839373785627877
1.4660766        0.0755683106928769
1.6755161        0.0610888133870998
1.8849556        0.0416396783046324
2.0943951        0.019519065516734
2.3038346        -0.00226168631502244
2.5132741        -0.0206005920207564
2.7227136        -0.0329220995968774
2.9321531        -0.0376176957124688
3.1415927        -0.0342741951426507
3.3510322        -0.0237389947743921
3.5604717        -0.00803680859667298
3.7699112        0.00991370306354871
3.9793507        0.026808362307764
4.1887902        0.039582546783106
4.3982297        0.0459562014148314
4.6076692        0.0447762626323625
4.8171087        0.0361645432606872
5.0265482        0.0215109045350385
5.2359878        0.00329051787674805
5.4454273        -0.0153456621450357
5.6548668        -0.0311675068120555
5.8643063        -0.0414765431671204
6.0737458        -0.0445541890163028
6.2831853        -0.0398986264011418


5楼2011-11-06 21:43:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026国自然函评费到账 +19 羊腰板 2026-08-21 22/1100 2026-08-25 10:38 by popular289
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 今天基金会出结果吗?20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 范进中举一文的中心思想 +7 炎黄贵胄 2026-08-22 8/400 2026-08-25 08:48 by ZJTJZ
[基金申请] 只有每年这种时候来逛逛小木虫 +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[基金申请] 明天应该可查了!? +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见