% estudo de uma possivel função de resposta dos inclinómetros

b = 0:90;
x = -90:90;
a = polyfit([-90 90], [-90 90], 3)

a(3) = -.2
plot(-x, -a(1)*x.^3*2-a(2)*x.^2+a(3)*x)
grid on


x = -90:90;
a = polyfit([-90 90], [-90 90], 3)

a = [ 0 0 0 0 .1 0];
a(1) = (90-a(5)*90)/90^5;
y = a(1).*x.^5+a(2).*x.^4+a(3).*x.^3+a(4).*x.^2+a(5).*x+a(6);
plot(x, y)
grid on


x = 0:90;
a=polyfit([0 5 90], [0 10 210], 1);

y = polyval(a, x); %a(1).*x+a(2);
plot(x, y)
grid on
