行处理法求解三对角线性方程组的C语言实现.pdfVIP

维普资讯 2 28 行处理法求解三对角线性方程组的c语言实现 第 13卷 行处理法求解三对角线性 方程组 的 C语言实现 n I． 、一， L／ 塑垦生 壹—兰 张 玲 李安志 杨本立 ． ／ (中国工程物理研究院职工工学院 四川绵阳621900) -rf~。lt。 摘 要 本文给出利用线性代数方程组的行处理法求解 三对角线性代数方 程组的 c语 言程序实现方法 关键词 三对角线性方翟组 盟竺矍鎏 压缩存储 排序加速技术 1 概 述 =，《秀雩序 在样条插值和微分方程的差分解法等问题中出现的线性代数方程组都是三对角线性代 数方程组 。用行处理法求解三对角线性方程组 ，明显地兼有直接法和迭代法二者 的优点 。 用(‘，)‘表示两向量的内积 ，用 lf．II2表示向量的欧式范数。 定义 i．I 记NxN阶线性代数方程组 r( ，)=b {IIa。II2=I (I—1) 【i：l2，3，… ，n ， 为 AX =b (I一2) 定义 1．2 设有 N阶线性代数方程组 (1—1)或 (1—2) 记A=(alj)… l≤i、j≤n ． 若 当l j—il≤l时 atjC-R而 当lj—il1时 a“=O，则称 (1一1)或 (1—2)为三对角线性方程组。 定义 1．3 设有 N阶三对角线性代数方程组fl 1一)或 (1—2)。称下面算法 P为 (1 1)或 (1—2)的排序加速技术。 Mgoriflun P： t：=O： ForK：=0to(n—I) Dobe西ns：=3kmo~(n)； Ifs=0 thent：= t+1 S： s+t： P【k+1J：=s · 中国工程物理研究院预研基金资助课题 收稿 日期 ：1999—03—20 维普资讯 2O0O年 3月 第 1期 教学与科技 end； 定义 1．4 设有 N阶三对角线性代数方程组 (1一1)或 (1—2)，记其超平面序号序列为 I = (1，2，…，n)。称下面算法A为求解此方程组的堑熊里鲨。 Al dⅡn A： P：，_+P[，] 排序加速 任取 ∈ 妒 “： +( 一( ， ))a = P[i]