| 查看: 1020 | 回复: 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的取值，这是什么原因呢？解区间的开始的值会导致发散的情况，这个解区间的开始值应怎么取呢？ |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有65人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有10人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
找到一些相关的精华帖子，希望有用哦~
求高手解答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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
