matlab求解微分方程的通解问题
题目
matlab求解微分方程的通解问题
解某微分方程:x^2*D2y+x*Dy+(x^2-1/2)*y=0,初值:y(pi/2)=2,Dy(pi/2)=-2/pi,书本标准答案为:
ans=2^(1/2)*pi*^(1/2)/x^(1/2)*sin(x);
但我运行:dsolve('x^2*D2y+x*Dy+(x^2-1/2)*y=0','y(pi/2)=2,Dy(pi/2)=-2/pi','x')
ans =
(2*bessely(-2^(1/2)/2,x)*(2^(1/2)*besselj(-2^(1/2)/2,pi/2) - besselj(-2^(1/2)/2,pi/2) + pi*besselj(1 - 2^(1/2)/2,pi/2)))/(pi*(besselj(1 - 2^(1/2)/2,pi/2)*bessely(-2^(1/2)/2,pi/2) - bessely(1 - 2^(1/2)/2,pi/2)*besselj(-2^(1/2)/2,pi/2))) - (2*besselj(-2^(1/2)/2,x)*(2^(1/2)*bessely(-2^(1/2)/2,pi/2) - bessely(-2^(1/2)/2,pi/2) + pi*bessely(1 - 2^(1/2)/2,pi/2)))/(pi*(besselj(1 - 2^(1/2)/2,pi/2)*bessely(-2^(1/2)/2,pi/2) - bessely(1 - 2^(1/2)/2,pi/2)*besselj(-2^(1/2)/2,pi/2)))
即使我simpify(ans),j结果还是一样,怎么回事啊?
答案
syms x y
y=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)
结果:
h =
-1/4/x*2^(1/2)*sin(x)*pi^(1/2)/(1/x)^(1/2)
可见不等于0,=y,说明原答案有问题
举一反三
已知函数f(x)=x,g(x)=alnx,a∈R.若曲线y=f(x)与曲线y=g(x)相交,且在交点处有相同的切线,求a的值和该切线方程.
我想写一篇关于奥巴马的演讲的文章,写哪一篇好呢?为什么好
最新试题
热门考点
- 瓶子高度为20厘米,下部呈圆柱形.瓶内装300克水,水面高12厘米.若将其倒立,则水面高14厘米.这个瓶子可装水多少克?
- 五彩缤纷怎么造句?
- "再别康桥'英文译稿是什么?拜托了各位 谢谢
- 小斌,英文直译的名字怎么说,就是比如“金”英文直译是“king”
- 从甲地到乙地,汽车第一天行驶全程的25%,第二天行驶全程的40%,这时超过中点175千米,甲乙两地的距离是
- 烽火戏诸侯的主要内容~要简短的!
- 如图所示的电路中,闭合开关S后,比较a、b、c、d四处电流大小,其中不正确的是( ) A.Ia>Ib B.Ia=Id C.Ia>Id D.Id>Ic
- 3,8,10,8用加减乘除和一个括号做出一个等式,结果为24.
- 英语翻译
- 要使浑浊的水变澄清且无色透明,可加入( )这道题选明矾还是活性炭