Diferencia entre revisiones de «PrInf17: Algunos ejemplos con funciones»

Revisión del 13:13 9 ago 2013

Este artículo está en versión beta. El autor de este artículo no lo ha terminado todavía, por favor no lo edites hasta que elimine este mensaje.

En esta práctica vamos a escribir dos funciones y usarlas dentro de un programa, para ilustrar la metodología de diseño top-down. Tenemos que realizar un programa que, dado un polígono, pregunte un punto al usuario, y encuentre el vértice del polígono más próximo al punto introducido por el usuario. El programa debe dibujar el polígono, el punto introducido por el usuario y el vértice encontrado. Como es un programa algo complejo, vamos a escribir dos funciones, que usaremos desde el programa principal.

1 Requisitos previos

Antes de realizar esta práctica, es imprescindible haber realizado antes la siguiente práctica:

• PrInf16: Introducción a las funciones

Además, una de las funciones que vamos a programar tiene varios argumentos de salida, por lo que puede ser útil visualizar este vídeo (aunque en nuestro caso, los argumentos de salida no son opcionales):

• Prog24: Funciones con varios argumentos de salida

2 Contenido de la práctica

El programa necesita encontrar el vértice más cercano, dado un punto que introduce el usuario. El vértice más cercano será el que tiene la menor distancia al punto en cuestión. No existe ninguna función en Octave UPM para calcular la distancia entre dos puntos, así que tendremos que escribirla nosotros. Sabemos que la distancia entre dos puntos [math](x_0,y_0)[/math] y [math](x_1, y_1)[/math] viene dada por [math]\displaystyle d = \sqrt{(x_1-x_0)^2+(y_1-y_0)^2}[/math]

Tarea:

Crea una función de nombre distancia, que devuelva el valor de d dados cuatro argumentos de entrada: x0, y0, x1 e y1.

Una vez que tenemos una función que calcula la distancia entre dos puntos en el plano, vamos a empezar a programar la función que encuentra el vértice más cercano de un polígono. Esta función acepta como argumentos de entrada las coordenadas [math](x_p,y_p)[/math] del punto, y dos vectores x e y, que contienen los vértices del polígono. Estos vectores contienen los vértices en el mismo formato que usan las funciones fill e inpolygon. La función devolverá las coordenadas [math](x_v, y_v)[/math] del vértice más cercano al punto.

Por tanto, la cabecera de la función será:

function [xv, yv] = encuentraVerticeMasCercano(xp, yp, x, y)

Tendremos que guardar esta función en un fichero de nombre encuentraVerticeMasCercano.m.

El algoritmo para buscar el vértice más cercano consistirá en ir comprobando la distancia del punto a cada uno de los vértices, y buscar el valor más pequeño de la distancia. Vamos a ir viéndolo paso a paso. Lo primero que necesitamos conocer es el número de vértices del polígono, ya que nuestra función debería funcionar con cualquier tipo de polígono. Como los vectores x e y contienen las coordenadas de los vértices, el número de vértices será igual al número de elementos en cualquiera de esos dos vectores:

  nVertices = length(x);

  xv = x(1);
  yv = y(1);
  d = distancia(xp, yp, xv, yv);

Hemos usado la función distancia, que hemos programado anteriormente.

Ahora será necesario recorrer todos los vértices restantes, para ver si hay algún vértice más cercano al punto [math](x_p, y_p)[/math]. Para realizar esta tarea, usaremos un bucle for. Ignoraremos el primer vértice, puesto que ya hemos calculado su distancia:

  for k=2:nVertices
    x1 = x(k);
    y1 = y(k);
    dNueva = distancia(xp, yp, x1, y1);
    if dNueva < d
      d = dNueva;
      xv = x1;
      yv = y1;
    end
  end

Con esto, hemos terminado la función. El código completo quedaría así:

function [xv, yv] = encuentraVerticeMasCercano(xp, yp, x, y)
% Función que encuentra el vértice más cercano a un punto dado
%
%   xp: Coordenada x del punto
%   yp: Coordenada y del punto
%   x:  Vector con las coordenadas x de los vértices del polígono
%   y:  Vector con las coordenadas y de los vértices del polígono
%
%   xv: Coordenada x del vértice más cercano
%   yv: Coordenada y del vértice más cercano

  nVertices = length(x);
  xv = x(1);
  yv = y(1);
  d = distancia(xp, yp, xv, yv);

  for k=2:nVertices
    x1 = x(k);
    y1 = y(k);
    dNueva = distancia(xp, yp, x1, y1);

    if dNueva < d
      d = dNueva;
      xv = x1;
      yv = y1;
    end
  end
end

3 Ejercicio post-práctica