%Łolsz Nikodem Klimala
%Labolatorium teori sygnałów
%Torre Verona

function jus ()
t = -2:0.01:2;
%y = xpi(t, -1/2, 3);
%y2 = xpi(t, 0, 1/2);

y3 = xpi(t, -1/2, 3)-xpi(t, 0, 1/2);

w1=xpi(t,-1,2) - xpi(t,1,2);
w2=xpi(t,-3/2,1)-xpi(t,0,2)+xpi(t,3/2,1);
w3=xpi(t,-3/2,1)-xpi(t,-1/2,1)+xpi(t,1/2,1)-xpi(t,3/2,1);

w1_calka = calka(w1.*w1,0.001),
w2_calka= calka(w2.*w2, 0.001),
w3_calka= calka(w3.*w3, 0.001),

a1 = calka(y3.*w1,0.001) ./ w1_calka,
a2 = calka(y3.*w2,0.001) ./ w2_calka,
a3 = calka(y3.*w3,0.001) ./ w3_calka,
a = [a1,a2,a3],

xa(1)=calka(a1.*xpi(t,-1,2),0.001) + calka(a2.*xpi(t,-3/2,1),0.001) + calka(a3.* xpi(t, -3/2 , 1), 0.001);
xa(2) = calka(a1.* xpi(t, -1, 2), 0.001) + calka(a2.* xpi(t,-1/2, 1),0.001).*(-1) + calka(a3.*xpi(t, -1/2, 1),0.001).*(-1);
xa(3) = calka(a1.*xpi(t, 1, 2),0.001).*(-1) + calka(a2.*xpi(t, 1/2, 1),0.001).*(-1) + calka(a3.*xpi(t, 1/2,1),0.001);
xa(4) = calka(a1.*xpi(t, 1,2),0.001).*(-1) + calka(a2.* xpi(t, 3/2, 1),0.001) + calka(a3.*xpi(t,3/2,1),0.001).*(-1);

xa


grid on;
figure(1)
  subplot(5,1,1);plot(t,y3,'r');hold on; grid on;
  subplot(5,1,2);plot(t,w1);hold on; grid on;
  subplot(5,1,3);plot(t,w2);hold on; grid on;
  subplot(5,1,4);plot(t,w3);hold on; grid on;
  xapp=  xa(1).* xpi(t, -3/2,1) + xa(2).*xpi(t, -1/2,1) + xa(3).*xpi(t, 1/2 ,1) + xa(4).*xpi(t,3/2,1);

  subplot(5,1,5);plot(t, xapp);hold on; grid on;
 grid on;