givens.m

MATLAB科学计算与工程分析源代码源程序4

function [P,y]=Givens(x)
% 此函数用来将一个n维列向量化为:y=[norm(x) 0 ... 0]'
% 输出参数P为变换矩阵，即y=P*x
% 编写日期：2007-5-13

n=length(x);
P=eye(n);
for i=n:-1:2
    [G,x(i-1:i)]=planerot(x(i-1:i));
    P(i-1:i,:)=G*P(i-1:i,:);
end
y=x;