数值分析 追赶法求解三对角线性方程组 MATLAB实现

原创 2015年11月05日 11:48:53
函数主体部分编程算法  参考 数值分析 第四版 颜庆津 P27
运行结果截图:

%追赶法求解三对角线性方程组,Ax=b,A用一维数组a,c,d存储。
function [L,U,x]=crout(a,c,d,b)%数组a存储三角矩阵A的主对角线元素,c、d存储主对角线上边下边带宽为1的元素
    n=length(a);
    n1=length(c);
    n2=length(d);
    %错误检查
    if n1~=n2%存储矩阵的数组维数错误
        error('MATLAB:Crout:不是三对角矩阵,参数数组中元素个数错误.');
    elseif n~=n1+1
        error('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-1
        q(i)=c(i)/p(i);
        p(i+1)=a(i+1)-d(i)*q(i);%d的下标改为1到n-1
    end
    %正解y
    y(1)=b(1)/p(1);%用x存储y
    for i=2:n
        y(i)=(b(i)-d(i-1)*y(i-1))/p(i);
    end
    %倒解x
    x(n)=y(n);
    for i=(n-1):-1:1
        x(i)=y(i)-q(i)*x(i+1);
    end
    %L,U矩阵
    for i=1:n
        L(i,i)=p(i);
        U(i,i)=1;
    end
    for i=1:n-1
        L(i+1,i)=d(i);
        U(i,i+1)=q(i);
end %end of function

求解三对角方程组的追赶法(Matlab程序)

clear all;clc; fprintf('输入n:(10,20,30)\n'); n=input(''); n a=zeros(1,n);b=zeros(1,n);c=zeros(1,n...
  • zhangchao3322218
  • zhangchao3322218
  • 2012年03月30日 18:43
  • 10354

追赶法解三对角线性方程组(Matlab)

%追赶法解三对角线性方程组 clear;clc A=[1,2,0,0,0;2,3,1,0,0;   0,-3,4,2,0;0,0,4,7,1;0,0,0,-5,6]; d=[5 9 2 19 ...
  • zhangchao3322218
  • zhangchao3322218
  • 2011年09月23日 08:07
  • 5341

追赶法求解三对角方程组

在这篇文章里,我们介绍追赶法的基本原理,以及用追赶法求解三对角方程组的算法....
  • u010450214
  • u010450214
  • 2015年12月04日 09:37
  • 7658

追赶法求解三对角方程的MATLAB程序代码

  • 2013年05月14日 16:52
  • 695B
  • 下载

计算方法之用追赶法求线性方程组

/************************************* * 用追赶法求线性方程组 * * |- -| |- -| |- -| * | 2 -1...
  • wzhg0508
  • wzhg0508
  • 2013年06月11日 12:57
  • 1197

追赶法的MATLAB程序

  • 2009年11月24日 18:51
  • 106KB
  • 下载

数值分析 追赶法求解三对角线性方程组 MATLAB实现

函数主体部分编程算法  参考 数值分析 第四版 颜庆津 P27 运行结果截图: %追赶法求解三对角线性方程组,Ax=b,A用一维数组a,c,d存储。 function [L,U,x]=cr...
  • jingmiaa
  • jingmiaa
  • 2015年11月05日 11:48
  • 4890

Guass列主元、平方根法、追赶法求解方程组的C++实现

一,要解决的问题选用合适的算法,求解三种线性方程组:一般线性方程组,对称正定方程组,三对角线性方程组。 方程略。二,数值方法1,使用Guass列主元消去法求解一般线性方程组。Guass列主元是为了防...
  • qq_24421591
  • qq_24421591
  • 2016年03月28日 22:05
  • 1471

求解线性方程组--追赶法

求解线性方程组--追赶法 /* 求解线性方程组--追赶法 *///####################################//# 本分解法仅适用于当 #...
  • gouki04
  • gouki04
  • 2009年08月27日 15:41
  • 1157

MATLAB常用算法

  • 2010年04月05日 10:34
  • 129KB
  • 下载
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:数值分析 追赶法求解三对角线性方程组 MATLAB实现
举报原因:
原因补充:

(最多只允许输入30个字)