maple常微分方程组,作图作业

这个方程要解出来还要初始条件吧不妨设(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]再问：你好，请问能提供一下具体的求救过

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(

solve函数,查帮助有例子

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

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

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

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

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

代码供参考：n=1.5;%x1=φ,x2=Rdx=@(z,x)[1/(x(2)*(n+(n+1)*(2*cos(x(1))+cos(x(1))*cos(x(1)))));-1/tan(x(1))];[

就是先把方程组的系数写成矩阵的形式再解特征根~比如说方程组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')就一个式子啊

分享到：\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