解三对角线性方程组的追赶法

function [x,y]=after_method()
%数组a存储三角矩阵A的主对角线元素，c、d存储主对角线上边下边带宽为1的元素
%追赶法a=[2,2,2,2,2];c=[-1,-1,-1,-1];d=[-1,-1,-1,-1];b=[1,0,0,0,0];n=length(a);n1=length(c);n2=length(d);%错误检查if n1~=n2%存储矩阵的数组维数错误error('MATLAB:Crout:不是三对角矩阵，参数数组中元素个数错误.');elseif n~=n1+1error('MATLAB:Crout:不是三对角矩阵，参数数组中元素个数错误.');end%初始化L=zeros(n);%生成n*n的全零矩阵U=zeros(n);p=1:n;q=1:n-1;x=1:n;y=1:n;%追赶法程序主体p(1)=a(1);for i=1:n-1q(i)=c(i)/p(i);p(i+1)=a(i+1)-d(i)*q(i);%d的下标改为1到n-1end%正解yy(1)=b(1)/p(1);%用x存储yfor i=2:ny(i)=(b(i)-d(i-1)*y(i-1))/p(i);end%倒解xx(n)=y(n);for i=(n-1):-1:1x(i)=y(i)-q(i)*x(i+1);end%L,U矩阵for i=1:nL(i,i)=p(i);U(i,i)=1;endfor i=1:n-1L(i+1,i)=d(i);U(i,i+1)=q(i);disp([x;y])
end %end of function

解三对角线性方程组的追赶法相关推荐

c语言程序用追赶法求解方程组,编写用追赶法解三对角线性方程组的程序,并解下列方程组（3页）-原创力文档...
计算方法与实习上机实验(二) 实验名称: 编写用追赶法解三对角线性方程组的程序,并解下列方程组: (1) (2)Ax=b,其中 A10×10=, b10×1= 程序代码: #include using ...
c语言程序用追赶法求解方程组,编写用追赶法解三对角线性方程组的程序,并解下列方程组...
计算方法与实习上机实验(二) 实验名称: 编写用追赶法解三对角线性方程组的程序,并解下列方程组: 2x1 x2 5, x 2x x 12, 123(1) x 2x x 11,34 2 x3 2x4 1 ...
三对角线性方程组的追赶法
提出背景: 在一些实际问题中,例如解常微分方程边值问题.求热传导方程及三次样条插值函数等,都会遇到系数矩阵是三对角矩阵的方程组对于这种特殊的方程组,若还用原有的一般方法来求解,势必造成存储和计算的浪费 ...
Python解线性方程组的直接法(6)————求解三对角方程组的追赶法
求解三对角方程组的追赶法 import numpy as npdef zuiganfa(A, d):n = A.shape[0]l = np.mat(np.zeros(n, dtype=float
追赶法求解三对角线性方程组的MATLAB程序
function [M] = after_method(a,b,c,g,tolerance) %% 追赶法求解三对角矩阵 % a为三对角矩阵左下对角线上的值 % b为三对角矩阵中间对角线上的值 % c ...
求解线性方程组（SVD，QR，Gauss，LU）
曲线拟合过程中,需要求解线性方程组,下面谈谈线性方程组的求解方法: 1)svd求解对于齐次线性方程 A*X =0; 当A的行数大于列数时,就需要求解最小二乘解,在||X||=1的约束下,其最小二乘解 ...
用追赶法求方程组c语言程序,数值计算——追赶法求解三对角方程组（附代码）...
目录追赶法基础理论追赶法c++程序代码程序运行结果源码文件下载地址追赶法基础理论在数值计算中,对三次样条曲线插值和用差分方法求解常微分方程边值问题时,通常会遇到Ax=d三对角形式的方程组: ...
二、解线性方程组的直接方法
https://zhuanlan.zhihu.com/p/30485749 设 n n n阶线性方程组: { a 11 x 1 + a 12 x 2 + . . . + a 1 n x n = b 1 ...
现代数值分析 matlab,现代数值分析（MATLAB版）
商品描述: [编辑推荐]: 马昌凤编著的<现代数值分析>是普通高等教育十二五规划教材之一.本书共九章节,内容包括现代数值分析引论.非线性方程的求根方法.线性方程组的直接解法.线性方程组的迭 ...

解三对角线性方程组的追赶法

解三对角线性方程组的追赶法相关推荐

最新文章

热门文章