原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
的特解进行讨论并用Matlab程序实现绘制特解的变化图像[精选]
０、前言：
到目前为止已经有不少人对Bessel方程进行了讨论，但对Bessel方程的特定边界条件下的特解情况的讨论却很少。本文主要是针对一类零阶Bessel方程在如下边界条件
的特解进行讨论并用Matlab程序实现绘制特解的变化图像。
Matlab是类似于Fortran和C的一种语言，虽然很难简短描述其特点，但对于科学计算而言，其主要特色有：
易于编程
整数、实数、复数之间的统一性
高精度及扩充的数值范围
综合性的数学工具库
包含图形用户界面在内的功能完备的图形工具
与传统编程语言的接口
Matlab程序的可移植性
本文便是利用的强大图形工具绘制出方程特解的图像，以便形象地描述Bessel方程特解的振荡性、衰减性、周期性等。
１、数学模型的求解及分析：
数学模型：
1.1 方程通解的求解过程
首先对方程
(1)
进行化简。引入，则有
=
(2)
=
将方程(1)两边同时乘以后，经整理后得:
(3)
将方程(2)代入到方程(3)有:
(4)
对比Bessel标准方程
（这里是非负常数，不一定是正整数）知，方程（4）其实就是一个的Bessel方程；即零阶的Bessel方程。
根据参考文献⑴158页定理11:若方程
①
中系数具有这样的性质,即和均能展成的幂级数,且收敛区间为,则方程①有形如
即
②
的特解,这里,是一个特定的常数。级数②也以为收敛区间。
将方程（4）改写成
易见，它满足定理11的条件，且, , 按展开成幂级数,它的收敛区间为，由定理11，方程有形如
(5)
的解.这里。而和是待定常数。
将（5）代如（4）有
按t的同幂次项整理有：
令各项系数等于零，得一系列的代数方程：
(6)
因为，故从(6)的第一个方程解得：
考虑时方程(4)的一个特解.这时总可以从(6)中逐个地确定所有的系数。把代入（6）得到
k=2,3,……
对讨论奇偶两种情况，分别有
k=1,2,3,……..
从而求得
k=1,2,3,……..
……………………
k=1,2,3,……
将各代入方程(5)得到方程(4)的一个解
(7)
不妨令
2
则(7)变为
根据参考文献⑵P404 (7.9):
(n=0时,右边第二项不存在)
可以得到
(8)
当为实数时, 是一个衰减振荡函数, 在有奇性。 、都是Bessel方程的解，且其中任意两个都是线性无关的。所以Bessel方程（4）的通解可写成：
(9)
1.2 在左右边界条件下方程特解的求解过程
1.2.1左边界条件
根据Bessel函数性质有
即有
文档评论（0）