MATLAB幂法原点平移法,幂法及其MATLAB程序.doc
幂法及其MATLAB程序
5.2 幂法及其MATLAB程序
5.2.2 幂法的MATLAB程序
用幂法计算矩阵的主特征值和对应的特征向量的MATLAB主程序
function [k,lambda,Vk,Wc]=mifa(A,V0,jd,max1)
lambda=0;k=1;Wc =1; ,jd=jd*0.1;state=1; V=V0;
while((k<=max1)&(state==1))
Vk=A*V; [m j]=max(abs(Vk)); mk=m;
tzw=abs(lambda-mk); Vk=(1/mk)*Vk;
Txw=norm(V-Vk); Wc=max(Txw,tzw); V=Vk;lambda=mk;state=0;
if(Wc>jd)
state=1;
end
k=k+1;Wc=Wc;
end
if(Wc<=jd)
disp('请注意:迭代次数k,主特征值的近似值lambda,主特征向量的近似向量Vk,相邻两次迭代的误差Wc如下:')
else
disp('请注意:迭代次数k已经达到最大迭代次数max1,主特征值的迭代值lambda,主特征向量的迭代向量Vk,相邻两次迭代的误差Wc如下:')
end
Vk=V;k=k-1;Wc;
例5.2.2 用幂法计算下列矩阵的主特征值和对应的特征向量的近似向量,精度.并把(1)和(2)输出的结果与例5.1.1中的结果进行比较.
(1); (2);(3);(4).
解 (1)输入MATLAB程序
>>A=[1 -1;2 4]; V0=[1,1]';
[k,lambda,Vk,Wc]=mifa(A,V0,0.00001,100),
[V,D] = eig (A), Dzd=max(diag(D)), wuD= abs(Dzd- lambda), wuV=V(:,2)./Vk,
运行后屏幕显示结果
请注意:迭代次数k,主特征值的近似值lambda,主特征向量的近似向量Vk,相邻两次迭代的误差Wc如下:
k = lambda = Wc =
33 3.00000173836804 8.691862856124999e-007
Vk = V = wuV =
-0.49999942054432 -0.70710678118655 0.44721359549996 -0.89442822756294
1.00000000000000 0.70710678118655 -0.89442719099992 -0.89442719099992
Dzd = wuD =
3 1.738368038406435e-006
由输出结果可看出,迭代33次,相邻两次迭代的误差Wc 8.69 19e-007,矩阵的主特征值的近似值lambda3.000 00和对应的特征向量的近似向量Vk (-0.500 00,1.000 00, lambda与例5.1.1中的最大特征值近似相等,绝对误差约为1.738 37e-006,Vk与特征向量 的第1个分量的绝对误差约等于0,第2个分量的绝对值相同.由wuV可以看出,的特征向量V(:,2) 与Vk的对应分量的比值近似相等.因此,用程序mifa.m计算的结果达到预先给定的精度.
(2) 输入MATLAB程序
>>B=[1 2 3;2 1 3;3 3 6]; V0=[1,1,1]';
[k,lambda,Vk,Wc]=mifa(B,V0,0.00001,100), [V,D] = eig (B),
Dzd=max(diag(D)), wuD= abs(Dzd- lambda), wuV=V(:,3)./Vk,
运行后屏幕显示结果
请注意:迭代次数k,主特征值的近似值lambda,主特征向量的近似向量Vk,相邻两次迭代的误差Wc如下:
k = lambda = Wc = Dzd = wuD =
3 9 0 9 0
Vk = wuV =
0.50000000000000 0.81649658092773
0.50000000000000 0.81649658092773
1.00000000
MATLAB幂法原点平移法,幂法及其MATLAB程序.doc相关推荐
- 基于matlab的车牌识别系统程序,基于matlab的车牌识别系统的设计(附程序).doc
基于matlab的车牌识别系统的设计(附程序).doc 1车牌识别系统的设计1.摘要:汽车牌照自动识别系统是制约道路交通智能化的重要因素,包括车牌定位.字符分割和字符识别三个主要部分.本文首先确定车辆 ...
- matlab 计算变异系数,变异系数法求权重matlab 代码
利用matlab编程,很容易根据变异系数法,求得多指标综合评价模型的权重. 代码如果有不懂的地方,可以联系我. 变异系数法求权重matlab 代码 clear;clc; [data1,header1] ...
- matlab 计算变异系数,变异系数法求权重matlab代码
<变异系数法求权重matlab代码>由会员分享,可在线阅读,更多相关<变异系数法求权重matlab代码(1页珍藏版)>请在读根文库上搜索. 1.变异系数 法求权重 matlab ...
- 基于AR模型的功率谱matlab,基于AR模型法功率谱估计的Matlab实现
2005iF 12月 笫 21卷第 6期 武警工程学 院学报 jOURNAL OF ENGG COLLEGE OF ARMED POLICE FORCE DeC.2005 Vo1.21 No.6 [计 ...
- matlab符号函数绘图法_《MATLAB符号运算及其应用》(黄忠霖著).pdf
<MATLAB符号运算及其应用>(黄忠霖著).pdf [General Information] 书名=MATLAB符号运算及其应用 作者=黄忠霖,黄京编著 页数=425 出版社=北京市: ...
- 用matlab编程节点电压法求电路,MATLAB在电路中的应用
<MATLAB在电路中的应用>由会员分享,可在线阅读,更多相关<MATLAB在电路中的应用(59页珍藏版)>请在人人文库网上搜索. 1.MATLAB应用(三) Matlab在电 ...
- matlab转换为部分分式,【MATLAB用部分分式展开法资讯】MATLAB用部分分式展开法足球知识与常识 - 足球百科 - 599比分...
MATLAB中对多项式进行部分分式展开. 用MATLAB进行部分分式展开 MATLAB有1个命令用于求B(s)/A(s)的部分分式展开式. 设s的有理分式为 式中 (i=)和(j=)的某些值可能为零. ...
- matlab中值滤波法算法,基于MATLAB中值滤波算法的优化与实现
总第238期2014年第4期 舰 船 电 子 工 程 Ship Electronic Engineering Vol.34No.437 基于 MATLAB中值滤波算法的优化与实现* 赵建春 刘力源 ( ...
- MATLAB之牛顿下山法
MATLAB之牛顿下山法 算法原理 matlab程序 算法原理 上一篇博客,我介绍了牛顿法迭代法,接下来我就们接着讲解一下什么是牛顿下山法. 一.迭代公式 在牛顿迭代过程中,若满足单调性|f(x(k+ ...
最新文章
- 损失函数(Loss Function) -1
- 配置ubuntu - tftp server服务器步骤
- 3. std::string::size_type
- GPS-nmealib学习
- 【FPGA】相关介绍
- winform改变启动界面
- electron 渲染进程调用主进程_万物皆可快速上手之Electron(第一弹)
- 【JY】基于OpenSees和SAP2000静力动力计算案例分析
- python之class
- Bert预训练新法则!
- SecureCRT教程——如何使用SecureCRT命令窗口将命令发送到多个会话
- 7-5 盲盒包装流水线
- 共有24款C/C++ 地理信息系统GIS开源软件
- 教你做Ghost系统盘--Windows2000、XP、2003系统万能Ghost全攻略New!!!!!!
- 字符串逆序输出----Python
- fairuse4wm源码_使用FairUse Wizard 2.9将DVD电影直接转换为AVI
- 傻白入门芯片设计,典型的2D/2D+/2.5D/3D封装技术(六)
- 弘辽科技:拼多多销量排名是怎么算的?影响因素有哪些?
- 吃鸡游戏跑毒和倍镜开发
- window XP安装
热门文章
- 【RFID】射频滤波器的设计
- 高数_第6章无穷级数_幂级数_欧拉公式
- 学习Linux命令(20)
- 2022-2-24 MIT 6.268 lab:Booting PC Part 3: The Kernel —— exercise7-8
- pandoc实现LaTeX转word+公式格式批量转换
- P3353 在你窗外闪耀的星星————树状数组,前缀和
- leetcode系列-104.二叉树的最大深度
- 一个早起的程序员简介
- minimum required php,今天部署zabbix所遇到问题
- xnawindowsphCocos2d-x for XNA游戏引擎全面支持WP7