10楼: Originally posted by 吕春莹 at 2014-07-20 19:54:58
好吧，一开始我们老师也是教我们maple，后来就用matlab了，加上要用神经网络什么的，后来我看到另外一个老师用mathmatic的积分和求导等式子的形式漂亮多了，但是找不到安装包...

6楼: Originally posted by 月只蓝 at 2014-07-20 16:21:44
结果见附图。

附图1.png

7楼: Originally posted by 吕春莹 at 2014-07-20 19:08:10
我在3楼写的是用matlab算
syms a s y x
Y=dsolve('a*D2y-Dy/x-s*y=0','x')

解析解是
Y =
C3*x^((a + 1)/(2*a))*besselk(-(a + 1)/(2*a), -(s^(1/2)*x)/a^(1/2)) + (C2*x^((a + 1)/(2*a))*besseli(-(a + 1) ...

11楼: Originally posted by 月只蓝 at 2014-07-20 20:00:59
安装包小木虫里就有的吧。...

14楼: Originally posted by 孙小辉 at 2014-07-21 08:35:30
我刚才表达的不准确，师兄，为什么正常解法和laplace解法得到的解析解形式不一样，还有，第一种解析解的求解过程，师兄您有相关的资料吗...

6楼: Originally posted by 月只蓝 at 2014-07-20 16:21:44
结果见附图。

附图1.png

4楼: Originally posted by 吕春莹 at 2014-07-20 11:09:10
但是，反过来 ans=a*diff(Y,x,2)-diff(Y,x)/x-s*Y

ans =

a*(C3*x^((a + 1)/(2*a))*((s^(1/2)*((besselk(1 - (a + 1)/(2*a), -(s^(1/2)*x)/a^(1/2))*((a + 1)/(2*a) - 1))/x + (s^(1/2)*besselk(-(a + 1)/( ...

16楼: Originally posted by 孙小辉 at 2014-07-21 16:28:09
这位师兄，送佛送到西啊。得到Y（x)的解析解之后，用Maple软件进行：
（1）比如由Y(0)=1;Y(1)=2确定出Y（x)的C1,C2
（2）画出0<x<10之间的Y（x）曲线。
谢谢啦。。。...

16楼: Originally posted by 孙小辉 at 2014-07-21 16:28:09
这位师兄，送佛送到西啊。得到Y（x)的解析解之后，用Maple软件进行：
（1）比如由Y(0)=1;Y(1)=2确定出Y（x)的C1,C2
（2）画出0<x<10之间的Y（x）曲线。
谢谢啦。。。...

17楼: Originally posted by 孙小辉 at 2014-07-21 16:29:55
这位师姐，得到Y（x)的解析解之后，用Matlab进行：
（1）比如由Y(0)=1;Y(1)=2确定出Y（x)的C1,C2
（2）画出0<x<10之间的Y（x）曲线。
谢谢啦。。。
                                   求截图...