2021-01-07 matlab数值分析 矩阵特征值与特征向量的计算 改进乘幂法 反幂法
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.矩阵特征值和特征向量定义 A为n阶矩阵,若数λ和n维非0列向量x满足Ax=λx,那么数λ称为A的特征值,x称为A的对应于特征值λ的特征向量.式Ax=λx也可写成( A-λE)x=0,并且|λE-A ...
- 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 ...
- QR分解求矩阵特征值、特征向量 C语言
最近在看一个高光谱图像压缩算法,其中涉及到正交变换,计算正交变换时,需要对普通矩阵求其特征向量.想要在网上找一个现成的程序,可能是我百度的能力不强吧,居然真的没找见.好了废话不多说,下面进入正题. 计 ...
- python numpy逆_Python使用numpy计算矩阵特征值、特征向量与逆矩阵
原标题:Python使用numpy计算矩阵特征值.特征向量与逆矩阵 Python扩展库numpy.linalg的eig()函数可以用来计算矩阵的特征值与特征向量,而numpy.linalg.inv() ...
- Python使用numpy计算矩阵特征值、特征向量与逆矩阵
Python扩展库numpy.linalg的eig()函数可以用来计算矩阵的特征值与特征向量,而numpy.linalg.inv()函数用来计算可逆矩阵的逆矩阵. >>> impor ...
- 使用MTL库求解矩阵特征值和特征向量
关于矩阵的特征值和特征向量求解,大部分的数学运算库都进行了提供,下面是使用MTL库的接口进行封装. #include <mtl/matrix.h> #include <mtl/mtl ...
- 求矩阵特征值和特征向量
求矩阵特征值和特征向量的一个小程序 代码较长,如果不能执行,就是要建立结构体,大家试试吧,希望能用. // // 实对称三对角阵的全部特征值与特征向量的计算 // // 参数: // 1. doubl ...
- pythonnumpy库求特征向量_Python使用numpy计算矩阵特征值、特征向量与逆矩阵
Python扩展库numpy.linalg的eig()函数可以用来计算矩阵的特征值与特征向量,而numpy.linalg.inv()函数用来计算可逆矩阵的逆矩阵. >>> impor ...
- ktt算法 约化_矩阵特征与特征向量的计算
矩阵特征与特征向量的计算 第三章第三章 矩阵特征与特征向量的计算矩阵特征与特征向量的计算3.1 引言引言在科学技术的应用领域中,许多问题都归为求解一个特征系统.如动力学系统和结构 系统中的振动问题,求 ...
最新文章
- linux 修改java版本_Linux 有问必答:如何在 Linux 中改变默认的 Java 版本
- 中文Python:中文编程不是梦
- Chrome原生工具实现长截图
- tab css html,纯css的tab 切换
- Vue warn Failed to mount component: template or render function not defined
- re.compile
- 【转载】架构师需要了解的Paxos原理、历程及实战
- 第一个冲刺周期(第四天)
- 【Elasticsearch】es fielddata 字段数据
- php忘记密码啦怎么办,后台账号密码忘记了怎么办.
- Graph Embedding及其在知乎的实践(附pdf下载链接)
- IT群侠传第一回内外兼修
- idea project settings在哪_IntelliJ IDEA如何对project的目录进行筛选显示?
- 基于bs4的python爬虫-链家新房(广州页面)
- 怎么用计算机表达爱意,不知道如何向心仪的她表达爱意的你 可能需要Zooids这帮小助攻...
- JS基础-特效篇-随机点名册
- 解决win10/Office2016/Onedrive/上载中心经常冲突的问题
- 巴菲特08年致股东信
- python 字典的值可以为集合吗_Python 集合(Set)、字典(Dictionary)
- PyPI 推送自己的