关于二阶二阶线性常微分方程程第10题 求大神解答

二阶常微分方程解_百度文库
两大类热门资源免费畅读
续费一年阅读会员,立省24元!
二阶常微分方程解
上传于||暂无简介
阅读已结束,如果下载本文需要使用
想免费下载本文?
下载文档到电脑,查找使用更方便
还剩22页未读,继续阅读
你可能喜欢百度文库-信息提示
两大类热门资源免费畅读
续费一年阅读会员,立省24元!
对不起,该文档已被删除,无法查看
4秒后,自动返回首页小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
师兄师姐们,帮小弟看看这个 Laplace变换后的二阶常微分方程怎么解啊
师兄师姐们,看看这个变系数的二阶常微分方程怎么求解析解啊
但是,反过来 ans=a*diff(Y,x,2)-diff(Y,x)/x-s*Y
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)/(2*a), -(s^(1/2)*x)/a^(1/2)))/a^(1/2)))/a^(1/2) + (besselk(-(a + 1)/(2*a), -(s^(1/2)*x)/a^(1/2))*(a + 1))/(2*a*x^2) - ((a + 1)*((s^(1/2)*besselk(1 - (a + 1)/(2*a), -(s^(1/2)*x)/a^(1/2)))/a^(1/2) - (besselk(-(a + 1)/(2*a), -(s^(1/2)*x)/a^(1/2))*(a + 1))/(2*a*x)))/(2*a*x)) + (C2*x^((a + 1)/(2*a))*((s^(1/2)*((besseli(1 - (a + 1)/(2*a), (s^(1/2)*x)/a^(1/2))*((a + 1)/(2*a) - 1))/x + (s^(1/2)*besseli(-(a + 1)/(2*a), (s^(1/2)*x)/a^(1/2)))/a^(1/2)))/a^(1/2) + (besseli(-(a + 1)/(2*a), (s^(1/2)*x)/a^(1/2))*(a + 1))/(2*a*x^2) - ((a + 1)*((s^(1/2)*besseli(1 - (a + 1)/(2*a), (s^(1/2)*x)/a^(1/2)))/a^(1/2) - (besseli(-(a + 1)/(2*a), (s^(1/2)*x)/a^(1/2))*(a + 1))/(2*a*x)))/(2*a*x)))/(-1)^((a + 1)/(2*a)) + (C3*x^((a + 1)/(2*a) - 1)*(a + 1)*((s^(1/2)*besselk(1 - (a + 1)/(2*a), -(s^(1/2)*x)/a^(1/2)))/a^(1/2) - (besselk(-(a + 1)/(2*a), -(s^(1/2)*x)/a^(1/2))*(a + 1))/(2*a*x)))/a + (C3*x^((a + 1)/(2*a) - 2)*((a + 1)/(2*a) - 1)*besselk(-(a + 1)/(2*a), -(s^(1/2)*x)/a^(1/2))*(a + 1))/(2*a) + (C2*x^((a + 1)/(2*a) - 1)*(a + 1)*((s^(1/2)*besseli(1 - (a + 1)/(2*a), (s^(1/2)*x)/a^(1/2)))/a^(1/2) - (besseli(-(a + 1)/(2*a), (s^(1/2)*x)/a^(1/2))*(a + 1))/(2*a*x)))/((-1)^((a + 1)/(2*a))*a) + (C2*x^((a + 1)/(2*a) - 2)*((a + 1)/(2*a) - 1)*besseli(-(a + 1)/(2*a), (s^(1/2)*x)/a^(1/2))*(a + 1))/(2*(-1)^((a + 1)/(2*a))*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) - (besselk(-(a + 1)/(2*a), -(s^(1/2)*x)/a^(1/2))*(a + 1))/(2*a*x)) + (C2*x^((a + 1)/(2*a))*((s^(1/2)*besseli(1 - (a + 1)/(2*a), (s^(1/2)*x)/a^(1/2)))/a^(1/2) - (besseli(-(a + 1)/(2*a), (s^(1/2)*x)/a^(1/2))*(a + 1))/(2*a*x)))/(-1)^((a + 1)/(2*a)) + (C3*x^((a + 1)/(2*a) - 1)*besselk(-(a + 1)/(2*a), -(s^(1/2)*x)/a^(1/2))*(a + 1))/(2*a) + (C2*x^((a + 1)/(2*a) - 1)*besseli(-(a + 1)/(2*a), (s^(1/2)*x)/a^(1/2))*(a + 1))/(2*(-1)^((a + 1)/(2*a))*a))/x - s*(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)/(2*a), (s^(1/2)*x)/a^(1/2)))/(-1)^((a + 1)/(2*a)))
这回我不明白了
您好,能不能发个截图上去,后面这个是它的解析解吗?
我在3楼写的是用matlab算
syms a s y x
Y=dsolve('a*D2y-Dy/x-s*y=0','x')
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)/(2*a), (s^(1/2)*x)/a^(1/2)))/(-1)^((a + 1)/(2*a))
而我在4楼计算的是想通过得到的解析解,反过来计算微分方程右边是不是等于0
请问你用的是软件是mathematica吗
不是哈,是Maple。
好吧,一开始我们老师也是教我们maple,后来就用matlab了,加上要用神经网络什么的,后来我看到另外一个老师用mathmatic的积分和求导等式子的形式漂亮多了,但是找不到安装包
安装包小木虫里就有的吧。
第一个解的时候,这种明显不是贝塞尔方程的形式(一次项分母有x),这位师兄,能不能和小弟我说一下,第一种解形式是怎么得来的……在线等……
嗷嗷,太谢谢了,不过是什么原因呢,为什么重新代入就不行了,,,,:sweat:
我刚才表达的不准确,师兄,为什么正常解法和laplace解法得到的解析解形式不一样,还有,第一种解析解的求解过程,师兄您有相关的资料吗
这是软件解出来的。具体解的过程涉及的专业知识,我并不是很了解。
这位师兄,送佛送到西啊。得到Y(x)的解析解之后,用Maple软件进行:
(1)比如由Y(0)=1;Y(1)=2确定出Y(x)的C1,C2
(2)画出0&x&10之间的Y(x)曲线。
谢谢啦。。。
这位师姐,得到Y(x)的解析解之后,用Matlab进行:
(1)比如由Y(0)=1;Y(1)=2确定出Y(x)的C1,C2
(2)画出0&x&10之间的Y(x)曲线。
谢谢啦。。。
& && && && && && && && && && && &&&求截图
Y(0)=1,Y(1)=2?
有没有Y'(0)=的数据?
Y(0)=1;Y(1)=2的条件的话,就是边值问题,哪有(或者很难)有解析解。
Y(0)=,Y'(0)= 的条件是初值问题,可以试着找解析解。
你的方程x在分母位置,x不可以为0,所以你的条件Y(0)=1需要修改
用matlab解方程组用
x=1;y1=subs(Y);
x=2;y2=subs(Y);
= solve(‘y1=1’,‘y2=2’)
解出C1,C2,你再代回去
给我红花吧^^
师兄啊,三个困惑(相互没有联系):
1.既然解析解里只有两个未知量C1,C2,按理说不是随便代入两个条件就可以确定吗?
2.再就是边值问题和初值问题有啥区别啊,但就像师兄所说的,小弟我确实遇到了Y(0)=1;Y(1)=2的条件不能解,而Y(0)=,Y'(0)= 的条件能解的情况。& &&&
不是有三类边界条件吗,Y(0)=1;Y(1)=2为第一类边界条件& & 而Y'(0)=为第二类边界条件,
好像都是边界(边值)条件啊,咋又分出的初值和边值问题啊
3.第三个问题是这样的。我遇到的模型条件是
边界条件:第一类边界条件
& && && && &&&Y(x=0) =45
& && && && &&&Y(x=正无穷大)=30
& && && && && &第二类边界条件
& && && && &&&Y'(x=正无穷大)&&=0
初始条件:&&Y(s&0) =30& && && && &
如果照师兄你说的,这里面利用哪些条件能得到确定的解析解。
& && && && & 非常希望留下师兄您的QQ或邮箱,方便小弟向您请教。
必须红花啊,这位师姐,方便留下你的QQ,可以向你继续请教吗
求帮助。。。
研究生必备与500万研究生在线互动!
扫描下载送金币 上传我的文档
 下载
 收藏
该文档贡献者很忙,什么也没留下。
 下载此文档
正在努力加载中...
.二阶常微分方程的幂级数解法(精)
下载积分:1500
内容提示:.二阶常微分方程的幂级数解法(精)
文档格式:PPT|
浏览次数:15|
上传日期: 17:37:33|
文档星级:
该用户还上传了这些文档
.二阶常微分方程的幂级数解法(精)
官方公共微信

我要回帖

更多关于 一阶常微分方程解法 的文章

 

随机推荐