急求matlab复数矩阵QR分解代码我知道matlab本身有qr函数 但是我想知道利用householder变换递归实现的原理。
E称为置换矩阵，满足 A*E = Q*R
完全可以，是一样的调用QR。
>> A=rand(5)
A =
0.5828 0.2259 0.2091 0.5678 0.4154
0.4235 0.5798 0.3798 0.7942 0.3050
0.5155 0.7604 0.7833 0.0592 0.8744
0.3340 0.5298 0.6808 0.6029 0.0150
0.4329 0.6405 0.4611 0.0503 0.7680
>> [Q,R]=qr(A)
Q =
-0.5598 0.8235 -0.0747 0.0193 0.0492
-0.4068 -0.1939 0.5250 -0.6391 -0.3359
-0.4952 -0.3400 -0.2889 0.4805 -0.5699
-0.3208 -0.2913 -0.6630 -0.4618 0.3993
-0.4159 -0.2894 0.4425 0.3833 0.6329
R =
-1.0410 -1.1752 -1.0696 -0.8846 -1.1138
0 -0.5246 -0.4996 0.1034 -0.2410
0 0 -0.2899 -0.0200 0.2063
0 0 0 -0.7274 0.5207
0 0 0 0 -0.0882
至于递归实现，只要对第一列进行消去后再递归就行了
function [Q,R]=qrhs(A)
n=size(A,1);
R=A;
Q=eye(n);
for i=1:n-1
x=R(i:n,i);
y=[1;zeros(n-i,1)];
Ht=householder(x,y);
H=blkdiag(eye(i-1),Ht);
Q=Q*H;
R=H*R;
end
求教matlab画复数函数图象的问题
复数(a+bj)可以写成如下幅相形式的 y = A*exp(j*oumige)其中A = sqrt(a^2+b^2);oumige = tan^(-1)(b\/a)你这个题目要求画幅(amplitude)相(phase)特性 所以需要首先将复数化成幅相形式 t = linspace(0,4,50);f = (1+0.25j).*t-2;amp = abs(f);pha = angle(f)*180.\/...
怎样用matlab进行复数运算
例如想用matlab将一个带变量的复数式(5+i*b)\/(3-2i*a)整理为实部+虚部的形式。则可以用如下指令：syms a b real z=(5+i*b)\/(3-2i*a);simple([real(z),imag(z)])ans = [ (15-2*b*a)\/(9+4*a^2), (10*a+3*b)\/(9+4*a^2)]...
matlab钟输入一个复数,如何输出复数的模和相角?
你好，这是最简单的求复数模和相角的程序。> x=1+1*i x = 1.0000 + 1.0000i >> abs(x)ans = 1.4142 >> angle(x)ans = 0.7854
matlab求解复数方程
>> solve('x^2=-1')ans i -i 第二次 >> solve('x^2=-1')ans = sqrt(-1)-sqrt(-1)注意ans是<2*1 sym>型 我的版本是：>> ver --- MATLAB Version 7.6.0.324 (R2008a)
下图中求复数的辐角运算在matlab中怎么编
设gamma(d)是一个复数值，用 f_ML = -angle(gamma(d))\/(2*pi);
用MATLAB编程解复数方程组z1+2*z2=1+i 3*z1+I*z2=2-3*i
用MATLAB编程解复数方程组，可以用solve（）函数命令。求解过程如下：>> syms z1 z2 >> [z1,z2]=solve(z1+2*z2==1+i,3*z1+i*z2==2-3*i)运行结果 z1 =1 - i z2 =i
Matlab求复数所有根
p(r)=r^5-c p=[1,0,0,0,0,-c]; %从前到后依次为r的5次方，四次方，三次方，二次方，一次方，零次方系数 r=roots(p) %求根
matlab中求解含复数方程的程序怎么写
复数也是一样的 先把格式转化下 format rat 再用solve函数，可以实现的
matlab中复数j怎么输入
按照复数形式输入。matlab中用i或j来表示虚部，创建复数时，可以直接按照复数形式输入。复数是指把形如a+bi（a，b均为实数）的数称为复数，其中a称为实部，b称为虚部，i称为虚数单位。
MATLAB的solvepde函数能求解复数方程么
solvepde Solve PDE specified in a PDEModel collapse all in page Syntax result = solvepde(model)result = solvepde(model,tlist)Description example result = solvepde(model) returns the solution to the stationary PDE represented in model. A stationary PDE has the property model.IsTime...
唱环心灵： 从你提问的水平来判断,你可能根本不知道QR分解的实现机理,那么比较务实的方法是先去了解一下镜像变换(Householder变换)的原理和实现方法,然后可以参考最简单的level 2 QR的实现 给你一些fortran代码,自己想办法看懂 QR分解 http:...
阜阳市13943672377： 急求matlab复数矩阵QR分解代码我知道matlab本身有qr函数 但是我想知道利用householder变换递归实现的原理. - ？
唱环心灵： 实Householder变换和复Householder变换没有本质区别,只不过是把H=I-2ww^T改成H=I-2ww^H 至于递归实现,只要对第一列进行消去后再递归就行了
阜阳市13943672377： 急求,matlab中, 已知矩阵A,已完成对A的QR分解,下一步求A的特征值和特征向量,程序怎么编写? - ？
唱环心灵： 楼主的问题是自己写程序完成矩阵的QR分解,既然是迭代实现QR分解,就与矩阵论中说的计算特征值和特征向量的方法有些区别了.大体的步骤应该是首先将矩阵化成双对角矩阵,然后追赶计算特征值和特征向量,程序代码可以参考 徐士良编的 常用数值算法 c语言描述
阜阳市13943672377： matlab QR分解用什么算法实现的 - ？
唱环心灵： QR分解法是目前求一般矩阵全部特征值的最有效并广泛应用的方法,一般矩阵先经过正交相似变化成为Hessenberg矩阵,然后再应用QR方法求特征值和特征向量.它是将矩阵分解成一个正规正交矩阵Q与上三角形矩阵R,所以称为QR分解法,与此正规正交矩阵的通用符号Q有关.
阜阳市13943672377： 题目具体意思是:对一个矩阵A进行QR分解,只有唯一的一种情况吗?通过Matlab自带程序 qr分解,即 [Q R]=qr(A)上式Q和R分别为正交矩阵和上三角矩阵,... - ？
唱环心灵：[答案] 是唯一的! 矩阵A进行QR分解后,R中对角元素一定是实数,由于这个条件,使得它唯一!
阜阳市13943672377： 如何用matlab对矩阵进行正交分解 - ？
唱环心灵： 矩阵分解 (decomposition,factorization)是多半将矩阵拆解为数个三角形矩阵(triangular matrix). 依使用目的的不同 ,可分为三种矩阵分解法:1)三角分解法 (Triangular Factorization),2)QR 分解法 (QR Factorization),3)奇异值分 解...
阜阳市13943672377： matlab中的QR分解都能分解什么样的矩阵?? - ？
唱环心灵： % 正交分解(qr) 对于矩阵A(n*n),如果A非奇异,则存在正交矩阵Q和上三角矩阵R,使得A满足关系式 A=Q*R,并且当R的对交元都为正时,qr分解是唯一的.
阜阳市13943672377： 对矩阵x进行QR分解和LU分解,QR分解和LU分解是什么意思呢 - ？
唱环心灵： 为了求解线性方程组,我们通常需要一定的解法.其中一种解法就是通过矩阵的三角分解来实现的,属于求解线性方程组的直接法.在不考虑舍入误差下,直接法可以用有限的运算得到精确解,因此主要适用于求解中小型稠密的线性方程组. ...
阜阳市13943672377： matlab中如何用qr函数求特征值和特征向量,矩阵是mxn - ？
唱环心灵： 直到A收敛到对角块不超过2阶的分块上三角阵, 但实际上也并不该qr这个函数来实现.当然先不要考虑matlab了; A=Q'*A*Q, 对每个特征值各解一次方程组就行了.就讲这些, 即使你看不明白, 先回去复习一下线性代数, 单个的矩阵但不是方阵何谈特征值 即使是方阵, QR分解也不是直接用来求特征值和特征向量的.尽管求所有特征值和特征向量最重要的算法是QR算法, 数学上可以解释为反复做QR分解.至于求特征向量, 如果你一定想用qr, 那么可以反复迭代 [Q,R]=qr(A)
阜阳市13943672377： 用Matlab中的qr分解含字母的矩阵 - ？
唱环心灵： lu分解牵涉到选主元,符号本身不带大小,所以lu就不行了,除非你自己知道某些未知量的性质然后自己写一个.qr分解的计算当中倒是没有比较,不过我估计直接调qr还是不行的,如果不能直接用的话就自己写一段代码,反正块运算对符号也没用.我不知道你为什么有这些需求,我估计你本身的问题可能有别的办法来解决,符号计算大多数时候是没用的.
