现在时间是:
当前位置:首 页 >> 矩阵计算>> 教学区>> 文章列表

130213012李文萍 上机作业2

作者:李文萍   发布时间:2015-06-23 21:24:44   浏览次数:167

文件1:

% 书本P37 T1
%A*x=b    A=L*U<=>L*U*x=b   (1)用前代法解L*y=b得y   (2)用回代法解U*x=y得x
function x=jiex(A,b)
B=ones(1,83);  %1行83列全为1数组
C=6*ones(1,84);  %1行84列全为6数组
D=8*ones(1,83);
A=diag(B,1)+diag(C,0)+diag(D,-1);  %表示出A  diag(B,1):上对角线的值 diag(C,0):主对角线的值
b=[7;15*ones(82,1);14];
[L,U]=gauss(A);   %高斯分解出 L U
y=qiandaiy(L,b); 
x=huidaix(U,y);

文件2:

function [L,U]=gauss(A)
n=size(A,1);
tol=10^(-16);
for k=1:n-1
   if abs(A(k,k))>tol
       A(k+1:n,k)=A(k+1:n,k)/A(k,k);
        A(k+1:n,k+1:n)=A(k+1:n,k+1:n)-A(k+1:n,k)*A(k,k+1:n);
   end
end
L=eye(n)+tril(A,-1);  %严格下三角阵
U=triu(A,0);          %上三角阵
 

文件3:

function y=qiandaiy(L,b)
n=size(L,1);
y=zeros(n,1);   %n行1列全为0
y(1)=b(1)/L(1,1);
for(j=2:n)
    y(j)=(b(j)-L(j,1:j-1)*y(1:j-1))/L(j,j);
end

文件4:

function x=huidaix(U,y)
n=size(U,1);
x=zeros(n,1);
x(n)=y(n)/U(n,n);
for(j=n-1:-1:1)
    x(j)=(y(j)-U(j,j+1:n)*x(j+1:n))/U(j,j);
end

结果:

ans =

  1.0e+008 *

    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
   -0.0000
    0.0000
   -0.0000
    0.0000
   -0.0000
    0.0000
   -0.0000
    0.0000
   -0.0000
    0.0000
   -0.0000
    0.0000
   -0.0001
    0.0002
   -0.0003
    0.0007
   -0.0013
    0.0026
   -0.0052
    0.0105
   -0.0209
    0.0419
   -0.0836
    0.1665
   -0.3303
    0.6501
   -1.2582
    2.3487
   -4.0263
    5.3684

 

 

 

 







上一篇:没有了    下一篇:没有了

Copyright ©2019    计算数学达人 All Right Reserved.

技术支持:自助建站 | 领地网站建设 |短信接口 |燕窝 版权所有 © 2005-2019 lingw.net.粤ICP备16125321号 -5