| 查看: 2786 | 回复: 9 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
owxiao木虫 (小有名气)
|
[求助]
求救，谁会解这个二阶微分方程，重金答谢
|
|||
求救。 之后还有金币。。。 |
回复此楼» 收录本帖的淘帖专辑推荐
 感兴趣的问题 |
» 猜你喜欢
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有113人回复
-
课题组招1-2名硕士
已经有2人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高金求助matlab解微分方程组
已经有12人回复
- 那位大侠，帮帮忙，求助一个二阶微分方程的解，谢谢！
已经有10人回复
- 数据如何根据二阶微分方程进行拟合
已经有5人回复
- 用MATLAB ode45求解2阶微分方程
已经有4人回复
- 【求助】怎么计算二阶的偏微分方程
已经有7人回复
- 【求助】一个二阶偏倒方程变换成微分方程，看不懂，求高手指点
已经有5人回复
- 【求助】求解一个二阶非线性微分方程x''+Ax-B/(x^3)=C
已经有6人回复
- 【求助】一个二阶常微分方程的解法（重奖金币！）
已经有16人回复
- 【求助】有限元法求解二阶偏微分方程组
已经有6人回复
- 【求助】如何解这个常微分方程【已完结】
已经有9人回复
- 【求助】如何用Runge-Kutta迭代求解二阶常微分方程组【已解决】
已经有9人回复
刘伟利
木虫 (文坛精英)
- 应助: 14 (小学生)
- 金币: 1881.1
- 红花: 8
- 帖子: 14472
- 在线: 767.5小时
- 虫号: 1339791
- 注册: 2011-07-07
- 性别: GG
- 专业: 化工系统工程
【答案】应助回帖
★ ★
dbb627(金币+1): 有问题，可以名言 2011-11-14 23:38:04
dbb627(金币-3): 这个属于违规复制回答 2011-11-14 23:52:47
dbb627(金币+1): 有问题，可以名言 2011-11-14 23:38:04
dbb627(金币-3): 这个属于违规复制回答 2011-11-14 23:52:47
|
function Exercise110 clear all clc a = 0; b = 1; solinit = bvpinit(linspace(a,b,11),[0 0]); S = diag([0 -1]); % 奇异非线性微分方程组y' = S*y/x+f(x,y)的S矩阵,则ODEfun = f(x,y) options = bvpset('SingularTerm',S); % 设置奇异性参数S矩阵选项 sol = bvp4c(@ODEfun,@BCfun,solinit,options); x = 0:0.02:1; y = deval(sol,x); plot(x,y(1,,'ro-') xlabel('x') ylabel('y') function dy=ODEfun(x,y) dy1=y(2); dy2=2*y(1)/(1+y(1)); dy=[dy1;dy2]; function bc = BCfun(ya,yb) bc = [ya(2);yb(1)-1]; 如果一窍不通,请先看书学习回来再开口 2011-11-14 22:45:37 2 TA的主页 加收听 加好友 只看回帖 发短消息 送红花 引用回复 编辑 通知版主处理 dbb627 实习版主 (小有名气) 计算强帖: 2 听众: 9 红花: 5 贵宾: 0.012 金币: 2942.6 散金: 73 帖子: 231 在线: 194.3小时 虫号: 149791 注册: 2005-12-29 活跃: 2011-11-14 状态: 在线 性别: GG 专业: 污染控制化学 管辖: 计算模拟 【答案】应助回帖 注意，这里换了一个条件，dy/dx（x=0）=0 因为x=0 无意义 solinit = bvpinit(linspace(0,1,10),[5;1]); sol = bvp4c(odefun,bcfun,solinit); xint = linspace(0,1,30); yint = deval(sol,xint); hold on; plot(xint,yint(1,,'rp'); ??? Error using ==> bvp4c at 252 Unable to solve the collocation equations -- a singular Jacobian encountered CODE: [Copy to clipboard] odefun=@(x,y)[y(2);2*y(1)/(1+y(1))-y(2)/x]; bcfun=@(ya,yb)[ya(2);yb(1)-1]; solinit = bvpinit(linspace(1e-7,1,10),[5;1]); sol = bvp4c(odefun,bcfun,solinit); xint = linspace(1e-7,1,30); yint = deval(sol,xint); hold on; plot(xint,yint(1,  ,'-rp');结果如图  |
4楼2011-11-14 23:04:56
change0618
铁杆木虫 (著名写手)
方丈大师
- 计算强帖: 1
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
2楼2011-11-14 22:45:37
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 计算强帖: 12
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
【答案】应助回帖
★ ★ ★
cenwanglai(金币+3): 谢谢回复！ 2011-11-15 07:11:26
owxiao(金币+30): 2011-11-15 12:27:16
cenwanglai(金币+3): 谢谢回复！ 2011-11-15 07:11:26
owxiao(金币+30): 2011-11-15 12:27:16
|
注意，这里换了一个条件，dy/dx（x=0）=0 因为x=0 solinit = bvpinit(linspace(0,1,10),[5;1]); sol = bvp4c(odefun,bcfun,solinit); xint = linspace(0,1,30); yint = deval(sol,xint); hold on; plot(xint,yint(1, ,'rp');??? Error using ==> bvp4c at 252 Unable to solve the collocation equations -- a singular Jacobian encountered 结果如图  |
3楼2011-11-14 22:56:01
owxiao
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2160.5
- 红花: 1
- 帖子: 180
- 在线: 104.3小时
- 虫号: 1196613
- 注册: 2011-01-26
- 性别: GG
- 专业: 环境工程
5楼2011-11-15 00:18:53
