5点差分格式的Matlab程序.docVIP

三对角与块三对角方程组课程设计 一、基于高斯消元法的三对角方程组求解 三对角矩阵是一类重要的特殊矩阵，在数学计算和工程计算中有广泛应用。例如，二阶常微分方程边值问题数值求解，一维热传导方程数值求解，以及三次样条函数计算等都会涉及到三对角方程组求解。由于三对角矩阵的稀疏性质，用直接法求解三对角方程组的算法效率较高，很有实用价值。 考虑n阶三对角矩阵和n维向量 A =， 求解方程组 Ax = f 的高斯消元法的程序如下 function f=triGauss(gama,alpha,bata,f) %Solving TriDiag(gama,alpha,bata)systems by Gauss method n=length(alpha); for k=1:n-1 m=gama(k)/alpha(k); alpha(k+1)=alpha(k+1)-m*bata(k); f(k+1)=f(k+1)-m*f(k); end f(n)=f(n)/alpha(n); for k=n-1:-1:1 f(k)=(f(k)-bata(k)*f(k+1))/alpha(k); end 由程序知，对于n阶三对角方程组，高斯消元法只用到 5n –4 次乘法和除法。 例1．求二阶常微分方程边值问题 数值解，并与解析解：作对比。 解：对正整数n，取h= 1/(n + 1)，令