Advertisement
Coriic

Untitled

Dec 5th, 2016
152
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.30 KB | None | 0 0
  1. %Zadanie 2
  2. imie='Kamil';
  3. nazwisko='Sobolewski';
  4. A=[0,1; -length(imie), -length(nazwisko)];
  5. B=[0;1];
  6. C=[1,0;0,1];
  7. h=1;
  8. t=sym('t');
  9. Aplus=expm(h*A);
  10. Bplus=double( int( expm( h*A )*B, t, 0, h ) );
  11. Cplus=C;
  12. figure();
  13. sim('model2');
  14. i=1:h:11;
  15. hold on;
  16. plot(i, x_d_2(:,2), 'o');
  17. plot(x_c_2.Time, x_c_2.Data(:,2));
  18. hold off;
  19.  
  20. %Zadanie 3
  21. K_0=1000; %kapitał początkowy
  22. n=5; %czas trwania inwestycji przy założeniu że jest wielkrotnością okresów kapitalizacji
  23. r=0.05; %nominalna stopa procentowa
  24. figure();
  25. for m=1:12
  26. hold on;
  27. K=K_0*(1+r/m)^(m*n);
  28. plot(m, K, 'o');
  29. end;
  30. hold off;
  31. xlabel('Liczba kapitalizacji odsetek w roku(m)');
  32. ylabel('Kapitał końcowy (K)');
  33.  
  34. %Zadanie 4
  35. liczba_iteracji=15; %maksymalna liczba iteracji
  36. epsilon_1=10^-6; %dokładność z jaką funkcja jest bliska zeru w wyznaczonym punkcie
  37. epsilon_2=10^-6; %dokładność określająca minimalną długość kroku w kolejnej iteracji
  38. x0=-4; %punkt startowy
  39. f=@(x)x^3-2*x^2-11*x+12; %funkcja f(x)
  40. g=@(x)3*x^2-4*x-11; %pochodna funkcji f(x)
  41. i=0; %licznik iteracji
  42. while i<liczba_iteracji
  43. i=i+1;
  44. x=x0-(f(x0)/g(x0));
  45. if(abs(f(x))<epsilon_1)
  46. break;
  47. end;
  48. if(abs(x-x0)<epsilon_2)
  49. break;
  50. end;
  51. x0=x;
  52. end;
  53. disp(x); %wypisuje znalezione miejsce zerowe na ekran
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement