| 查看: 2170 | 回复: 6 | |||
youxinda金虫 (小有名气)
|
[交流]
关于打靶法解二阶微分方程的数值解的问题 已有4人参与
|
|
原题 c''+2/r*c'=5c 0≤r≤1 c'(0)=0,c(1)=1 转化为二元方程组求解：设c1=c,c2=c'=f1（r)；则c2'=-2/r*c2+5*c1=f2(r)；c2(0)=0,令c1(0)=M,解出c1(1)即c(1)使其满足值为1. 也就是转化成 c1'=c2=f1 c2'=-2/r*c2+5c1=f2 c1(0)=m,c2(0)=0,目标是c（1）=1也就是c1（1）=1 求取采用四阶R-K（龙格库塔）并迭代m。 但是问题是用R-K过程中，需要f2（0），但是c2'=-2/r*c2+5*c1=f2（r)中的r不能为零，于是出现问题。不知道该如何解决，老师说可以用罗比达法则求极限来求f2（0）所以来请教.谢谢大家~刚刚用没金币。。。。听师兄说小木虫有很多大神，没什么金币希望能帮忙~ |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有14人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有15人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有17人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有21人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有21人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请问这个偏微分方程组有一般的解法么？
已经有12人回复
- 微分方程的数值解求法
已经有3人回复
- 如何把一个方程的时域数值解变换到频域数值解？
已经有5人回复
- 怎么解这个复数方程
已经有3人回复
- 2阶偏微分方程的解法
已经有11人回复
- 如何用 数值解法 求 二阶微分方程的本征值
已经有18人回复
- 求解一个二阶偏微分方程~~有结果，但是不知道中间过程（完美解决者答谢100￥）~
已经有4人回复
- 抛物型偏微分方程的数值解求解问题，谢谢！
已经有14人回复
- 如何用matlab求解矩阵系数的二阶微分方程
已经有13人回复
- 求助一阶非线性微分方程的解
已经有7人回复
- 微分方程解的个数与方程的阶数关系
已经有5人回复
- 如何求解描述振动的二阶微分方程
已经有21人回复
- 【求助】怎么计算二阶的偏微分方程
已经有7人回复
- 常微分方程(ODEs)的MATLAB数值解法
已经有119人回复
- 【求助】请教二阶微分方程的特解（含初始条件）
已经有8人回复
- 【求助】常微分方程组求解中系数与某变量值关联的问题
已经有12人回复
- 【求助】有限元法求解二阶偏微分方程组
已经有6人回复
- 【求助】先解常微分方程 然后画图【已完结】
已经有11人回复
2楼2014-05-13 18:58:46
peterflyer 
木虫之王 (文学泰斗)
peterflyer
- 数学EPI: 10
- 应助: 20282 (院士)
- 金币: 146755
- 红花: 1374
- 帖子: 93091
- 在线: 7695.5小时
- 虫号: 1482829
- 注册: 2011-11-08
- 性别: GG
- 专业: 功能陶瓷
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
c''+2/r*c'=5c 0≤r≤1 (1) c'(0)=0,c(1)=1 (2) (1)两边同乘r: r*c''+2*c'=5*r*c (3) 设U(s)为c的拉氏变换，（3）两边求拉氏变换，得到： -d[s^2*U-s*c(0)]/ds+2*[s*U-c(0)]=-5*dU/ds U(s)=A+c(0)/[2*sqrt(5)]*Ln{[s+sqrt(5)]/[s-sqrt(5)]} c(r)=A*δ(r)+c(0)/[2*sqrt(5)]*L^(-1){Ln{[s+sqrt(5)]/[s-sqrt(5)]}} L^(-1){Ln{[s+sqrt(5)]/[s-sqrt(5)]}}是Ln{[s+sqrt(5)]/[s-sqrt(5)]}的拉氏逆变换，δ(r)为阶跃函数。 |
3楼2014-05-13 21:08:21
youxinda
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1102.4
- 散金: 164
- 红花: 1
- 帖子: 53
- 在线: 101.7小时
- 虫号: 2947843
- 注册: 2014-01-25
- 性别: GG
- 专业: 分离过程
4楼2014-05-14 16:34:23
youxinda
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1102.4
- 散金: 164
- 红花: 1
- 帖子: 53
- 在线: 101.7小时
- 虫号: 2947843
- 注册: 2014-01-25
- 性别: GG
- 专业: 分离过程
5楼2014-05-14 16:42:36
6楼2015-01-30 10:20:11
pippi6
铁杆木虫 (著名写手)
工程和科学数值计算咨询
- 数学EPI: 6
- 应助: 413 (硕士)
- 贵宾: 0.002
- 金币: 7116.5
- 散金: 15
- 红花: 63
- 帖子: 1639
- 在线: 798.9小时
- 虫号: 2469437
- 注册: 2013-05-14
- 专业: 计算数学与科学工程计算
7楼2015-01-30 22:16:32
20