tadejpetric

Untitled

Jan 10th, 2021
1,116
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. function x = potencna(A, zac, tol)
  2.   rayleigh = @(x) x'*A(x)/(x' * x);
  3.   x.vrednosti = [rayleigh(zac)];
  4.   x.vektorji = [zac];
  5.   k = 1;
  6.   while norm(A(x.vektorji(:, k)) - x.vrednosti(k)*x.vektorji(:, k)) >= tol
  7.     k = k+1;
  8.     x.vektorji(:,k) = A(x.vektorji(:, k-1));
  9.     x.vektorji(:,k) = x.vektorji(:, k)/norm(x.vektorji(:, k));
  10.     x.vrednosti(k) = rayleigh(x.vektorji(:, k));
  11.     if k == 501
  12.       break
  13.     endif
  14.      
  15.   endwhile
  16.  
  17. endfunction
  18.  
  19.  
  20.  
  21. function x = hotelling(A, zac, hotel_steps)
  22.   u = potencna(A, zac, eps);
  23.   x.vrednosti(1) = u.vrednosti(end);
  24.   x.vektorji(:,1) = u.vektorji(:,end);
  25.   for i = 2:hotel_steps+1
  26.     A = @(y) A(y) - x.vrednosti(i-1) * x.vektorji(:,i-1)*x.vektorji(:,i-1)'*y;
  27.     u = potencna(A, zac, eps);
  28.     x.vrednosti(i) = u.vrednosti(501);
  29.     x.vektorji(:,i) = u.vektorji(:,501);
  30.   endfor
  31. endfunction
  32.  
RAW Paste Data

Adblocker detected! Please consider disabling it...

We've detected AdBlock Plus or some other adblocking software preventing Pastebin.com from fully loading.

We don't have any obnoxious sound, or popup ads, we actively block these annoying types of ads!

Please add Pastebin.com to your ad blocker whitelist or disable your adblocking software.

×