1 disp("Wspolrzedne wezlow");
2
3 w0=[-5,5]
4 w1=[-3,1]
5 w2=[1,2]
6 w3=[2,-1]
7 w4=[-2,4]
8 w5=[4,-1]
9 X=[w0(1,1)^0, w0(1,1)^1, w0(1,1)^2, w0(1,1)^3, w0(1,1)^4, w0(1,1)^5;
10 w1(1,1)^0, w1(1,1)^1, w1(1,1)^2, w1(1,1)^3, w1(1,1)^4, w1(1,1)^5;
11 w2(1,1)^0, w2(1,1)^1, w2(1,1)^2, w2(1,1)^3, w2(1,1)^4, w2(1,1)^5;
12 w3(1,1)^0, w3(1,1)^1, w3(1,1)^2, w3(1,1)^3, w3(1,1)^4, w3(1,1)^5;
13 w4(1,1)^0, w4(1,1)^1, w4(1,1)^2, w4(1,1)^3, w4(1,1)^4, w4(1,1)^5;
14 w5(1,1)^0, w5(1,1)^1, w5(1,1)^2, w5(1,1)^3, w5(1,1)^4, w5(1,1)^5];
15
16 disp("Macierz X");
17 disp(X);
18 Y=[w0(1,2); w1(1,2); w2(1,2); w3(1,2); w4(1,2); w5(1,2)];
19 disp("Macierz Y");
20 disp(Y);
21 A=inv(X)*Y;
22 disp("Macierz wspolczynnikow");
23 disp(A);
24 fprintf(' W(x)=%f*x^5 + %f*x^4 + %f*x^3 + %f*x^2 + %f*x + %f \n',A(6,1), A(5,1), A(4,1),
A(3,1), A(2,1), A(1,1));
25
26 x=[-5:1:5];
27 y=zeros(1,11);
28 for i=1:11
29 y(i)=A(1,1)+A(2,1)*x(i)+A(3,1)*x(i)^2+A(4,1)*x(i)^3+A(5,1)*x(i)^4+A(6,1)*x(i)^5;
30 end
31
32 plot(x, y);
33 grid on;
34 hold on;
35 vector_x=[w0(1,1), w1(1,1), w2(1,1), w3(1,1), w4(1,1), w5(1,1)];
36 vector_y=[0,0,0,0,0,0];
37 for i=1:6
38 vector_y(1,1)=vector_y(1,1)+A(i,1)*w0(1,1)^(i-1);
39 vector_y(1,2)=vector_y(1,2)+A(i,1)*w1(1,1)^(i-1);
40 vector_y(1,3)=vector_y(1,3)+A(i,1)*w2(1,1)^(i-1);
41 vector_y(1,4)=vector_y(1,4)+A(i,1)*w3(1,1)^(i-1);
42 vector_y(1,5)=vector_y(1,5)+A(i,1)*w4(1,1)^(i-1);
43 vector_y(1,6)=vector_y(1,6)+A(i,1)*w5(1,1)^(i-1);
44 end
45 plot(vector_x,vector_y,'rx');
46 plot(w0(1,1),w0(1,2),'--gs');
47 plot(w1(1,1),w1(1,2),'--gs');
48 plot(w2(1,1),w2(1,2),'--gs');
49 plot(w3(1,1),w3(1,2),'--gs');
50 plot(w4(1,1),w4(1,2),'--gs');
51 plot(w5(1,1),w5(1,2),'--gs');
52 title('Interpolacja wielomianowa');
53 ylabel('Os y');
54 xlabel('Os x');