% Aproksymacja - liniowa regresja danych
clear
clc
close all

% Dane pomiarowe
x=[0:0.1:10]';
[m n]=size(x);
 for i=1:m
     y(i,1)=20+20*randn(1,1)+i;
 end;
  plot(x,y,'.'); hold on;
 
%  Funkcja y=ax przyblizajaca najlepiej w sensie sumy kwadratow
% a=x\y;
a=([x, ones(m,1)]\y)';
y1a=a(1)*x+a(2);
plot(x,y1a);
