matlab2015使用dsolve错误,用matlab解常微分方程,调用dsolve是老出错 y=dsolve,用matlab解微分方程y=dsolve(’D2y+Dy2...

导航:网站首页 >

用matlab解常微分方程,调用dsolve是老出错 y=dsolve,用matlab解微分方程y=dsolve(’D2y+Dy2

用matlab解常微分方程,调用dsolve是老出错 y=dsolve('Dy=x*sin(x)/cos(y...,用matlab解微分方程y=dsolve('D2y+Dy-2y=2x','y(0)=0,Dy(0)=1','x')为什...

匿名网友:

大哥,你是中英文不分吗?你的程序没错,错的是你把里面很多的括号打成了中文dsolve('200000*D2x+4.5*(Dx)^2=0','200000*D2y=-2000000+4.5*(Dy)^2+30*((Dy)^2+(Dx)^2)','Dx(0)=240','x(0)=0','Dy(0)=0','y(0)=10668','t')不过貌似解不了

匿名网友:

> syms x y a;>> dsolve('Dy=-2*x*y/(x^2+2*a)','x')ans =C1/(x^2+2*a)我这里运行没错,你检查一下你的命令中括号等符号是否有问题,是否有输成全角的情况。

另外注意:语句中要加上‘x',因为matlab默认自变量是t,你这样的式子算出来了也不对。

匿名网友:

>> syms x y0>> y=dsolve('Dy=y+1/y','y(0)=y0','x')y = (-1+exp(2*x)*(1+y0^2))^(1/2) -(-1+exp(2*x)*(1+y0^2))^(1/2)>> help dsolve DSOLVE Symbolic solution of ordinary differential equations. DSOLVE('eqn1','eqn2', ...) accepts symbolic equations representing ordinary differential equations and initial conditions. Several equations or initial conditions may be grouped together, separated by commas, in a single input argument. By default, the independent variable is 't'. The independent variable may be changed from 't' to some other symbolic variable by including that variable as the last input argument. The letter 'D' denotes differentiation with respect to the independent variable, i.e. usually d/dt. A "D" followed by a digit denotes repeated differentiation; e.g., D2 is d^2/dt^2. Any characters immediately following these differentiation operators are taken to be the dependent variables; e.g., D3y denotes the third derivative of y(t). Note that the names of symbolic variables should not contain the letter "D". Initial conditions are specified by equations like 'y(a)=b' or 'Dy(a) = b' where y is one of the dependent variables and a and b are constants. If the number of initial conditions given is less than the number of dependent variables, the resulting solutions will obtain arbitrary constants, C1, C2, etc. Three different types of output are possible. For one equation and one output, the resulting solution is returned, with multiple solutions to a nonlinear equation in a symbolic vector. For several equations and an equal number of outputs, the results are sorted in lexicographic order and assigned to the outputs. For several equations and a single output, a structure containing the solutions is returned. If no closed-form (explicit) solution is found, an implicit solution is attempted. When an implicit solution is returned, a warning is given. If neither an explicit nor implicit solution can be computed, then a warning is given and the empty sym is returned. In some cases involving nonlinear equations, the output will be an equivalent lower order differential equation or an integral. Examples: dsolve('Dx = -a*x') returns ans = C1*exp(-a*t) x = dsolve('Dx = -a*x','x(0) = 1','s') returns x = exp(-a*s) y = dsolve('(Dy)^2 + y^2 = 1','y(0) = 0') returns y = sin(t) -sin(t) S = dsolve('Df = f + g','Dg = -f + g','f(0) = 1','g(0) = 2') returns a structure S with fields S.f = exp(t)*cos(t)+2*exp(t)*sin(t) S.g = -exp(t)*sin(t)+2*exp(t)*cos(t) dsolve('Dy = y^2*(1-y^2)') returns Warning:Explicit solution could not be found; implicit solution returned. ans = t+1/2*log(y-1)-1/2*log(y+1)+1/y+C1=0 dsolve('Df = f + sin(t)', 'f(pi/2) = 0') dsolve('D2y = -a^2*y', 'y(0) = 1, Dy(pi/a) = 0') S = dsolve('Dx = y', 'Dy = -x', 'x(0)=0', 'y(0)=1') S = dsolve('Du=v, Dv=w, Dw=-u','u(0)=0, v(0)=0, w(0)=1') w = dsolve('D3w = -w','w(0)=1, Dw(0)=0, D2w(0)=0') y = dsolve('D2y = sin(y)'); pretty(y) See also solve, subs.Reference page in Help browser doc dsolve

匿名网友:

>>clear>>syms a b c d e;>>y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')y =(2*a*e + b*d + d*(b^2 - 4*a*c)^(1/2))/(2*exp((t*(b - (b^2 - 4*a*c)^(1/2)))/(2*a))*(b^2 - 4*a*c)^(1/2)) - (2*a*e + b*d - d*(b^2 - 4*a*c)^(1/2))/(2*exp((t*(b + (b^2 - 4*a*c)^(1/2)))/(2*a))*(b^2 - 4*a*c)^(1/2))>>a=1;b=1;c=1;d=1;e=1;%若常数已知>>t=1;%desolve中没指定自变量x,这里默认为t,带入任意一个x值>>y=eval(y) y =1.1932 + 0.0000i这是解方程的方法,画图也差不多,我在命令窗写的改麻烦,下面是我在editor写的。

syms a b c d e;y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');x=200:0.1:400;%取值范围,步长a=1;b=1;c=1;d=1;e=1;%带入参数y=eval(y);plot(x,y)你补充问题的方程解得也是一个空的,并报错,或许就是无解的,我也很疑问。

问题推荐

  • 0
    点赞
  • 0
    评论
  • 0
    收藏
  • 一键三连
    一键三连
  • 扫一扫,分享海报

表情包
插入表情
评论将由博主筛选后显示,对所有人可见 | 还能输入1000个字符
相关推荐
©️2020 CSDN 皮肤主题: 数字20 设计师:CSDN官方博客 返回首页
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、C币套餐、付费专栏及课程。

余额充值