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

Diff - Derivada Simbólica

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 5

SOLUCIONES DE PROBLEMAS CON MATLAB.

pdf
DERIVADA
Usamos la funcin diff para determinar la derivada simblica de una expresin simblica.
Hay cuatro formas de usar la funcin diff para realizar una derivacin simblica:
diff (f ) : Devuelve la derivada de la expresin f respecto a la variable independiente por omisin.
diff (f , t) : Devuelve la derivada de la expresin f respecto a la variable t.
diff (f , n) : Devuelve la n-sima derivada de la expresin f respecto a la variable independiente por
omisin.
diff ( f, t, n) : Devuelve la n-sima derivada de la expresin f respecto a la variable t.

Puesto que la funcin diff tambin se usa para la derivacin numrica, usted podra preguntarse cmo sabe la
funcin si debe calcular diferencias numricas o realizar derivacin simblica. La funcin puede determinar
qu es lo que se desea analizando los argumentos de entrada: si el argumento es un vector, calcula
diferencias numricas; si el argumento es una expresin simblica, realiza derivacin simblica.
Ahora presentamos varios ejemplos del uso de la funcin [diff] para derivar simblicamente. Suponga que
se han definido las siguientes expresiones:

La siguiente lista muestra referencias de funcin y sus valores correspondientes:


Referencia Valor de la funcin

diff (S1) 18*x^2-8*x+b


diff (S1,2) 36*x-8
diff (SI, 'b') x
diff (S2) cos (a)
diff (S3) -3*t^2/(l+t^4)-4*(1-t^3)/(1+t^4)A2*t^3
simplify(diff(S3)) t^2*(-3+t^4-4*t)/(l+t^4)^2

INTEGRAL:
Usamos la funcin int para integrar una expresin simblica f. Esta funcin intenta encontrar la expresin
simblica F tal que diff (F) = f. Es posible que la integral (o antiderivada) no exista en forma cerrada o que
MATLAB pueda obtener la integral.
En estos casos, la funcin devolver la expresin sin evaluarla. La funcin [int] puede usarse en las
siguientes formas:
int (f) Devuelve la integral de la expresin f respecto a la variable independiente por omisin.
int (f,'t') Devuelve la integral de la expresin f respecto a la variable t.
int(f, a, b) Devuelve la integral de la expresin f respecto a la variable independiente por omisin
evaluada en el intervalo [a, b], donde a y b son expresiones numricas.
int (f, 't', a, b) Devuelve la integral de la expresin f respecto a la variable t evaluada en el intervalo [a, b],
donde a y b son expresiones numricas.
int (f. 'm , n) Devuelve la integral de la expresin f respecto a la variable independiente por omisin
evaluada en el intervalo [m, n], donde m y n son expresiones simblicas.

Para evitar posibles problemas, es recomendable especificar la variable independiente en la derivacin


simblica y en la integracin simblica.
Ahora presentamos varios ejemplos del uso de la funcin int para la integracin simblica. Suponga que se
han definido las siguientes expresiones:
S1 = '6*X^3-4*X^2+b*x-5';
S2 = 'sin(a)';
S3 = 'sqrt (x)' ;
La siguiente lista muestra referencias de funcin y sus valores correspondientes:

Referencia Valor de la funcin


int (SI) 3/2*x^4-4/3*x^3+1/2*b*x^2 - 5*x
int(S2) -cos (a)
int(S3) 2/3*X^(3/2)
int(S3, 'a', 'b') 2/3*b^(3/2) - 2/3*a^(3/2)
int(S3, 0.5, 0.6) 2/2 5 *1 5^(1/2) - 1/6 * 2^(1/2)
numeric (int (S3, 0.5, 0.6) ) = 0. 0741

Resolucin aplicada de problemas: Globos meteorolgicos


Se usan globos meteorolgicos para obtener datos de temperatura y presin a diferentes alturas en la
atmsfera. El globo se eleva porque la densidad del helio en su interior es menor que la del aire que rodea al
globo. Al subir el globo, el aire circundante se vuelve menos denso, y el ascenso se va frenando hasta que el
globo alcanza un punto de equilibrio. Durante el da, la luz del Sol calienta el helio atrapado dentro del
globo; el helio se expande y se vuelve menos denso, y el globo sube ms. Durante la noche, en cambio, el
helio del globo se enfra y se vuelve ms denso, y el globo desciende a una altura menor. El da siguiente, el
Sol calienta el helio otra vez, y el globo sube. Este proceso genera una serie de mediciones de altura con el
transcurso del tiempo que se pueden aproximar con una ecuacin polinmica.
Suponga que el siguiente polinomio representa la altura en metros durante las primeras 48 horas despus del
lanzamiento de un globo meteorolgico:

Donde las unidades de t son horas. Genere curvas para la altura, velocidad y aceleracin de este globo
usando unidades de metros, m/s y m/s2. Adems, determine y exhiba la altura mxima y su hora
correspondiente.
PROBLEMAS:
a) Flujo de agua: Suponga que se bombea agua a un tanque que inicialmente est vaco. Se sabe que la
tasa de flujo de agua hacia el tanque en el instante t (en segundos) es de 50 litros por segundo. Puede
demostrarse que la cantidad de agua Q que fluye hacia el tanque durante los primeros X segundos es igual a
la integral de la expresin 50 - t evaluada de 0 a x segundos.
1. Determine una ecuacin simblica que represente la cantidad de agua en el tanque despus de x segundos.
2. Determine la cantidad de agua en el tanque despus de 30 segundos.
3. Determine la cantidad de agua que fluy al tanque entre los 10 y los 15 segundos despus de iniciarse
el flujo.
b) Resorte elstico. Considere un resorte con el extremo izquierdo fijo y el derecho libre para moverse
en el eje x. Suponemos que el extremo derecho est en el origen x = 0 cuando el resorte est en reposo.
Si el resorte se estira, su extremo derecho est en algn nuevo valor de x mayor que 0; cuando el resorte se
comprime, su extremo derecho est en algn valor de x menor que 0. Suponga que el resorte tiene una
longitud natural de 1 ft y que se requiere una fuerza de 10 lb para comprimir el resorte a una longitud de 0.5
ft. Entonces, puede demostrarse que el trabajo (en ft/lb) realizado para estirar este resorte desde su longitud
natural hasta una longitud total de n ft es igual a la integral de 20 x en el intervalo de 0 a n - 1.

4. Use MATLAB para determinar una expresin simblica que represente el trabajo necesario para estirar el
resorte a una longitud total de n ft.
5. Calcule el trabajo realizado al estirar este resorte a una longitud de 2 ft.
6. Si la cantidad de trabajo ejercido es de 25 ft/lb, qu longitud tiene el resorte estirado?

Fx = input('Ingrese la funcin (F(x))= ');


xo = input('Ingrese el valor inicial(Xo)= ');
syms x;
dFx = diff(Fx, x);
Fx = inline(Fx);
%dFx=inline(dFx);
F = feval(Fx, xo);
dF = subs(dFx, xo);
disp (dFx); disp ('F(xo)= '); disp (F); disp ('dF(xo)= '); disp (dF);

También podría gustarte