用Ode15s函数编写化学反应模式
放牛娃,2007年3月26日
去年使用烟雾箱模拟了几种物质在光化学烟雾中的臭氧活性，实验完成后打算独立编写一个反应模式对实验进行数值模拟。初步打算使用C语言编写，总共涉及一百多个反应方程，最艰难的部分是解决算法问题和编写60多个严重刚性的微分方程，要编写这样一个程序实在是一个浩大的工程，所以迟迟下不了决心。
在高教看了很多C的算法，正在一筹莫展的时候发现Matlab提供解决常微分方程问题的程序包Ode45，Ode15s。Ode45适合解决非刚性问题，Ode15s适合解决刚性问题（s为stiff刚性之意）。柳暗花明，根据情况选择了Ode15s，整个问题的重心由编写合适的算法转到编写方程组，工作量大大降低了，整个程序简洁友好，一个很复杂的问题瞬间简单化了，将我从浩瀚的代码中解放出来。 如果想学习Matlab我推荐大家看《Matlab 数值计算》这本书，由Matlab创始人之一Cleve B. Moler 编写，深入浅出有很多活泼有趣的例子，如密码学、Google网页分级和大气科学等。
分类: 网络资源 | 评论:0 | 浏览:
相关文章:
