当前位置：
&MATLAB 求解常微分方程组
MATLAB 求解常微分方程组
作者 MJ2016
本人MATLAB菜鸟，想用MATLAB求解一个常微分方程组，如下是编写的M文件，但是求解的曲线图不收敛，不知问题出在什么地方，急求各路大神帮忙解决！初值改变后经常半天没有反应，是死循环了还是要运算很久没反应啊，只好终止运算！
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function dy = vdp3000(t,y)
dy = zeros(6,1);& & % a column vector
dy(1) = y(2);
dy(2) = (1/0.00207)*((3.14e-6)*y(6)-8.15*((4.9e-3)+y(1))-(3.028e-3)*y(1)*y(6));
dy(3) = y(4);
dy(4) = (1/0.031)*((1.206e-3)*y(5)-3.02*((8.02e-3)+y(3))-(1.256e-3)*y(6)-0.0173*y(3)*y(5));
dy(5) = (0.6e7)*((6.15e-4)-0.135*((((0.5e-3)^8)/(5e-3))^(1/3))*((y(5)-y(6))^0.667)-0.0126*y(3)*((2*y(5)/)-(1.206e-3)*y(4));
dy(6) = (0.662e8)*(0.135*((((0.5e-3)^8)/(5e-3))^(1/3))*((y(5)-y(6))^0.667)-(1.611e-3)*y(1)*((2*y(6)/)-(3.14e-6)*y(2)+(1.256e-3)*y(4));
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[t,Y] = ode15s(@vdp],[0 0 0 0 40 0]);
plot(t,Y(:,3))
xlabel('t');
ylabel('y');
微分方程拟合用1stOpt吧，简单方便
具体怎么用呢，能给个例子吗？
从你的方程看，应该是有二阶微分，然后被你替换了，所以你最好把你的原始方程发上来，如果你用
dt/dz = d2y/dz2
那么你需要确定的知道 t 的初值，也就是 dy/dz 的初值。而且有的时候这种方程会有一些约束，最后形成的不一定完全是微分方程，还可能带有代数方程。
所以还是拿出原方程才好判断你哪里写的有问题。
我完全复制你的程序做出来的图是这样的，稍微改动了一下，采用了ode45，tspan参数应该是[0:10],不知道能不能帮到你。。。。
QQ截图39.png
function sol = testyours()
sol = ode23s(@vdp],[0 0 0 0 40 0]);
t = sol.x;
Y = sol.y;
plot(t,Y(3,)
xlabel('t');
ylabel('y');
function dy = vdp3000(t,y)
dy = zeros(6,1);& & % a column vector
dy(1) = y(2);
dy(2) = (1/0.00207)*((3.14e-6)*y(6)-8.15*((4.9e-3)+y(1))-(3.028e-3)*y(1)*y(6));
dy(3) = y(4);
dy(4) = (1/0.031)*((1.206e-3)*y(5)-3.02*((8.02e-3)+y(3))-(1.256e-3)*y(6)-0.0173*y(3)*y(5));
dy(5) = (0.6e7)*((6.15e-4)-0.135*((((0.5e-3)^8)/(5e-3))^(1/3))*((y(5)-y(6))^0.667)-0.0126*y(3)*((2*y(5)/)-(1.206e-3)*y(4));
dy(6) = (0.662e8)*(0.135*((((0.5e-3)^8)/(5e-3))^(1/3))*((y(5)-y(6))^0.667)-(1.611e-3)*y(1)*((2*y(6)/)-(3.14e-6)*y(2)+(1.256e-3)*y(4));
running results:
Warning: Imaginary parts of complex X and/or Y arguments ignored
图片是一个那样的东西，我不知道怎么发图，你把上面代码复制把，如果出错那就是你matlab的问题了，ode45，15s，23s都能得到结果，
24小时热帖
下载小木虫APP
与700万科研达人随时交流常微分方程的求解_百度文库
赠送免券下载特权
10W篇文档免费专享
部分付费文档8折起
每天抽奖多种福利
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
常微分方程的求解
阅读已结束，下载本文需要
想免费下载更多文档？
定制HR最喜欢的简历
下载文档到电脑，同时保存到云知识，更方便管理
加入VIP
还剩7页未读，
定制HR最喜欢的简历
你可能喜欢求解常微分方程 怎么做啊_百度知道
求解常微分方程 怎么做啊
求解常微分方程怎么做啊
我有更好的答案
答案是这个图
只不过现在因变量变成P了，理解下
我是按着这个公式做的啊。。可是后面好像没法积分了
不会的，亲
你就按照公式来，把因变量改过来嘛
采纳率：64%
来自团队：
为您推荐：
其他类似问题
常微分方程的相关知识
换一换
回答问题，赢新手礼包
个人、企业类
违法有害信息,请在下方选择后提交
色情、暴力
我们会通过消息、邮箱等方式尽快将举报结果通知您。【图文】第十章 常微分方程(组)的求解-0541_百度文库
赠送免券下载特权
10W篇文档免费专享
部分付费文档8折起
每天抽奖多种福利
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
第十章 常微分方程(组)的求解-0541
阅读已结束，下载本文到电脑
想免费下载本文？
登录百度文库，专享文档复制特权，积分每天免费拿！
你可能喜欢当前位置：
&怎样的常微分方程的解是正弦或者是余弦， 急， 谢谢！
怎样的常微分方程的解是正弦或者是余弦， 急， 谢谢！
请构造个方程，&&怎样的常微分方程的解是正弦或者是余弦， 急， 谢谢！
如果定义域是变动的更好， 特别感谢
二阶特征值问题的解都是正弦或者余弦形式的；另外，定义域是变动的我无法理解，因为在计算的时候不需要固定定义域，只需要是定义域是有界的即可（端点的位置影响的是特征值序列）。
二阶常系数常微分方程的特征方程的共轭复根对应的解即为正弦或余弦函数。如y''+y=0 ,特征方程r^2+1=0,r=±i,其解为：
y=A*Cosx+B*Sinx
d2y/dx2=﹣kx/m，这是标准的简谐振动的表达式结果是y=Acos(ωx+φ)，一般形如d2y/dx2=﹣ax（a为常数）的解析解都是三角函数，一定要有负号，像d2y/dx2=bx的解析解就可以是exp(√b*x)，其他也还有不少埃
求两次微分后函数形式不变，所以是二阶微分的方程，更高阶也行
Y''+a^2*Y=0& & (a为实数)，该形式的方程的通解为Y=A*sin(a*x)+B*cos(a*x)
常系数线性微分方程如果特征根全为虚数
24小时热帖
下载小木虫APP
与700万科研达人随时交流}