power.m

来自「提供一个直接搜索算法程序」· M 代码 · 共 23 行

function [lambda,V]= power(A,X,epsilon,max1)

%  NUMERICAL METHODS: Matlab Programs
% (c) 2004 by John H. Mathews and Kurtis D. Fink
%  Complementary Software to accompany the textbook:
%  NUMERICAL METHODS: Using Matlab, Fourth Edition
%  ISBN: 0-13-065248-2
%  Prentice-Hall Pub. Inc.
%  One Lake Street
%  Upper Saddle River, NJ 07458

lambda=0;cnt=0;err=1;iterating=1;state=iterating;

while ((cnt<=max1)&(state==iterating))
   Y=A*X;
   [m j]=max(abs(Y));
   c1=Y(j);
   cnt=cnt+1;
end

V=Y;