regression lineaire sur deux subplot

Essaie simplement ceci :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
function test
 
x = 0:4;
nx = numel(x);
y1 = x+rand(1,nx);
y2 = x+rand(1,nx);
 
M = [ones(nx,1) x(:)] ;
K1 = M\y1(:);
K2 = M\y2(:);
 
equ{1} = sprintf('y = %.2f*x + %.2f',K1(2), K1(1));
equ{2} = sprintf('y = %.2f*x + %.2f',K2(2), K2(1));
 
figure(1)
clf
 
subplot(211)
plot(x,y1,'b*',x,K1(2)*x+K1(1),'r-');
hold on
text(0.5,3.5,equ{1})
axis equal tight
 
subplot(212)
plot(x,y2,'b*',x,K2(2)*x+K2(1),'r-');
hold on
text(0.5,3.5,equ{2})
axis equal tight

Voir : Ajustement aux moindres carrés de courbes et de surfaces