Ajuste Linear SImples

Prévia do material em texto

//Ajuste Linear Simples Lucas Gomes dos Ramos 
 
clear 
clc 
disp('Ajuste Linear Simples Lucas Gomes dos Ramos') 
x=[1.3 3.4 5.1 6.8 8.0]; 
y=[2.0 5.2 3.8 6.1 5.8]; 
 
n=length(x); 
s1=0; 
s2=0; 
s3=0; 
s4=0; 
j=0; 
 
for i=1:n, 
 s1=s1+x(i); 
 s2=s2+((x(i)^2)); 
 s3=s3+(y(i)); 
 s4=s4+((x(i))*(y(i))); 
 j=j+1; 
end 
printf(" \n\n O valor de A11 é: %0.2f\n O valor de A12 é: %0.2f\n O valor de A21 é: %0.2f\n O valor de 
A22 é: %0.2f\n\n",s1,s2,s3,s4); 
b1=((n*s4)-(s1*s3))/((n*s2)-(s1)^2); 
mprintf(" O valor de Alpha 1 é:% 0.4f\n",b1); 
 
b0=((s3)-(s1*b1))/n; 
mprintf(" O valor de Alpha 2 é: %0.4f",b0); 
 
deff ("y=f(x)","y=b0+b1*x"); 
h=1:0.05:8.5; 
z=1:0.05:8.5; 
plot(x,y,'o',h,z)