matlab画分段常微分方程组相图

这个方程要解出来还要初始条件吧不妨设(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

在Matlab下输入：edit,然后将下面两行百分号之间的内容,复制进去,保存%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

告诉我邮箱,发你代码

注意,12.19t这种写法不对,应该写成12.19*t.结果如下>>symsS1S2t>>[S2,S1]=dsolve('DS2=12.19*t-24209-0.3*S2-0.01*S2','DS1=

用MATLABode45命令求解,对于分段函数可以用if语句,初值可以任意设定,只是不能全部为0,例如可以设为(0,1,-1,0).

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

functionhhX0=[17,195,1500,4,62,200];TSPAN=[04];[t,X]=ode45(@myfun,TSPAN,X0)plot(t,X)legend('x','y','

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

dsolve('m*D2y-a-b*f*Dx-sqrt(Dx*Dx+Dy*Dy)*Dy','m*D2x-b*f*Dy-c*sqrt(Dx*Dx+Dy*Dy)*Dx')正常情况使用上面的命令就可以解决问

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

先编一个M函数文件 function x=rr(t,u) %下面的参数值自己设置合适的,否则不一定有解,画不出图来. b1=2; d1=-3;&nbs

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

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

求系数矩阵的特征值,特征向量；特征向量求出后,构造基解矩阵,就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(

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

这个好像是常微分方程组把,先建立一个调用式子的函数：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