[Supuestamente resuelto] Ej. de Matlab nº3

Dice. Supongamos datos \(({X}_{i},{Y}_{i})\), generados según la siguiente ley:

(ie, los Xi están distribuidos de manera aleatoria en el intervalo [-1,1]).

Para cada valor j entre 5 y 100:

a) Generar una muestra de \(({X}_{i}.{Y}_{i})\)
b) Para la muestra dada, obtener los coeficientes A, B y C de cuadrados mínimos suponiendo una regresión \({Y}_{i}=Acos({X}_{i})+B{{X}_i}}^{2}+C\)
c)Guardar el SSR obtenido

Finalmente, graficar todos los SSR obtenidos, en función de j.

Nota: \(SSR=\sum_{i=1}^{n}{({Y}_{i}-F({X}_{i}))}^{2}\)

Mi programa:

E=zeros(1,100)
for j=5:100;
X=zeros(j,1);
Y=zeros(j,1);
for i=1:j;
X(i)=2rand-1;
Y(i)=3cos(X(i))+5*(X(i))^2+2*(X(i))+5*(rand-1);
end
A=ones(j,3);
A(:,1)=cos(X);
A(:,2)=X.^2;
B=A'A;
K=A'Y;
r=B\K;
a=r(1);
b=r(2);
c=r(3);
F=acos(X)+bX.^2+c;
e=0;
for i=1:j;
e=e+(Y(i)-F(i))^2;
end
E(j)=e;
end
E
N=1:100
plot(N,E,'.')

Al menos corre.

Qué sueño.

Esto es de un parcial??

A la noche, Lucas lo trajo como ejercicio

claro, yo lo hice para practicar para el parcial porque lo trajo lucas.. no creo que lo hayan tomado en un parcial, dada la dificultad despreciable del parcial quetomaron frente a este ej