matlab求解偏微分方程组

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

这个方程要解出来还要初始条件吧不妨设(x(0),y(0))=(A,B);还有用ode45只能得到数值解,不能给出解析式表达那么先建立一个函数文件test_fun.mfunctiondx=tsst_fu

常微分方程的解有标准格式,需要根据议程形式假设其通解,比如假设y=c1x+c2x'+c3,代入方程解通解,再加边界条件x(0)=y(0)=1解特解.

dsolve命令,具体用法看Helpdsovle哥们,加油哦,用这个命令能搞定的,没问题的

functionhhX0=[17,195,1500,4,62,200];TSPAN=[04];[t,X]=ode45(@myfun,TSPAN,X0)plot(t,X)legend('x','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

这个问题很难,

lumda=1;mu=0.3;fun=@(t,is)[lumda*is(2)*is(1)-mu*is(1);-lumda*is(2)*is(1)];[T,Y]=ode45(fun,[0100],[0.

解方程很简单,调用ode45即可关键是求运动时间,假设z方向坐标变为0即为落至地面运动结束,可求得运动时间的近似值及落地点坐标代码如下：clearallclcc=0.07;u0=25;a=pi/6;f

lz需要给出y_0的初值~这个方程满足李普希兹条件,因此,解存在唯一并且可以唯一延拓到边界,应用lax等价定理可以知道,向前欧拉法具有二阶的收敛速度……所以可以尝试用向前欧拉法编写：我刚学matlab

dsolve()函数,请自己在matlab里边搜索用法,很简单的,看看下边的例子就会了~

clearall;clc;load15.mat;t=celiang(:,1);u1=celiang(:,2);symsu2i5iixu5%定义u1,L1,L2,R1,C1,C3,ii等vU1=celi

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

你这个问题问了两次啊,上次给你回另一个问题了,这边来收分,不知道你改程序的参数没有?

试了一下,发现在一定的边界条件下,可以如：①设定每个变量的初值边界条件M文件functionres=ivbc(ya,yb)q=5;res=[ya(1)-qya(2)-2ya(3)-3ya(4)];②设

1、后两个方程都是dH?2、最后一个dH的表达式有误（括号不匹配）.3、st的初值?

clearclc[t,x]=ode45(@(t,x)[-x(1)^2-x(2);x(1)-x(2)^3],[05],[10.5]);plot(t,x(:,1),'r',t,x(:,2)

使用符号运算工具.首先定义变量符号symssxt;使用字符串定义公式eq1='diff(x,t)=0.11*s*x/(0.06*x+s)';eq2='diff(s,t)=(-1)*0.029*x-0.

把第三个式子也微分.>>[x,y,z]=dsolve('Dx=(-2)*A*x+u*y','Dy=2*A*x+2*u*z-(A+u)*y','Dx+Dy+Dz=0','x(0)=1,y(0)=0,z(

[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',