1楼：电灯剑客
你说的是把m阶复矩阵嵌入到2m阶实矩阵再算特征值的方法吧
这种情况最好是通过特征向量来判断
当然，如果没有特殊需求的话不如直接对复矩阵用qr算法，没必要嵌入到实矩阵
2楼：由江桐安荷
query取得iframe中元素的几种方法在iframe子页面获取父页面元素
**如下:$(
用qr方法怎样求矩阵的特征值？
3楼：匿名用户
对称矩阵用雅可比迭代法求特征值；一般矩阵 (非对称矩阵) 用qr迭代法求特征值。看下面例题。
4楼：匿名用户
a=q r
a1=r q=q1 r1
a2=r1 q1=q2 r2
...注意所有的a,a1,...相似(a1=rq=q^t a q =q^ a q),具有相同的特征值。
在一定条件下最终收敛到一个上三角阵，把对角线上的元拿出来就是特征值。事实上，因为a是对称矩阵，a1=q1^t a q 所以a1是对称阵(显然a1^t=a1)，以此类推，a2,a3...都是对称阵。
所以当a是对称阵时an收敛于对角阵（既是上三角阵又是对称阵，所以是对角阵），对角线上的元素就是特征值。
如何用qr算法求矩阵特征值？？
5楼：匿名用户
function l = rqrtz(a,m)%瑞利商位移的qr算法求矩阵全部特征
值%已知矩阵:a
%迭代步数:m
%求得的矩阵特征值:la = hess(a);
for(i=1:m)
n = size(a);
n = n(1,1);
u = a(n,n);
[q,r]=qr(a-u*eye(n,n));
a = r*q+u*eye(n,n);
l = diag(a);
end4.4 qr算 法 qr算法也是一种迭代算法,是目前计算任意实的非奇异矩阵全部特征值问题的最有效的方法之一.该方法的基础是构造矩阵序列 ,并对它进行qr分解.
由线性代数知识知道,若a为非奇异方阵,则a可以分解为正交矩阵q与上三角形矩阵r的乘积,即a=qr,而且当r的对角线元素符号取定时,分解式是唯一的. 若a为奇异方阵,则零为a的特征值.任取一数p不是a的特征值,则a-pi为非奇异方阵.
只要求出a-pi的特征值,就很容易求出a的特征值,所以假设a为非奇异方阵,并不妨碍讨论的一般性. 设a为非奇异方阵,令 ,对 进行qr分解,即把 分解为正交矩阵 与上三角形矩阵 的乘积 = 做矩阵 继续对 进行qr分解 并定义 一般地,递推公式为 qr算法就是利用矩阵的qr分解,按上述递推公式构造矩阵序列 .只要a为非奇异方阵,则由qr算法就完全确定 .
这个矩阵序列 具有下列性质. 性质1 所有 都相似,它们具有相同的特征值. 证明 因为 若令 ,则 为正交阵,且有 因此 与a相似,它们具有相同的特征值.
性质2 的qr分解式为 其中 证明 用归纳法.显然当k=1时,有 假设 有分解式 于是 因为 ,所以 因为 都是正交阵,所以 也是正交阵,同样 也是上三角形阵,从而 的qr分解式为 由前面的讨论知 .这说明qr算法的收敛性有正交矩阵序列 的性质决定.
定理1 如果 收敛于非奇异矩阵 为上三角形矩阵,则 存在并且是上三角形矩阵. 证明 因为 收敛,故下面极限存在 由于 为上三角形矩阵,所以 为上三角形矩阵.又因为 所以 存在,并且是上三角形矩阵.
定理2 (qr算法的收敛性)设a为n 阶实矩阵,且1) a的特征值满足: 2) ,其中 且设 有三角分解式 =lu(l为单位下三角阵,u为上三角阵),则由qr算法得到的矩阵序列 本质上收敛于上三角形矩阵.即 满足 当 当 的极限不一定存在 证明 因为 ,矩阵 决定 的收敛性.
又 我们利用 求 ,然后讨论 的收敛性. 由定理条件 得 令 其中 的(i,j)元素 为 于是 由假设,当i>j时, 故 设方阵x的qr分解式为 由 由 知,对充分大的 非奇异,它应有唯一的qr分解式 ,并且 于是 但上三角阵 的对角线元素不一定大于零.为此,引入对角矩阵 以便保证( )的对角线元素都是正数,从而得到 的qr分解式 由 的qr分解式的唯一性得到 从而 由于 ,所以 从而 其中 于是 因为 为上三角阵, 为对角阵,且元素为1或-1,所以 当 当 的极限不一定存在 例 用qr算法求矩阵 的特征值.
a的特征值为-1,4,1+2i,1-2i. 解 令 ,用施密特正交化过程将 分解为 将 与 逆序相乘,求出 用 代替a重复上面过程,计算11次得 由 不难看出,矩阵a的一个特征值是4,另一个特征值是-1,其他两个特征值是方程 的根.求得为
matlab中如何用qr函数求特征值和特征向量，矩阵是mxn
6楼：匿名用户
1.矩阵qr分解直接用函数qr就可以了。qr函数适用于不是方针的矩阵分解。
2.用法[q,r]=qr(a)得到q是mm矩阵，r是mn.
3.排列大小的可以采用sort函数。
具体情况建议打开matlab 帮助浏览器详细看qr函数的用法。
7楼：电灯剑客
先不要考虑matlab了, 先回去复习一下线性代数, 单个的矩阵但不是方阵何谈特征值
即使是方阵, qr分解也不是直接用来求特征值和特征向量的.
尽管求所有特征值和特征向量最重要的算法是qr算法, 数学上可以解释为反复做qr分解, 但实际上也并不该qr这个函数来实现.
当然, 如果你一定想用qr, 那么可以反复迭代[q,r]=qr(a); a=q'*a*q;
直到a收敛到对角块不超过2阶的分块上三角阵.
至于求特征向量, 对每个特征值各解一次方程组就行了.
就讲这些, 即使你看不明白, 我也不会继续回答了, 这纯粹是浪费时间.
matlab用qr方法怎么求特征值,把程序写出来，谢谢
8楼：匿名用户
function l = rqrtz(a,m)%qr算法求矩阵全部特征值
%已知矩阵:a
%迭代步数:m
%求得的矩阵特征值:l
a = hess(a);
for i=1:m
n = size(a);
n = n(1,1);
u = a(n,n);
[q,r]=qr(a-u*eye(n,n));
a = r*q+u*eye(n,n);
l = diag(a);
end------------------------------------
a=[0 5 0 0 0 0;1 0 4 0 0 0;0 1 0 3 0 0;0 0 1 0 2 0;0 0 0 1 0 1;0 0 0 0 1 0]
a =0 5 0 0 0 01 0 4 0 0 00 1 0 3 0 00 0 1 0 2 00 0 0 1 0 10 0 0 0 1 0>> rqrtz(a,50)
ans =
-3.2030
3.2030
-1.8837
1.8837
-0.6167
0.6167
>> eig(a)
ans =
-3.3243
3.3243
-1.8892
-0.6167
1.8892
0.6167
9楼：匿名用户
看matlab帮助文件，里边都有例程，电脑没有安装matlab，懒得安装了，但是我知道帮助文件里有这个例子
10楼：言清韵柯北
前面的函数是求特征值的函数，后面的是调用的例子和验证functionl=
rqrtz(a,m)
%qr算法求矩阵全部特征值
%已知矩阵:a
%迭代步数:m
%求得的矩阵特征值:la=
hess(a);
fori=1:mn=
size(a);n=
n(1,1);u=
a(n,n);
[q,r]=qr(a-u*eye(n,n));a=r*q+u*eye(n,n);l=
diag(a);
enda=[050
000;1040
00;010
300;0010
20;000
101;0000
10]a=
0500
0010
4000
0103
0000
1020
0001
0100
0010
>>rqrtz(a,50)
ans=
-3.2030
3.2030
-1.8837
1.8837
-0.6167
0.6167
>>eig(a)
ans=
-3.3243
3.3243
-1.8892
-0.6167
1.8892
0.6167
qr分解怎么求特征向量，求矩阵e的特征值和特征向量
11楼：馥馥幽襟披
楼主的问题是自己写程序完成矩阵的qr分解，既然是迭代实现qr分解，就与矩阵论中说的计算特征值和特征向量的方法有些区别了。大体的步骤应该是首先将矩阵化成双对角矩阵，然后追赶计算特征值和特征向量，程序**可以参考 徐士良编的 常用数值算...
(QR分解相似变换RQ迭代)反复循环，求出矩阵特征值。这
1楼 匿名用户 你这样问好像说这个矩阵肯定求不出，所以希望大家告诉你为什么，但是大家都不知道你怎么求的，求得过程是否有问题，所以需要你把计算过程列出来，让大家看看是不是你求措了 2楼 匿名用户 你问的这题太难了解，答不了 你好！简单矩阵 7，2 ， 2，3 用qr分解迭代，为啥求不出特征值？谢谢！ ...
怎么用matlab求矩阵的特征值和特征向量
1楼 天云一号 eig函数直接可以求特征值和特征向量 在matlab中，计算矩阵a的特征值和特征向量的函数是eig a ，常用的调用格式有5种 e eig a 求矩阵a的全部特征值，构成向量e。 v d eig a 求矩阵a的全部特征值，构成对角阵d，并求a的特征向量构成v的列向量。 v d eig...
怎样用Maple求矩阵的特征值和特征向量
1楼 匿名用户 with student linearalgebra b matrix 3 3 eigenvectors b eigenvectors b output list 红色字体表示特征根的重数 这个矩阵怎样用maple求出特征值和特征向量 2楼 匿名用户 求特征值 eigenvals a...