用dsolve求不出解析解,用数值解法吧,最好能提供初始条件,这里假设t=0时,x0=1,y0=2,z0=3. function hh [t,x]=ode23t(@fun,0:1,[1 2 3])。
syms A B C x y;>> dsolve('Dy=A-B*y-C*x*y','x')ans =-1/2*A*erf(1/2*(-2*C)^(1/2)。
微分方程的话可以用dsolve函数,格式可以总结为dsolve('equation','condition',自变量)的形式。自变量默认为t。这种方法解微分方程还好。如。
这个方程相当于e*sin(t*w)*(b*w^4 - a*w^6 - c*w^2 + d)=0即e*sin(t*w)*(b*w^4 - a*w^6 - c*w^2 + d)则有sin(t*w)=0或b。
syms x y t;S=dsolve('Dx=-10*x+10*y','Dy=28*x-y-x*z','Dz=-8/3*z+x*y','x(0。
【1】用解析法r1=0.01,r2=0.05,K1=1000,K2=1000,m=9e-6,n=4e-5[N1,N2]=dsolve('DN1=(r1*(1-N1/K1)-m*N2)*N1。
我只会用常规分离变量法,龙格库塔法是什么哦~没学到额~
你这个是二阶微分方程 应该加上二阶的特解才可能解出来另外dsolve只能达到解析解我估计你这个函数是没有解析解的,只有数值解所以解不出来 你这个是。
dsolve('Dy-2*y/(x+1)=(x+1)^5/2')ans =(C2*exp((2*t)/(x + 1)))/4 - (3*x)/2 - (15*x^2)/4 - 5*x^。
==>sinydy/cosy=sinxdx/cosx ==>d(cosy)/cosy=d(cosx)/cosx ==>ln│cosy│=ln│cosx│+ln│C│ (C是积分常数)。
回顶部 |