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

原创 2012年03月30日 18:43:49

clear all;clc;
fprintf('输入n:(10,20,30)\n');
n=input('');
n
a=zeros(1,n);b=zeros(1,n);c=zeros(1,n-1);
A=zeros(n,n);
for i=1:n
    for j=1:n
        if j==i
            A(i,j)=2;
        elseif (j-i==-1)||(j-i==1)
            A(i,j)=1;
        else
            A(i,j)=0;
        end
    end
end
fprintf('系数矩阵:');A
fprintf('主对角线元素:');b=diag(A)'
fprintf('第1条对角线元素:');c=diag(A,1)'
fprintf('第-1条对角线元素:');a(2:n)=c
d=zeros(1,n);d(1)=1;d(n)=((-1)^(n+1));
fprintf('给定的系数矩阵(转置后):');d
y=zeros(n,1);x=y;
u=zeros(1,n);l=u;%u矩阵U的主对角线元素,l矩阵L的第-1条对角线元素
%---------追过程--------
u(1)=b(1);y(1)=d(1);
for i=2:n
    l(i)=a(i)/u(i-1);
    u(i)=b(i)-l(i)*c(i-1);
    y(i)=d(i)-l(i)*y(i-1);
end
fprintf('中间解向量:');y
fprintf('单位下三角阵L的第-1条对角线元素:');l
fprintf('上三角阵U的主对角线元素:');u  
U=zeros(n);
L=eye(n);
for i=1:n-1
    L(i+1,i)=l(i+1);
end
fprintf('单位下三角阵:');L    %单位下三角阵

for i=1:n-1
    U(i,i)=u(i);
    U(i,i+1)=c(i);
end
U(n,n)=u(n);
fprintf('上三角阵:');U    %上三角阵
%---------赶过程--------
x(n)=y(n)/u(n);
for i=n-1:-1:1
    x(i)=(y(i)-c(i)*x(i+1))/u(i);
end
fprintf('解向量:');x                                                     

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

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

追赶法解三对角线性方程组(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
  • 5343

数值计算-线性方程组求解(2)-追赶法解三对角矩阵-MATLAB实现

问题描述  在用差分法求解二阶常微分方程的边值问题、热传导问题及三次样条插值函数的求解问题中,都会遇到下面形式的阶数较高的三对角方程组AX=fAX=f,即 ⎡⎣⎢⎢⎢⎢⎢⎢⎢b1a2⋱c1b2an−...
  • baidu_33351337
  • baidu_33351337
  • 2017年10月26日 01:06
  • 230

数值分析 反幂法求矩阵按模最小特征值 MATLAB实现

值分析第四版 颜庆津 计算实习题P238 结果: %function [lam]=jingmi(ep) %反幂法求矩阵按模最小特征值 %2015.11.8 密密编写  (*^__^...
  • jingmiaa
  • jingmiaa
  • 2015年11月11日 18:29
  • 1840

二分法求解对称三对角矩阵的特征值

  • 2011年06月20日 10:58
  • 508KB
  • 下载

数值分析 jacobi迭代法求解线性方程组 MATLAB程序实现

数值分析 jacobi迭代法MATLAB程序实现   函数内容(保存为jacobi.m文件,在MATLAB中运行): %function [G,d,x,N]=jacobi(A,b) %Jacob...
  • jingmiaa
  • jingmiaa
  • 2015年11月06日 19:01
  • 5922

追赶法求解三对角方程组

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

利用追赶法来求解方程Ax=b的C++程序

利用追赶法来求解矩阵Ax=b。
  • datouniao1
  • datouniao1
  • 2015年11月28日 18:56
  • 1079

追赶法的MATLAB程序

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

求解三对角方程组的追赶法(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
  • 10360
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:求解三对角方程组的追赶法(Matlab程序)
举报原因:
原因补充:

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