function [ D,P ] = Rayleigh( A,m )

%----------------------------------------------------------------

% m为多步瑞利迭代过程的次数

% A1表示A矩阵的实数部分

% A2表示A矩阵的虚数部分

% n表示A的阶数

%----------------------------------------------------------------

A1 = real( A );

A2 = A - A1;

n = length( A );

[x,D] = eig( A1 );

lamada = zeros( n,1 );

sign = 0;

for i = 1:n

for j =1:n

if A(i,j) == A1(i,j)

sign = 1;

break;

end

end

if sign == 1

break;

end

end

if sign == 0

P = x;

else

for i = 1:n

%         lamada(i) = D(i,i);

lamada(i) = 1;

end

for j = 1:m

A3 = j / m * A2;

A = A1 + A3;

y = zeros( n );

for k = 1:n

temp = 0;

while  abs(lamada(k)-temp)/abs(lamada(k)) > 0.001

lamada(k) = temp;

temp = x(:,k)' * A * x(:,k)/(x(:,k)'* x(:,k));

y(:,k) = x(:,k) / max(abs(x(:,k)));

if det(A-temp*eye(n)) <10e-5

break;

else

x(:,k) = ( A-temp*eye(n) ) \ y(:,k);

end

end

lamada(k) = temp;

end

end

P = y;

for i = 1:n

D(i,i) = lamada(i);

end

end

matlab迭代实现矩阵运算,用matlab实现Rayleigh迭代计算矩阵特征值的程序相关推荐

  1. 矩阵在Matlab中应用规则,1.3 MATLAB中的矩阵运算免费阅读_详解MATLAB在最优化计算中的应用免费全文_百度阅读...

    4.6 线性规划问题的MATLAB求解方法 与一般的线性规划理论一样,在MATLAB中有线性规划的标准型,和前面介绍的一般标准型有类似之处,也有不同之处.在调用MATLAB线性规划函数linprog时 ...

  2. matlab 自定义对象,面向对象:MATLAB的自定义类 [MATLAB]

    这几天刚刚开始学习MATLAB的面向对象编程.以前做的事情都是用MATLAB写一些简单的脚本或者函数,这方面MATLAB成熟的函数和直截了当的矩阵运算方法和语法都很容易上手,方便人专注于算法本身.前几 ...

  3. matlab的实用与仿真,MATLAB建模与仿真实用教程

    前言 第1章MATLAB建模基础 1.1 MATLAB简介 1.1.1 MATLAB的安装 1.1.2 MATLAB通用命令 1.1.3 MATLAB应用窗口简介 1.2 MATLAB数学建模概述 1 ...

  4. 线性代数 matlab,实用大众线性代数(MATLAB版)

    实用大众线性代数(MATLAB版) 作者:陈怀琛 著 出版时间:2014年版 内容简介 传统的线性代数源于数学家,教理论不教应用.工科需要的是应用,却学不到,2009年由本书作者牵头在19个大学实施了 ...

  5. matlab基础与实例教程,MATLAB基础与实例教程

    系统全面,实例丰富 考虑到Matlab进行仿真和运算分析时的基础知识和实践操作,讲解从基础的变量.函数.数据类型等入手,涉及到数学分析.图形可视化.Simulink仿真.文件读写等,全面地介绍了Mat ...

  6. bag of words matlab,Bag of words(matlab实现)

    Bag of Word主要思想:将训练样本特征Kmeans聚类,对测试样本的每个特征,计算与其最近的类心,相应类别计数count加1,这样每个测试样本可以生成ncenter维的直方图. 比如:训练样本 ...

  7. MATLAB学习笔记2:MATLAB基础知识(下)

    阅读前请注意: 1. 该学习笔记是华中师范大学HelloWorld程序设计协会2021年寒假MATLAB培训的学习记录,是基于培训课堂内容的总结归纳.拓展阅读.博客内容由 @K2SO4钾 撰写.编辑, ...

  8. 如何将excel表格导入matlab,将Excel数据导入MATLAB中的方法

    在使用MATLAB对矩阵进行数据处理时,为了方便编辑与修改,常常需要先将数据录入到Excel中,然后再将其导入到MATLAB中参与矩阵运算.那么下面小编教你怎么将Excel数据导入MATLAB中. 将 ...

  9. matlab牛顿法解非线性方程组,matlab实现牛顿迭代法求解非线性方程组.pdf

    资源描述 matlab 实现牛顿迭代法求解非线性方程组实现牛顿迭代法求解非线性方程组 已知非线性方程组如下 3*x1-cos(x2*x3)-1/2=0 x1^2-81*(x2+0.1)^2+sin(x ...

最新文章

  1. javascript图片懒加载与预加载的分析
  2. 基于ACE Proactor框架下高并发、大容量吞吐程序设计既最近的一个产品开发总结
  3. 《系统集成项目管理工程师》必背100个知识点-27产品范围和项目范围
  4. C++中的STRING数据 重载 = 和数组下标
  5. 设置代理_如何防止用户更改Windows 10上的代理设置
  6. JDBC中Statement与PreparedStatement的区别
  7. android使用handler记录
  8. scatter python_python数据可视化(matplotlib、scatter)
  9. 职场提醒:面试失败n次以后
  10. 面向IT专业人员的8个新兴AI工作
  11. 怎么在html页面和js里判断是否是IE浏览器
  12. vivado 一次性设置多线程编译
  13. 《CLR via C#》读书笔记 之 计算限制的异步操作
  14. 一款基于jQuery的图片分组切换焦点图插件
  15. SpringBoot实现文件在线预览
  16. 腾讯qq群推广“一键加群”的一个细节
  17. nw.js html5,用nw.js构建项目
  18. 体脂秤里的测脂模块方案,测量体脂全靠它?
  19. SSO(Single Sign On)系列(二)--SSO原理
  20. ssh连接提示IT IS POSSIBLE THAT SOMEONE IS DOING SOMETHING NASTY

热门文章

  1. 关于阿里云个人用户选择阿里云服务器配置说明
  2. mysql加减时间-函数-时间加减
  3. 【摸鱼吃瓜工作录】刚到公司,如何让项目经理对你刮目相看
  4. Mac 下有哪些能极大地提高工作效率的软件?
  5. 蝉知门户系统迁移到SAE平台-File模块扩展
  6. 基于android的有声听书系统
  7. 海面电磁散射MATLAB程序,matlab 电磁散射特性计算
  8. 10步搞定App内测发布(蒲公英内测平台)
  9. 微信小程序拼团头像叠加css 【亲测有效】
  10. Abbirb120型工业机器人_工控产品-ABB IRB120工业机器人-IRB120