PROGRAMACION
PROBLEMÁTICA 1
EQUIPO 4
REGLA DEL TRAPECIO
UNIVERSIDAD AUTONOMA DE COAHUILA
FACULTAD DE CIENCIAS QUIMICA
INTEGRANTES:
CLARISA MARICRUZ VILLALPANDO PEREZ
DIEGO ALEJANDRO MAGALLANES ACEVEDO
GUSTAVO ENOC FUENTES HERNANDEZ
JAZMIN ESMERALDA GUERRERO FLORES
JUAN PABLO TORRES MENDOZA
XIMENA GUADALUPE MORENO PECINA
PROBLEMÁTICA
Elabórese un programa que pida al usuario los limites de integración (a y b) y el numero
de subdivisiones (N) y que devuelva la aproximación mediante la suma de trapecios del
𝑏
𝑥𝛿ⅇ𝑛𝑥
valor de la integral de la función: න ⅆ𝑥
𝑥+ⅇ𝑥
𝑎
Los limites de integración a y serán introducidos por el usuario
Ejecútese el programa para a=0° y b=90°, comprobando que:
➢ Para 100 subdivisiones (N=10) el valor de la integral es: 0.264258
➢ Para 100 subdivisiones (N=100) el valor de la integral es: 0.264766
ANTECEDENTES
La integración numérica es una técnica que se puede usar para aproximar el valor de la
integral de una función que no sea posible anti diferenciar (integrar).Con el objeto de
integrar numéricamente la integral comprendida en el intervalo cerrado [a, b], lo
podemos hacer a través de dos métodos de integración numérica: la Regla del trapecio y
la Regla de Simpson.
➢ El objetivo de este trabajo es, por un lado, el de analizar y caracterizar algunas
formulas de cuadratura conocidas en la literatura, como método de aproximación de
integrales definidas en intervalos de la recta real, desde un punto de vista teórico, y
por otro lado, el de considerar ciertos aspectos computacionales para un cálculo
eficiente de estas reglas.
REGLA DEL TRAPECIO
Para aproximar una integral definida consiste en utilizar n trapecios, como se
muestra en la figura, en la formulación de este método, se supone que f es continua
y positiva en el intervalo [a , b] de tal modo, la integral definida
𝑏
න 𝑓 𝑥 ⅆ𝑥
𝑎
Representa el área de la región delimitada por la grafica de f y el eje x, desde x=a hasta x=b.
primero , se divide el intervalo [a, b] en n subintervalos, cada uno de ancho dx=(b-a)/n, de
modo que
Luego se forma el trapecio para cada subintervalo como el la figura, El área del i-
esimo del trapecio es :
HIPOTESIS
PROBLEMÁTICA 1
1.- 2.- 3.-
El ejercicio nos plantea una Para este ejercicio nos indica Para realizar el ejercicio
problemática que se puede que nos tiene que dar un tenemos que saber que el
resolver usando la regla del resultado de 0.264258, cuando límite inferior en ambos
trapecio, la cual es ´´ancho= tenemos 10 segmentos, y un casos es 0, y el límite
(b-a) /N´´ resultado de 0.264766 cuando superior también en ambos
tenemos 100 segmentos. casos es de 90*pi/180 o
pi/2, que representan lo
mismo.
4.- 5.-
. Esto se tiene que representar Al plasmar estos datos en el
así porque las funciones de código en Matlab,
Matlab funcionan en radianes, y conseguiremos sacar los
eso hubo que pasarlo de grado
resultados con los distintos
a radianes antes de poder
segmentos que se nos dan.
evaluarlo.
ALGORITMO
v
1. Inicio
2. Introducir la función sin evaluar la entrada como expresión y guardarlo como F
3. Construir el identificador de función a partir del vector de caracteres con la variable F, con el
siguiente comando:
f = str2func(['@(x) ‘, F])
3. Definir Límite Inferior y guardarlo como a, aceptando la entrada como solo y exclusivamente un
número
Definir Límite Superior y guardarlo como b, aceptando la entrada como solo y exclusivamente un
número
Definir Numero de segmentos y guardarlo como N, aceptando la entrada como solo y exclusivamente
un número
ALGORITMO
v
4. Si N<=0
disp “Valor de ingreso no valido”
Fin
5. Definir método de trapecio como:
• Delta:
h = (b-a) /N
• Definir función
Y = f(X)
Y(1,2:(N)) = Y(1,2:(N)) *2
• Formula de Método de Trapecios
suma = sum(Y)*h/2
Disp(suma)
INICIO
DIAGRAMA DE FLUJO
Función (F)
Str2func
(identificador de
función)
N<=0
a,b y N
Valor de
ingreso no
h=(b-a) /N invalido.
X=linspace(a,b, N+1)
Y=(1.2:(N))= Y(1.2:(N))*2
suma= sum(Y)*h/2
suma FIN
CODIGO
1. clc [Link] N<=0
2. clear all
3. close all
18. fprintf('Valor de ingreso no
4. F = input('f(x): ','s'); valido\n')
5. f = str2func(['@(x) ',F]); [Link]
% integrador
6. a = '0';
7. b = '0'; 20.h = (b-a)/N;
8. N = '0';
21.X = linspace(a,b,N+1);
9. while (isnumeric(a) == 0 && isnumeric(b) == 0 && isnumeric(N) == 0)
22.Y = f(X);
10. a = input('Limite inferior: ','s'); 23.Y(1,2:(N)) = Y(1,2:(N))*2;
11. b = input('Limite superior: ','s'); [Link] = sum(Y)*h/2;
12. N = input('Numero de segmentos: ','s');
[Link]('La integral vale %f\n',suma)
13. a = str2num(a);
14. b = str2num(b);
15. N = str2num(N);
16. end
PROBLEMÁTICA 2
TIRO AL BLANCO
Tenemos un cañón que lanza proyectiles con una velocidad inicial de disparo
𝑣0 y que nuestro objetivo es acertar un blanco que se encuentra localizado
en las coordenadas (𝑥𝐵 , 𝑦𝐵 ). Para acertar, tendremos que ajustar el ángulo
de tiro. Sabemos que los componentes de la velocidad inicial se calculan así:
dado que en la coordenada x el movimiento es uniforme ,en la coordenada y
es uniformemente acelerado, esperado y utilizando las ecuaciones
correspondientes a esos movimientos , Siguiente sistema de 2 ecuaciones
con 2 incógnitas ,t y a (𝑥𝐵 e 𝑦𝐵 son las coordenadas conocidas del blanco):
PROBLEMÁTICA 2
Resolviendo esta ecuación para la tan(a), encontramos los ángulos de disparo que
alcanzan el blanco.
Se pide hacer un script en Matlab que pida al usuario las coordenadas del blanco y la
velocidad de disparo del cañón y calcule :
➢ Los 2 ángulos que dan en el blanco
➢ El tiempo de vuelo del proyectil para cada ángulo hasta que impacta en el blanco .
➢ Una gráfica en la que se muestren en diferentes colores los posibles efectos del
proyectil desde que sale del cañón hasta que se impacta en el blanco. Esta gráfica
deberá estar formateada adecuadamente, con las leyendas y los títulos de los ejes
correspondientes y contener una malla o grid
Se supondrá que el cañón está situado en el origen de las coordenadas. Tómese el
valor de g= 9.8 m/s2
ANTECEDENTES
➢ Cualquier cuerpo puede realizar un movimiento parabólico: una pelota de futbol, de
tenis, un dardo, un misil. Cualquiera de ellos se denominan proyectil.
➢ En este tipo de movimiento el cuerpo se dispara con respecto a la horizontal y cae
bajo la acción de la gravedad. El vector de aceleración solo tiene componente Y por
lo cual la variable X se mantiene constante.
HIPOTESIS
PROBLEMÁTICA 2
1.- Punto de Partida:
Al analizar la problemática (tiro parabólico), se dedujo que, la ecuación
para obtener los ángulos era resultado de una ecuación cuadrática.
v
2.- Calculo de ángulos:
Utilizando la formula general, con un discriminante, obtendríamos las
condiciones necesarias para obtener los dos ángulos considerando a
tan(alfa) como la variable a obtener.
3.- Calculo del tiempo:
Sustituyendo las variables obtenidas de nuestros ángulos en la formula del
movimiento uniforme en “X” podríamos obtener nuestros derivados tiempos al
sustituir para “t”.
4.-Comparacion de graficas:
Al crear dos vectores que fueran desde 0 hasta el primer y segundo tiempo
obtendríamos lasa áreas a graficar en una comparativa entre el eje “x” y el eje ”y” y
mediante fórmulas, que multiplicarían a sus componentes con sus respectivos vectores,
obtendríamos las diferentes curvas del tiro hasta el blanco.
ALGORITMO
v
ALGORITMO
v
DIAGRAMA DE FLUJO
CODIGO