第2章:误差分析 三对角线性方程组的存贮方式 对于N个变元利用三对角线性方程组我们可以说明 4个长度为N的一维数组AN],BN],cN],DN来存 贮所有的非零元素注意到方程组中的a1和cN没有 定义,我们可以简单地把数组中相应的元素设为零 由于c语言数组的下标是从零开始的,所以我们编程 时也要作相应的调整。 解的递推形式 为了讨论,计算,编程的方便,我们再定义两个数组 UN]和VN]计算格式为: a,v d.-akvk-l Vk=1,2,…,N b 结论:对于上面定义的向量U,,我们有 Xk=uk-k°Xk+1,K=1,2,…,N-1第 2 章：误差分析 - 6 - 6/11 三对角线性方程组的存贮方式 对于 N 个变元利用三对角线性方程组，我们可以说明 4 个长度为 N 的一维数组 A[N],B[N],C[N],D[N]来存 贮所有的非零元素。注意到方程组中的 a1 和 cN 没有 定义，我们可以简单地把数组中相应的元素设为零。 由于 C 语言数组的下标是从零开始的，所以我们编程 时也要作相应的调整。 解的递推形式 为了讨论，计算，编程的方便，我们再定义两个数组 U[N]和 V[N],计算格式为： k N b a v d a v v b a v d a u u b c v b d u k k k k k k k k k k k k k k , 1,2, , 1 1 1 1 1 1 1 1 1 1  =                − − = − − =               = = − − − − 结论：对于上面定义的向量 U,V，我们有 xk=uk-vk·xk+1,K=1,2,…,N-1