24小时热门版块排行榜    

查看: 1790  |  回复: 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)所以来请教.谢谢大家~刚刚用没金币。。。。听师兄说小木虫有很多大神,没什么金币希望能帮忙~
回复此楼

» 猜你喜欢

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

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

xcenxcen

银虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
lim(r->0)c'(r)/r=c''(0)
f2(0)=lim(r->0)(-2/r*c'(r))+5c(0)=-2c''(0)+5c(0)=-2c2'+5c1
待我长发及腰,遮住一身肥膘。纵然虎背熊腰,也要高冷傲娇。
2楼2014-05-13 18:58:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer



小木虫: 金币+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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

youxinda

金虫 (小有名气)

引用回帖:
2楼: Originally posted by xcenxcen at 2014-05-13 18:58:46
lim(r->0)c'(r)/r=c''(0)
f2(0)=lim(r->0)(-2/r*c'(r))+5c(0)=-2c''(0)+5c(0)=-2c2'+5c1

c2'=-2/r*c2+5c1,那么c2‘(0)还是不能算呀。我也有试过这么做。。还是谢了~
4楼2014-05-14 16:34:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

youxinda

金虫 (小有名气)

引用回帖:
3楼: Originally posted by peterflyer at 2014-05-13 21:08:21
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/ds+2*=-5*dU/ds
U(s)=A+c(0)/* ...

不好意思啊。要求的不是解析解是数值解哦。
5楼2014-05-14 16:42:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

嘻哈学士

银虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
楼主我也困扰,你解决了吗?
6楼2015-01-30 10:20:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询


小木虫: 金币+0.5, 给个红包,谢谢回帖
是这样,在r=0附近不能直接采用数值积分,要先做级数展开,得到r=0附近的展开式,然后在稍远离r=0的地方开始使用数值方法。
7楼2015-01-30 22:16:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 youxinda 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 317求调剂 +7 申子申申 2026-03-19 12/600 2026-03-20 22:45 by JourneyLucky
[考研] 287求调剂 +7 晨昏线与星海 2026-03-19 8/400 2026-03-20 22:19 by JourneyLucky
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 5/250 2026-03-20 22:11 by 云游重阳
[考研] 一志愿南昌大学,327分,材料与化工085600 +8 Ncdx123456 2026-03-19 8/400 2026-03-20 22:06 by JourneyLucky
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 316求调剂 +5 梁茜雯 2026-03-19 5/250 2026-03-20 21:26 by 你好你好666
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 286求调剂 +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 085601专硕,总分342求调剂,地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 机械专硕325,寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
信息提示
请填处理意见