以下伪码来自某《数值线性代数》教科书,非原创

乘幂法

function [a,k,er]=Power(A)
%定义乘幂法函数文件
%a:表示该方法下矩阵A的最大特征值
%k:表示停机时实际的迭代次数
%er:表示停机时实际的绝对误差
tol=1e-6;%绝对误差限
[n,m]=size(A);
if n~=m %判断输入的合法性error('Wrong  Input');
end
u=ones(n,1);q=max(abs(u));
k=0;
while truey=A*u;q0=q;q=max(abs(y));if abs(q-q0)<tol|k>1e4 %设置停机条件break;endu=y/q;k=k+1;
end
er=abs(q0-q);
a=q;

反幂法

function [a,k,er]=InversePower(A)
%定义反幂法函数文件
%a:表示该方法下矩阵A的最小特征值
%k:表示停机时实际的迭代次数
%er:表示停机时实际的绝对误差
tol=1e-6;%绝对误差限
[n,m]=size(A);
if n~=m %判断输入的合法性error('Wrong  Input');
end
u=ones(n,1);q=max(abs(u));
k=0;
while truey=A\u;%解出yq0=q;q=max(abs(y));if abs(q-q0)<tol|k>1e4 %设置停机条件break;endu=y/q;k=k+1;
end
er=abs(q0-q);
a=1/q;

注:反幂法在计算量较大的情况下,貌似计算不出结果,考考你,找出问题,请留言。个人猜测:y=A\u处会出问题,考虑到矩阵A的奇异情况。

Matlab实现 乘幂法反幂法相关推荐

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

    matlab数值分析  矩阵特征值与特征向量的计算 1改进乘幂法 function [t,y]=eigIPower(A,v0,ep) [tv,ti]=max(abs(v0)); lam0=v0(ti) ...

  2. 乘幂法计算矩阵主特征值和特征向量-Matlab实现

    文章目录 1.前言 2.方法介绍 3.算法步骤 4.数值实验 5.总结 6.Matlab代码 1.前言 乘幂法主要用于求实矩阵按模最大的特征值(主特征值)和相应特征向量.本文通过Matlab解决实际例 ...

  3. Matlab:序列分析法MATLAB代码

    Matlab:序列分析法MATLAB代码 目录 输出结果 设计代码 输出结果 更新-- 设计代码 ###下面所有带代码中的n值需要以自己输入的数据为准###1.简单一次滑动平均法预测MATLAB程序代 ...

  4. matlab窗函数带通滤波器,Matlab结合窗函数法设计数字带通FIR滤波器

    Matlab结合窗函数法设计数字带通FIR滤波器 课程设计任务书学生姓名: 专业班级: 通信工程 指导教师: 工作单位: 信息工程学院 题 目:利用 Matlab 仿真软件系统结合窗函数法设计一个数字 ...

  5. 频域串联滞后校正matlab,自动控制原理实验七基于MATLAB控制系统频域法串联校正设计.doc...

    自动控制原理实验七基于MATLAB控制系统频域法串联校正设计.doc 实验七 基于MATLAB控制系统频域法串联校正设计 一.实验目的 (1)对给定系统设计满足频域或时域指标的串联校正装置: (2)掌 ...

  6. 乘幂法、复化梯形以及二分法求特征值【相关例题python实现】

    目录 写在前面的话 乘幂法 例题 代码 结果 二分法求对称三对角矩阵特征值 代码 结果 复化梯形 例题 代码 结果 写在前面的话 考试结束,<数值计算方法>编程作业发出来给大家看看~ 乘幂 ...

  7. matlab求刚度,matlab直接刚度法计算结构频率

    matlab直接刚度法计算结构频率 syms E I K r l x cja1=-E*I*K*r^3*(cos(x)*sinh(x)+sin(x)*cosh(x))jc1=-E*I*K*r^2*sin ...

  8. matlab 龙格-库塔 法求解常微分方程

    最近学习分室模型,里面碰到了用matlab 龙格-库塔 法求解常微分方程 研究了一阵子终于明白到底怎么实现了: 1. matlab 新建.m文件,编写龙格-库塔法求解函数 function [x,y] ...

  9. MATLAB差影法的人体姿态识别

    GUI框架源码: https://download.csdn.net/download/weixin_44748303/12682348 该课题为基于MATLAB差影法的人体姿态识别.带有一个GUI可 ...

最新文章

  1. Linux I2C工具查看配置I2C设备【转】
  2. javascript迭代器_JavaScript符号,迭代器,生成器,异步/等待和异步迭代器-全部简单解释...
  3. 使用OpenCV与sklearn实现基于词袋模型的图像分类预测与搜索
  4. 使用WakeLock将Android应用程序保持后台唤醒
  5. request.getParameterValues与request.getParameter的区别 想搞清楚为什么前者返回的是数组...
  6. (JavaWeb)Filter过滤器
  7. 美团笔试1--螺旋矩阵
  8. sql每一个join都要加on
  9. 重大变动:OpenCV拟更改授权协议!
  10. Java中的编译、反编译和反编译工具全家桶分享
  11. linux服务器ftp连接失败的原因,错误:无法与 SFTP 服务器建立 FTP 连接
  12. 渗透测试专业术语——防守篇
  13. Office办公软件测试题
  14. Android 面试题(答案最全) 转:http://www.jobui.com/mianshiti/it/android/2682/
  15. Direct3D 的进化史
  16. 使用计算机管理文件教后反思,《管理计算机中的文件》教学设计
  17. Composer 简单介绍
  18. SIMD、SIMD、SIMT、MISD、MIMD详解与比较
  19. java把分数区分ABCD四个等级_java编程输入成绩的英文等级abcd输出对应的中文等级优秀良好及格不及格...
  20. WebLogic 10.3.6.0 部署安装问题(简装版)

热门文章

  1. ie浏览器下载门户网站的excel文档出现乱码现象的解决办法
  2. 每日新闻:华为发布5大AI战略;支付宝建议苹果用户调低免密支付额度;胡润百富榜:马云居首马化腾第三;微软重新推送新版Win10...
  3. 热电偶和热电阻的区别
  4. CodeForces - 3A Shortest path of the king【水题】
  5. 有没有在线抠图网页版?这两个方法轻松制作
  6. 联合 EMQ 发布云原生物联网消息服务联合解决方案,云上轻松构建 IoT 应用
  7. 什么是STN黄绿模、什么是蓝模、什么是灰模
  8. kubernetes——安全认证
  9. Aspose.word java 实现word转pdf
  10. 乔布斯初心与终点:疯子一样的不知妥协