matlab 三阶微分方程的解法 含参数

请参考pdepe函数的用法

dsolve函数默认自变量是t,试试soln_1=dsolve('Dy=3*x^2','y(2)=0.5',’x‘)如果还是不行就不知道了可以了我运行过了>>soln_1=dsolve('Dy=3*x

可以用simulink做啊,很简单的,只需要托模块就可以了

这里面fo(t)是什么假设fo（t）=f[c1,c2]=dsolve('Dc1=-(k12+k13)*c1+V2/V1*k21*c+f/V1','Dc2=V1/V2*k12*c1-k12*c2','t

当然不是了,首先解齐议程对应的特征方程r^2-r+1=0r=(1±√3i)/2所以齐次通解是y=e^(1/2x)(C1cos√3x+C2sin√3x)特解可能观察得得y=a因此非齐次通解为y=e^(1

你可以用mathematic试试,不懂的可以通过帮助函数查找,但是势都有边界条件,比其编程加入,个人感觉还是手算方便的多.

看了看,运行了一下确实出问题 原因是在用ode数值求解时,x并不是1:0.01:3均匀分散的 解决方法：1.在画解析解和欧拉解时横轴用x的转置；在画数值解时横轴用x,已运行成功2.

symsxyy=sqrt(2*pi/x)*sin(x);y1=diff(y,'x');y2=diff(y,'x',2);h=x^2*y2+x*y1+(x^2-1/2)*y;h=simplify(h)结

我运行的>>symsaknNzz1>>z=dsolve('Dx=a*x*(N-x)','t')z1=dsolve('Dx=a*x*(n-x)','t')结果：z=N0N/(exp(-N*(C3+a*t

在matlab中默认的自变量是t,所以你的结果中在自变量是t,x被当做一个参数了如果你的自变量是x,你在写方程的时候把x写出t就行了!例如y=dsolve('Dy=t')y=t^2/2+C18

symsxa=5;b=a*x;w=dsolve('Dw=b','x')subs(w)%将a,b代入w的表达式.

这个你没给初值条件啊我假设是x0=[0,0];x(1)是A,x(2）是R程序如下首先定义这个方程functiondx=funct(t,x)dx=zeros(2,1);dx(1)=0.00875+7.5

matlab里面常使用龙格库塔方法求解常微分方程组,命令是ode45,还有其他一些函数,但是最常用的是ode45,lz可以help一下,很简单的,另外给你一个文档,讲的还是比较详细,希望可以帮到你ht

新的matlab版本好像不鼓励采用global了.你的全局变量有点多了,哈哈.简单例子:m=2;[t,y]=ode45(@(t,x)f1(t,x,m),[0,10],[2])functiondy=f1

不知道LZ的教材是哪一本,是不是高教出的王高雄,周之铭版的,常微的教材很多,你可以多参考几本：吉林大学出版社,王怀柔,伍卓群版的不错或者东北师大版的也不错但个人认为包络和奇解这个方向好像没什么深入的必

用幂级数法：设y=c0+c1x+c2x^2+...+cnx^n+...则y'=c1+2c2x+3c3x^2+...+ncnx^(n-1)y"=2c2+6c3x+12c4x^2+...+n(n-1)cn

因为你x=0时2/x是无穷大呀,然后y'又是0,然后(2/x)y'就是nan了,所以后面算的全是nan了.

见图

解肯定是存在且唯一的（在三个常数意义下唯一）,但是解析解不一定能用初等函数表达出来mathematica算是现在最好的符号计算软件之一,算出来了一个很复杂的结果：这个解析解就是写出来也没有太大意义（因

dy/dt+y+y²=2dy/dt=-y²-y+2dy/(-y²-y+2)=dt两边积分得1/3*㏑(y+2)/(y-1)=t+C