傅里叶积分、傅里叶变换的malab实现.doc

《傅里叶积分、傅里叶变换的malab实现.doc》由会员分享，可在线阅读，更多相关《傅里叶积分、傅里叶变换的malab实现.doc（10页珍藏版）》请在沃文网上搜索。

1、 目 录1. 引言2. 理论推导2.1傅里叶级数 2.2傅里叶积分及傅里叶变换 2.3傅里叶积分、傅里叶变换的应用 2.3.1对无限长的细杆导热问题的研究 2.3.2对长度为的细杆导热问题的研究2.3.3波动方程的定解条件 3. matlab模拟结果4. 总结5. 参考文献9傅里叶积分、傅里叶变换及其应用的matlab实现摘要：根据傅里叶积分、傅里叶变换理论，计算了若干例题，并利用此理论模拟了无限长细竿、有限长细竿的导热问题及波动方程的定解条件问题，做出了细竿导热情况的图像。关键词：傅里叶积分 傅里叶变换 热传导 定解问题1. 引言计算物理学是以计算机及计算机技术为工具和手段，运用计算数学的方

2、法，解决复杂问题的一门学科。傅里叶积分及傅里叶变换在物理学中有着重要的应用，而其运算相对繁琐，利用计算机技术可以方便地帮助我们解决这一问题，大大节省时间，提高研究效率。傅里叶积分及傅里叶变换作为重要的计算方法被应用在物理学中的各个领域。如量子力学、电动力学等等。我们选择用matlab解决傅里叶变换的计算问题；绘制出有限长和无限长细竿热传导温度分布图像，并对其作深入分析；解决波动方程定解条件的问题。2.理论推导2.1傅里叶级数若函数以为周期，即则，将展开为级数其中 若是定义在上的非周期函数，则可以采取延拓的方法，使其成为某种周期函数，而在上，。然后再对作傅里叶级数展开，使级数和在区间上代表。2.

3、2傅里叶积分及傅里叶变换傅里叶积分实际上是把定义在上的非周期函数进行积分形式的展开。即把展开为如下形式：其中第一个式子是傅里叶积分表达式，第二组式子为傅里叶变换式。把傅里叶积分写成复数形式就为傅里叶变换为下面举两道例题。例1 求矩形函数的傅里叶变换，其中解例2 求的傅里叶变换，其中，定义在上。解2.3傅里叶积分、傅里叶变换的应用基于maltab在数学物理方法中利用分离变数（傅里叶级数）法求解一维（线性）热传导方程问题的研究，在一维细杆热传导问题的研究将细杆分为有限长度与无限长度两方面来求解问题。2.3.1对无限长的细杆导热问题的研究无限长细杆的热传导的定解问题：细杆上任意一点的温度是时间t和位

4、置x的函数u(x,t)泛定方程 初始条件 利用傅里叶级数求得细杆上任意一点的温度为：若取初始温度分布设为 一个高度为一得矩形脉冲波；则得到2.3.2对长度为的细杆导热问题的研究讨论有限长度的细杆，在一端为第一类齐次边界条件，另一端为第二类边界条件下的热传导问题的研究。有限长细杆热传导定解问题就是将上述无限细杆的长度有限化，对取一确定有限值：泛定方程 边界条件 初始条件 当，时，解得将上述问题具体化为，初始时刻杆的一端温度为零度，另一端温度为,杆上的温度均匀（），零度温度一端保持温度不变（），另一端跟外界温度绝热（），这细杆上温度随时间与空间变化的函数关系设为。细杆上温度综合上述条件：泛定方程

5、边界条件 初始条件 由齐次方程的定解问题的求解方法求得将上述参数具体化，设定，则可化为2.3.3波动方程的定解条件一根长为两端固定的弦，用手把它的中点横向拉开距离 ，然后放手任其自由振动，写出它的初始条件。时各点的位移由图中折线确定，所以研究两端固定均匀弦的自由振动,即定解问题是：它的解是：其中对于有限长的弦，如果在讨论的时间范围内，边界的影响还没有到达，则产生的现象与无限长的弦是一样的。3. matlab模拟结果图1为例题1傅里叶变换的函数图像图1图2为例题2傅里叶变换的函数图像图2图3为无限长杆温度随时间和空间变化的瀑布图图3从图3中可以看出，在开始时刻，温度分布在原点附近定义为一个脉冲函

6、数，在沿着细杆的方向上，温度逐渐降低形成一个平缓的波包，并向周围传导，如果时间足够长，最终细杆上的温度为零。在前面的程序上加上以下程序，则图4表示杆上温度暂停0.1s时刻的传导情况：图4图5为有限长杆温度传导函数的图像图5由于初始条件相同，有限杆与无限杆的温度分布是一样的，无限长的杆热传导现象只是边界条件还没有产生影响的有限杆上热传导现象的一种近似。由于在理论的计算中，n的叠加到无穷，而以上程序中n只取到50，在图像中可以看到，在x=10到11的两端，温度出现较小幅度的波动，没有无限杆的热传导温度在两端都减小而不增大的现象，而当杆长趋近于无穷时，使得两图可以近似为同一图形。图6为有限长细杆上温

7、度随时间和空间变化的三维曲面图图6从图中可以看出，在沿杆的方向上，温度是随条件定义的线性传导变化，在同样的泛定函数下，给定不同的初始条件下，图四可近似看图三的一个部分。 如果考虑先前时刻即考虑t0,则随k的增大而急剧减小，从而级数收敛的很快。t越大，级数收敛越快。在t0.18时，可以只保留k=0的项，略去k0的项，从而简化，在MATLAB的程序中也可不使用forend的循环语句，画出的图像也大致相同，其误差不超过1%。4. 总结基于理论分析，我们展示了傅里叶变换的函数图像，无限长及有限长细竿热传导温度的分布规律并作出其图像，波动方程的定解问题。其中，我们利用图像形象直观地表现了所研究的问题，且理论分析也较为透彻，有待改进之处在与所研究问题关联性不强。5. 参考文献【1】 梁坤淼.数学物理方法.第三版.北京：高等教育出版社，1998.【2】 张志涌等.MATLAB教程.北京：航空航天大学出版社，2006.【3】 彭芳麟.数学物理方程的MATLAB解法与可视化.【4】 郑阿奇.MATLAB实用教程.第二版.电子工业出版社，2007.【5】 周晓阳.数学实验与MATLAB.华中科技大学出版社，2002.