二阶矩阵特征值与特征向量的计算matlab程序(可指定计算精度)_求二次型的特征值的matlab算法-程序员宅基地

技术标签: 信号与图像处理  

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

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

 

 

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/u014230360/article/details/107998584

智能推荐

数据库MySQL-查询练习2_查询选修“3105”课程的成绩高于“109”号同学成绩的所有同学的记录。-程序员宅基地

文章浏览阅读3.6k次。11.查询每门课的平均成绩关键字 avg()单门课: SELECT c_no,AVG(sc_degree) FROM SCORE where c_no = '3-105';每门课:SELECT c_no,AVG(sc_degree) FROM SCORE GROUP BY c_no;_查询选修“3105”课程的成绩高于“109”号同学成绩的所有同学的记录。

windows10安装MongoDB基础详细教程_mongodb安装教程win10-程序员宅基地

文章浏览阅读1.5w次,点赞13次,收藏62次。MongoDB 是一个基于分布式文件存储的数据库。由 C++ 语言编写。旨在为 WEB 应用提供可扩展的高性能数据存储解决方案。MongoDB 是一个介于关系数据库和非关系数据库之间的产品,是非关系数据库当中功能最丰富,最像关系数据库的。......_mongodb安装教程win10

mysql:排错技巧(1)general_log的使用_general_log=on-程序员宅基地

文章浏览阅读1.3k次。来源:<<MySql排错指南>>**作用:**从mysql.general_log(mysql自带的) 表里找到出错时执行的完整sql语句,提取出来手动执行,看看到底发生了什么。使用方法:注意事项:1.建议不要在生产环境使用,因为每条通过该数据库的语句都会被记录到general_log表里。2.使用前用 truncate 清理下 general_log。步骤一:..._general_log=on

keytool生成keystore、truststore、证书以及SSL单向认证在服务端tomcat和客户端的配置-程序员宅基地

文章浏览阅读1.3k次。1 SSL单向认证概念  当客户端(服务请求方)向服务端(服务提供方)发起请求时,服务器端需要向客户端提供认证。服务端需要生成一个keystore和一个服务器密钥对儿(公钥和私钥),客户端需要生成一个truststore,然后导入服务端的公钥证书。2 keystore以及服务器密钥对儿的生成keytool -genkeypair -alias certificatekey _truststore

tf.contrib.crf.crf_log_likelihood()的用法-程序员宅基地

文章浏览阅读7.5k次,点赞2次,收藏10次。函数:tf.contrib.crf.crf_log_likelihood(inputs, tag_indices, sequence_lengths, transition_params=None)在一个条件随机场里面计算标签序列的log-likelihood,函数的目的是使用crf来计算损失,里面用到最大似然估计的优化方法参数:inputs: 一个形状为[batch_size,..._crf_log_likelihood

Ipopt安装_ipopt安装matlab-程序员宅基地

文章浏览阅读1.3k次。https://blog.csdn.net/qq_24649627/article/details/103084849_ipopt安装matlab

随便推点

matlab添加新工具出现与原matlab函数冲突_警告: 函数 labindex 与某个 matlab 内置函数同名。建议您重命名该函数以避免潜在-程序员宅基地

文章浏览阅读8.5k次。matlab添加新工具导致函数冲突_警告: 函数 labindex 与某个 matlab 内置函数同名。建议您重命名该函数以避免潜在

uniapp详细介绍-程序员宅基地

文章浏览阅读2.3w次,点赞21次,收藏202次。uniapp的详细介绍_uniapp

【转】ora-00031:session marked for kill(标记要终止的会话)解决方法-程序员宅基地

文章浏览阅读147次。今天碰到一个问题,有一张表不能操作,很可能是被锁了,首先想到的是kill session,于是执行了下列的脚本找到是哪个session有问题:selectobject_name, machine, s.sid,s.serial#fromv$locked_object l, dba_objects o, v$session swherel.object_id = o.obj...

激光雷达:905与1550的战争-程序员宅基地

文章浏览阅读1.5k次,点赞2次,收藏11次。交流群 |进“传感器群/滑板底盘群”请加微信号:xsh041388交流群 |进“汽车基础软件群”请加微信号:ckc1087备注信息:传感器/滑板底盘/汽车基础软件+真实姓名、公司、岗位在去年8月底发布的《什么是激光雷达的“发动机技术”?一文讲透行业技术壁垒》一文中,九章智驾曾提出:激光雷达真正的壁垒在于激光收发系统,而不是“机械旋转”“半固态”“转镜”“棱镜”“MEMS”这些词语所代表的扫描系..._128个盖格探测器

python-给定任意起始点-角度后同等切分圆-circle-segementation_python将一个圆等分-程序员宅基地

文章浏览阅读742次。简介文中主要帮助一位小哥实现对圆的任意切分,切分需要从给定点,和给定角度进行切分。现把算法记录如下,需要使用的可以参考,当然比较简单。算法和效果如下。代码import mathimport matplotlib.pyplot as pltimport numpy as np# def return_xy_list(alfa,):# passr = 11.599xa = ..._python将一个圆等分

html5开发之ios屏幕适配,iOS开发屏幕尺寸以及屏幕适配等问题(转载内容)-程序员宅基地

文章浏览阅读383次。原帖地址:http://blog.csdn.net/phunxm/article/details/42174937/仅供我个人收藏学习,原博主如不同意请联系qq651263878进行删除,在此表示感谢以及歉意。1.iPhone尺寸规格后续上市的iPhone7以及iPhone7plus 与六代相同1 inch = 2.54cm = 25.4mm上表中的宽高(width/height)为手机的物理尺..._ios html注入屏幕尺寸适配