MATLAB在数理方程的应用.docx

MATLAB在数理方程的应用欧峥 2011326660117 11应数一班一、MATLAB在解偏微分方程中的应用 应用MATLAB 求解数学物理方程， 可通过编程或直接利用偏微分方程工具箱求解，直接利用偏微分方程工具箱更为简单、方便。在数理方程课上我们学习解矩形域方程的问题： 例如：在矩形域-0.5x0.5,-0.5y0.5 上求解Δu＝－2，且u 在边界上的值为零。利用偏微分方程工具箱求解过程如下：(1)启动偏微分方程求解界面。(2)在MATLAB 命令窗口中键入pderect ([-0.5 0.5 -0.5 0.5]), 选择Boundary→Remove All SubdomainBorders 菜单项，得出偏微分方程的求解区域。（3）单击偏微分方程界面工具栏中的PDE 图标，选择其中的Parabolic 选项，将给定的偏微分方程的参数输入到该对话框中。（4） 边界条件由Boundary 菜单下的Specify Boundary conditions 确定，输入边界条件h=1，r=0。（5）单击工具栏等号按钮，得到偏微分的解，单击图形设置按钮，得出如图1所示的结果。 MATLAB程序如下：clear all pderect([-0.5 0.5 -0.5 0.5 ]); %微分方程的求解区域h1=1,r=0; %初始条件二、MATLAB在解特殊方程中的应用：例如：画出第一类贝塞尔函数J15 （x）的曲线分布。MATLAB 语句为：x=0:0.01:50 %求解范围plot(x,BESSELJ(15,x)) %函数语句，输出图形运行结果如图3所示。图3贝塞尔函数（15 阶）J15 （x）的曲线分布三、MATLAB在傅里叶变换的应用在MATLAB语言中使用fourier函数来实现Fourier变换。首先定义符号变量t描述时域表达式f,直接调用fourier函数就可求出所需的时域函数的Fourier变换式。该函数的调用格式为F=fourier(f)采用默认的t为时域变量F=fourier( f，u，v),用户指定时域变量u和频域变量v使用ifourier函数来实现Fourier逆变换。它的使用格式如下:f=ifourier(F),按默认变量进行Fourier逆变换;f=ifourier(F，u，v),将u的函数变换成v的函数。例如:已知f(t)=e^(－t)sin2tu(t)，求f(t)的Fourier变换，并对结果进行Fourier逆变换。输入MATLAB 语句如下:syms t w %定义变量 f=exp(－t)*sin(2*t)*sym(‘Heaviside(t)’);%时域表达式f F=fourier( f) %求Fourier变换运行结果F= 2/(1+2*i+i*w)/(1－2*i+i*w) f=ifourier(F，w，t) %求Fourier 逆变换运行结果f=exp(－t)*sin(2*t)*Heaviside(t)返回原函数。 上面证明了在MATLAB语言中可以使用fourier函数来实现Fourier变换。下面我们来探讨能否在MATLAB语言中使用laplace函数来实现Laplace变换。四、MATLAB在Laplace变换的应用在MATLAB语言中使用laplace函数来实现Laplace变换。首先定义符号变量t，描述时域表达式f，直接调用laplace函数就可求出所需的时域函数的Laplace变换式。该函数的调用格式为:F=laplace(f)，采用默认的t为时域变量;F=laplace(f，u，v)，用户指定时域变量u和频域变量v。使用ilaplace函数来实现Laplace变换。它的使用格式如下:f=ilaplace(F)，按默认变量进行Laplace逆变换;f=ilaplace(F，u，v)，将u的函数变换成v的函数。例如:设函数f(t)=e^（－t）t^2.sin(t+π)，求其Laplace变换，并对结果进行Laplace逆变换。输入MATLAB 语句如下:syms t s %定义变量f=exp(－t)*t^2*sin(t+pi);%时域表达式f F=laplace(f，t，s) %求Laplace变换。运行结果F=－8/((s+1)^2+1)^3*(s+1)^2+