数值代数第二章实验报告

最新推荐文章于 2023-05-24 16:30:00 发布

扛着豆芽菜

最新推荐文章于 2023-05-24 16:30:00 发布

阅读量1.4k

点赞数

分类专栏：计算数学

本文链接：https://blog.csdn.net/baidu_41855511/article/details/88956623

版权

该实验报告详细记录了使用算法2.5.1对5到20阶Hilbert矩阵进行∞范数条件数的估计过程。随着阶数增加，条件数呈现指数增长，表明高阶Hilbert矩阵的求解可能面临数值稳定性问题。

摘要由CSDN通过智能技术生成

数值代数第二章实验报告

吕甜心 201511200135

算法2.5.1

function [v] = EsNorm1(A)
%求A的逆的的转置的一范数估计
k=1;
n=length(A);
x=(1/n)*ones(n,1);
while k==1
    w=ColMain(A',x);
    v=sign(w);
    z=ColMain(A,v);
    if norm(z,inf) <= z'*x
        v=norm(w,1);%w的1范数
        k=0;
    else
        [~,j]=max(abs(z));
        x=zeros(n,1);
        x(j)=1;
        k=1;
    end
end
end