UBACS Q&A Foro WikiCS
Fecha actual 18 Dic 2018, 22:49




 Página 1 de 1 [ 3 mensajes ] 
Autor Mensaje
 Asunto: [Supuestamente resuelto] Ej. de Matlab nº3
NotaPublicado: 21 Jun 2009, 04:53 
Doctor
Avatar de Usuario

Registrado: 09 Ago 2008, 19:30
Mensajes: 380
Dice. Supongamos datos ({X}_{i},{Y}_{i}), generados según la siguiente ley:

{X}_{i}=2rand-1
{Y}_{i}=3cos({X}_{i})+5{{X}_{i}}^{2}+2{X}_{i}+5(rand-1)

(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:

Resolución
E=zeros(1,100)
for j=5:100;
X=zeros(j,1);
Y=zeros(j,1);
for i=1:j;
X(i)=2*rand-1;
Y(i)=3*cos(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=a*cos(X)+b*X.^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.



_________________
Measure your life in love.
("Vete, pues. Existen otros mundos aparte de estos.")
Desconectado
 Perfil  
 
 Asunto: Re: [Supuestamente resuelto] Ej. de Matlab nº3
NotaPublicado: 01 Jul 2009, 18:49 
Profesor

Registrado: 08 Ago 2008, 21:57
Mensajes: 299
Esto es de un parcial??

A la noche, Lucas lo trajo como ejercicio


Desconectado
 Perfil  
 
 Asunto: Re: [Supuestamente resuelto] Ej. de Matlab nº3
NotaPublicado: 01 Jul 2009, 19:36 
Doctor
Avatar de Usuario

Registrado: 09 Ago 2008, 19:30
Mensajes: 380
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



_________________
Measure your life in love.
("Vete, pues. Existen otros mundos aparte de estos.")
Desconectado
 Perfil  
 
Mostrar mensajes previos:  Ordenar por  
 Página 1 de 1 [ 3 mensajes ] 


¿Quién está conectado?

Usuarios navegando por este Foro: No hay usuarios registrados visitando el Foro y 1 invitado


No puede abrir nuevos temas en este Foro
No puede responder a temas en este Foro
No puede editar sus mensajes en este Foro
No puede borrar sus mensajes en este Foro
No puede enviar adjuntos en este Foro

Buscar:
Saltar a:  

cron