9/14/24, 1:01 PM MATLAB Grader
Relatividad My Solutions
Supongamos que estamos evaluando la dilatación del tiempo Δt de una partícula que se mueve a una velocida
embargo, en este caso, la partícula no solo se mueve a alta velocidad, sino que también experimenta una acele
externas. La dilatación del tiempo y la velocidad de la partícula se describe mediante la siguiente ecuación no li
(c3)
− α sin ( )
Δt0 v3 v
Δt(v) = +k
c
√
v2
1−
c2
Para este modelo:
Δt0 = 1s representa el tiempo propio (tiempo medido en el marco de referencia de la partícula)
c = 3 × 108m/s es la velocidad de la luz en el medio.
k = 0.02 es un factor de corrección relacionado con efectos relativistas no lineales.
α = 0.5s es un término que representa la influencia de fuerzas externas en la dilatación del tiempo.
v es la velocidad de la partícula.
El objetivo es encontrar la velocidad v a la cual la dilatación del tiempo se duplica en comparación con el tiemp
valor inicial de v = 2.5 × 108 m/s para comenzar con las iteraciones del método de Newton. Para ello:
Defina una función f (v) cuya raíz resuelva dicho problema. Dar como respuesta ∣ f (2.8 × 108)∣ y alma
Aproxime la velocidad v con una tolerancia de 1e − 3 utilizando el método de Newton - Raphson, alma
Almacene la cantidad de iteraciones empleadas en la variable totaliter.
Script
Save Reset MATLAB Documentation ([Link]
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
[Link] 1/3
9/14/24, 1:01 PM MATLAB Grader
17
18 % Parámetros del problema
19 format long
20 c = 3e8; % Velocidad de la luz en m/s
21 dT0 = 1; % Tiempo propio en segundos
22 k = 0.02; % Factor de corrección relacionado con efectos relativistas no lineales
23 alpha = 0.5; % Influencia de fuerzas externas en segundos
24 v0 = 2.5e8; % Velocidad inicial estimada en m/s
25
26 % Definir la función de la dilatación del tiempo
27
28 syms v
29 f=dT0/sqrt(1-(v^2/c^2))+ k*((v^3)/(c^3))-alpha*sin(v/c)-2;
30 f1=matlabFunction(f);
31
32 % Evaluar la función en un valor cercano para verificar
33
34 Fp=2.8*10^8;
35 nvalue=f1(Fp)
36
37 % Mostrar los resultados
38
39 Tol=1e-3;
40 x0=2.5*10^8;
41 z=newton(f1,x0,Tol)
42 vapprox=z(end)
43 totaliter=length(z)-1
Run Script
Assessment: All Tests Passed (100%)
Run Pretest Submit (max. reached)
La variable nvalue almacena el valor de la funcion en un punto de
34% (34%)
prueba?
vapprox almacena la velocidad requerida? (Pretest) 33% (33%)
Encuentra la cantidad de iteraciones empleadas? 33% (33%)
Total: 100%
[Link] 2/3
9/14/24, 1:01 PM MATLAB Grader
Output
nvalue =
0.399886681549311
z =
1.0e+08 *
2.500000000000000
2.855300380725175
2.775591141457275
2.730338323146660
2.722285862264590
2.722093846144570
vapprox =
2.722093846144570e+08
totaliter =
© 2024 The MathWorks, Inc.
[Link] 3/3