matlab数值分析程序--高等数学,数值代数的matlab实现-文字版, matlab电子书, 和matlab 有关的电子书:

5.5.1 豪斯霍尔德方法及其MATLAB程序

5.5.1  豪斯霍尔德方法及其MATLAB程序

求初等反射矩阵,使得的第一个分量以外的其余的分量都为零的MATLAB主程序

function [xigema,rou,miou,P,PX]=Householder(X)

n=size(X);nX=norm(X,2);

xigema=nX*sign(X(1));

rou=xigema*(xigema+X(1));

miou=[xigema,zeros(1,n-1)]'+X,

E=eye(n,n); C=2*miou*(miou)';

P=E-C/(norm(miou,2)^2); PX=P*X;

 

例5.5.1  设向量,确定一个初等反射矩阵,使得的后两个分量为零.

解  输入MATLAB程序

>> X=[2 2 1]'; [xigema,rou,miou,P,PX]=Householder(X)

运行后屏幕显示结果

P =                                       PX =

   -0.6667   -0.6667   -0.3333                 -3.0000

   -0.6667    0.7333   -0.1333                 0.0000

   -0.3333   -0.1333    0.9333                 0.0000

欢迎转载,转载请注明来自一手册:http://yishouce.com/book/3/3050501.html
友情链接It题库(ittiku.com)| 版权归yishouce.com所有| 友链等可联系 admin#yishouce.com|粤ICP备16001685号-1