matlab求解微分方程dy dx 2*y=2*x^2 2*x
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/10 02:44:22
请参考pdepe函数的用法
x=dsolve('Dx=r*(1-x/xm)*x','x(0)=x0','t')x=xm/(1+exp(-r*t)*(xm-x0)/x0)
结果:代码:clearallclcf=@(x,y)([y(2); 0.357*y(1)-0.1905*y(1)*y(2)]);[x,Y]=ode45(f,[0100]
朋友,要根据初值积分对dM/dt积分求得M,才能求解M(t)=0或t(M)=0.solver(积分函数,积分时间,初值,设置)也是这样的数学方法.[时间,解]=solver(积分函数,积分时间,初值)
用dsolve求不出解析解,用数值解法吧,最好能提供初始条件,这里假设t=0时,x0=1,y0=2,z0=3. functionhh [t,x]=ode23t(@fun,0:1,[123]) p
你可以用mathematic试试,不懂的可以通过帮助函数查找,但是势都有边界条件,比其编程加入,个人感觉还是手算方便的多.
symstv=dsolve('Dv=(190.708-90.64*v^2)/47.27','v(0)=0','t');t=0:0.00001:0.002;v=eval(v);plot(t,v)使用这样
看了看,运行了一下确实出问题 原因是在用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
functiontest()[t,y]=ode45(@func1,[0,1],[0;0;1;2;2;2]);figure(1);clf;plot(t,y);legend('x','y','z','dx
Matlab偏微分方程工具箱应用简介1.概述本文只给出该工具箱的函数列表,读者应先具备偏微分方程的基本知识,然后根据本文列出的函数查阅Matlab的帮助,便可掌握该工具箱的使用.2.偏微分方程算法函数
由于方程比较复杂,解析解不能用初等函数表示只是要获得图像的话,用数值计算的办法可能更方便fun=@(t,y)1.44*(10^9)*(1-y).*exp(-109170./(8.314*t));[T,
是可以求解的,参看pdepe函数
i是虚数单位?那个1/3γ(u*v)中*是什么符号
解一阶偏微分方程,这是一道难题哦.你连初始条件都没有啊,没法做再问:x=0u=5000K,x=Lu=2000K,L=10m,t=0u=5000K.
最常用的就是广义特征向量基础矩阵解方法.你要一个思路,我给一个2维情况的例子,其中特解x(t0)=x0的理解和如何使用都有,你看看是否够用.. Matlab下二维的例子:再问:嗯,这个不错,
改动:functionS=history(t)S=ones(3,1);改为S=ones(4,1);OVER!PS.多翻翻help!
y=dsolve('2000*Dy-(0.08-y*0.08)','y(0)=0','t')y=1-exp(-1/25000*t)即:C(t)=1-exp(-1/25000*t)