谱方法解偏微分程.docVIP

谱方法解偏微分方程 学生：石幸媛，数学与计算机科学学院 指导老师：陈慧琴 ，江汉大学数学与计算机科学学院 学号：200808101125  摘要 本论文分析的是偏微分方程的谱方法解。在此，我借用向新民编的《谱方法的数值分析》中第67页例2.1方程进行计算。根据例2.1的谱方法计算方式，给该方程具体的函数进行计算，求解其值，并绘图。最后研究比较一阶波动方程的Fourier谱方法与Fourier配点逼近有什么不同与相近之处，做出结论。 关键词：Fourier配点逼近，截断函数，插值函数,Fourier谱方法 Abstract This paper analyses the partial differential equations of the spectral method. Here, I use the Xiang Xinmin series numerical analysis of spectral method on page sixty-seventh example 2.1equation. According to the case of 2.1spectral methods for computing method, give the specific function for calculating equation, solving its value, and drawing. The final study comparing a first-order wave equation in Fourier spectral method and Fourier collocation approximation of what is the difference and similarities, make a conclusion. Key words: Fourier collocation approximation, truncated function, interpolation function, Fourier spectral method  目录 绪论 4 论文主题 5 §1定义引用： 5 §2论文内容： 5 2.1：Fourier配点法 5 结论 11 致谢 12 参考文献 13  绪论 谱方法是70年代发展起来的一种数值求解偏微分方程的方法，它具有“无穷阶”收敛性，可采用快速算法，现已被广泛用于气象、物理、力学等诸多领域，成为继差分法和有限元法之后又一种重要的数值方法。π)=o,改变a(x),f(x)或u(x)函数，然后在计算机中运用Matlab程序 进行计算。Matlab程序设置或步骤如下： lamal=0;%可以对应修改的常数 N=100; h=2*pi/(2*N); for k=-N:N-1 for j=1:2*N x=j*h; C(k+N+1,j)=1/(2*N)*exp(-i*k*x); end end k1(1)=0; for k=-N+1:N-1 k1(k+N+1)=k*i; end K=diag(k1); for j=1:2*N x=j*h; a(j)=1; %可以对应修改的函数 %a(j)=x; f(j)=-sin(x);%可以对应修改的函数 %f(j)=x; end A=diag(a); F=f; I=ones(2*N,2*N); D=inv(C)*K*C; AA=D*A*D-lamal*I; for j=1:2*N-2 AA1(:,j)=AA(:,j+1); end