常微分方程组dx dt=4x y dy dt=-2x y 的

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

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

x(t)=e^(a*t);y(t)=b/(a+c)*[e^(a*t)-e^(-c*t)]z(t)=b*c/(a+c)*[e^(a*t)/a+e^(-c*t)/c]再问：你好，请问能提供一下具体的求救过

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

调用方法：[t,x]=ode45(@eqx,[00.2],[011.11.25])x的初始值写在一个矩阵里．另外你的函数定义时K少一维K(4),我给你加了一个：functionxdot=eqx(t,x

你让李雅普诺夫指数1恒为正.再使得它对t的导数为定号即可.你会画出它的轨线关于奇异点的逼近方式（可以远离.）事实上,要掌握的是图像,而不仅仅是具体方法.再问：谢谢你！我是新手，能不能说的具体点？谢谢啊

格式上的错误是,大括号有一对就可以了,不需要把初始条件和微分方程分割开来.然后,即使改了也是解不了的……要不是这个没有解析解,要不就是mathematica还不会解这个方程的解析解.（mathemat

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

clcYd=inline('[(YY(2)-1)*2*pi*50;(1-(1.47*1/2.82)*sin(YY(1)))/11.28]','t','YY');[t,YY]=ode45(Yd,0.3,

提供一点信息供楼主参考：1、使用dsolve求微分方程的解析解,只在一些比较简单的情况下能得到结果,很多情况下难以求出显式解.2、在MATLAB6.可以求出,用时约一分钟,但表达式过于复杂,根本没有实

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

dslove('D2y2+0.01-2y/5/(1+(1-x)^2)=0','x','y(0)=0','Dy(1)=0')就一个式子啊

eqns={x'[t]+y[t]==Cos[t],y'[t]==-x[t]+Sin[t]};sol=DSolve[eqns,{x,y},t]

分享到：\x0d收藏推荐0引言在科学应用中,常需要建立实际问题的数学模型,建立各种变量的常微分方程组及对其求解.微分方程是指方程中未知的是一个变量或几个变量的函数,并且在方程中不仅有函数本身而且有它们

你有初始值吗?解常微分方程方法有很多,比如：欧拉方法、龙格库塔法、多步法等等,龙格库塔方法用得比较多.数值解的原理基本是利用前面已知的点求后面的点,比如：dx可以近似的写为x(t0+h)-x(t0),

求系数矩阵的特征值,特征向量；特征向量求出后,构造基解矩阵,就ok.看书,照猫画虎做一遍,就会了.再问：能给出具体过程吗？再问：这个我真不会

第二个方程右边的dy1是什么意思?再问：右边的d应该是一个系数，可以是0.9再答：clcf=@(t,y)([y(1)*(1-y(1))-.5*y(1)*(y(2)+y(3)+y(4))/(1+2*y(

这个好像是常微分方程组把,先建立一个调用式子的函数：functionf=hanshu(y,x)f=[a(x(3)-x(2));b(x(1)-0.5*sqrt(x(2)-c));d(2.5-x(3))]

这个题目并不难,只是你的题目写太晦涩难懂.我也是看了好长时间才猜出来的,还不知道对不对呢.题意:xm(t),xp(t),r(t)是要求解的未知函数,而um(t),up(t)是已知的关于t的函数.程序如

用例子吧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