Seguimos buscando a Arshak. Ayudanos compartiendo!
Encuesta no oficial de docentes
Resultados de la encuesta no oficial de docentes
Probaste el SIGA Helper?

Donar $100 Donar $200 Donar $500 Donar mensualmente


Enviar respuesta 
 
Calificación:
  • 0 votos - 0 Media
  • 1
  • 2
  • 3
  • 4
  • 5
Buscar en el tema
[Pedido][Analisis Numerico] Tp de Optimizacion de castro 2011
Autor Mensaje
charrly Sin conexión
Empleado del buffet
esperando el milagro
*

Ing. Industrial
Facultad Regional Buenos Aires

Mensajes: 1
Agradecimientos dados: 0
Agradecimientos: 0 en 0 posts
Registro en: Oct 2011
Mensaje: #1
[Pedido][Analisis Numerico] Tp de Optimizacion de castro 2011 Trabajo practico Análisis Numérico y Cálculo Avanzado
Se solicita con suma urgencia dadores del metodo del gradiente en matlab de la materia analisis numerico o bien el tp de optimizacion de castro de este año.

Si alguno lo tiene que lo suba porque estoy al horno!!!
(Este mensaje fue modificado por última vez en: 24-11-2011 14:09 por gonnza.)
24-11-2011 10:11
Encuentra todos sus mensajes Agregar agradecimiento Cita este mensaje en tu respuesta
UES Sin conexión
Empleado del buffet
Sin estado :(
*

Ing. Industrial
Facultad Regional Buenos Aires

Mensajes: 20
Agradecimientos dados: 0
Agradecimientos: 0 en 0 posts
Registro en: Mar 2010
Mensaje: #2
RE: [Pedido][Analisis Numerico] Tp de Optimizacion de castro 2011
Aca te mado la codificacion que use el año pasado para hacer el TP, espero que te sirva


Funciones utilizadas:
Parametrización:
function param=parametrizacion (t, x, y,derparcx, derparcy)
param=2*(1-(x+derparcx*t))^2 + 200* ((y+derparcy*t)-(x+derparcx*t)^2)^2;
end

Derivadas parciales:
function derx = vx( x,y )
derx=4*(1-x)+800*x*(y-x^2);
end

function dery = vy( x,y )
dery=-(400*(y-x^2));
end

Función a Minimizar:
function z = funcion( x,y )
z=2*(1-x)^2 + 200*(y-x^2)^2;
end







Codificacion
format long;
clear;
clc;
it=0;
x=-1;
y=1;
it=0; b=6; c=0; t=0; d=0;
r2= (sqrt (5.0) - 1.0)/2.0;
r1= 1.0 - r2;
derparcx=vx(x,y);
derparcy=vy(x,y);
p=parametrizacion(t,x,y,derparcx,derparcy);
display(derparcx);
display(derparcy);
error2d=10^6;
N=sqrt(derparcx^2 + derparcy^2);
tol1d=10^(-7);
tol2d=10^(-6);

while (error2d > tol2d)
t=0;
b=6;
derparcx=derparcx/N;
derparcy=derparcy/N;
derparcx=vx(x,y);
derparcy=vy(x,y);
it=it + 1;
error1d=10^6;
c=t + r1 * (b-t);
d=t + r2 * (b-t);
fc= parametrizacion (c, x, y, derparcx, derparcy);
fd= parametrizacion (d, x, y, derparcx, derparcy);
while (error1d > tol1d)
c=t + r1 * (b-t);
d=t + r2 * (b-t);
fc= parametrizacion (c, x, y, derparcx, derparcy);
fd= parametrizacion (d, x, y, derparcx, derparcy);
if (fc <= fd)
b=d;
d=c;
c=t + r1 * (b-t);
fd=fc;
fc=parametrizacion (c, x, y, derparcx, derparcy);
end;
if (fc>fd)
t=c;
c=d;
d=t + r2 * (b-t);
fc=fd;
fd=parametrizacion (d, x, y, derparcx, derparcy);
end;
tmin= (c+d) /2;
error1d= abs (d-c);
error1d = sqrt(((d-c)*derparcx)^2+((d-c)*derparcy)^2);
end;
dx=tmin * derparcx;
dy=tmin * derparcy;

error2d=sqrt (dx^2 + dy^2);
x=x+dx;
y=y+dy;
end;
display(it);
24-11-2011 20:13
Encuentra todos sus mensajes Agregar agradecimiento Cita este mensaje en tu respuesta
romanpiacquadio Sin conexión
Empleado del buffet
Sin estado :(
*

Ing. Industrial
Facultad Regional Buenos Aires

Mensajes: 5
Agradecimientos dados: 4
Agradecimientos: 0 en 0 posts
Registro en: Jul 2010
Mensaje: #3
RE: [Pedido][Analisis Numerico] Tp de Optimizacion de castro 2011
Capo me darias una mano con este quilombito que tengo? estoy cursando con castro ya aprobe los parciales pero me faltan los tp 2, 3 y 4 para poder firmala. Con cualquiera de ellos que me puedas ayudar me viene barbaro. Te lo recontra agredecería!!
06-12-2011 15:35
Envíale un email Encuentra todos sus mensajes Agregar agradecimiento Cita este mensaje en tu respuesta
Buscar en el tema
Enviar respuesta 




Usuario(s) navegando en este tema: 1 invitado(s)