Advertisement
Guest User

Untitled

a guest
Jan 25th, 2015
165
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.73 KB | None | 0 0
  1. clear all;
  2. close all;
  3. clc;
  4.  
  5. %przedział
  6. a=-2;
  7. b=2;
  8.  
  9. punkty=a:0.1:b; %w tych punktach szukamy co dzieje się z funkcją, zwiększajac zwiekszamy dokladnosc
  10.  
  11. %ilosc wezlow
  12. ilosc=16;
  13. m=9;
  14.  
  15.  
  16. %wezly rownolegle
  17. n=ilosc;
  18. for i=1:n
  19. x(i)=a+((b-a)/(ilosc))*(i-1);
  20. end
  21.  
  22. %wezly czebyszewa
  23. % n=ilosc;
  24. % for i=1:ilosc
  25. % x(i)=((a-b)*cos(((2*i-1)*pi)/(2*n)))/2+(a+b)/2;
  26. % end
  27.  
  28. %dla funkcji:
  29. %y=exp(x);
  30. y=1./(1+x.^2);
  31. %tworznie macierzy
  32. for i=1:1:ilosc
  33. for j=1:1:m
  34. D(i,j)=x(i)^(j-1);
  35. end
  36. end
  37.  
  38. %Ze wzoru
  39. Lewa=D'*D;
  40. Prawa=D'*y';
  41.  
  42. %Eliminacja Gaussa
  43. U=[Lewa Prawa]; %Tworzymy macierz rozszerzoną
  44. [W K]=size(U); %Pobieramy rozmiar macierzy rozszerzonej
  45.  
  46. for i=1:W
  47. for j=i:W
  48. if abs(U(j,i))==max(abs(U(i:W,i)))
  49. tmp1=U(i,:);
  50. U(i,:)=U(j,:);
  51. U(j,:)=tmp1;
  52. end
  53. end
  54. end
  55.  
  56. for j=1:W-1
  57. for i=j+1:W
  58. U(i,j:K)=U(i,j:K)-U(i,j)/U(j,j)*U(j,j:K); %U(i,j)/U(j,j) oblicza nam wpsółczynik
  59. end %prze jaki należy pomnożyć elementy z wiersza
  60. end %np. pierwszego by później odjąć je od np.wiersza drugiego
  61. %w celu wyzerowania danego elementu
  62. a=zeros(W,1);
  63. for i=W:-1:1
  64. a(i)=(U(i,K)-U(i,i+1:W)*a(i+1:W))/U(i,i); %Obliczanie niewiadomych
  65. end
  66.  
  67. a=a';
  68. %układamy w odpowiedniej kolejności
  69. for i=1:9
  70. wsp(1,i)=a(1,10-i);
  71. end
  72.  
  73. plot(punkty, polyval(wsp,punkty), 'o', 'LineWidth', 2);
  74. hold on;
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement