【原创】矩阵Given变换Matlab代码_风雪夜归人_新浪博客

本帖由坛原创,原帖参见,转载请注明Matlab技术论坛原创

 

如果有任何问题、建议,或者更多资源、代码、视频,欢迎您访问专业的Matlab技术交流平台——Matlab技术论坛

 

Wish my dear girl,healthy pretty and oveliness

 

2010-1-16 23:42:17

 

function [R,y]=givens(x,i,j)
% 求解标准正交的Given变换矩阵R,使用Rx=y,其中y(j)=0,y(i)=sqrt(x(i)^2+x(j)^2)
%
% 参数说明
% x:需要进行Givens变换的列向量
% i:变为sqrt(x(i)^2+x(j)^2)的元素下标
% j:变为0的元素的下标
% R:Givens变换矩阵
% y:Givens变换结果
%
% 实例说明
% x=[1 3 5 9 6]'; % 将3等效到9上
% [R,y]=givens(x,4,2) % 注意3的下标为2,9的下标为4
% R*x-y % 验证Rx=y
%
% 关于Givens变换说明
% 1.Givens矩阵是标准正交矩阵,也叫平面旋转矩阵,它是通过坐标旋转的原理将元素j的数值等效到元素i上
% 2.Givens变换每次只能将一个元素变为0,而Householder变换则一次可以将任意个元素变为0
% 3.Givens变换常用于将矩阵A变为对角阵
%
% by dynamic of Matlab技术论坛
% see also

% contact me
% 2010-01-16 23:31:05
%
xi=x(i);
xj=x(j);
r=sqrt(xi^2+xj^2);
cost=xi/r;
sint=xj/r;
R=eye(length(x));
R(i,i)=cost;
R(i,j)=sint;
R(j,i)=-sint;
R(j,j)=cost;
y=x(:);
y([i,j])=[r,0];

已投稿到:
郑重声明:资讯 【【原创】矩阵Given变换Matlab代码_风雪夜归人_新浪博客】由 发布,版权归原作者及其所在单位,其原创性以及文中陈述文字和内容未经(企业库qiyeku.com)证实,请读者仅作参考,并请自行核实相关内容。若本文有侵犯到您的版权, 请你提供相关证明及申请并与我们联系(qiyeku # qq.com)或【在线投诉】,我们审核后将会尽快处理。
—— 相关资讯 ——