二阶矩阵特征值与特征向量的计算matlab程序（可指定计算精度）

最新推荐文章于 2024-07-22 09:13:01 发布

光电学子

最新推荐文章于 2024-07-22 09:13:01 发布

阅读量2.5k

点赞数 1

分类专栏：信号与图像处理

本文链接：https://blog.csdn.net/u014230360/article/details/107998584

版权

信号与图像处理专栏收录该内容

45 篇文章 71 订阅

订阅专栏

先前由于项目需要，自己实现了二阶矩阵特征值与特征向量的求取代码，其实很简单，主要部分就是求解一个二次方程的根，现在分享给大家作为参考（因为组长说求平方根可能比较费资源，为方便硬件实现，甚至连平方根求取的代码都是自己写的\笑哭）：

function [V,lamda]=myEig(A,err);
% 求特征值与特征向量的程序
% err为指定的计算误差
a11=A(1,1);
a12=A(1,2);
a21=A(2,1);
a22=A(2,2);
delta=(a11-a22)^2+4*a12*a21;
if delta<0
    disp('the input matrix has no eigen values and featrue vector!');
end
sqrtDelta=sqrtByBisection(delta,err);
eig1=((a11+a22)-sqrtDelta)/2;
eig2=((a11+a22)+sqrtDelta)/2;
lamda=[eig1,eig2];
v1=[a12,eig1-a11];
v2=[a12,eig2-a11];
v1=v1/norm(v1);
v2=v2/norm(v2);
V=[v1;v2];

function mid=sqrtByBisection(x,err)
% 二分法开平方根程序
% x为输入值
% err为计算精度
if err^2>x
    mid=0;
    return;
end
up=max(1,x);
low=0;
mid=(up+low)/2;
while ~( (mid+err)^2>x && (max(0,mid-err))^2<x)
    if mid^2<x
        low=mid;
    else
        up=mid;
    end
    mid=(up+low)/2;
end