请教大家：非线性方程的求解，非常感谢！

cornfield · 发表于 2006-6-9 14:34

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

方程如下： dyc/dt=a-by 其中y是t的函数，y在初始时刻的值已知。c是常数。a是已知的关于t的序列（离散的）。
求出y在任一时刻t的表达式。

谢谢大家帮忙
新手发帖，如有不当，请多多包涵

[ 本帖最后由 lxq 于 2007-5-3 00:02 编辑 ]

cornfield · 发表于 2006-6-9 14:36

不好意思的是，上午发帖的位置不对，被删了。 我很想知道这个方程怎样求解，最好是解析解，数值解也好，请大家帮忙，谢了

feifeifool · 发表于 2006-6-9 17:12

可以方程可以改成 cy^(c-1)dy/dt=a-by,可以考虑一下ode算法求解如果y^(c-1)的倒数存在，可以写成dy/dt=(a-by)*y^(1-c)/c.

cornfield · 发表于 2006-6-9 23:20

感谢feifei的建议，变成你说的形式后再怎样求解呢 好像还是没办法

feifeifool · 发表于 2006-6-10 09:19

Matlab里面对于Mdy/dt=f(y,t)的方程都可以用ode解决，常微分的初值问题、边值问题，都可以。

cornfield · 发表于 2006-6-13 08:52

非常感谢你的解答，我对Matlab不熟悉，matlab里有解这个方程的源代码吗？我想把解这个方程的代码嵌套到我的程序中，因为a虽然已知，但不同的时间段里值是不同的。 有源代码的话，能不能提供一下，非常感谢

cornfield · 发表于 2006-6-13 18:53

谁有源代码，能不能提供一下啊，十分感激

feifeifool · 发表于 2006-6-14 11:09

a这个可以用if语句实现，你先看看相关的书吧！ 我现在忙着准备毕业答辩，没那么多时间，还有不到两周。

cornfield · 发表于 2006-6-14 16:55

恩，我在看书，想用五阶变步长的runge-kutta法求解，据说在《计算机数值计算及程序设计》这本书中有描述，不知道哪里有这本书 非常感谢你，也祝愿你顺利毕业。

cornfield · 发表于 2006-6-14 16:58

我的qq号是37854947，希望大家能给我点帮助，谢谢

feifeifool · 发表于 2006-6-14 17:19

不是让你看这个，是看matlab中odefile的用法。不需要自己进行龙个库塔的编程，：）

xuebx · 发表于 2006-6-14 19:19

看书是第一位的，所以要重视基础知识的积累

bainhome · 发表于 2006-6-14 19:48

In this level,maybe you can get something useful in here. <a href="http://www.simwe.com/forum/thread-541373-1-1.html" target="_blank" >http://www.simwe.com/forum/thread-541373-1-1.html</A> And I do believe there is no such a wall between two forums,so I give the link directly instead of moving that stuff from one place into another.

账号		自动登录	找回密码
密码			我要加入

[编程技巧] 请教大家：非线性方程的求解，非常感谢！

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（cornfield）感谢feifei的建议，变成你说的形式...

回复：（cornfield）非常感谢你的解答，我对Matlab不熟...

回复：（cornfield）恩，我在看书，想用五阶变步长的r...