autumn20080101的专栏

function x=MyTDMA %首先说明：追赶法是适用于三对角矩阵的线性方程组求解的方法，并不适用于其他类型矩阵。 %定义三对角矩阵A的各组成单元。方程为Ax=f % b为A的对角线元素(1~n)，a为-1对角线元素(2~n)，c为+1对角线元素(1~n-1)。 % A=[2 -1 0 0 % -1 3 -2 0 % 0 -2 4

//追赶法 #include #include #define N 4 //定义一个4*4的矩阵，改变N可以改变矩阵的大小 void TDMA(float a[N],float b[N],float c[N-1],float f[N],float x[N]) { float u[N],l[N],y[N]; u[0]=b[0]; y[0]=f[0