| 查看: 966 | 回复: 0 | |||
prince1000银虫 (小有名气)
|
[求助]
再次恳求高手帮忙ode45解微分方程组的问题
|
|
function dy=diff(r,y) dy=zeros(3,1) dy(1)=-((371.04*cos(123.68*r)/(20+sin(123.68*r))+1/r+140/(20+sin(123.68*r))^2/r)/(1+140/(20+sin(123.68*r))^2))*y(1)-((137.2848*cos(123.68*r)/(20+sin(123.68*r))/r-1/r^2)/(1+140/(20+sin(123.68*r))^2))*y(2)+(140/(20+sin(123.68*r))^2/r^2/(1+140/(20+sin(123.68*r))^2))*y(2) dy(2)=y(1) 用ode45（'diff',[r0,r],[y1,y2]）命令解以上一阶常微分方程组，这个解区间[r0,r]应该为[0 100],可是零是要出错的，因为方程中有1/r,那这个r0怎么设定呢？是不是取一个很小的数，如1.D-50，可是我试了用[1.D-50 100]和[1.D-40 100]，算出来的y（1）差好几个数量级，也就是说y（1）强烈的依赖于r0的取值，这是什么原因呢？解区间的开始的值会导致发散的情况，这个解区间的开始值应怎么取呢？ |
回复此楼» 猜你喜欢
职称评审没过，求安慰
已经有15人回复
-
垃圾破二本职称评审标准
已经有11人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有16人回复
-
谈谈两天一夜的“延安行”
已经有14人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
EST投稿状态问题
已经有4人回复
-
聘U V热熔胶研究人员
已经有10人回复
-
求助文献
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
-
三无产品还有机会吗
已经有6人回复
找到一些相关的精华帖子，希望有用哦~
求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复- 一阶常微分方程组，matlab进行数值求解
已经有3人回复
- 如何用matlab求解矩阵系数的二阶微分方程
已经有13人回复
- matlab解微分方程
已经有10人回复
- lsqnonlin函数拟合微分方程组参数拟合问题
已经有10人回复
- matlab解微分方程组
已经有15人回复
- 一个微分方程组，求解
已经有18人回复
- 求助利用matlab的ODE45求解微分方程
已经有16人回复
- 【紧急求助】 带有微分平方项的ODE求解
已经有11人回复
- 用MATLAB ode45求解2阶微分方程
已经有4人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】急求各位大侠，关于用ode45解微分方程
已经有8人回复
- 【求助】请教非齐次常微分方程组的解析解法
已经有4人回复
- 【求助】ode45数值求解常微分方程总是错，找不出原因
已经有3人回复
- 【求助完毕】用ode逆向求解微分方程问题
已经有4人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】常微分方程组求解中系数与某变量值关联的问题
已经有12人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】解常微分方程组
已经有7人回复
- 【求助】一个时变微分方程组的matlab实现
已经有8人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人