24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1892  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 评审有感 +10 popular289 2026-05-18 15/750 2026-05-19 09:18 by jurkat.1640
[考博] 26/27申博自荐-锂/钠电池方向 5+3 狗头军师. 2026-05-15 4/200 2026-05-19 09:10 by moonboat
[基金申请] 别被青基扩招骗了!26年科研内卷才刚刚开始 +3 国自然国社科中 2026-05-14 4/200 2026-05-19 08:48 by archvillain
[考博] 2026博士还有哪些学校有名额 +7 小王求读研 2026-05-15 8/400 2026-05-19 08:27 by zhyzzh
[基金申请] 面上本子正文33页,违规吗?会被低分嘛? +8 1234567wang 2026-05-17 10/500 2026-05-18 18:52 by zzahkj
[硕博家园] 考博自荐 +5 科研狗111 2026-05-13 6/300 2026-05-18 11:22 by 糊糊涂涂好
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 ky2p12rrjj 2026-05-15 4/200 2026-05-17 19:47 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +6 l7k6xnh0yc 2026-05-14 7/350 2026-05-17 19:42 by Equinoxhua
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 5/250 2026-05-17 18:39 by Equinoxhua
[考博] 光量子物理方向 博士招生 1人(2026.09) +3 sandyworld 2026-05-15 4/200 2026-05-17 14:38 by sandyworld
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +6 l7k6xnh0yc 2026-05-14 6/300 2026-05-17 07:16 by 11n4dfd8yn
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +6 l7k6xnh0yc 2026-05-14 6/300 2026-05-17 07:11 by 11n4dfd8yn
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 cjf4bx70cj 2026-05-14 7/350 2026-05-17 06:55 by 11n4dfd8yn
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 3/150 2026-05-17 00:45 by ue3ir18jc3
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:35 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:25 by ue3ir18jc3
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
[有机交流] 如何实现卤原子转化 +3 BT20230424 2026-05-15 5/250 2026-05-16 16:20 by czyzsu
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 6/300 2026-05-16 11:29 by h3oerqvkv9
[教师之家] 教学课件你会给同学吗 +8 硕士研究生吗 2026-05-13 8/400 2026-05-14 22:23 by 常规沥青
信息提示
请填处理意见