Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 11 vistas

    Integración

    Cargado por

    Juan Solis
    Tarea de métodos numéricos, en el método de integración.

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd

    Integración

    Cargado por

    Juan Solis
    11 vistas11 páginas
    Tarea de métodos numéricos, en el método de integración.

    Título original

    Integración.

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    DOCX, PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Tarea de métodos numéricos, en el método de integración.

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd
    Descargar como docx, pdf o txt
    11 vistas11 páginas

    Integración

    Cargado por

    Juan Solis
    Tarea de métodos numéricos, en el método de integración.

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd
    Descargar como docx, pdf o txt
    de 11

    INSTITUTO TECNOLGICO DE LOS MOCHIS

    ING. MECATRNICA
    Mtodo numricos

    Trabajo:
    A5-3. Integral.
    Profesor:
    Marco Antonio Rodrguez Rodrguez
    Alumno:
    Juan Manuel Solis Robles
    Nmero de control:
    15440167
    Carrera:
    Ing. Mecatrnica
    3er Semestre.
    Fecha: 04-12-2016
    1. Evale la integral siguiente:
    f1=(1-e^(-2*x))dx de 0 a 4
    a) En forma analtica; b) con aplicacin mltiple de la regla del trapecio,
    con n=4; c) con una sola aplicacin de la regla de Simpson 1/3, con n=4; d)
    Compare sus resultados usando el comando integral de Matlab u Octave.
    a) En forma analtica
    x+(e^(-2x))/2 de 0 a 4=
    VA1=3.500167731

    VA1 =
    3.5002
    b) Con aplicacin mltiple de la regla del trapecio
    Funcin regla del trapecio:
    ft=@(a,b,f)(b-a)*((f(a)+f(b))/2)

    ft =
    @(a,b,f)(b-a)*((f(a)+f(b))/2)
    Funcin:
    f1=@(x)(1-exp(-2*x))

    f1 =
    @(x)(1-exp(-2*x))

    a1=0

    a1 =
    0
    b1=4

    b1 =
    4
    n1=4

    n1 =
    4

    h1=(b1-a1)/n1

    h1 =
    1
    Primer trapecio
    F11=ft(a1,h1,f1)

    F11 =
    0.4323
    Segundo trapecio
    F12=ft(a1+h1,2*h1,f1)

    F12 =
    0.9232
    Tercer trapecio
    F13=ft(a1+2*h1,3*h1,f1)

    F13 =
    0.9896
    Cuarto trapecio
    F14=ft(a1+3*h1,4*h1,f1)

    F14 =
    0.9986
    Suma de todos los trapecios
    F1=F11+F12+F13+F14

    F1 =
    3.3437
    Error absoluto:
    ea1=abs(((VA1-F1)/VA1)*100)

    ea1 =
    4.4702
    c) Con una sola aplicacin de la regla de Simpson 1/3, con n=4;
    Funcin regla de Simpson de 1/3
    frs=@(a,b,f,xi)(b-a)*((f(a)+4*f(xi)+f(b))/6)

    frs =
    @(a,b,f,xi)(b-a)*((f(a)+4*f(xi)+f(b))/6)
    Xi es el punto medio entre a y b

    xi1=(b1-a1)/2

    xi1 =
    2

    Aplicacin de la funcin
    f1c=frs(a1,b1,f1,xi)

    f1c =
    3.2843

    Error absoluto:
    ea1c=abs(((VA1-f1c)/VA1)*100)
    ea1c =
    6.1683
    d) Compare sus resultados usando el comando integral de Matlab u
    Octave.
    f1d=integral(f1,a1,b1)
    f1d =
    3.5002
    Error absoluto:
    ea1d=abs(((VA1-f1d)/VA1)*100)

    ea1d =
    8.9696e-09
    2. Evale la integral siguiente:
    f2=(1-x-4*x^3+2x^5)dx de -2 a 4
    a) Con la regla del trapecio compuesta, con n=4; b) con la regla de
    Simpson 1/3, con n=4; c) Compare los resultados usando el comando
    integral.
    Funcin:
    f2=@(x)(1-x-4*x.^3+2*x.^5)

    f2 =
    @(x)(1-x-4*x.^3+2*x.^5)
    Datos
    a2=-2

    a2 =
    -2
    b2=4

    b2 =
    4
    n2=4

    n2 =
    4

    a) Con la regla del trapecio compuesta


    Funcin del trapecio compuesta:
    ftc=@(a,b,n,f)(b-a)*((f(a)+2*(sum(f(a2+((b2-a2)/n2):((b2-a2)/n2):b2-
    ((b2-a2)/n2))))+f(b))/(2*n))

    ftc =
    @(a,b,n,f)(b-a)*((f(a)+2*(sum(f(a2+((b2-a2)/n2):((b2-a2)/n2):b2-
    ((b2-a2)/n2))))+f(b))/(2*n))

    F2a=ftc(a2,b2,n2,f2)

    F2a =
    1.5169e+03
    b) Con la regla de Simpson 1/3, con n=4
    frs

    frs =
    @(a,b,f,xi)(b-a)*((f(a)+4*f(xi)+f(b))/6)
    h2=(b2-a2)/n2

    h2 =
    1.5000
    Primer intervalo:
    F2b1=frs(a2,a2+h2,f2,a2+((a2+h2)-(a2))/2)

    F2b1 =
    -2.8066
    Segundo intervalo:
    F2b2=frs(a2+h2,a2+2*h2,f2,(a2+h2)+((a2+2*h2)-(a2+h2))/2)

    F2b2 =
    0.6738
    Tercer intervalo:
    F2b3=frs(a2+2*h2,a2+3*h2,f2,(a2+2*h2)+((a2+3*h2)-(a2+2*h2))/2)

    F2b3 =
    42.9668
    Cuarto intervalo:
    F2b4=frs(a2+3*h2,a2+4*h2,f2,(a2+3*h2)+((a2+4*h2)-(a2+3*h2))/2)

    F2b4 =
    1.0657e+03
    Sumatoria de los intervalos:
    F2b=F2b1+F2b2+F2b3+F2b4

    F2b =
    1.1065e+03
    c) Compare los resultados usando el comando integral.
    F2c=integral(f2,-2,4)

    F2c =
    1104
    Errores absolutos:
    F2=[F2a F2b F2c]

    F2 =
    1.0e+03 *
    1.5169 1.1065 1.1040
    ea=abs(((F2(3)-F2)./F2(3))*100)

    ea =
    37.3981 0.2293 0
    3. Evala la integral de los datos tabulados a continuacin, con a) La regla
    del trapecio compuesta, y b) La regla de Simpson de 1/3.
    Nota: Se puede usar trapz para verificar sus resultados
    x3=[0:0.1:0.5]

    x3 =
    0 0.1000 0.2000 0.3000 0.4000 0.5000
    f3=[1 8 4 3.5 5 1]

    f3 =
    1.0000 8.0000 4.0000 3.5000 5.0000 1.0000
    n3=6

    n3 =
    6
    a) La regla del trapecio
    a3=0

    a3 =
    0
    b3=0.5

    b3 =
    0.5000
    h3=(b3-a3)/5

    h3 =
    0.1000

    F3a1=(h3/2)*(f3(1)+f3(2))

    F3a1 =
    0.4500

    F3a2=(h3/2)*(f3(2)+f3(3))

    F3a2 =
    0.6000

    F3a3=(h3/2)*(f3(3)+f3(4))

    F3a3 =
    0.3750

    F3a4=(h3/2)*(f3(4)+f3(5))

    F3a4 =
    0.4250

    F3a5=(h3/2)*(f3(5)+f3(6))

    F3a5 =
    0.3000

    F3a=F3a1+F3a2+F3a3+F3a4+F3a5

    F3a =
    2.1500
    b) La regla de Simpson de 1/3.
    x3

    x3 =
    0 0.1000 0.2000 0.3000 0.4000 0.5000
    f3

    f3 =
    1.0000 8.0000 4.0000 3.5000 5.0000 1.0000
    Funcin modificada
    frs2=@(fx0,fx2)(0.1)*(fx0+4*(fx0+((fx2-fx0)/2))+fx2)*(1/(6))

    frs2 =
    @(fx0,fx2)(0.1)*(fx0+4*(fx0+((fx2-fx0)/2))+fx2)*(1/(6))

    F3b1=frs2(f3(1),f3(2))

    F3b1 =
    0.4500

    F3b2=frs2(f3(2),f3(3))

    F3b2 =
    0.6000

    F3b3=frs2(f3(3),f3(4))

    F3b3 =
    0.3750

    F3b4=frs2(f3(4),f3(5))

    F3b4 =
    0.4250

    F3b5=frs2(f3(5),f3(6))

    F3b5 =
    0.3000

    F3b=F3b1+F3b2+F3b3+F3b4+F3b5

    F3b =
    2.1500
    Nota: Se puede usar trapz para verificar sus resultados
    trapz(x3,f3)

    ans =
    2.1500
    4. La Ley de Faraday caracteriza la cada de voltaje a travs de un conductor,
    as:
    VL=L(di/dt)
    Con base en la ley de Faraday, use los datos siguientes de voltaje para estimar
    la inductancia en henrios si se pasa durante 400 milisegundos una corriente de
    2A por el inductor.
    T,s 0 10 20 40 60 80 120 180 280 400
    V,volts 0 18 29 44 49 46 35 26 15 7
    ts=[0 10 20 40 60 80 120 180 280 400]

    ts =
    0 10 20 40 60 80 120 180 280 400
    v=[0 18 29 44 49 46 35 26 15 7]

    v =
    0 18 29 44 49 46 35 26 15 7
    F3=(400-0)*(0+2*(18+29+44+49+46+35+26+15)+7)*(1/(2*9))

    F3 =
    11800
    Por Regla del trapecio
    f31=((10-0)/2)*(0+18)

    f31 =
    90
    f32=((20-10)/2)*(18+29)

    f32 =
    235
    f33=((40-20)/2)*(29+44)

    f33 =
    730
    f34=((60-40)/2)*(44+49)

    f34 =
    930
    f35=((80-60)/2)*(49+46)

    f35 =
    950
    f36=((120-80)/2)*(46+35)

    f36 =
    1620
    f37=((180-120)/2)*(35+26)

    f37 =
    1830
    f38=((280-180)/2)*(26+15)

    f38 =
    2050
    f39=((400-280)/2)*(15+7)

    f39 =
    1320
    F3a=f31+f32+f33+f34+f35+f36+f37+f38+f39

    F3a =
    9755
    L en t= 400
    F3a/(7v*1000)=1.4Henrys
    Usando la frmuloa:
    VL=L(di/dt)
    VL(dt/di)=L
    7V(400ms)/2A=1400mH
    La inductancia es la integral del V/A con respecto al tiempo.

    5. si inicialmente un capacitor no tiene carga, el voltaje atravs de l como


    funcin del tiempo se calcula por medio de:
    V(t)=1/Ci(t)dt
    Si C=10^-5 faradios, use los datos de corriente que siguen para elaborar
    una grfica del voltaje versus el tiempo:
    T,s 0 0.2 0.4 0.6 0.8 1 1.2
    I,10^3A 0.2 0.3683 0.3819 0.2282 0.0486 0.0082 0.1441

    x5=[0:0.2:1.2]

    x5 =
    0 0.2000 0.4000 0.6000 0.8000 1.0000 1.2000
    f5=[0.0002 0.0003683 0.0003819 0.0002282 0.0000486 0.0000082 0.0001441]

    f5 =
    1.0e-03 *
    0.2000 0.3683 0.3819 0.2282 0.0486 0.0082 0.1441

    Funcin regla del trapecio:


    ft5=@(a,b,fa,fb)(b-a)*(1/2)*(fa+fb)*(1/(1e-5))

    ft5 =
    @(a,b,fa,fb)(b-a)*(1/2)*(fa+fb)*(1/(1e-5))

    Lo que se har en este problema, es obtener los voltajes en cada intervalo, esto
    mediante la integracin de cada intervalo, usando la frmula dada para el
    voltaje.
    v1=ft5(x5(1),x5(2),f5(1),f5(2))

    v1 =
    5.6830

    v2=ft5(x5(2),x5(3),f5(2),f5(3))
    v2 =
    7.5020

    v3=ft5(x5(3),x5(4),f5(3),f5(4))

    v3 =
    6.1010

    v4=ft5(x5(4),x5(5),f5(4),f5(5))

    v4 =
    2.7680

    v5=ft5(x5(5),x5(6),f5(5),f5(6))

    v5 =
    0.5680

    v6=ft5(x5(6),x5(7),f5(6),f5(7))

    v6 =
    1.5230

    VT=v1+v2+v3+v4+v5+v6

    VT =
    24.1450
    Ahora para graficar crear una nueva matriz para ordenar los valores del voltaje
    acumulado, y junto con los ya existentes del tiempo graficar.
    V=[0 v1 v1+v2 v1+v2+v3 v1+v2+v3+v4 v1+v2+v3+v4+v5 v1+v2+v3+v4+v5+v6]

    V =
    0 5.6830 13.1850 19.2860 22.0540 22.6220 24.1450
    x5

    x5 =
    0 0.2000 0.4000 0.6000 0.8000 1.0000 1.2000
    Ahora grafcar con splines:
    X=0:0.001:1.2;
    Y=spline(x5,V,X);
    plot(x5,V,'pr',X,Y,'b')

    También podría gustarte