Advertisement
nAIMx

Octave deberes

Nov 2nd, 2020 (edited)
2,318
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Octave 1.78 KB | None | 0 0
  1. BOLETIN 2 DE OCTAVE EJERCICIO 3(SANGRE)
  2.  
  3. %Definimos variables usando libreria symbolic
  4.  
  5. syms k
  6. syms t
  7.  
  8. % Definimos la ecuacion
  9.  
  10. f = t * (e ^ ((3 - t) / 3)) - 0.25
  11.  
  12. % Resolvemos para asi encontrar y escoger un intervalo cercano con el fin de que
  13. % el número de iteraciones que realicemos sea la adecuada
  14. double(solve(f, t))
  15.  
  16. a1 = -0.4
  17. b1 = 0.5
  18.  
  19. % Calcular el mínimo de iteraciones con el intervalo anterior(para que el error<2min)
  20.  
  21. solve((((b1-a1)/(2^k)) < (2/60)), k)
  22. log(27) / log(2)
  23.  
  24. % El resultado es 4.7549, se deben hacer 5 iteraciones
  25.  
  26. % PRIMERA ITERACION
  27.  
  28. x1 = (a1 + b1) / 2
  29.  
  30. %x1 = 0.05000
  31. double(subs (f, t, x1))
  32. % ans =  -0.11633 es diferente de 0, lo que quiere decir que x1 no es solución
  33. double(subs (f, t, x1) * subs (f, t, a1))
  34. % ans =  0.17361 > 0, el intervalo de la iteración  es (x1, b1)
  35.  
  36. % SEGUNDA ITERACION
  37.  
  38. x2 = (x1 + b1) / 2
  39.  
  40. % x2= 0.27500
  41. double(subs (f, t, x2))
  42. % ans =  0.43205, no es 0 osea que x2 no es solución
  43. double(subs (f, t, x2) * subs (f, t, x1))
  44. % ans = -0.050262 < 0, siguiendo la misma metodologia, el intervalo de la iteración es (x1, x2)
  45.  
  46. % TERCERA ITERACION
  47.  
  48. x3 = (x1 + x2) / 2
  49. % x3= 0.16250
  50. double(subs (f, t, x3))
  51. % ans =  0.16843 es diferente de 0, lo que quiere decir que x3 no es solución
  52. double(subs (f, t, x3)*subs (f, t, x1))
  53. % ans = -0.019594 < 0 es decir, el intervalo de la iteración es (x1, x3)
  54.  
  55. % CUARTA ITERACION
  56.  
  57. x4 = (x1 + x3) / 2
  58. % x4=  0.10625
  59. double(subs (f, t, x4))
  60. % ans =  0.028768 sigue siendo diferente de 0 por lo que x4 tampoco es solución
  61. double(subs (f, t, x4) * subs (f, t, x1))
  62. % ans = -0.0033466 < 0, el intervalo de la iteración es (x1, x4)
  63.  
  64. % QUINTA ITERACION
  65.  
  66. Se realiza la ultima iteracion para comprobar el ejercicio
  67.  
  68. x5 = (x1 + x4) / 2
  69. % x5=  0.078125
  70.  
  71.  
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement