matlab求微分方程初值问题的符号解,并与数值解进行比较
来源:学生作业帮助网 编辑:作业帮 时间:2024/05/14 12:25:19
x=dsolve('Dx=r*(1-x/xm)*x','x(0)=x0','t')x=xm/(1+exp(-r*t)*(xm-x0)/x0)
建议看看这个:利用有限差分和MATLAB矩阵运算直接求解二维泊松.和你要的原理是完全一样的.MATLAB程序...自己写吧.
你这个方程应该是在初条件v(0)=0时解的咯,没有初条件matlab解不出来.语句如下(我现在没条件在软件中试,如果不行你告诉我再改~)G=527.436*0.4536*9.8;F=470.327*0
请参考ode45函数的用法将所求值设为终值即可
%给你参考一下我当年做的一个类似的程序%cleardt=0.004;t=0:dt:5;a=-pi:0.2:pi;l=100;r=8;d=10;%导弹尺度参数:长度、半径、火焰尺度L=l;D=l/3;c
一维导热问题,应该采用有限差分法求解.具体方法随便找个数值计算方法书都有
ode45可以求解
matlab中导数是diff.help一下,就什么都有了
朋友,要根据初值积分对dM/dt积分求得M,才能求解M(t)=0或t(M)=0.solver(积分函数,积分时间,初值,设置)也是这样的数学方法.[时间,解]=solver(积分函数,积分时间,初值)
f(x,y)=x-y^2|f(x,y1)-f(x,y2)|<|y1^2-y2^2|
y=dsolve('D3y=D2y-Dy-y+t^2','y(0)=0','Dy(0)=1','D2y(0)=-1');tt=linspace(0,10,1001);%表示t从0到10fori=1:1
参考答案:少而好学,如日出之阳;壮而好学,如日中之光;志而好学,如炳烛之光.——刘向
y=dsolve('x*(D2y)+(1-2)*(Dy)+y=0','y(0)=0','Dy(0)=0','x')y=C6*x*besselj(2,2*x^(1/2))
前两天刚回答了一个类似的问题:再问:非常感谢你,那个问题也是我发的,但我处理二阶就不会了,失败了。听同学说没有解析解,想问下怎么写二阶形式的求解,只要ODE45的方法就行了,别的不麻烦您了再答:>>o
方程为:x*y''=k(1+y'^2)^0.5(y''表示y的二阶导数,y'表示y的一阶导数,y'^2表示y的一阶导数的平方,(1+y'^2)^0.5表示(1+y'^2)开根号)其中k=a/b初始条件
解一阶偏微分方程,这是一道难题哦.你连初始条件都没有啊,没法做再问:x=0u=5000K,x=Lu=2000K,L=10m,t=0u=5000K.
用matlab的dsolve命令求解失败,结果如下:>>symsf;a=dsolve('Dx=(sqrt((f-x)^2+y^2)-f+x)/y')Warning:Explicitsolutionco
x1=dsolve('(x^2+2*x*y-y^2)+(y^2+2*x*y-x^2)*Dy=0','y(1)=1','x')x1=1/2+1/2*(1+4*x-4*x^2)^(1/2)x2=dsolv
1(7)cosxsinydx+sinxcosydy=0,-cosxsinydx=sinxcosydy,-cosxdx/sinx=cosydy/siny,ln(siny)=-ln(sinx)+lnCsi
这是我在matlab里面运行得到的结果,a、b都是常数>>y=dsolve('D2y=-a^2*y+x','y(0)=b','Dy(pi/a)=0')y=cos(a*t)*(-x+b*a^2)/a^2