Advertisement
Guest User

Octave Nauka

a guest
Jun 24th, 2018
87
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Octave 1.58 KB | None | 0 0
  1. % Zadanie 1
  2. x1 = -5 : 0.01 : -0.01;
  3. x2 = 0 : 0.01 : 5;
  4.  
  5. f1 = @(x) cos(x.^2);
  6. f2 = @(x) 1 ./(x+1);
  7.  
  8. plot(x1, f1(x1), "r", x2, f2(x2), "r")
  9.  
  10.  
  11. % Zadanie 2
  12. C = randi(5, 5)
  13. C(C != 1)
  14.  
  15.  
  16. % Zadanie 3
  17. A = randi(5, 5)
  18. n = size(A, 1);
  19. B = triu(ones(n)) & fliplr(triu(ones(n)));
  20. sum(A(B))
  21.  
  22.  
  23. % Zadanie 4
  24.  x(1)=3;
  25.  y(1)=4;
  26.  z(1)=5;
  27.  for n=1:3
  28.    x(n+1)=3*x(n)+2*z(n)+1;
  29.    y(n+1)=3*x(n)+2*z(n)+2;
  30.    z(n+1)=4*x(n)+3*z(n)+2;
  31.  endfor
  32.  for n=1:4
  33.    disp([num2str(x(n)^2)," + ",num2str(y(n)^2)," = ",num2str(z(n)^2)])
  34.  endfor
  35.  
  36. ---------------------------------------------------------------
  37. % Zadanie 1
  38. x1 = -2 : 0.01 : 0;
  39. x2 = 0 : 0.01 : 2;
  40.  
  41. f1 = @(x) sin(x.^5);
  42. f2 = @(x) x.^3 - sqrt(x);
  43.  
  44. plot(x1, f1(x1), "r", x2, f2(x2), "r")
  45.  
  46.  
  47. % Zadanie 2
  48. C = rand(5) * 1000
  49. C(C < 10)
  50. C(C >= 10 & C < 100)
  51. C(C >= 100)
  52.  
  53.  
  54. % Zadanie 3
  55. A = randi(5, 5)
  56. n = size(A, 1);
  57. B = eye(n) | fliplr(eye(n))
  58. sum(A(B))
  59.  
  60.  
  61. % Zadanie 4
  62. function rek(c)
  63.     if (mod(c(end), 2) == 0)
  64.         c(end+1) = 0.5*c(end);
  65.     else
  66.         c(end+1) = 3*c(end)+1;
  67.     endif
  68.     c(end)
  69.     rek(c);
  70. endfunction
  71.  
  72. -----------------------------------------------------------
  73. n = 5; %liczba linii
  74. s = 1; %dlugosc podstawowa
  75. k = 4; %co ktora linia zmienionej dlugosci
  76. f = 0.5; %wspolczynnik (o ile ma byc dluzszy co k-ty)
  77. h = 2*pi/n; %krok
  78.  
  79. theta = linspace(0, 2*pi - h, n);
  80. r = s*ones(1,n);
  81. D = ones(1,n); %wektor zmieniajacy rozmiar
  82.  
  83. for i = 1:n
  84.   if mod(i,k) == 0
  85.     D(i) = f*D(i);
  86.   endif
  87. endfor
  88.  
  89. r = r.*D;
  90. [x,y] = pol2cart(theta,r);
  91.  
  92. hold on
  93. for i= 1:n
  94.   plot([x(i) 0], [y(i) 0]);
  95. endfor
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement