% File matlab per la risoluzione dell'Esercizio 3.8

clear
syms x y

% pdf congiunta
fxy=-3/2*x^2+2*x*y;
xl=0; xu=1;
yl=1; yu=2;

% verifica che sia una pdf!
int(int(fxy,x,xl,xu),y,yl,yu)

% pdf di x e di y
fx=int(fxy,y,yl,yu);
fy=int(fxy,x,xl,xu);

% pdf di x condizionata a y
fxcy=fxy/fy;

% stima MEQM
xmeqm=int(x*fxcy,x,xl,xu)

% medie e covarianze
mx=int(x*fx,x,xl,xu);
my=int(y*fy,y,yl,yu);
Rxy=int(int((x-mx)*(y-my)*fxy,y,yl,yu),x,xl,xu);
Ry=int((y-my)^2*fy,y,yl,yu);

%stima LMEQM
xlmeqm=mx+Rxy*inv(Ry)*(y-my)

figure
ya=[yl:0.05:yu];
plot(ya,subs(xmeqm,ya),'r')
hold on
plot(ya,subs(xlmeqm,ya),'g')
xlabel('y')
ylabel('stime di X')
legend('MEQM','LMEQM')