用Matlab求解一阶线性微分方程y 2xy=xe^(-x)^2

1定义函数：function y=fun(t,x)y=zeros(3,1);x1=x(1);x2=x(2);x3=x(3);y(1)=x1*(1-x1/150000-0.5*x2/30000

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

这个问题很难,

这个方程应该可以用特征线法去求解,但是还要给出u,v的边界条件才能给出具体的表达式.建议你看本数学物理方程的书都有讲特征线法的.英文书你可以看Evans的PDE再问：能帮我解一下吗？本人不是数学系的，

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

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

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

就是先把方程组的系数写成矩阵的形式再解特征根~比如说方程组dx/dt=3x-5ydy/dt=5x+3y那么该矩阵A就是[3-5](不会打大的括号,凑合看吧)[53]下面算det(A-λE)=|3-λ-

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

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

解法一：（全微分法）∵y'-2y/x=x^3==>xy'-2y=x^4==>xdy-2ydx=x^4dx==>x²dy-2xydx=x^5dx==>x²dy-yd(x²)

i是虚数单位?那个1/3γ（u*v）中*是什么符号

我来解答你的关键问题：将之前主程序计算出来的系数用globalcof定义,再在ode45的子函数里同样定义globalcof,这样就能把值计算出来并且传递进去了．

不太理解你说的公式法是什么……是指的一阶微分方程的y'+P(x)y=Q(x)么?结果当然是一样的再问：第二行那个就是公式法不是可以直接带入求出来的么再答：我刚修改了……你看过程发的图片，你在电

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)

把第三个式子也微分.>>[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(

再问：这个我明白，但是积分符号怎么处理的？再答：这两步之间哪涉及积分符号？你是问哪个积分符号？这两步间只有一个积分符号，照抄下来，没变化啊