Universidad Nacional de Ingeniería Curso MB536-C
DACBHCC - FIM-UNI PA-2021-1
Test Laboratorio 2: Programación en MATLAB y teoría de errores
P1
Desarrolle una función llamada “nt”, que retorne el número de términos necesarios para
aproximar el número π hasta n cifras decimales exactos, usando la siguiente serie:
%n=cifras significativas
function s=nt(n)
y=1;
error=1;
t=4;
s=t;
while ((error)>(5*10^-n))
ts=4*((-1)^y)/(2*y+1);
s=s+ts;
y=y+1;
error=abs(s-pi)/pi*100;
end
s=vpa(pi,n);
end
P2
El cuadrado inclinado
Se determinará el área A de un cuadrado.
El rango h es 50.0 cm±1.0 cm, el rango s es 170.0 cm ±2.0 cm de largo.
a) ¿Qué tan grande es el área A?
A=(51^2+172^2)=32185cm2
b) ¿Qué error absoluto tiene A según la propagación del error?
Error abs.=(Amax-Amin)/2
Error abs.=(32185-30625)/2=780cm
P3
¿Cuál es valor del número en notación interna de la máquina?
0 1 0 0 0 0 0 1 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 =+8*1.75=14
2^(130-127)=8
-1^0=+1
Universidad Nacional de Ingeniería Curso MB536-C
DACBHCC - FIM-UNI PA-2021-1
P4
Como se escribe el número decimal -15.6875 en notación de punto flotante interna
para una máquina de simple precisión.
-15.6875= 1 | 1 0 0 0 0 0 1 0 | 1 1 1 1 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
P5
¿Cuál es el epsilón de la máquina hipotética F(2,9,-2,3)?
Eps=1+eps-1
1+eps= 0 | 0 1 1 | 0 0 0 0 0 0 0 0 1
1= 0|011|000000000
Eps= 0 | 0 0 0 | 0 0 0 0 0 0 0 0 1
P6 ¿Cuál es la salida del siguiente código?
Meansqr(NaN)
a) 0
b) NaN
c) Error
d) 1
P7
1 r2
Si V =
2 h ( )
+h y el error en V no es mayor al 0.4%, encontrar el porcentaje de error
permisible en r y h donde r ≈ 5.1 cm y h ≈ 5.8 cm
Rpsta:desarrollando por principio de igual efecto.
%εr=0.2293%
%εh=1.5635%
P8
Sea el polinomio de segundo grado x 2+ 2 px−q , con √ p 2+ q ≥ 0, en los problemas mal
condicionados cuando p ≫≫ q, se produce error de cancelación.
Evalúe el error relativo porcentual de los algoritmos a continuación a probar:
q
1) x=− p+ √ p2 +q 2) x=
p + √ p 2+ q
Evalué y reporte cuál de los casos usaría 1 ó 2, si p=4.999999999995*10 4 y q=0.01
Se evalúa el manejo de las cifras significativas exactas y de la teoría de errores.
Utilice Matlab para sus cálculos.
Caso 1 o 2 : 2
Raíz = 9.999999999999998e-08
Error relativo = 3.8072e-04
p=4.999999999995*10^4;
q=0.01;
format long
Universidad Nacional de Ingeniería Curso MB536-C
DACBHCC - FIM-UNI PA-2021-1
x1=-p+(p^2+q)^0.5
x2=q/(p+(p^2+q)^0.5)
era=(abs(x1-x2))/2
num=(x1+x2)/2
err=(era/num)*100
P9
a) δA ≤ 2∗δh +δR
b) δA ≤ δh+ δR
c) δA ≤ δh+ 2∗δR
d) δA ≤ 4∗δh +δR
P10
Grafique la función si usar lazos de control
Universidad Nacional de Ingeniería Curso MB536-C
DACBHCC - FIM-UNI PA-2021-1
Solamente vectores y direcciones o posiciones.
x=-2:0.2:2
x1=[zeros(1,11) (0.2:0.2:2)];
x2=[zeros(1,11) (0.2:0.2:0.8) zeros(1,6)];
x3=[(-2:0.2:0.8) zeros(1,6)];
y1=-1*x.^2-(-1*x1.^2);
y2=x2;
y3=x.^2+1-(x3.^2+1);
y=y1+y2+y3;
plot(x,y)
grid on
La Profesora.
R.M.G.J.