Advertisement
Guest User

Untitled

a guest
May 31st, 2016
69
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.67 KB | None | 0 0
  1. format short;
  2.  
  3. fun1 = @(x) cos(x)*cos(x)*cos(x) + 0.7*cos(x*x);
  4. n = 10;
  5. N = 2^n;
  6.  
  7. % uzmemo neki gusti interval
  8. interval = 0:2*pi/(3*N):2*pi-2*pi/N;
  9.  
  10. for k = 1:N
  11. x(k) = (2*(k-1)*pi) / N;
  12. f(k) = fun1( x(k) );
  13. end
  14.  
  15. s = prir_kub_splajn(x, f);
  16.  
  17. for i = 1 : length(interval)
  18. [y, ~, ~] = vrij_kub_splajna( interval(i), x, f, s);
  19. vrij(i) = y;
  20. end
  21.  
  22. % racunamo gresku
  23. for j = 1 : N
  24. [ y, ~, ~ ] = vrij_kub_splajna( x(j), x, f, s );
  25. greska( j ) = abs( f(j) - y );
  26. end
  27.  
  28. maks_greska = max( greska );
  29.  
  30. whitebg('black');
  31. fplot( fun1, [0 max(interval)], 'r');
  32. hold on
  33. plot( interval, vrij, 'w');
  34. hold off
  35. legend('funkcija', 'aproksimacija splajnom');
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement