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

《傅里叶积分、傅里叶变换的matlab实现.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、院校：物理与电子科学学院班级：0801班姓名：9目录1.引言………………………………………………………………………………2.理论推导…………………………………………………………………………2.1傅里叶级数……………………………………………………………………2.2傅里叶积分及傅里叶变换……………………………………………………2.3傅里叶积分、傅里叶变换的应用……………………………………………2.3.1对无限长的细杆导热问题的研究…………………………………………2.3.2对长度为的细杆导热问题的研究…………………………………………2.3.3波动方程的定解条件…………………………………

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

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

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

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

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

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