用Mathematica实现经典龙格库塔法(RK4)的教程没人写过?
mathematica吧
全部回复
仅看楼主
吧务
level 15
xzcyr 楼主
用“龙格库塔 mathematica”在百度翻了10页,只看到2个基于Mathematica的经典龙格库塔法实现。一个是一篇论文,用的是For循环。(顺便摘要里还说经典龙格库塔是常微分方程近似解法里“最常用的一种”,这句话即使放在论文发表的06年也是大有问题。)还有一个是我在知道上的回答,引用的是Stackexchange上使用MapThread和Rule的代码,这个答案如今看来写得也不好。
自编RK4也算是个常见的作业题,如果目前中文网络上的相关实现真的只有上面两个,那是不是有必要写一个相关的教程?(帖子《NDSolve解不出来的东西,你也别指望经典龙格库塔能解!》( tieba.baidu.com/p/8054474975 )里面提及的primaryRK`程序包不太适合入门。)以此为切入点,可写的东西还挺多的其实。谁有兴趣写写吗?(我?我现在累了……)
大家手里要是有现成教程的话,也不妨分享一下。
2022年11月05日 08点11分 1
吧务
level 10
主要是tutorial/NDSolveExplicitRungeKutta不够“入门”
不过最基本版的就NestWhile套公式就完事了吧
2022年11月05日 11点11分 2
我觉得至少可以考虑把算术函数的Listable属性放进来讲,也算非常好用也常用的功能了。
2022年11月05日 11点11分
此外方程右侧显含时的话,又会略复杂一点,甚至可以上升到哲学之争:在含时的情形下,用Nest/Fold系函数真的比Do更简洁吗?
2022年11月05日 12点11分
@xzcyr Nest的目的主要在于直接用计算的栈而不用符号表。
2022年11月06日 03点11分
@asdasd1dsadsa ……啥栈?
2022年12月03日 03点12分
吧务
level 15
xzcyr 楼主
顶楼所引帖子《NDSolve解不出来的东西,你也别指望经典龙格库塔能解!》的链接更新为
tieba.baidu.com/p/8443197888 (先前那版被我误操作删除了——百度贴吧现在居然允许直接删除精华帖……)
2023年06月03日 03点06分 4
level 3
狠狠赞了,RK4解定态薛定谔方程就跟写不出来作文一样难受
2024年03月17日 08点03分 5
1