matlab数值分析  矩阵特征值与特征向量的计算

1改进乘幂法

function [t,y]=eigIPower(A,v0,ep)
[tv,ti]=max(abs(v0));
lam0=v0(ti);
u0=v0/lam0;
err=ep*10;             %为第一步循环做准备,此处不考虑0次循环的情况
while(err>ep)v1=A*u0;[tv,ti]=max(abs(v1));lam1=v1(ti);err=abs(lam0-lam1);u0=v1/lam1;lam0=lam1;
end
t=lam1;
y=u0;

调用程序

A=[12,6,-6;6,16,2;-6,2,16];xinit=[1;0.5;-0.5];[t,y]=eigIPower(A,xinit,1e-4)

2 反幂法

function [t,y]=eigIPower_inv(A,v0,ep)
[tv,ti]=max(abs(v0));
lam0=v0(ti);
u0=v0/lam0;
err=ep*10;
while(err>ep)v1=A\u0;[tv,ti]=max(abs(v1));lam1=v1(ti);err=abs(1/lam0-1/lam1);      %反幂法在误差计算时用的是特征值的倒数u0=v1/lam1;lam0=lam1;
end
t=1/lam1;y=u0;

调用程序

​​​​​​​

A=[12,6,-6;6,16,2;-6,2,16];xinit=[1;0.5;-0.5];[t,y]=eigIPower_inv (A,xinit,1e-4)

2021-01-07 matlab数值分析  矩阵特征值与特征向量的计算 改进乘幂法 反幂法相关推荐

  1. 矩阵特征值和特征向量详细计算过程(转载)

    1.矩阵特征值和特征向量定义 A为n阶矩阵,若数λ和n维非0列向量x满足Ax=λx,那么数λ称为A的特征值,x称为A的对应于特征值λ的特征向量.式Ax=λx也可写成( A-λE)x=0,并且|λE-A ...

  2. matlab求矩阵特征值和特征向量、行列式

    如果A为方阵,满足AX=λX的λ称为A的特征值,X称为A的特征向量. 计算A的特征值用eig(A). 例:A=[1 2 3;4 5 6;7  8 9]; Z=eig(A) Z =    16.1168 ...

  3. QR分解求矩阵特征值、特征向量 C语言

    最近在看一个高光谱图像压缩算法,其中涉及到正交变换,计算正交变换时,需要对普通矩阵求其特征向量.想要在网上找一个现成的程序,可能是我百度的能力不强吧,居然真的没找见.好了废话不多说,下面进入正题. 计 ...

  4. python numpy逆_Python使用numpy计算矩阵特征值、特征向量与逆矩阵

    原标题:Python使用numpy计算矩阵特征值.特征向量与逆矩阵 Python扩展库numpy.linalg的eig()函数可以用来计算矩阵的特征值与特征向量,而numpy.linalg.inv() ...

  5. Python使用numpy计算矩阵特征值、特征向量与逆矩阵

    Python扩展库numpy.linalg的eig()函数可以用来计算矩阵的特征值与特征向量,而numpy.linalg.inv()函数用来计算可逆矩阵的逆矩阵. >>> impor ...

  6. 使用MTL库求解矩阵特征值和特征向量

    关于矩阵的特征值和特征向量求解,大部分的数学运算库都进行了提供,下面是使用MTL库的接口进行封装. #include <mtl/matrix.h> #include <mtl/mtl ...

  7. 求矩阵特征值和特征向量

    求矩阵特征值和特征向量的一个小程序 代码较长,如果不能执行,就是要建立结构体,大家试试吧,希望能用. // // 实对称三对角阵的全部特征值与特征向量的计算 // // 参数: // 1. doubl ...

  8. pythonnumpy库求特征向量_Python使用numpy计算矩阵特征值、特征向量与逆矩阵

    Python扩展库numpy.linalg的eig()函数可以用来计算矩阵的特征值与特征向量,而numpy.linalg.inv()函数用来计算可逆矩阵的逆矩阵. >>> impor ...

  9. ktt算法 约化_矩阵特征与特征向量的计算

    矩阵特征与特征向量的计算 第三章第三章 矩阵特征与特征向量的计算矩阵特征与特征向量的计算3.1 引言引言在科学技术的应用领域中,许多问题都归为求解一个特征系统.如动力学系统和结构 系统中的振动问题,求 ...

最新文章

  1. linux 修改java版本_Linux 有问必答:如何在 Linux 中改变默认的 Java 版本
  2. 中文Python:中文编程不是梦
  3. Chrome原生工具实现长截图
  4. tab css html,纯css的tab 切换
  5. Vue warn Failed to mount component: template or render function not defined
  6. re.compile
  7. 【转载】架构师需要了解的Paxos原理、历程及实战
  8. 第一个冲刺周期(第四天)
  9. 【Elasticsearch】es fielddata 字段数据
  10. php忘记密码啦怎么办,后台账号密码忘记了怎么办.
  11. Graph Embedding及其在知乎的实践(附pdf下载链接)
  12. IT群侠传第一回内外兼修
  13. idea project settings在哪_IntelliJ IDEA如何对project的目录进行筛选显示?
  14. 基于bs4的python爬虫-链家新房(广州页面)
  15. 怎么用计算机表达爱意,不知道如何向心仪的她表达爱意的你 可能需要Zooids这帮小助攻...
  16. JS基础-特效篇-随机点名册
  17. 解决win10/Office2016/Onedrive/上载中心经常冲突的问题
  18. 巴菲特08年致股东信
  19. python 字典的值可以为集合吗_Python 集合(Set)、字典(Dictionary)
  20. PyPI 推送自己的

热门文章

  1. 里面使用轮播_小程序ColorUI框架初步使用教程及个人项目实战
  2. clodeblocks debug断点调试_Intellij IDEA高阶DEBUG大杀器
  3. LeetCode-剑指 Offer 50. 第一个只出现一次的字符
  4. 搜索引擎solr系列---多字段匹配的实现方法
  5. 社交产品后端架构设计--转载
  6. solr4.2增量索引之同步(修改,删除,新增)--转载
  7. Java Math 类中的新功能--浮点数
  8. 数模学习笔记——微分方程(传染病模型)
  9. 机器学习算法基础——k-means原理
  10. vue实战第一课:天气查询