close all,clear all
x=2:45:722; y=[0.0175902,0.30468,0.64733,0.94591,0.9588,0.80693,0.66372,0.55691,0.47833,0.41902,0.37288,0.33603,0.30592,0.28085,0.25965,0.24147,0.22571];
xi=2:20:722;
xii=2:0.1:722;
yi=interp1(x,y,xi,'spline');
m=[1 1010 1.01*1e6];%współczynniki wielomianu w mianowniku
s=j*2*pi*xii;
hs=1000*(s+10)./polyval(m,s);
A=abs(hs);
plot(x,y,'or',xi,yi,'o',xii,A),grid,legend("Pomiar","Interpolacja","Teoria"),xlabel("f [hz]"),ylabel("A");
%aproksymacja
a=polyfit(x,y,8);
y2=polyval(a,x);
figure(2),plot(x,y,'o',xii,A,x,y2),grid,legend("Pomiar","Aproksymacja","Teoria"),xlabel("f [hz]"),ylabel("A");