Deber 1 Metodos Numericos

Escuela Superior Politécnica de Chimborazo
Ingeniería Automotriz
Deber No. 1
Métodos Numéricos
Edy Ulloa
1314
29 de octubre del 2015
1. Calcule el error absoluto y error relativo en las aproximaciones de p mediante 𝑝∗ :
p p* E absoluto E relativo
a) 3,141592654 3 1/7 ,126448926735E-02 ,4025E-03
b) 2,718281828 2,718 ,281828459045E-03 ,1037E-03
c) 3,141592654 3,1416 ,734641020683E-05 ,2338E-05
d) 1,414213562 1,414 ,213562373095E-03 ,1510E-03
e) 22026,46579 22000 ,264657948067E+02 ,1202E-02
2. Determine el mayor intervalo en que debe estar p* para aproximar p con un error relativo de a
lo sumo 10-4 para cada valor de p.
p E relativo p*
a) 3,14159265 ,10E-04 ,314156E+01
b) 2,71828183 ,10E-04 ,271825E+01
c) 1,41421356 ,10E-04 ,141420E+01
3. Realice los siguientes cálculos
a) En forma exacta
b) Mediante una aritmética de truncamiento a tres cifras
c) Mediante una aritmética de redondeo a tres cifras
d) Calcule los errores relativos en los incisos b y c
Forma exacta At(3) Ar(3) E rel(t) E relr

1) 1 2/15 ,113E+01 ,113E+01 ,2941176E-02 ,294118E-03
2) 4/15 ,266E+00 ,267E+00 ,250000E-02 ,125000E-02
3) 139/660 ,210E+00 ,211E+00 ,287770E-02 ,187050E-02
Datos
4/5 3/11 8,00E-01 2,727E-01
1/3 3/20 3,33E-01 1,50E-01
47 Use una aritmética de redondeo a tres cifras para los siguientes cálculos. Calcule el error absoluto
y el error relativo con el valor exacto determinado a por lo menos cinco cifras.
a) 133 + 0.921
b) 133 - 0,499
3
c) -10π+6℮-
62
22
𝜋− 7
d) 1
17
Datos truncados 3 Datos truncados 4 Datos redondeados 3 Datos redondeados 4

,133E+03 ,1330E+03 ,133E+03 ,1330E+03
A
,921E+00 ,9210E+00 ,921E+00 ,9210E+00
,133E+03 ,1330E+03 ,133E+03 ,1330E+03
B
,499E+00 ,4990E+00 ,499E+00 ,4990E+00
,314E+01 ,314E+02 ,3141E+01 ,3141E+02 ,314E+01 ,314E+02 ,3142E+01 ,3142E+02
C ,271E+01 ,163E+02 ,2718E+01 ,1630E+02 ,271E+01 ,163E+02 ,2718E+01 ,1631E+02
,483E-01 ,4838E-01 ,484E-01 ,4839E-01
,314E+01 ,3141E+01 ,314E+01 ,3142E+01
D ,314E+01 ,3142E+01 ,314E+01 ,3143E+01
,588E-01 ,5882E-01 ,588E-01 ,5882E-01
Datos reales
133 0,499 3,141592654 3/62 3 1/7
0,921 10 2,718281828 6 1/17
C Calculado V Real E Abs E rel

4 ,134E+03 ,13392E+03 ,790E-01 ,590E-03
5 ,1339E+03 ,13392E+03 ,2100E-01 ,1568E-03
A
6 ,133E+03 ,13392E+03 ,921E+00 ,687E-02
7 ,1339E+03 ,13392E+03 ,2100E-01 ,1568E-03
4 ,133E+03 ,13250E+03 ,499E+00 ,377E-02
5 ,1325E+03 ,13250E+03 ,1000E-02 ,7547E-05
B
6 ,132E+03 ,13250E+03 ,501E+00 ,378E-02
7 ,1325E+03 ,13250E+03 ,1000E-02 ,7547E-05
4 -,152E+02 -,15155E+02 ,454E-01 ,299E-02
5 -,1519E+02 -,15155E+02 ,3538E-01 ,2334E-02
C
6 -,151E+02 -,15155E+02 ,546E-01 ,360E-02
7 -,1518E+02 -,15155E+02 ,2538E-01 ,1674E-02
4 -,503E+02 -,21496E-01 ,503E+02 ,234E+04
5 -,5026E+02 -,21496E-01 ,5024E+02 ,2337E+04
D
6 -,502E+02 -,21496E-01 ,502E+02 ,233E+04
7 -,5026E+02 -,21496E-01 ,5024E+02 ,2337E+04
a b
,100E+01 ,100E+01
,100E+01 ,100E+01 Sumatoria e E abs E rel
,500E+00 ,500E+00 a ,27167E+01 2,71828183 ,16152E-02 ,59418E-03
,167E+00 ,167E+00 b ,27183E+01 2,71828183 ,27313E-07 ,10048E-07
,417E-01 ,417E-01
,833E-02 ,833E-02
,139E-02
,198E-03
,248E-04
,276E-05
,276E-06
a)
𝑥𝑐𝑜𝑠(𝑥) − 𝑥𝑠𝑒𝑛(𝑥)
lim
𝑥→0 𝑥 − 𝑠𝑒𝑛𝑥
L´Hospital
𝑐𝑜𝑠𝑥 − 𝑥𝑠𝑒𝑛𝑥 − (𝑠𝑒𝑛𝑥 + 𝑥𝑐𝑜𝑠𝑥)
lim
𝑥→0 1 − 𝑐𝑜𝑠𝑥
L´Hospital
−𝑠𝑒𝑛𝑥 − 𝑥𝑐𝑜𝑠𝑥 − 𝑠𝑒𝑛𝑥 − 𝑐𝑜𝑠𝑥 − 𝑐𝑜𝑠𝑥 + 𝑥𝑠𝑒𝑛𝑥
lim
𝑥→0 𝑠𝑒𝑛𝑥
−2𝑠𝑒𝑛𝑥 − 𝑥𝑐𝑜𝑠𝑥 − 2𝑐𝑜𝑠𝑥 + 𝑥𝑠𝑒𝑛𝑥
lim
𝑥→0 𝑠𝑒𝑛𝑥
L´Hospital
𝑥𝑠𝑒𝑛𝑥 − 𝑐𝑜𝑠𝑥 − 2𝑐𝑜𝑠𝑥 + 𝑥𝑐𝑜𝑠𝑥 + 𝑥𝑠𝑒𝑛𝑥 + 2𝑠𝑒𝑛𝑥
lim
𝑥→0 −𝑐𝑜𝑠𝑥
𝑥𝑠𝑒𝑛𝑥 − 3𝑐𝑜𝑠𝑥 + 𝑥𝑐𝑜𝑠𝑥 + 3𝑠𝑒𝑛𝑥
lim
𝑥→0 −𝑐𝑜𝑠𝑥
b)
f(0.1) =2,58896=0,258896E+01
c)
V calculado V real E relativo

,2588960E+01 ,1999000E+01 ,2951276E+00
At(3) Ar(3) At(4) Ar(4)

,142857E+00 ,142E+00 ,143E+00 ,1428E+00 ,1420E+00
,204082E-01 ,204E-01 ,204E-01 ,2040E-01 ,2040E-01
,291545E-02 ,291E-02 ,292E-02 ,2915E-02 ,2910E-02
,416493E-03 ,416E-03 ,416E-03 ,4164E-03 ,4160E-03
p* ,165E+00 ,167E+00 ,1665E+00 ,1657E+00
p 0,16659725 0,16659725 0,16659725 0,16659725
E abs ,159E-02 ,139E-03 ,658E-04 ,871E-03
E rel ,954E-02 ,834E-03 ,395E-03 ,523E-02
11. Use el formato real largo de 64 bits para determinar el decimal que equivale a los siguientes
números de máquina de punto flotante
a) 0 10000001010 1001001100000000000000000000000000000000000000000000
>> ((-1)^0)*(2^( bin2dec('10000001010')-1023))*(1+(1/2)^1+(1/2)^4+(1/2)^7+(1/2)^8)
ans =
3224
b) 1 10000001010 1001001100000000000000000000000000000000000000000000
>> ((-1)^1)*(2^( bin2dec('10000001010')-1023))*(1+(1/2)^1+(1/2)^4+(1/2)^7+(1/2)^8)
ans =
-3224
c) 0 01111111111 0101001100000000000000000000000000000000000000000000
>> ((-1)^0)*(2^( bin2dec('01111111111')-1023))*(1+(1/2)^2+(1/2)^4+(1/2)^7+(1/2)^8)
ans =
1.3242
Factor comun
L´Hospital
c E relativo
a 2
b 2,00940167 2,003335 0,00302828
13. Use el sistema de ayudas de MATLAB para encontrar información sobre las siguientes constantes
especiales: realmin, realmax, inf, NaN, pi, eps.
>> realmin
ans =
2.2251e-308
>> realmax
ans =
1.7977e+308
>> pi
ans =
3.1416
>> eps
ans =
2.2204e-16
>> inf
ans =
Inf
>> NaN
ans =
NaN
14. Obtenga información sobre comandos utilizados en Matlab, por medio de los sistemas de ayuda de
MATLAB. En particular, consulte sobre for, while, if, end, tres puntos (...), hold y line.
For
Declaraciones programáticas, extremo ejecuta repetidamente una o más sentencias de
MATLAB en un bucle. Valores tiene una de las siguientes formas:
initval:endval incrementos de la variable índice de valor

inicial al valor final por 1 , y repite la
ejecución de sentencias de programa
hasta que el índice es mayor que endval
initval:step:endval incrementa el índice por el paso de valor
en cada iteración , o decrementos cuando
paso es negativo.
valArray crea un índice de vector columna de las
columnas posteriores de matriz valarray
en cada iteración. Por ejemplo, en la
primera iteración , índice = valarray (: , 1
) . El bucle se ejecuta por un máximo de
n veces , donde n es el número de
columnas de valarray , dada por Numel (
valarray , 1 , :) . La entrada valarray
puede ser de cualquier tipo de datos
MATLAB , incluyendo una cadena,
matriz celular, o struct .
while
mientras que la expresión, declaraciones, extremo ejecuta repetidamente una o más

sentencias de programa MATLAB en un bucle mientras una expresión sigue siendo
cierto.
Una expresión evaluada es verdadera cuando el resultado es no vacío y contiene todos

los elementos no nulos (lógicos o numéricos real). De lo contrario, la
expresión es falsa.
if,
si la expresión , declaraciones , extremo evalúa una expresión y ejecuta un grupo de
sentencias cuando la expresión es verdadera .
elseif y else son opcionales , y ejecutar sentencias sólo cuando las expresiones
anteriores en el bloque si son falsas . Un bloque
if puede incluir múltiples sentencias elseif .
Una expresión evaluada es verdadera cuando el resultado es no vacío y contiene todos

los elementos no nulos ( lógicos o numéricos real) . De lo contrario , la expresión es
falsa
End
extremo termina por, mientras que, conmutadores, intente, si, y parfor declaraciones.
Sin una declaración final, porque, mientras que, switch, tratar, si, y parfor esperar a
nuevas aportaciones. Cada extremo está emparejado con el más cercano anterior no
pareado para, mientras, switch, tratar, si, o parfor y sirve para delimitar su alcance.
final también marca la terminación de una función, aunque en muchos casos es

opcional. Si su función contiene una o más funciones anidadas, entonces usted debe
terminar cada función en el archivo, si anidado o no, con el extremo. Esto incluye
funciones primarias, anidados, privadas y locales.
Hold
Controla la función de retención si MATLAB® borra el gráfico actual cuando usted
hace llamadas posteriores a tablas de funciones (por defecto), o agrega un nuevo gráfico
para el gráfico actual.
aferrarse conserva el gráfico actual y agrega otro gráfico a la misma. MATLAB ajusta
los límites de ejes, marcas de verificación, y garrapatas etiquetas como sea necesario
para mostrar la gama completa de la gráfica añadido.
contener a reajustes tienen estado al comportamiento predeterminado, en el que

MATLAB borra el gráfico actual y restablece ejes propiedades a sus valores
predeterminados antes de sacar nuevas parcelas.
mantenga todas las bodegas del gráfico y el color de la línea actual y estilo de línea para
que los comandos de trazado subsiguientes no restablecen los valores de propiedad
ColorOrder y LineStyleOrder al principio de la lista. Trazado de comandos siguen
ciclismo a través de los colores predefinidos y estilos de línea desde donde el último
gráfico se detuvo en la lista.
bodega invierte el estado de espera actual. Si el estado de retención se encuentra

actualmente en, a continuación, un comando de detención establece el estado de
apagado. Del mismo modo, si el estado de espera es actualmente, luego un comando de
detención establece el estado activado.
Line
línea crea un objeto de línea en los ejes actuales con los valores por defecto x = [0 1] ey
= [0 1]. Se puede especificar el color, ancho, estilo de línea y tipo de marcador, así
como otras características.
15. Desarrolle, depure y pruebe un programa en MATLAB, a .n de calcular la velocidad del paracaídas
que cae como se explicó en el ejemplo desarrollado en clase. Diseñe el programa de modo que permita
al usuario introducir valores para el coeficiente de arrastre y la masa. Pruebe el programa con la
reproducción de los resultados del ejemplo de clase. Repita el cálculo pero utilice tamaños de paso de
0.2 y 0.1 s (Grafique sus resultados en Matlab). Compare sus resultados con la solución analítica que se
obtuvo previamente. Un tamaño de paso más pequeño, ¿hace que los resultados sean mejores o
peores?
Explique sus resultados.
Datos
Masa=70kg
Constate=12.5
Gravedad=9.8m/s2
Paso= 0.1 y 0.2
% ESCUELA DE INGENIERIA AUTOMOTRIZ

% Materia:Metodos Numericos
% Semestre: Quinto A
% Autor: Pedro Puyol
v=0
t=0
clc;
disp('Calculo de velocidad');
m=input('ingrese el valor de la masa m = ');
c=input('ingrese el valor de la constante c = ');
g=input('ingrese el valor de la gravedad g = ');
p=input('ingrese el valor del paso del tiempo p = ');
v(1)=0;
A=10/p;
for i= 1:A
v(i+1) =v(i)+(g-((c/m)*v(i)))*p;
end
t(1)=0;
for x= 1:A
t(x+1)=x*p;
end
t
v
plot(v,t)

Deber 1 Metodos Numericos

Cargado por

Copyright:

Formatos disponibles

Deber 1 Metodos Numericos

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Deber 1 Metodos Numericos

Cargado por

Copyright:

Formatos disponibles

Escuela Superior Politécnica de Chimborazo

Forma exacta At(3) Ar(3) E rel(t) E relr

Datos truncados 3 Datos truncados 4 Datos redondeados 3 Datos redondeados 4

C Calculado V Real E Abs E rel

V calculado V real E relativo

At(3) Ar(3) At(4) Ar(4)

>> ((-1)^0)(2^( bin2dec('10000001010')-1023))(1+(1/2)^1+(1/2)^4+(1/2)^7+(1/2)^8)

>> ((-1)^1)(2^( bin2dec('10000001010')-1023))(1+(1/2)^1+(1/2)^4+(1/2)^7+(1/2)^8)

>> ((-1)^0)(2^( bin2dec('01111111111')-1023))(1+(1/2)^2+(1/2)^4+(1/2)^7+(1/2)^8)

initval:endval incrementos de la variable índice de valor

mientras que la expresión, declaraciones, extremo ejecuta repetidamente una o más

Una expresión evaluada es verdadera cuando el resultado es no vacío y contiene todos

Una expresión evaluada es verdadera cuando el resultado es no vacío y contiene todos

final también marca la terminación de una función, aunque en muchos casos es

contener a reajustes tienen estado al comportamiento predeterminado, en el que

bodega invierte el estado de espera actual. Si el estado de retención se encuentra

% ESCUELA DE INGENIERIA AUTOMOTRIZ

También podría gustarte