Skriv först in de uppmätta värdena för längd,l, och
periodtid, t:
l=[ 0.95 0.70 0.40 0.20 ];
t=[1.99 1.71 1.31 0.92];
% Rita t som funktion av l
plot (l,t,'o')
axis([0 1 0 2])
Visst liknar det en "rot-kurva"
Prova att rita t som funktion av roten ur l:
plot (l.^0.5,t,'o')
axis([0 1 0 2])
% Rita in de exakta värdena
fplot('2*pi*(x/9.81^0.5)',[0 1])
% Beräkna g i varje punkt
g = (l./t.^2) *(2*pi)^2
gmedel=sum(g)/4
Prova ocksä en log-log plot
plot (log10(l),log10(t),'o'); title "log-log"
Minsta kvadrat-anpassning
(Kap 7.3 i matlab-manualen) ax=b, loeses genom x=a\b
Vi har t=x * l.^2, dvs:
x= l'.^0.5\t';
% Om vi vet att t=2*pi*(l/g)^0.5 sa kan vi bestamma g:
glsf=(2*pi/x)^2
tfit=2*pi*(l/glsf).^0.5)
plot (tfit,l,'-')