用matlab求微分方程组的数值解

[x,y]=dsolve('Dx=x-y-x*(x^2+y^2)','Dy=x+y-y*(x^2+y^2)','x(0)=2','y(0)=1')得到的结果是解析解没有找到.用数值解.在Matlab下

functionmyCalcclear;clc;closeall;opt=odeset('RelTol',1E-5,'AbsTol',[1E-61E-6]);[TY]=ode45(@getdY,[01

我改了你的一些内容,不知道有没有被我改错.在Matlab下输入：edit,然后将下面两行百分号之间的内容,复制进去,保存%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

functiondx=appollo(t,x)mu=1/82.45;mustar=1-mu;r1=sqrt((x(1)+mu)^2+x(3)^2);r2=sqrt((x(1)-mustar)^2+x(

这题的符号解很复杂,maple可以求出来我可以用matlab帮你求数值解,但是你要告诉我求解区间就是y(t),t的变化范围再问：t变化范围是【0，1】，y(t)的变化范围是【1，0】，猜测实际结果中y

dsolve('Dx1=-8/3*x1+x2*x3','Dx2=-10*x2+10*x3','Dx3=-x2*x1+28*x2-x3')ans=[x3(t)=0,x2(t)=0,x1(t)=C1*ex

这个问题很难,

clc;clear;[x,y]=dsolve('Dx=y','Dy=x+y')%test检验s1=diff(x)-ys11=simplify(s1)s2=diff(y)-x-ys22=simplify

先在一个空白的M文件里创建一个M函数,文件名要与函数名相同,例如：文件名：differential.m函数内容：function[dxdt]=differential(t,x)%x:二维列向量,x(1

你是要画图还是求a=?时可以使两式相等?另外你说两式相等是说x(t)=y(t)还是dx/dt=dy/dt?dsolve可以解微分方程组,ode45可以数值地解微分方程组,你需要数值解还是解析解?再问：

一般解非线性微分方程组可以用ode45,ode23等命令,编程时将p,v,th,x,y记为x(5),x(1),x(2),x(3),x(4)(方便方程组函数的编写就这么记号）下面便方程组函数,文件命名为

假设x(1)=x1,x(2)=x1',x(3)=x2;把高阶导数换算成低阶倒数；得到3个方程联立的方程组1：x(1)'=x(2)2：x(2)'=-2*x(2)-2*x(3)3：x(3)'=-x(2)-

dsolve('Dx=-y-z','Dy=x+0.1*y','Dz=2+z*(x-4)')Errorusing==>mapleError,wrongnumber(ortype)ofparameters

ode45的帮助例子那里有一阶常方程组的求解方法你可以把上面的方程组改写成一阶常方程组令x1=xx2=dx1/dtx3=dx2/dty1=yy2=dy1/dty3=dy2/dtx3=-u*x1/r3y

建立myfun.m文件,把以下代码复制进去function dx=myfun(t,x)dx=zeros(2,1);dx(1)=4*x(1)+x(2).^2;dx(2)=8*x(1)+2*x(

a=1;b=.1;[t,x]=ode45(@(t,x)[ab;-ba]*x,[0,10*pi],[1;1]);plot(x(:,1),x(:,2));xlabel('x1');ylabel('x2')

问的是微分方程!

[xy]=dsolve('D2x=(x*x+3*x)*Dy','D2y=-(x*x+3*x)*Dx','x(0)=0','y(0)=0','Dx(0)=a','Dy(0)=a','D2x(0)=a',

没解析解你肯定得用数值解.以下是步骤.第一步,在Matlab里写m-file,命名微分方程对应的函数：functionf=myfun(t,x)f=[-x(1)+x(2)^2;-2*x(2)+x(1)^

用例子吧dy/dt+2x=0dx/dt+3y=0>>[x,y]=dsolve('Dy+2*x=0','Dx+3*y=0')x=C1*exp(6^(1/2)*t)+C2*exp(-6^(1/2)*t)y