Parámetro Característico de A Catenaria
Parámetro Característico de A Catenaria
Parámetro Característico de A Catenaria
A B
80
70
f
60
50
40
30
α
20
10
0
- 50 - 45 - 40 - 35 - 30 - 25 - 20 - 15 - 10 -5 0 5 10 15 20 25 30 35 40 45 50
Cuando los puntos de los que pende están situados a la misma altura y se escogen los
ejes de tal manera que, con x = 0 , y ′ = 0 , entonces la ecuación de la catenaria es:
x
y = α cosh (1)
α
donde a es el parámetro característico que se define como el cociente entre la tensión horizon-
tal en cualquier punto de la catenaria T0 y el peso por unidad de longitud q:
T0
α= (2)
q
CÁLCULO DEL PARÁMETRO CARACTERÍSTICO DE LA CATENARIA 5
2. Resolución
De la integral de un diferencial de longitud de la catenaria entre 0 y una abscisa x, se obtiene
la siguiente relación:
x
∫0
ds → α 2 + s 2 = y 2 (3)
donde α es el parámetro característico, s la longitud e y la ordenada del punto x extremo de la
integral. Por tanto, si se tiene en cuenta que la ordenada de los puntos A y B de la Figura 2-1
es α + f , la longitud total del cable entre esos puntos será:
l = 2 (α + f ) − α 2
2
(4)
donde L valdrá 100 m. Además, conociendo la altura de la que cuelgan los cables y teniendo en
cuenta que la normativa española exige que la distancia mínima al suelo sea:
U
3.3 + m (6)
100
donde U es la tensión de la línea en KV, el valor de f resulta:
380
f = 45 − 3.3 − = 37.9 m (7)
100
En definitiva la ecuación que hay que resolver es:
100
α cosh − 1 − 37.9 = 0 (8)
α
Dado el programa que resuelve la ecuación por el método de la bisección, se pide pro-
gramar el método de la falsa posición para la resolución de ecuaciones no lineales y aplicarlo a
la ecuación (8) con el fin de determinar el valor de α y posteriormente el de la longitud de la
catenaria.
3. El método de la bisección
Para saber cuánto vale la función catenaria con los sucesivos valores aproximados de α que se
vayan obteniendo, se ha programado la función catenar que se encuentra en el directorio
publico\numerico\catenar.m:
function cate = catenaria(x)
cate=x.*(cosh(100./x)-1)-37.9;
6 PRÁCTICAS MÉTODOS NUMÉRICOS / MÉTODOS SIMBÓLICOS Y NUMÉRICOS
Así pues, el primer argumento de entrada a la función bisect será ‘catenar’. La fun-
ción bisect se encuentra también en el directorio publico\numerico\bisect.m y se escribe
a continuación.
function [c,ite,elapsed,oper,steps] = bisect(fun, x, tol)
%
% [xsol,ite,elapsed,oper,steps] = bisect(fun, x, tol)
%
%BISECT Se busca el cero de una función de una variable por el
% método de la bisección.
% fun es el nombre de la función (debe estar en un m-file). El
% nombre de la función tiene que ir escrito entre comillas.
% x es el punto central del intervalo en el que se va a buscar
% la raíz. Debe ser cercano a un punto donde la función cambie
% signo.
% tol es un argumento opcional que dice cuánto vale el error
% máximo admisible (la tolerancia para considerar una solución
% válida).
%
% salida
% xsol: valor de la raíz,
% ite: es el número de veces que ha habido que dividir el intervalo
% para encontrar la raíz.
% elapsed: tiempo invertido en encontrar la raíz.
% oper: número de operaciones necesarias para encontrar la raiz.
% steps: que da la historia del proceso (las sucesivas raíces y los
% valores que toma la función en ellas).
oper=flops; t0=clock;
if x ~= 0, dx = x/20; else, dx = 1/20; end
a = x - dx; fa = feval(fun,a);
b = x + dx; fb = feval(fun,b);
% Búsqueda del cambio de signo (no se hace el producto fa*fb porque es más
% "cara" una multiplicación que una comparación).
% Bucle principal. Obsérvese que, para evitar goto’s (se ha sustituido por un while),
% El criterio de convergencia se pone al principio y por tanto cambia el orden de los
% pasos que se han visto en clase.
ite=0;
while abs((b - a)/(b + a)) > tol
c = (a + b)/2; fc = feval(fun,c);
if trace, steps = [steps; [c fc]]; end
if (fb > 0) == (fc > 0)
CÁLCULO DEL PARÁMETRO CARACTERÍSTICO DE LA CATENARIA 7
b = c; fb = fc;
else
a = c; fa = fc;
end
ite=ite+1;
end
oper=flops-oper;
elapsed=etime(clock,t0);
a f (b ) − b f ( a ) f ( b )( b − a )
xm = =b−
f (b ) − f ( a ) f (b ) − f ( a )
xmi +1 − xmi
Era = <ε
xmi +1
5. Newton-Raphson y bisección
Dado el programa que resuelve la ecuación por el método de la bisección, se pide programar el
método de la iteración de punto fijo, el de Newton-Raphson y el de la Secante.
6. Aceleración de la convergencia
Probar los métodos de la ∆2 de Aitken y de Steffensen para acelerar la convergencia.