解常微分函数方程组 matlab

[x,y]=dsolve('Dx=x-y-x*(x^2+y^2)','Dy=x+y-y*(x^2+y^2)','x(0)=2','y(0)=1')得到的结果是解析解没有找到.用数值解.在Matlab下

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

在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).

把方程和初值条件都列出来这一大堆东西几个人看得懂?再问：公式（1）（2）就是方程啊，我写的那两个就是对应的方程。倒数第二张图片有实验数据。还需要什么参数吗？非常抱歉，我学的专业不太扎实，问题可能有点弱

怎么说呢...可能是因为所求系数(整型)太大了,也可能是你初始量设置的有问题.具体这会儿我没参考书,要不然就告诉你怎么改进.那句"Error,(incombine/radical)integertoo

这题的符号解很复杂,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

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

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

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

大概这样写：首先先建立2个editor第一个用来描述方程,建立M文件；第二个用来算方程第一个editor：%%%%设y(1)=y,y(2)=dyfunctiondy=fun(t,y)globalMBC

人家问的就是解析解,不知道上边几个在说什么你这个方程高版本的matlab可以解,只是算出来之后带积分项目：symsabQSTzZmcAPMeESS=dsolve('Dx=a*x','Dy=(b*Q*x

1、你贴出来的报错信息和代码对不上号：前面显示错误的那行代码和你贴出来的完全不一样；而后面的错误（Inputargument'u1'isundefined）也不可能是目前的代码所导致的

第二个方程右边的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))]

角度用sind、cosd,例如>>sind(30)ans=0.5000%真有意思,用不同ode,得到的结果不同～!图形见：http://hi.baidu.com/chemical%5Fengineer

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