11 votos

Solución numérica de x = tan (x)

Necesitaba encontrar, usando el método de bisección, el primer valor positivo que satisfaga $x = \tan(x)$ . Así que fui a Scilab, escribí el método de bisección y obtuve $1.5707903$ . Pero después de razonar un poco he llegado a la conclusión de que este valor es erróneo:

  1. $\tan(1.5707903) \approx 1.6x10^5$ . Ni siquiera cerca de $1.5707903$ .
  2. Olvida por un momento lo anterior. $x = \tan(x)$ es en realidad encontrar puntos fijos de $f(x) = \tan(x)$ ; $(x, f(x))$ debe estar en la línea $y = x$ . Aquí es la trama: plot(tan(x), x)

En $(0, \frac{3}{2}\pi)$ Sólo puedo ver un punto fijo a la derecha de $x = 4$ Por lo tanto $1.5707903$ se equivoca.

Aquí viene la parte interesante. Si usted vaya a Wolfram Alpha y escriba $x = \tan(x)$ Verás que $1.5708$ en la sección de Parcela: x = tan(x)

Sin embargo, no hay $1.5708$ en la sección Soluciones numéricas. Wolfram Alpha encontró $0, \pm 4.49340945790906, \ldots$ .

Pero si escribe $\tan(x) = x$ no verás $1.5708$ en el apartado de Parcela!: tan(x) = x

Para resumir:

  1. Es $4.49340945790906$ el primer valor positivo que satisfaga $x = \tan(x)$ ?
  2. ¿Sabes por qué Wolfram Alpha muestra $1.5708$ como solución al escribir $x = \tan(x)$ pero no cuando se escribe $\tan(x) = x$ ?

Gracias.

0 votos

0 votos

(1) $\tan x$ mapas $(-\pi/2,\pi/2)$ de forma biyectiva sobre $\mathbb{R}$ . (2) Además, $\tan$ es periódica con periodo $\pi$ . De (1) y (2) se deduce que el $x$ que buscamos es en $(\pi/2,3\pi/2)$ .

0 votos

BTW la bisección no es muy rápida en absoluto, yo sugeriría que usted intentó con Newton-Raphson con un punto de partida como $\pi$ .

5voto

Halfgaar Puntos 2866

Como se ve en el gráfico de $\tan x$ , estás interceptando la asíntota, lo que no es realmente el comportamiento deseado. La bisección no es el mejor método a utilizar.

Sin embargo, si se requiere el uso de la bisección, entonces en lugar de ello, tenga en cuenta que $\tan x = \frac{\sin x}{\cos x}$ por lo que, para valores relevantes de $x$ ,

$$x = \tan x \implies x\cos x - \sin x = 0$$

Esta última función es continua, y debería obtener la solución deseada de $x \approx 4.49$ .

0 votos

Me gusta tu idea. Ahora me pregunto: ¿Qué función usarías como $g(x)$ si se requiere usar el método de iteración de punto fijo?

0 votos

@DavidRobertJones De entrada, sospecho que $\sin^{-1}(x\cos x) = x$ debería funcionar, siempre que se inicie la iteración en $(0,2\pi)$ .

0 votos

No estoy seguro ya que $\sin^{-1}(x\cos x)$ no tiene puntos fijos en el intervalo. Tampoco olvides las condiciones de convergencia.

4voto

jlupolt Puntos 369

La razón por la que obtienes esta "solución" es porque el método de bisección asume que la función es continua en el rango, que no es . Dado que la función a ambos lados de $x=\pi/2$ es $\pm \infty$ el método de bisección siempre convergerá a esta "solución".

0 votos

Tampoco es continuo en $[\pi/2,3\pi/2]$ Sin embargo, la solución está ahí.

0 votos

@DavidRobertJones - es continuo en el intervalo abierto.

-1voto

merletta Puntos 41

Tengo un código de MathLab, que se basa en el método de bisección; Aquí está:

clc borrar % Ceros de la ecuación x = tan(x); % Las raíces de la ecuación x=tanx mediante el método de bisección % z representa la raíz % epsilon es el error % NR es el número de raíces que se desea; número de raíces = 2*NR+1 % Max es el número de bisecciones

epsilon = 0,0001; NR = 2; Max = 100;

g=@(x) x-tan(x);

para k = -NR:NR a = (2*k-1)*pi/2 + épsilon; b = (2*k+1)*pi/2 - epsilon; si feval(g,a)*feval(g,b) < 0 para j=1:Max z = (a+b)/2; si (feval(g,a)*feval(g,z) < 0) b=z; elseif (feval(g,a)*feval(g,z) > 0) a=z; fin fin
si feval(g,z) < epsilon %error
z fin end end

% de salidas para NR=2 son : -7.7253 , -4.4934 , 0 , 4.4934 , 7.7253

-1voto

Darsh Puntos 1

1.57 es incorrecto porque tan π/2 no está definido 3.14/2 es 1.57 por lo tanto no está definido y no puede ser denominado como el menor valor positivo de tan X = x..por lo tanto el menor valor positivo es 4.49 para tanx = x

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X