Ecuaciones Matlab

TEORA DE ECUACIONES CON MATLAB
SESIN 01
1. ENTORNO DE TRABAJO DEL MATLAB
Los nmeros reales se pueden insertar tambin en notacin cientfica, muy adecuados si se traza de
nmeros grandes o muy pequeos (en valor absoluto). As, se tiene la siguiente regla de
construccin para representar nmeros reales:
m.10 r entonces su equivalente es x
Si x
0.007
7e 3
453.1012
7 x10
m.e.r as como por ejemplo:
453e12
1201200000
0.00031415
1.2012e009
3.1415e 004
Existen adems dos nmeros especiales: Inf y NaN . El primer signo representa la cantidad infinita
(). La segunda, es una abreviatura de no es un nmero (Not a Number) y es el resultado que se
devuelve ante una operacin indefinida como 0/0.
2. OPERADORES ARITMTICOS
Los operadores aritmticos que se posee MATLAB se ven a continuacin:
Operacin
Descripcin
Elementos a y b
a+b
: Realiza la suma de a y b
Escalares, vectores o matrices
a-b
: Realiza la resta de a y b
Escalares, vectores o matrices
a*b
: Realiza la multiplicacin de a y b
Escalares o matrices
.*
a .* b : Realiza la multiplicacin de a y b
Escalares o de vectores
a / b : Realiza la divisin de a entre b
Matrices
./
a / b : Realiza la divisin de a entre b
Vectores
a \ b : Realiza la divisin de b entre a
Escalares o de matrices
.\
a .\ b : Realiza la divisin de b entre a
a ^ b : Eleva la base a al exponente b
.^
a .^ b : Eleva la base a al exponente b
Vectores
Escalares o escalar de matrices ( M p )
vectores
3. OPERADORES RELACIONALES
Los operadores relacionales que se posee MATLAB se ven a continuacin:
Operacin
Descripcin
~=
a ~= b : Establece la condicin de a
>
a>b
: Establece la condicin de a mayor que b
<
a<b
: Establece la condicin de a menor que b
>=
a >= b : Establece la condicin de a mayor o igual que b
<=
a <= b : Establece la condicin de a menor o igual que b
Dr. SORIA QUIJAITE JUAN JESS
a|b
: Establece la disyuncin de a o b
==
a==b
: Establece la condicin de a igual a b
4. FORMATOS
MATLAB, presenta los valores de acuerdo a un formato, el que por defecto es
format short ; existen otros
Valor de pi
Comandos
long
3.14159265358979
short e
3.1416e+000
long e
3.141592653589793e+000
hex
400921fb54442d18
bank
3.14
rat
355/113
short
3.1416
5. OPERADORES LGICOS
Los operadores relacionales que se posee MATLAB se ven a continuacin:
Descripcin
Operacin
Negacin lgica (NOT) o complementario de A
~A
A&B
Conjuncin lgica (AND) o interseccin de A y B
A|B
Disyuncin lgica inclusiva (OR) o unin de A y B
xor(A,B)
OR exclusivo (XOR) o diferencia simtrica de A y B
6. VARIABLES
En Matlab como en cualquier otro lenguaje de programacin se utilizan variables, estas deben tener
un nombre segn ciertas reglas las cuales son:
No pueden comenzar con un nmero, aunque si pueden tener nmeros en su estructura:
variable1 es un nombre vlido.
Las maysculas y minsculas se diferencian en los nombres de variables: A y a son variables
diferentes.
Los nombres de variables no pueden contener operadores ni puntos. No es vlido usar / , * , ,+,. , ; , : ,^.
Para el uso de una variable no es necesario declarar sus nombres , en la siguiente tabla se presenta
las variables predefinidas que posee Matlab.
Nombre de la
Significado
variable
Pi
i y j
Unidad imaginaria =
Inf
Infinito =
eps
psilon de la mquina =1.0000e-006
NaN
No es un nmero
realmin
Menor nmero 2
relamas
Mayor nmero (2 e)21023
Date
Fecha
flops
Contador de operaciones de punto flotante
nargin
Nmero de argumentos de entrada de una funcin
nargout
Nmero de argumentos de salida de una funcin
1022
7. COMANDOS DE LECTURA Y ESCRITURA

Matlab provee una forma sencilla de leer variables desde el teclado y visualizar mensajes en la
pantalla de la computadora a travs de las siguientes funciones:
input .- Permite el ingreso de datos al programa a travs del teclado asignndolo a una
variable, esta orden puede usarse con un mensaje en la lnea de comandos. Despus de
imprimir el mensaje, la orden espera que el usuario digite el valor numrico, un vector, una
matriz o una expresin vlida del matlab.
Ejemplo:
z= input(ingrese un nmero, s)
Asigna a la variable z la cadena ingresada
s : indica que la entrada que se har por el teclado es una cadena
fprintf.- Permite la visualizacin de un valor numrico o el resultado de una expresin
guardada por el usuario.
Ejemplo:
>> vol=25
>>fprintf (el volumen de la esfera es:, %12.0f \ n, vol )
\ n
: indica que la impresin de la variable vol ser en la siguiente lnea
%12.0f
: formato de un nmero entero
%12.5f
: formato de un nmero real con 5 decimales.
disp.- Permite visualizar en pantalla un mensaje de texto o el valor de una matriz, pero sin
imprimir su nombre. En realidad, disp siempre imprime vectores y/o matrices, las cadenas de
caracteres se consideran un caso particular de vectores.
Ejemplos:
>> disp(Esta es una prueba)
>> disp(pi)
>> disp(El programa ha terminado)
>> A=rand(4,4)
>> disp(A)
clear : Borra las variables usadas de la memoria
clc
: Limpia la informacin de la ventana de comandos
8. FUNCIONES MATEMTICAS EN MATLAB

Matlab ofrece un sinnmero de funciones las que aceptan como argumento variables reales
y/o complejas sin discriminacin, as como con argumentos matriciales.
FUNCIONES TRIGONOMTRICAS
FUNCIN- Matlab
DESCRIPCIN
FUNCIN INVERSA
sin(x)
Seno de x
asin(x)
cos(x)
Coseno de x
acos(x)
tan(x)
Tangente de x
atan(x) y atan2(x)
cot(x)
Cotangente de x
acot(x)
sec(x)
Secante de x
asec(x)
csc(x)
Cosecante de x
acsc(x)
FUNCIONES HIPERBLICAS
FUNCIN- Matlab
DESCRIPCIN
FUNCIN INVERSA
sinh(x)
Seno hiperblico de x
asinh(x)
cosh(x)
Coseno hiperblico de x
acosh(x)
tanh(x)
Tangente hiperblica de x
atanh(x)
coth(x)
Cotangente hiperblica x
acoth(x)
sech(x)
Secante hiperblica de x
asech(x)
csch(x)
Cosecante hiperblica x
acsch(x)
FUNCIONES EXPONENCIALES Y LOGARTMICAS
FUNCIN- Matlab
DESCRIPCIN
exp(x)
Funcin exponencial en base e (e^x)
f ( x)
ex
log(x)
Funcin logaritmo en base e de x
f ( x)
Ln( x)
log10(x)
Funcin logaritmo en base 10 de x
Lenguaje matemtico
f ( x)
Log ( x)
log2(x)
Funcin logaritmo en base 2 de x
f ( x)
Log2 ( x)
pow2(x)
Funcin potencia de base 2 de x
f ( x)
2x
Funcin raz cuadrada de x
f ( x)
sqrt(x)
FUNCIONES ESPECFICAS DE VARIABLE NUMRICA

DESCRIPCIN
FUNCIN- Matlab
abs(x)
Valor absoluto del real x
floor(x)
El mayor valor entero o igual que el
Lenguaje matemtico
f ( x)
f ( x)
Sing ( x)
real x
ceill(x)
El menor entero mayor o igual que el

real x
round(x)
El entero ms prximo al real x
rem(a,b)
Da el resto de la divisin entre los

reales a y b
sing(x)
fix(x)
Signo del real x (1 si x>0 ; -1 si x<0)

Elimina la parte decimal del real x
format rat
pi
%ans =355/113
format short
pi
%ans =3.1416
format short e
pi
%ans =3.1416e+000
format bank
pi
%ans = 3.14
format long
pi
%ans =3.14159265358979
format long e
pi
%ans =3.141592653589793e+000
str2num('15/14')
%ans=1.071428571428571e+000
num2str(pi)
%ans=3.1416
str2mat('IEP "DE LA CRUZ"','MATEMATICA','LOS TRIANGULOS','PROPIEDADES','EJERCICIOS')

%ans =IEP "DE LA CRUZ"
MATEMATICA
LOS TRIANGULOS
PROPIEDADES
EJERCICIOS
SISTEMAS DE NUMERACIN
Matlab permite trabajar con sistemas de numeracin de base cualquiera.
dec2base(decimal, n_base)
Convierte el nmero decimal (base 10) especificado a la
nueva base n base dada
base2dec(numero,B)
Convierte el nmero dado en base B a decimal
dec2bin(decimal)
Convierte el nmero decimal especificado a base 2(binario)
dec2hex(decimal)
Convierte el nmero decimal especificado a base 16
(hexadecimal)
bin2dec(binario)
Convierte el nmero binario especificado a base decimal
hex2dec(hexadecimal)
Convierte el nmero base 16 especificado a base decimal
Transformar el numeral 324 al sistema
binario
dec2base(324,2)
% ans =101000100
4.
Transformar el numeral 3121(4) al

sistema decimal
base2dec('3121',4)
ans = 217
Transformar el numeral 3242 al
sistema binario
dec2bin(3242)
ans =110010101010
Transformar el numeral 125 al sistema
hexadecimal
dec2hex(125)
ans =7D
5.
6.
123 a base 6 :
254 a base 7: .
3.
Convertir:
432(5) a base 7 :
Colocar V o F segn corresponda:

I. 27 = 102(5) (
)
II. 57 = 321(6) (
)
III. 10 = 1010(2) (
)
IV. 22 = 113(4) (
)
Colocar > ; < = segn corresponda:
16(7)
15(8)
23(5)
23(6)
28(9)
121(4)
Si:
mnp (8)
= 312(7)
Hallar: m + n + p
a) 7
b) 8
d) 10
e) 11
7.
Transformar el numeral 101101(2) al

sistema decimal
bin2dec('101101')
ans = 45
EJERCICIOS DE APLICACIN
1. Relacionar
ambas
columnas
adecuadamente:
I)
23(5)
(
) 15
II) 15(7)
(
) 13
III) 33(4)
(
) 12
2. Convertir:
202(3) a base 8 :
8.
9.
Si: abc (9) = 175

Hallar: a + b + c
a) 3
b) 5
d) 9
e) 11
Hallar x si:
xxx = 4210(5)
a) 3
d) 6
b) 4
e) 7
Convertir:
A. 1023(5) a base 25
a) 513(25) b) 5(13) (25)
d) 512(25) e) 5(12) (25)
B. 11102(3) a base 9
a) 442(9) b) 142(9)
d) 342(9) e) 742(9)
c) 9
c) 7
c) 5
c) 6(13) (25)
c) 332(9)
10. Si: N = 73 x 5 + 72 x 4 + 7 x 3 + 9
Convertir N a base 7
a) 5439(7) b) 5432(7)
c) 5442(7)
d) 5437(7) e) 5449(7)
a) 7542(8) b) 5472(8)
d) 7564(8) e) 8564(8)
c) 754(20)(8)
11. Si: N = 83 x 7 + 82 x 5 + 8x 4 + 20
Convertir N a base 8.
SESIN 02
MATLAB es una abreviatura de Matriz Laboratory (Laboratorio de Matrices). La implementacin
original de Matlab la realiz Cleve Moler a finales de los aos 70, las aplicaciones de Matlab se
fueron extendiendo a otras ramas del clculo cientfico y de las ciencias aplicadas en general. Dichas
extensiones se consiguieron en gran parte mediante la implementacin de toolboxes, libreras que
utilizan Matlab para ampliar el rango de problemas que puede cubrir.
1. ENTORNO DE TRABAJO DEL MATLAB

Entre las primeras caractersticas de Matlab destacamos las siguientes:
El prompt de Matlab es >>. El usuario escribe a continuacin
Para ejecutar se pulsa la tecla Enter
Se puede recuperar comandos anteriores navegando con las flechas y
Cuando se trabaje con Matlab, debemos tener en cuenta que:
Se distinguen las maysculas de las minsculas
El carcter % se utiliza para insertar comentarios. Todo lo que sigue es ignorado por Matlab
Si se teclea al final de una orden ; sta se ejecuta pero el resultado no se visualiza por
pantalla.
RESOLVIENDO ECUACIONES POLINMICAS CON MATLAB
I.- El comando roots,
determina las
races(soluciones) de un polinomio de
cualquier grado n.
Su sintaxis es: r=roots(p)
Descripcin: r=roots(p) retorna un vector
columna cuyos elementos son las races del
polinomio p
E-1) Resolver la ecuacin cuadrtica
2 x 2 13 x 15
Resolucin
%Digitar en MATLAB
p1=[2 -13 15];
roots(p1)
%resultado
ans =
5.0000
1.5000
El conjunto solucin es CS {5; 32 }
E-2) Resolver la ecuacin cuadrtica
2 x2
3x 6
Resolucin
%Digitar en MATLAB
P2=[sqrt(2) -sqrt(3) 6];
roots(p2)
%resultado
ans =
0.6124 + 1.9666i
0.6124 - 1.9666i
E-3) Resolver la ecuacin polinmica
2x4
3x3 12 x 2
7x 6
Resolucin
%Digitar en MATLAB
p3=[2 -3 -12 7 6];
roots(p3)
%Resultado
ans = 3.0000
-2.0000
1.0000
-0.5000
E-4) Resolver la ecuacin polinmica
12 x5 8 x 4
45 x3
45 x 2
8 x 12
Resolucin
%Digitar en MATLAB
p3=[12 -8 -45 45 8 -12];
roots(p3)
%Resultados
ans = -2.0000
-0.5000
1.5000
1.0000
0.6667
II.- El comando solve,
determina las
soluciones de una ecuacin algebraica en una
variable as como resuelve sistemas de
ecuaciones lineales y no lineales.
Su sintaxis es: g=solve(eq)
g = solve(eq,var)
g = solve(eq1,eq2,...,eqn)
g= solve(eq1,eq2,...,eqn,var1,var2,...,varn)
Descripcin:
g=solve(eq) retorna las
soluciones de una ecuacin algebraica.
E-1) Resolver la ecuacin algebraica
3x 2
x 3
x 3 1 0
%Digitar en MATLAB
syms x
solve(sqrt(3*x-2)-sqrt(x+3)-1)
%resultado
ans = 6
x2
x2
21x 90
3 x 54
x2
21x 90
3x 54
x 6
x 6
%Digitar en MATLAB
syms x
solve(sqrt(x.^2-21*x+90)-sqrt(x.^2+3*x-54)x+6)
%Resultado
ans = 6
1
x
1
x 10
1
12
1
x 10
1
12
%Digitar en MATLAB
syms x
solve((1./x)+(1./(x+10)-(1/12)))
%resultado
ans =
-6
20
E-4) Resolver el sistema de ecuaciones
x2
x
xy
4x 3
15
xy 6
(5;2), (5; 2), ( 5;2), ( 5; 2)

GUA DE EJERCICIOS PRCTICOS
1.
2.
3.
4.
5.
2 x 2 3x 2 0
3x3 6 x 2 5 x 2 0
12 x 4 4 x 3 41x 2 4 x 12 0
6 x 5 41x 4 97 x 3 97 x 2 41x 6
x 6 28 x 3 27 0
2
6. 2 x
2x 3 x2 x 3
7. (2 x 7)( x 2 9)( 2 x 5)
8.
1 x
x
1 x
x
24
x
4x
5
4
1
4
y 1
7
y 1
3
91
13
6
6
x
1
x
10.
4
x
1
x
11.
x
1
y
y
20
xy
1
2
ACTIVIDADES DOMICILIARIAS
I. Resolver las siguientes ecuaciones y sistemas:
1.- 2 x 4 3x 3 13 x 2 3 x 2 0
5
4
3
2
2.- 2 x 3 x 5 x 5 x 3 x 2 0
3.- x 7 2 x 5 3 x 4 3x 2 3x 1 0
2
4x 5 0
4.- 3 x
3
2
5.- 3 x 2 x 3 x 4 0 3
Resolucin
%Digitar en MATLAB
syms x y
[x,y]=solve('x^2 + x*y + y=3','x ^2 - 4*x + 3 = 0')
%resultados
x=1
3
y=1
-3/2
xy
I. Resolver las siguientes ecuaciones y sistemas:
9.
Resolucin
Lo llevamos a su equivalente
1
x
CS
Resolucin
x2
x2
Resolucin
syms x y
[x,y]=solve('x^2-x*y-15=0','y^2 -x*y+6=0')
x= 5
-5
y=2
-2
Luego el conjunto solucin es:
Resolucin
3x 2
E-5) Resolver el sistema de ecuaciones
6.
7.
2 x (1 2 y )
2
x 1
yx
8.
x2 y
x3
(1 y )
x
20
8
5x4
x
2
xy2
yx
y3
5x2
x
9.
y 2 xy
x
36
x2
10.
12.
y2
z2
z
xy
( x2
y 2 )( x
14
x 1
11.
x
y)
y 2
x
y
y
x
15
2
65
y 13
SESIN 03
PLOTEANDO FUNCIONES REALES EN 2D

CON EL SOFTWARE MATLAB
Funcin Matemtica
En MATLAB
f ( x)
f ( x)
sqrt ( x)
f ( x)
f ( x)
abs ( x)
f ( x)
f ( x)
f ( x)
x
log( x)
Ln( x)
f ( x)
x.^ (1/n)
f ( x)
2x
f ( x)
f ( x)
f ( x)
log 10( x)
log( x)
pow2( x)
f ( x)
f ( x)
f ( x)
f ( x)
ax
sen ( x)
cos( x)
tan( x)
f ( x)
a. ^ ( x)
f ( x)
ex
f ( x)
f ( x)
f ( x)
f ( x)
sin( x)
cos( x)
tan( x)
exp( x)
f ( x)
f ( x)
x.^ (n)
13.
ylabel('EJE DE ORDENADAS')
title('PLOTEANDO CON SORIA')
gtext('mximo')
gtext('y=-3*x.^2-5*x+12')
%resultando el grfico
E-2)
I. El comando plot grafica el vector X versus el

vector Y. Si X e Y son matrices entonces el
vector es ploteado por filas y columnas de
matriz. La sintaxis es:
Graficar la funcin f ( x) 2 x
Resolucin
%Digitar en MATLAB
x=-6:0.05:6;
y=2*x.^3-2*x.^2-4*x+8;
plot(x,y,'r')
grid on
2x
plot(X1,Y1,...)
plot(X1,Y1,LineSpec,...)
plot(...,'PropertyName',PropertyValue,...)
plot(axes_handle,...)
h = plot(...)
hlines = plot('v6',...)
E-1) Graficar la funcin f ( x)
Resolucin
%Digitar en MATLAB
x=-9:0.05:9;
y=-3*x.^2-5*x+12;
plot(x,y,'r')
grid on
xlabel('EJE DE ABSCISAS')
3x 2
5 x 12
gtext('raiz')
gtext('y=2*x.^3-2*x.^2-4*x+8')
%Resultando
4x 8
E-3) Interceptar la funciones reales
f ( x)
g ( x)
cos( sen ( x)) sen (tan( x)) 2 x

e
2x
.sen ( x)

gtext('Total=72900')
gtext('f(x)=100*(3.^x)')
x 4
Resolucin
%Digitar en MATLAB
x =-pi:0.005:pi;
f=cos(sin(x))-sin(tan(x))-2*x.^2-3;
g=exp(-2*x).*sin(x)+x.^2-x+4;
plot(x,f,'r',x,g,'k')
grid on
gtext('f=cos(sin(x))-sin(tan(x))-2*x.^2-3')
gtext('g=exp(-2*x).*sin(x)+x.^2-x+4;')
%Resultando
E-4. Se dispone de 24 m de malla metlica

para construir un corral rectangular. Cules
deben ser las medidas del corral para que este
tenga la mayor superficie (rea)?
Resolucin
x
y
Permetro
2( x
x+y=12
y=12-x
Area
E-3) Supongamos que cierto tipo de bacterias

se reproduce al triple cada media hora. Qu
cantidad habr al cabo de tres horas si
inicialmente hay 100 bacterias?
Resolucin
Sea el siguiente cuadro
Tiempo transcurrido
Nmero de bacterias
0 hora
100= 30.100
1
hora
300= 31.100
2
1 hora
900= 32.100
1 12 hora
2700= 33.100
2 hora
4
8100= 3 .100
y ) , entonces 2(x+y)=24
Area
x. y
x 2 12 x
x.(12 x)
Luego para X=6 se tiene mayor superficie

%Ploteando en MATLAB
x=-3:0.05:15;
y=-x.^2+12*x;
plot(x,y,'b')
grid on
xlabel('EJE X')
ylabel('EJE Y ')
gtext('MXIMO=36')
gtext('y=-x.^2+12*x')
x horas
3 x.100
Luego la regla de correspondencia

es : p( x) 100.3x
Como en tres horas hay seis medias horas,
entonces p(6) 100.36 100.729 72 900
%ploteo de MATLAB
x=0:0.05:6;
y=100*(3.^x);
plot(x,y,'b')
grid on
xlabel('HORAS')
ylabel('POBLACIN EN MILES ')
10
GUA DE EJERCICIOS
I. Plotear las siguientes funciones reales
1. f ( x)
2.
( x)
4x
12 x 24 ; g ( x)
4 x3 5 x 2
2
x2 3
4
x 1
2 x 46
3.
( x)
4x 2
4x2
x2
4.
( x)
sen ( x 2 )
2x2
log x 2
5.
( x)
sen( x)
6. g ( x)
x2
3x 4
sen( x3
2x2
3x 1) x 2
x3
7. g ( x)
senx
3x
x2
8. g ( x) 5 sen( x 3x 1) ln( x
2) 4
9. Los registros de temperatura tomados entre
las
0h y las 24 h en una zona rural de Ica se
ajustan a la funcin T ( x)
Las
curvas
f ( x)
c x2
con
ecuaciones
3.7
4.1
12.3
Log 83.7
Log 2
4.4
4.2
17.1
6.8 .Log 7
4.8 .Log 6
(a b 3)3 3(a b 2)3 3(a b 1)3

Resolucin
%MATLAB
syms a b p
simplify(p)
%ans = 6
2. Simplificar
R
c 3 (a b)3 a 3 (b c)3 b3 (c a)3

(a b)3 (b c)3 (c a)3
Resolucin
syms a b c n d R
4.3
3.7
7.3
%MATLAB
%MATLAB
En Matlab el mximo entero del nmero real x
lo determina el comando floor(x).
Ejemplo
n=c.^3*(a-b).^3+a.^3*(b-c).^3+b.^3*(c-a).^3;
d=(a-b).^3+(b-c).^3+(c-a).^3;
R=n./d;
5.2
pretty(R)
7.1
simplify(R)
%Sentencia en Matlab
sqrt(floor(4.3)+floor(5.2))) / (floor(-3.7)+floor(-7.1))
%ans = -1/4
Simplificar las siguientes expresiones
11
3.3
p=(a+b+3).^3-3*(a+b+2).^3+3*(a+b+1).^3-(a+b).^3;
Definicin.- Si x IR, el mximo entero que se

x
denota por
se define por el smbolo
x
n
, donde n es el mayor entero menor
que, o igual a x. Es decir:
x
n
n x n 1 ,n Z
Ejemplos:
5.3
5
Calcular el valor de:
Log 2
2.2
PRODUCTOS NOTABLES
1. Reducir
se llaman curvas de nariz de
EL MXIMO ENTERO
5.7
D) M
0.6
19.1
B) M
C) M
10.8
6.7
bala. Grafique algunas para ver el porqu de

este nombre. Qu sucede al crecer c ?
5.3
6.5
1
( x 12) 2 10 ,
10
donde T es la temperatura en grados

centgrados y x es la hora del da.
a) Cul fue la temperatura mxima? A qu
hora se registr?
b) A que hora la temperatura fue de 0 C?
c) Qu temperatura se registro a las tres de
la tarde?
10. Trace la grfica de la funcin
f ( x) x n .2 x ; x 0 , para n=1,2,3,4,5 y 6.
Cmo cambia la grfica al crecer n?
11.
9.1
A) M
3
3
3
3 3
3
c (a - b) + a (b - c) + b (c - a)
--------------------------------------3
3
3
(a - b) + (b - c) + (c - a)
%ans =c*a*b
( a b) 3
3. Calcular ExF si las expresiones son
( x 3)( x 5)( x 2)( x 4) ( x 2
x 13)2 50
( x 2 5x 5)2 ( x 1)( x 2)( x 3)( x 4)
subs(f,2)+subs(f,-3)
%ans = 47
%realiza lo mismo
>> f='3*x^2-4*x+2';
Resolucin
subs(f,2)+subs(f,-3)
%MATLAB
%ans = 47
syms x E F P
E=(x-3).*(x-5).*(x+2).*(x+4)-(x.^2-x-13).^2+50;
x3
2. Sea f ( x)
4x2
f (a 1)
F=(x.^2+5*x+5).^2-(x+1).*(x+2).*(x+3).*(x+4);
P=E.*F;
Resolucin
pretty(P);
syms x a
simplify(P)
f=x^3+4*x^2-6*x+5;
2
2
% ((x - 3) (x - 5) (x + 2) (x + 4) - (x - x - 13) + 50)
2
2
% ((x + 5 x + 5) - (x + 1) (x + 2) (x + 3) (x + 4))
% ans =1
p=subs(f,a+1)
1) Simplificar
(2 x 3 y )
la
siguiente
(3 y 2 x)
6 xy
2) Simplificar
6 x 5 . Calcular
la
expresin
%p =(a+1)^3+4*(a+1)^2-6*a-1
3
2
%
(a + 1) + 4 (a + 1) - 6 a - 1
%ans =a^3+7*a^2+5*a+4
siguiente
la
simplify(p)
expresin
( x y 3)( x y 3) 9 ( x y )2 4 xy
3) Simplificar
pretty(p)
siguiente
expresin
4. Sea f ( x, y, z )
f ( 2;2;3)
x2
y2
z 2 . calcular
f ( 1; 2;1)
( x 9) ( x 13)( x 5)
( x 10)( x 9) ( x 16)( x 3)
4) Simplificar
la
siguiente
Resolucin
expresin
(a b c)3 (a b)3 3(a b c)(a b)c
5) Simplificar
a b
a b
la
siguiente
2
a b
a b
expresin
a b
a 2 b2
15
7) Sea x
Hallar el
valor
15
3
p=subs(f,{x,y,z},{-2 2 3})+subs(f,{x,y,z},{-1 -2 1})
1. Calcular el valor numrico de
B x3 3x 2 y 3xy 2 y 3 ; para x= -1 ;
y=1
P( x) 4 x 2 4 x 1 ,
2. Si
calcular
numrico
de
donde
x3 3x 23
VALOR NUMRICO
2
4 x 2 . Hallar f(2)+f(-3)
1. Sea f ( x ) 3 x
Resolucin
syms x
f=3*x^2-4*x+2;
12
f=x^2+y^2+z^2;
% p =23
6) Sea P( x) ( x 1)( x 1)( x 2 x 1)( x 2 x 1)

Hallar el valor numrico de P(x) para
syms x y z
P P
1
2
P( x) x 2 6 x 9 ,
calcular
P (5) P (4)
E
P (4)
4. Si P ( x) 2 x 2 53 x 3 13 , calcular P ( 21 )
3. Si
P( x) x 2 3x 1 ,
P( 2) P( 1)
P(4) P(3)
5. Si
6. Dados
los
A( x) 1 x x
7. Hallar P(0)+P(1) si
calcular
P( x)
x 2 3x 5
C ( x)
4 5x
EF
F 2 , donde E=x+3 ;
F=2-x+x2 .
polinomios
P( x 3) 5 x 7
B( x)
E2
8. Si
Adems P M ( x ) 3 15x 2
Hallar el valor P M (1)
Evaluar P(-1/7), si P(x)=A+B-C

SESIN 04 - POLINOMIOS
Polinomio.- Un polinomio es una funcin de una sola variable que se puede expresar en la
siguiente forma general, f ( x) a0 x N a1 x N 1 a2 x N 2 aN 2 x 2 aN 1 x aN donde:
a0 ; a1 ; a2 ;; aN = Coeficientes ; x= variable
Si a0
aN
; N=grado de f(x) .
1,entonces f(x) es mnico.

T .I . = trmino independiente ;
coeficientes
f (1)
Tr min o.Ind . f
f ( 0)
Grado del polinomio.- El grado del polinomio es igual al valor ms alto empleado como
exponente.
Evaluacin de polinomios.- Existen varias formas de evaluar un polinomio para un conjunto
de valores usando MATLAB . Los polinomios pueden evaluarse con el comando polyval.
E-1) Sea el polinomio p( x) 3x 4 7 x3 2 x 2 x 1 . Evaluar P(2. 5) en MATLAB
p=[3,-7,2,1,1]; %coefientes del polinomio
xi=2.5;
yi=polyval(p,xi) %evaluacin del polinomio en xi=2.5
yi = 23.812
Las raices del polinomio
p=[3,-7,2,1,1]; %coefientes del polinomio p(x)
r=roots(p)
%soluciones o races de p(x)
r=
r1= 1.8050
r2=1.0000
r3=-0.2358 + 0.3592i
r4=-0.2358 - 0.3592i
Ajuste de polinomios.- Un polinomio de orden N est determinado de forma nica si se dan
N+1 puntos. En otras palabras, el polinomio de orden N ajustados a N+1 puntos de datos,
( xi ; yi ) , i=1,2,3,, N+1, es nico. Los coeficientes del polinomio se pueden determinar
fcilmente con polyfit .
Supongamos que un conjunto de datos est dado por :
x=[1.1,2.3,3.9, 5.1];
y=[3.887,4.276,4.651,2.117];
entonces :a=polyfit(x,y,length(x)-1) % produce coeficientes del polinomio reajustado
produce
a=-0.2015 1.4385 -2.7477 5.4370
que es un vector de los coeficientes del polinomio.
El polinomio aqu determinado es :
y
0.2015 x3 1.4385 x 2 2.7477 x 5.4370 , cuya grafica es :
13
E-2) Dados lo polinomios p( x) 2 x5 5 x 4 4 x3 7 x 2 8 x 2 y q( x) 2 x3 4 x 2

El producto de los polinomios lo calculamos con el comando conv(p,q)
%calcula el producto de p(x).q(x)
p=[2,-5,4,-7,8,-2]
q=[2,-4,3,-5]
prod=conv(p,q)
p = 2 -5 4 -7 8 -2
q = 2 -4 3 -5
prod = 4 -18 34 -55 81 -77 67 -46 10
luego p( x).q( x) 4 x8 18 x 7 34 x6 55 x5 81x 4 77 x3 67 x 2 46 x 10
3x 5
La divisin de los polinomios lo calculamos con el comando [Q,R]=deconv(p,q)

donde Q y R representan los coeficientes del cociente y residuo respectivamente.
p=[2,-5,4,-7,8,-2]
q=[2,-4,3,-5]
[Q,R]=deconv(p,q)
p = 2 -5 4 -7 8 -2
q = 2 -4 3 -5
Q = 1.0000 -0.5000 -0.5000
R= 0
0
0 -2.5000 7.0000 -4.5000
Luego : Q( x) 1.0000 x 2 0.5000 x 0.5000 y R( x)
2.5000 x 2 7.0000 x 4.5000
PRCTICA DE LABORATORIO CON MATLAB
E-1) Suponga que se han dado los siguientes polinomios
f1 ( x)
x3 3x 2 x 3
f 4 ( x)
x3 5 x 2 7 x 3
f 2 ( x)
f 5 ( x)
x3
6 x 2 12 x 8
f 3 ( x)
x 3 8 x 2 20 x 16
x 2
Grafique cada una de las siguientes funciones en el intervalo [0; 4]. Use funciones MATLAB
con vectores de coeficientes de polinomios para evaluar las expresiones:
1)
f1 ( x)
2)
4)
f1 ( x) * f3 ( x)
5)
7)
f 34 ( x) * f 32 ( x)
8)
14
f 2 ( x) 2 f 4 ( x)
f 4 ( x)
x 1
3) 3 f5 ( x) f 2 ( x) 2 f3 ( x)
f1 ( x) * f 2 ( x)
6)
f 5 ( x)
f 54 ( x)
x2 1
8)
f 23 ( x) * f 54 ( x)
f 3 ( x)
E-2) Determine las races reales de los siguientes polinomios. Luego grafique cada polinomio
en un intervalo apropiado a fin de verificar que cruza el eje x en las posiciones de las races
reales.
a) g1 ( x) x3 5 x 2 2 x 8
b) g 2 ( x)
c) g 3 ( x)
x2 4x 4
x2 2x 2
e) g 5 ( x)
x5 3x 4 11x3 27 x 2 10 x 24
x 5 4 x 4 9 x 3 32 x 2 28 x 48
f) g 6 ( x)
x5
3x 4
4 x3
g) g 6 ( x)
x5
9x4
35 x 3
h) g 5 ( x)
2 x9
d) g 4 ( x)
5 x8
4 x7
26 x 2
65 x 2
x6
40 x 24
64 x 26
2 x5
7 x4
2 x3
2x2
8x 4
AJUSTE DE POLINOMIOS
function interpolacion
fprintf('\n')
x = input('ingrese los valores de x=');
y = input('ingrese los valores de y=');
n=length(x);
V=ones(n,n);
x=x';
y=y';
for j=2:n
V(:,j)=x.*V(:,j-1);
end
a=(V\y)'
%GRAFICA
z=x(1):0.2:x(n);
imagen=a(n);
for i = n-1: -1 :1
imagen=a(i)+z.*imagen;
end
plot(z,imagen,x,y,'o')
ylabel('g(x):puntos de datos')
title('INTERPOLACION DE POLINOMIOS')
grid on
E-3) En la tabla se lista el nivel promedio de dixido de carbono en la atmsfera, medido en partes por
milln (ppm) en el observador de Mauna Loa, desde 1972 hasta 1970.
a) Haga un diagrama de dispersin de los datos
b) Determine y grafique el polinomio interpolador
c) Use el polinomio interpolador de b) para estimar el nivel promedio de CO2 para 1987 y predecir el
nivel para los aos 2005 y 2008.
d) De acuerdo al modelo Cunto exceder 400 partes por milln el nivel de CO2?
Ao
Nivel de CO 2 en ppm
1972
327.3
1974
330.0
1976
332.0
1978
335.3
1980
338.5
1982
341.0
1984
344.3
15
1986
1988
1990
347.0
351.3
354.0
Resolucin
a) Haciendo el diagrama de dispersin para los puntos desde 1970 en el eje x
x=[2,4,6,8,10,12,14,16,18,20];
y=[327.3,330,332,335.3,338.5,341,344.3,347,351.3,354];
b)
Polinomio
interpolador
%compilacin
ingrese los valores de x=[2,4,6,8,10,12,14,16,18,20]
ingrese los valores de y=[327.3,330,332,335.3,338.5,341,344.3,347,351.3,354]
%respuesta
a = Columns 1 through 8
402.0000 -110.9756 65.6330 -20.4850
Columns 9 through 10
0.0000 -0.0000
p ( x)
3.7962 -0.4372
0.0316 -0.0014
402 110.9756 x 65.6330 x 2 20.4850 x3 3.7962 x 4 0.4372 x5 0.0316 x6 0.0014 x7 0 x8 0 x9
E-1) Los bilogos han observado que la cantidad de chirridos por minuto de los grillos de una
especie, est relacionada con la temperatura ambiente. La tabla muestra el nmero de
chirridos por minuto para varias temperaturas:
Temperaturas(F) Chirridos por minuto
50
20
55
46
60
79
65
91
70
113
75
140
80
173
85
198
90
211
a) Haga un diagrama de dispersin de los datos
b) Determine y grafique el polinomio interpolante
c) Use el polinomio para estimar el nmero de chirridos a los 100F
E-2) (Contaminacin del agua) Al depositarse en un lago, los desperdicios orgnicos disminuyen el
contenido de oxgeno del agua. Si t denota el tiempo en das despus que se deposita el desperdicio, se
encuentra experimentalmente en un caso que el contenido de oxgeno es y
t 3 30 t 2
6000 con
0 t 25 . Encuentre los valores mximo y mnimo de y durante los 25 das siguientes al vaciado del
desperdicio.
16
E-3) (Contaminacin atmosfrica) El ndice de contaminacin atmosfrica en cierta ciudad vara
durante el da de la siguiente manera: P (t )
2 4t
; si 0
6 2t
; si 2
14
; si 4
12
50 3 t ; si 12
16
Aqu t es el tiempo en horas, con t=0 correspondiente a 6 a.m. y t=16 a 10 p.m. Haga la grfica de esta
funcin. Cul son los niveles de contaminacin a las 8 a.m. , 12 del da, 6 p.m. y 8 p.m.?
E-4) Al depositarse en un lago, los desperdicios orgnicos disminuyen el contenido de oxgeno del agua.
Si t denota el tiempo en das despus que se deposita el desperdicio, se encuentra experimentalmente
3
en un caso que el contenido de oxgeno es y t

30 t 6000 con 0 t 25 . Encuentre los
valores mximo y mnimo de y durante los 25 das siguientes al vaciado del desperdicio.
SESIN 05- CLCULOS MATEMTICOS EN MATLAB

1) Clculos en una sola variable
end
2
%resultado
El rea de la esfera A 4 r
A = 50.2655
r=2;
El enunciado if se puede utilizar con
A=4*pi*r.^2
else o elseif
%resultado
r=2;
A = 50.2655
if r>3
b=1
Otra forma es utilizando el comando
elseif r==3
subs
b=2
syms r
else
A=4*pi*r.^2;
b=0
subs(A,2)
end
%resultado
%resultado b =0
ans = 50.2655
El enunciado if
r=2;
if r>0
A=4*pi*r.^2
end
%resultado
A = 50.2655
El operador igual(==)
r=2;
if r==2
A=4*pi*r.^2
end
%resultado
A = 50.2655
El operador no igual (~=)
r=2;
if r~=3
A=4*pi*r.^2
17
r=6;
if r>3
b=1
elseif r==3
b=2
else
b=0
end
%resultado
b=1
r=3;
if r>3
b=1
elseif r==3
b=2
else
b=0
end
%resultado
b =2
El comando disp exhibe un nmero,

vector, matriz o cadena en la ventana
de comandos
disp(pi)
%resultado
3.1416
disp('YO SOY UN EXCELENTE ALUMNO')
%resultado
YO SOY UN EXCELENTE ALUMNO
Los ciclos for/end y while/end
El rea de la esfera para r=1 hasta 5
for r=1:5
Area=4*pi*r.^2;
disp([r,Area])
end
%resultado
ri
Areas
1.0000
12.5664
2.0000
50.2655
3.0000
113.0973
4.0000
201.0619
5.0000
314.1593
r=0;
while r<5
r=r+1;
Area=4*pi*r.^2;
disp([r,Area])
end
%resultado
ri
1.0000
2.0000
3.0000
4.0000
5.0000
Areas
12.5664
50.2655
113.0973
201.0619
314.1593
El ndide de ciclo puede decrementarse

for r=5:-1:1
Area=4*pi*r.^2;
disp([r,Area])
end
%resultado
ri
Areas
5.0000
314.1593
4.0000
201.0619
3.0000
113.0973
2.0000
50.2655
1.0000
12.5664
Ciclos dobles y triples
for r=1:5
for s=1:r
Area=4*pi*r.^2;
18
disp([r,Area])
end
end
%resultados
ri
1.0000
2.0000
2.0000
3.0000
3.0000
3.0000
4.0000
4.0000
4.0000
4.0000
5.0000
5.0000
5.0000
5.0000
5.0000
Areas
12.5664
50.2655
50.2655
113.0973
113.0973
113.0973
201.0619
201.0619
201.0619
201.0619
314.1593
314.1593
314.1593
314.1593
314.1593
El comando break
for r=1:3
for j=1:10
Area=4*pi*r.^2;
disp([r,Area])
if j>2*r
break
end
end
end
%resultados
ri
Areas
1.0000
12.5664
1.0000
12.5664
1.0000
12.5664
2.0000
50.2655
2.0000
50.2655
2.0000
50.2655
2.0000
50.2655
2.0000
50.2655
3.0000
113.0973
3.0000
113.0973
3.0000
113.0973
3.0000
113.0973
3.0000
113.0973
3.0000
113.0973
3.0000
113.0973
VARIABLES Y ARREGLOS
x=[4,6,8,1,9];
x(2)%
x(5)%
%resultado
ans =6
ans =9
ADEMS
for i=1:6
x(i)=i*(i+1)/2
end
%RESULTADOS
x = 1 3 6 10
x = 1 3 6 10
x = 1 3 6 10
x = 1 3 6 10
x = 1 3 6 10
x = 1 3 6 10
15
15
15
15
15
15
21
21
21
21
21
21
Decrementa de dos en dos

x=14:-2:2
%resultado
x =14 12 10 8 6 4 2
Incrementa de dos en dos
x=2:2:12
%resultado
x = 2 4 6 8 10 12
Arreglos
x=[3,5,7,8];
y=[1,2,9,6];
for i=1:4
z(i)=x(i)+y(i)
end
%resultados
z=4
z=4 7
z = 4 7 16
z = 4 7 16 14
Posiciones de areglos
A=[4,6,7;9,2,5;-3,-5,9]
A(1,1)
A(2,3)
A(1,1)
%resultados
A =4 6 7
9 2 5
-3 -5 9
ans = 4
ans = 5
ans = 4
Posiciones
A=[4,6,7;9,2,5;-3,-5,9]
A(3,:)
%resultados
A =4 6 7
9 2 5
-3 -5 9
ans =-3 -5 9
19
A=[4,6,7;9,2,5;-3,-5,9]
A(:,3)
%resultados
A =4 6 7
9 2 5
-3 -5 9
ans = 7
5
9
Anidados
a=[3,4,6;3,8,3;4,3,6];
b=[1,2,3;7,8,9;4,6,6];
for i=1:3
for j=1:3
c(i,j)=a(i,j)+b(i,j)
end
end
%respuesta
c=4
c=4 6
c=4 6 9
c=4 6 9
10 0 0
c=4 6 9
10 16 0
c=4 6 9
10 16 12
c=4 6 9
10 16 12
8 0 0
c=4 6 9
10 16 12
8 9 0
c=4 6 9
10 16 12
8 9 12
sumatorias
x=[3,4,6,7,9];
sum(x)
%respuesta
ans = 29
x=[3,4,6;5,7,9];
sum(x)
%respuesta
ans = 8 11 15
x=[3,4,6;5,7,9;1,2,3];
sum(x)
%respuesta
ans = 9 13 18
x=[3,4,6,5,7,9];
MIN=min(x)
MAX=max(x)
%respuestas
MIN = 3
MAX = 9
3) La densidad D de un objeto cilndrico es
superior de la densidad del objeto donde:
Media aritmtica
function [media, desviacion]=medias(x)
n=length(x);
media=sum(x)/n;
desviacion=sqrt(sum(x.^2)/n-(media).^2);
%Compilando
x=[2,3,4,5,6,7];
[m,d]=medias(x)
%respuesta
m = 4.5000
d = 1.7078
2) El campo elctrico de un disco circular

con radio R y densidad de carga
con
distancia
Z
a
un
eje
es
3 2
20
15,431 0,012 cm
1
R
z2
12
0
c
10 5 m ;
2 10
16
3 109 =velocidad de la luz . Hallar su
velocidad.
4 3
r con r=2
3
1,2 10
4) Para un mesn 0 que decae en dos

fotones con una vida media
se quiere
calcular
su velocidad por la frmula
%compilando
fun('soria',2,4,-1)
%resultados
ans = 7.5000
LOBORATORIO DE MATLAB
1) Calcular el volumen de una esfera
para
8,20 0,01 cm
donde; L
function mp=fun(soria,a,b,c)
mp=((feval(soria,a)+2*feval(soria,b)+feval(soria,
c))/4);
.z 1
2 0 z
%Compilando
soria([3, 2;4,1])
%respuesta
ans = 5 3
9 3
Funciones de Usuario propias
EZ
0.029 0.005 gr
Creacin de archivos M
function y=soria(x)
y=x.^2-3*x+5;
%Compilando
soria(3)
%respuesta
ans = 5
M
, hallar la cota inferior y
R2 L
. Calcular EZ
8,85 10
12
5 4 3 e e=N neperiano
4)
(Contaminacin
del
agua)
Al
depositarse en un lago, los desperdicios
orgnicos disminuyen el contenido de
oxgeno del agua. Si t denota el tiempo en
das despus que se deposita el
desperdicio,
se
encuentra
experimentalmente en un caso que el
contenido
de
oxgeno
es
3
2
y t 30 t 6000
con 0 t 25 .
Encuentre los valores mximo y mnimo de
y durante los 25 das siguientes al vaciado
del desperdicio.
SESIN 06 - GRFICOS CON MARCAS Y PROGRAMACIN CON EL SOFTWARE MATLAB 7.0
GRFICOS CON MARCAS

TIPO DE MARCA
Punto
Ms
Estrella
Crculo
Marca x
TIPOS DE LNEA
TIPO DE LNEA
Continua
Guiones
Punteada
Guiones y puntos
SMBOLO
.
+
*
o
x
SMBOLO
-:
-.
COLORES DE LNEA
COLORES DE LNEA
SMBOLO
Rojo
r
Amarillo
y
Magenta
m
Turquesa
c
Verde
g
Azul
b
Blanco
w
Negro
k
E-1) Graficar la funcin f ( x )
ylabel('EJE Y')
text(4,5,'f1=10*sin(x)./(x)')
10 Sen( x )
x
Resolucin
%Formato para graficar
x=-15:0.05:15;
f1=10*sin(x)./(x);
plot(x,f1,'+')
21
grid on
title('G
RAFIC
AS
CON
MARC
AS')
xlabel(
'EJE
X')
El comando FontSize
%Este formato aumenta de tamao las

letras
x=-15:0.05:15;
f=sin(x)./(x);
plot(x,f,'*m')
grid on
title('GRAFICAS CON
MARCAS','FontSize',[20],'color','g')
xlabel('ABSCISAS','FontSize',[18],'color','r')
ylabel('ORDENADAS','FontSize',[18],'color','
r')
gtext('f=sin(x)./(x)','FontSize',[15],'color','b')
El
comando
superposicin
x=-10:0.05:10;
f1=sin(x)./(x+eps);
plot(x,f1,'*')
hold on
f2=sin(x);
f3=-cos(x);
plot(x,f2,'+g')
plot(x,f3,'.-r')
grid on
hold
on
en
la
title('GRAFICAS
CON
MARCAS','FontSize',[20],'color','b')
xlabel('ABSCISAS','FontSize',[18],'color','c')
ylabel('ORDENADAS','FontSize',[18],'color','c')
Ejemplo:
E-1)Graficar las funciones
subplot(2,2,1)
x=-5:0.05:5;
f1=cos(x);
plot(x,f1,'b')
grid on
title('subplot(2,2,1)')
subplot(2,2,2)
x=-5:0.05:5;
f2=sin(x);
plot(x,f2,'g')
grid on
subplot(2,2,3)
x=-2:0.05:2;
f3=-x.^2;
plot(x,f3,'r')
grid on
subplot(2,2,4)
x=-5:0.05:5;
f4=-x.*sin(x);
plot(x,f4,'k')
grid on
El comando subplot en grficos del

plano
Con subplot podemos graficar m por n
grficas en una sola figura; la sintaxis es:
subplot(m,n,k) donde k es el nmero
secuencial de la grfica.
%Este formato realiza 4 graficas de 2 filas y
2 %columnas
subplot(2,2,1)
gtext('figura1')
subplot(2,2,2)
gtext('figura2')
subplot(2,2,3)
gtext('figura3')
subplot(2,2,4)
gtext('figura4')
Prctica dirigida con MATLAB
1.2.3.4.5.6.7.-
Graficar cada una de las funciones reales
f ( x)
x2
f ( x)
f ( x)
x cos( x) 0
sen ( x) 2 cos( x)
f ( x)
f ( x)
f ( x)
f ( x)
8.- f ( x)
22
ex
cos( x) (1 x )
( x 2)
x
ln( x)
0
0
4 xsen( x) (2senx)2
log( x) 0.2 x
2
0
1
2 cos x 6
9. f ( x)
10. f ( x)
senx
5) f ( x)
x3 5
ex
cos x
Lnx
Utilice el comando subplot para graficar
sen 3 x
1. f ( x)
2 f ( x)
sen (tan( x 2 3))
x3 5
5
2 x 3 log( senx )
senx cos x x 3 2 x 6
3. f ( x)
4. f ( x) cos(tan( .x))
TAREA DOMICILIARIA
Graficar las funciones reales
1) f ( x) 10 2 x log x 1 x2
3) f ( x) Ln( x) 2 x 4
4)
( x)
x3
e senx
senx
cos x
log x
x3
cos x
7) f ( x) e
ln x x3 4
8) ( x) x 2 4 xsen( x) (2senx )2
9) f ( x)
x 1 x sen( x)
10) g ( x) x 2 4 x cos( x) (2 cos x)2
12) r ( x) 2 x e x 2 cos x 6
2x 1
SESIN 07 - ALGORITMOS Y PROGRAMACIN CON EL SOFTWARE MATLAB 7.0

I. Escriba un M-archivo de funcin que calcule la solucin de ax 2 bx c 0
Resolucin
Sabemos que la solucin la determina la formula general dada por
x1, 2
b2
2a
4ac
(1)
Esto tendramos que ingresarlo en MATLAB

PRIMERA FORMA CON MATLAB:
%programa soria1 para hallar la solucin de (1)
function raiz
a=input('INGRESE EL VALOR DE a=');
b=input('INGRESE EL VALOR DE b=');
c=input('INGRESE EL VALOR DE c=');
fprintf ('\n');
fprintf ('
raiz1
raiz2 \n');
r1=(-b+sqrt(b.^2-4*a*c))./(2*a);
r2=(-b-sqrt(b.^2-4*a*c))./(2*a);
disp([r1,r2])
E-1) Resolver la ecuacin 2 x 2 11x 15 0
Resolucin
%Compilacin
raiz
INGRESE EL VALOR DE a=2
INGRESE EL VALOR DE b=-11
INGRESE EL VALOR DE c=15
%Respuestas
raiz1
raiz2
3.0000 2.5000
E-2) Resolver la ecuacin x 2 x 1 0
Resolucin
%Compilacin
raiz
INGRESE EL VALOR DE b=1
%Respuestas
raiz1
raiz2
-0.5000 + 0.8660i
-0.5000 - 0.8660i
23
E-3) Resolver la ecuacin 9 x 2
12 x 4
0
Resolucin
%Compilacin
raiz
INGRESE EL VALOR DE b=-12
raiz1
raiz2
0.6667
0.6667
SEGUNDA FORMA CON MATLAB:
%programa soria2 para hallar la solucin de (1)
function [r1,r2]=solucion(a,b,c)
D=b.^2-4*a*c;
r1=(-b+sqrt(D))./(2*a);
r2=(-b-sqrt(D))./(2*a);
11x 15
0
Resolucin
%Compilacin
[r1,r2]=solucion(2,-11,15)
%Respuesta
r1 = 3
r2 = 2.5000
E-2) Resolver la ecuacin x 2
x 1
0
Resolucin
%Compilacin
[r1,r2]=solucion(1,1,1)
%Respuestas
r1 = -0.5000 + 0.8660i
r2 = -0.5000 - 0.8660i
12 x 4
0
Resolucin
%Compilacin
[r1,r2]=solucion(9,-12,4)
%Respuesta
r1 = 0.6667
r2 = 0.6667
EL COMANDO fplot del MATLAB 7.0
El comando fplot traza una funcin entre los lmites especificados. La funcin debe ser
y=f(x) donde x es un vector cuyo rango especifica los lmites y y es un vector del mismo
tamao como x y contiene el valor de la funcin a los puntos en x.
EJEMPLO1
Graficar la funcin f ( x)
300 Sen( x)
, en un dominio de [-15; 15]
x
Resolucin
fplot('300*sin(x)./x',[-15 15])
grid on
title('EL COMANDO fplot')
24
gtext('f(x)=300*sin(x)./x')
Cuya grfica es:
En el comando fplot si la funcin f(x) devuelve ms de un valor por un x dado, entonces y

es una matriz cuyas columnas contienen cada componente de f(x). Aqu se debe crear un
M-archivo especificando las funciones que se quieren graficar.
EJEMPLO2
Graficar la funcin f ( x)
2 e x .Cos(kx) , para cunado k varie de 1,2,3,4,5,6,7. con D(f)=[-5;1]

Resolucin
%Crear un M-archivo y guardarlo con grafico7.m
function y=grafico7(x)
y(:,1) = 2*exp(-x(:)).*cos(x(:));
y(:,2) = 2*exp(-x(:)).*cos(2*x(:));
y(:,3) = 2*exp(-x(:)).*cos(3*x(:));
y(:,4) = 2*exp(-x(:)).*cos(4*x(:));
y(:,5) = 2*exp(-x(:)).*cos(5*x(:));
y(:,6) = 2*exp(-x(:)).*cos(6*x(:));
y(:,7) = 2*exp(-x(:)).*cos(7*x(:));
%Para compilar los grficos crear un M-archivo y guardarlo con compfplot7 .m
function compfplot7
fh=@grafico7;
fplot(fh,[-5 1])
grid on
xlabel('EJE X')
ylabel('EJE Y')
title('EL COMANDO fplot DE MATLAB')
gtext('f1=2*exp(-x)*cos(x)')
gtext('f2=2*exp(-x)*cos(2*x)')
EJEMPLO3
Graficar la funcin para k=0,1,2,3,4,5,6. con dominio [-2 ; 1] siendo
f ( x)
50Sen(kx) 50Cos(( k 1) x)
Resolucin
%Crear un M-archivo y guardarlo con

grafico8.m
function y=grafico8(x)
y(:,1) = 50*sin(pi*x(:));
y(:,2) = 50*sin(2*pi*x(:))+50*cos(3*pi*x(:));
y(:,3) = 50*sin(4*pi*x(:))+50*cos(5*pi*x(:));
y(:,4) = 50*sin(6*pi*x(:))+50*cos(7*pi*x(:));
y(:,5) = 50*sin(8*pi*x(:))+50*cos(9*pi*x(:));
y(:,6) = 50*sin(10*pi*x(:))+50*cos(11*pi*x(:));
25
y(:,7) = 50*sin(12*pi*x(:))+50*cos(13*pi*x(:));
%Para compilar los grficos crear un M-archivo y guardarlo con compfplot8 .m
function compfplot8
fh=@grafico8;
fplot(fh,[-2 1])
grid on
xlabel('EJE X')
ylabel('EJE Y')
title('EL COMANDO fplot DE MATLAB')
Graficas de superficies
E-1) Graficar la superficie
f ( x, y)
Sen( xy) x 2
2Cos( y)
Resolucin
%Este formato grafica una superficie

x1=-1:0.05:1;
y1=-1:0.05:1;
[x,y]=meshgrid(x1,y1);
f1=sin(x.*y)-x.^2+2*cos(y)
mesh(x,y,f1)
grid on
title('f(x;y)=sin(x.*y)-x.^2+2*cos(y)')
xlabel('EJE X')
ylabel('EJE Y')
zlabel('EJE Z')
E-2) Graficar la superficie
f ( x, y)
x2
x. y 3 3x 2 . y
y3 5
Resolucin
x1=-1:0.05:1;
y1=-1:0.05:1;
f1=x.^2+x.*(y.^3)-3*(x.^2).*y+y.^3-5;
mesh(x,y,f1)
grid on
title('x.^2+x.*(y.^3)-3*(x.^2).*y+y.^3-5')
xlabel('EJE X')
ylabel('EJE Y')
zlabel('EJE Z')
%GRAFICA DE CONTORNO3
x1=-6:0.05:6;
y1=-6:0.05:6;
f1=x.^2+x.*(y.^3)-9;
contour(x,y,f1,[0.00 0.00],'k')
26
hold on
f2=3*(x.^2).*y-y.^3-4;
contour(x,y,f2,[0.00 0.00],'b')
grid on
title('Contorno de f1(x,y) y f2(x,y) ')
gtext('x.*exp(x.*y+0.8)+exp(y.^2)-3')
gtext('x.^2-y.^2-0.5*exp(x.*y)')
PRCTICA DIRIGIDA DE MATLAB
f ( x)
1. Cree un archivo M de funcin, fun_ex(x) que calcule la funcin

Pruebe su funcin tecleando fun_es(3) y fun_es([1 2 3]).
2.Escriba un M-archivo de funcin que calcule la solucin de ax 2

solucin la determine la formula general dada por
x1, 2
0.5 e
bx c
x
3
x 2 Senx .
0 , donde la
2c
b2
4ac
3. Utilice el programa anterior y resuelva:
x2
x2
x2
x2
1000.001x 1 0
10000.0001x 1 0
100000.00001x 1 0
1000000.000001x 1 0
4. El impulso de una partcula en un instante t0 es P(t )

Crear un M-archivo que al ingresar: k; m;
0
5. Para un mesn
k2
t.e
2m
t2
, en un tiempo determinado nos d el impulso P(t).
que decae en dos fotones con una vida media
velocidad por la frmula L
t2
(1 e
se quiere calcular su
Despeje la variable v y cree un archivo-M para calcular la velocidad ingresando valores para
2 10 16
L, v, c y . Pruebe con L 10 5 m ;
3 109 =velocidad de la luz.
GRAFICAR LOS SIGUIENTES CONTORNOS

1)
f1 ( x, y )
f 2 ( x, y )
3)
f1 ( x, y )
f 2 ( x, y )
27
x2
e
y2
xy x 2
xe
( x 4)
x
xe3 x
y2
x3
2)
( y 4)
4)
16
f1 ( x, y )
xexy
0.8
f 2 ( x, y )
f1 ( x, y )
x2
f 2 ( x, y )
3 cos( x)
ey
0.5e
y 1
y
3
xy
5)
7)
f1 ( x, y )
cos( x)
f 2 ( x, y )
x sen( y )
f1 ( x, y )
xe y 1
f 2 ( x, y )
x2
4 y2
6)
8)
f1 ( x, y )
x e
f 2 ( x, y )
f1 ( x, y )
x 2 sen ( y )
y2 1
f 2 ( x, y )
x 2 cos( y )
y2 1
PRCTICA DIRIGIDA
1
.e
2
1. La distribucin Normal lo determina la funcin real f ( x )
kx 2
2
Utilice el comando fplot para generar grficos donde k vare de 1,2,3,4,5,6,7,8 y 9.

2. La funcin de distribucin de Poisson es f ( x)
Genere un archivo M que al ingresar
e .
x!
;x
0,1,2,3
y x obtenga probabilidad de Poisson f(x).
3. La probabilidad de que un operador de telfonos reciba exactamente x llamadas durante

cierto periodo esta dado por f ( x)
e 3 . 3x
x!
Encuentre la probabilidad de que el operador reciba exactamente tres llamadas.

4. En estudios de redes elctricas, aparece la ecuacin S 2
R
.S
L
1
LC
.(1)
Escriba un M-archivo de funcin que calcule la solucin de (1), donde la solucin la determine
la formula general dada por
R
2L
S1, 2
R
2L
1
LC
5. Un termmetro con resistencia de platino, de ciertas especificaciones, opera de acuerdo con

la ecuacin.
10 000 (4.124 10 2 ) T
(1.779 10 5 ) T 2
Donde R es la resistencia (en ohms) del termmetro a la temperatura T (en grados celsius). Si
R=13.946, determine el valor de correspondiente de T.
6. Bajo ciertas condiciones, si dos padres con ojos de color caf tienen exactamente tres hijos,
la probabilidad P de que tengan exactamente r hijos con ojos azules est dada por la funcin
r
P(r )
3 r
1
3
3!
4
4
r !. (3 r ) !
;r
0,1,2,3.
Determine la probabilidad de que exactamente dos de los hijos tengan los ojos azules.
28
SESIN 08- PLOTEOS TRIDIMENSIONALES

En el cuadro siguiente se presentan los comandos de MATLAB ms comunes en la
representacin de grficos de lneas 3-D.
plot3(X,Y,Z)
Dibuja el conjunto de puntos (X,Y,Z), donde X,Y y Z son vectores fila. X,

Y, Z pueden ser coordenadas paramtricas o matrices de la misma
dimensin, en cuyo caso se hace una grfica por cada tripleta de filas y
sobre los mismos ejes. Para valores complejos de X, Y y Z se ignoran
las partes imaginarias.
plot3(X,Y,Z,S)
Grfica de plot(X,Y,Z) con las opciones definidas en S. Usualmente S se

compone de dos dgitos entre comillas simples, el primero de los cuales
fija al color de la lnea del grafico y el segundo el carcter a usar en el
graficado.
Los
valores
posibles
de
colores
caracteres
son,
respectivamente, los siguientes: y(amarillo), m(magenta), c(cyan), r(rojo),

g(verde), b(azul), w(blanco), k(negro), .(puntos), o(crculos), x(x-marcas)
, +(signo ms), -(slido), *(estrellas), :(dos puntos), -.(guiones y puntos) y
- (semislido).
plot3(X1,Y1,Z1,S1,
Combina, sobre los mismos ejes, los grficos definidos para las tripletas
X2,Y2,Z2,S2,
(Xi,Yi,Zi,Si). Se trata de una forma de representar varias funciones sobre
X3,Y3,Z3,S3, )
el mismo grfico.
fill3(X,Y,Z,C)
Dibuja el polgono compacto cuyos vrtices son las tripletas de

componentes (Xi,Yi,Zi) de los vectores columna X, Y y Z, que contiene
los colores Ci de cada punto (Xi,Yi,Zi). Los valores de Ci pueden ser y,
m, c, r, g, b, w, k. cuyos significados ya conocemos. Si C es un solo
caracter, se pintarn todos los puntos del polgono de color
correspondiente al carater.
fill3(X1,Y1,Z1,C1,
Dibuja el polgono compacto cuyos vrtices vienen dados por los puntos
X2,Y2,Z2,C2, )
(Xi,Yi,Zi,Ci)
29
Transforma el campo de definicin dado de las variables x e y de la

[X,Y]=meshgrid(x,y)
funcin a representar z=f(x,y) en argumentos matriciales utilizables por

los comandos surf y mesh para obtener grficos de superficie y malla,
respectivamente
surf(X,Y,Z,C)
Representa la superficie explcita z=f(x,y) o la paramtrica x=x(t,u),

y=y(t,u), z=z(t,u), realizando el dibujo con los colores especficos en C. El
argumento C se puede ignorar.
surfc(X,Y,Z,C)
Representa la superficie explcita z=f(x,y) o la paramtrica x=x(t,u),

y=y(t,u), z=z(t,u), junto con el grfico de contorno correspondiente(curvas
de nivel proyectadas sobre el plano XY).
surfl(X,Y,Z)
Representa la superficie explcita z=f(x;y) o la paramtrica x=x(t,u),

y=y(t,u), z=z(t,u), con el dibujo con sombreado.
mesh(X,Y,Z,C)
Representa la superficie explicita z=f(x,y) o la paramtrica x=x(t,u),

y=y(t,u), z=z(t,u), dibujando las lneas de la rejilla que componen la malla
con los colores especificadas en C(opcional).
meshz(X,Y,Z,C)

y=y(t,u), z=z(t,u), con una especie de cortina o teln en la parte inferior
meshc(X,Y,Z,C)

y=y(t,u), z=z(t,u), junto con el grfico de contorno correspondiente(curva
de nivel proyectadas sobre el plano XY )
contour(Z)
Dibuja el grfico de contorno (curvas de nivel) para la matriz Z. El

nmero de lneas de contorno a utilizar se elige automticamente.
contour(Z,n)
Dibuja el grfico de contorno (curvas de nivel) para la matriz Z usando n

lneas de contorno.
Dibuja el grfico de contorno (curvas de nivel) para la matriz Z usando en los
contour(x,y,Z,n)
ejes X e Y el escalado definido por los vectores x e y (n lneas de contorno).
contour3(Z)
contour3(Z,n)
Dibujan los grficos de contorno en 3 dimensiones.
contour3(x,y,Z,n)
contourf()
Dibuja un grfico de contorno y rellena las reas entre las isolneas.
pcolor(X,Y,Z)
Dibuja un grfico de contorno (curvas de nivel ) para la matriz (X,Y,Z)

utilizando una representacin basada en densidades de colores. Suele
denominarse grfico de densidad
comet3(z)
Grfico cometa relativo al vector z
comet3(x,y,z)
Grfico de cometa paramtrico x(t), y(t) , z(t)
comet3(x,y,z,p)
Grfico de cometa con cuerpo de longitud p*length(y)
[X,Y,Z]=cylinder
Da las coordenadas del cilindro unidad
30
[X,Y,Z]=cylinder(r)
Da las coordenadas del cilindro generado por la curva r
[X,Y,Z]=cylinder(r,n)
Da las coordenadas del cilindro generado por la curva r con n puntos en

la circunferencia seccin horizontal alineado con el eje Z (n=20 por
cylinder()
defecto)
Grafica los cilindros anteriores
sphere
Grafica la esfera unidad usando 20x20 caras
sphere(n)
Grafica la esfera unidad usando nxn caras
[X,Y,Z]=sphere(n)
Da las coordenadas de la esfera en tres matrices (n+1)x(n+1)
OPCIONES DE MANEJO DE GRFICOS 3D

colormap(M)
Sita la matriz M como el mapa corriente de colores. M debe tener

tres columnas y contener valores slo entre 0 y 1. Tambin puede
ser una matriz cuyas filas sean vectores RGB del tipo [r g b]. Existen
en MATLAB matrices M ya definidas, que son las siguientes:
bone(p), contrast(p), cool(p), copper(p), flag(p), gray(p), hsv(p),
hot(p), jet(p), pink(p), prism(p) y white(p). Todas las matrices tienen
3 columnas y p filas. Por ejemplo, la sintaxis colormap(hot(8)) sita
la matriz hot(8) como el mapa corriente de colores(sistema completo
de colores de la figura actual).
E-1) Grfica de una funcin paramtrica:
x(t )
sen(t )
y (t )
cos(t )
z (t )
Resolucin
%PLOTEO1 - Comando plot3
t=0:pi/50:10*pi;
x=sin(t);
y=cos(t);
z=t;
plot3(x,y,z,'k');
grid on
axis square
title('Hlice paramtrica x(t)=sin(t),y(t)=cos(t),z(t)=t');
xlabel('EJE DE ABSCISAS X');
ylabel('EJE DE ORDENADAS Y');
zlabel('EJE Z');
31
E-2) Grfica de la funcin
f ( x, y)
sen( x 2
y 2 ) cos 2 ( xy)
Resolucin
%PLOTEO2 - Comando mesh

xa=-1/2*pi:.1:1/2*pi;
xb=-1/2*pi:.1:1/2*pi;
[x,y] = meshgrid(xa,xb);
z = sin(x.^2+y.^2)-(cos(x.*y)).^2;
mesh(x,y,z)
colormap(cool)
grid on
title('z = sin(x.^2+y.^2)-(cos(x.*y)).^2');
xlabel('EJE X');
ylabel('EJE Y');
zlabel('EJE Z');
E-3) Grfica de la funcin f ( x, y)
sen( x 2
y 2 ) cos 2 ( xy)
Resolucin
%PLOTEO3 - Comando surf
xa=-1/2*pi:.1:1/2*pi;
xb=-1/2*pi:.1:1/2*pi;
[x,y]=meshgrid(xa,xb);
z=sin(x.^2+y.^2)-(cos(x.*y)).^2;
surf(x,y,z)
colormap(hsv)
grid on
title('z = sin(x.^2+y.^2)-(cos(x.*y)).^2');
32
xlabel('EJE X');
ylabel('EJE Y');
zlabel('EJE Z');
sen( x 2
y 2 ) cos 2 ( xy)
Resolucin
%PLOTEO4 - Comando surfc
xa=-1/2*pi:.1:1/2*pi;
xb=-1/2*pi:.1:1/2*pi;
[x,y]=meshgrid(xa,xb);
z=sin(x.^2+y.^2)-(cos(x.*y)).^2;
surfc(x,y,z)
colormap(gray)
grid on
title('z=sin(x.^2+y.^2)-(cos(x.*y)).^2');
xlabel('EJE X');
ylabel('EJE Y');
zlabel('EJE Z');
E-5) Grfica de la funcin
f ( x, y)
sen( x 2
y 2 ) cos 2 ( xy)
Resolucin
%PLOTEO5 - Comando contour

xa=-1/2*pi:.1:1/2*pi;
ya=-1/2*pi:.1:1/2*pi;
[x,y]=meshgrid(xa,ya);
z=sin(x.^2+y.^2)-(cos(x.*y)).^2;
contour(x,y,z,40)
colormap(cool)
33
grid on
title(' z=sin(x.^2+y.^2)-(cos(x.*y)).^2');
xlabel('EJE X');
ylabel('EJE Y');
zlabel('EJE Z');
sen( x 2
y 2 ) cos 2 ( xy)
Resolucin
%PLOTEO6 - Comando contour3
xa=-1/2*pi:.1:1/2*pi;
ya=-1/2*pi:.1:1/2*pi;
[x,y] = meshgrid(xa,ya);
z=sin(x.^2+y.^2)-(cos(x.*y)).^2;
contour3(x,y,z,40)
colormap(cool)
grid on
xlabel('EJE X');
ylabel('EJE Y');
zlabel('EJE Z');
sen( x 2
y 2 ) cos 2 ( xy)
Resolucin
%PLOTEO7 - Comando waterfall
xa=-1/2*pi:.1:1/2*pi;
ya=-1/2*pi:.1:1/2*pi;
z=sin(x.^2+y.^2)-(cos(x.*y)).^2
waterfall(x,y,z)
colormap(cool)
grid on
34
xlabel('EJE X');
ylabel('EJE Y');
zlabel('EJE Z');
sen( x 2
y 2 ) cos 2 ( xy)
Resolucin
%PLOTEO8 - Comando meshz
xa=-1/2*pi:.1:1/2*pi;
ya=-1/2*pi:.1:1/2*pi;
z=sin(x.^2+y.^2)-(cos(x.*y)).^2;
meshz(x,y,z)
colormap(cool)
grid on
xlabel('EJE X');
ylabel('EJE Y');
zlabel('EJE Z');
E-9) Grfica de la funcin f ( x, y )
sen x 2 y 2
x2 y 2
Resolucin
%PLOTEO9 - Comando surfl

xa=-7.5:.1:7.5;
ya=-7.5:.1:7.5;
z=sin(sqrt(x.^2+y.^2))./(sqrt(x.^2+y.^2));
surfl(x,y,z)
colormap(hot)
35
grid on
title('z=sin(sqrt(x.^2+y.^2))./(sqrt(x.^2+y.^2))');
xlabel('EJE X');
ylabel('EJE Y');
zlabel('EJE Z');
E-10) Grfica del contorno de las funciones f1 y f2, para resolver
f 1( x, y )
x.e( xy
0.8 )
f 2( x, y )
ey
0.5e
3
xy
Resolucin
%crear una ventana y guardarlo como f1.m
function f=f1(x,y)
f=x.*exp(x.*y+0.8)+exp(y.^2)-3;
------------------------------------------------------------------------------%crear una ventana y guardarlo como f2.m
function f=f2(x,y)
f=x.^2-y.^2-0.5*exp(x.*y);
------------------------------------------------------------------------------%crear una ventana y guardarlo como contorno1.m
%Grfica del contorno de f1 y f2
clear, clg ,clf, hold off
x1=-3:0.01:4;
y1=-3:0.01:4;
g1=f1(x,y);
g2=f2(x,y);
contour(x1,y1,g1,[0.00, 0.00],'k')
hold on
contour(x1,y1,g2,[0.00, 0.00],'b')
grid on
zoom on
title('f1=x*exp(x*y+0.8)+exp(y^2)-3;f2=x^2-y^2-0.5*exp(x*y)')
36
f 1( x, y )
f 2( x, y )
. Dados
E-11) Grfica del contorno de las funciones f1 y f2, para resolver
f 1( x, y )
sen( xy)
f 2( x, y )
x2
2 cos( y )
Resolucin
%crear una ventana y guardarlo como f1.m

function f=f1(x,y)
f=sin(x.*y);
-----------------------------------------------------------------------------%crear una ventana y guardarlo como f2.m
function f=f2(x,y)
f=x.^2+2*cos(y);
-----------------------------------------------------------------------------%crear una ventana y guardarlo como contorno2.m
clear, clg ,clf, hold off
x1=-2*pi:0.01:2*pi;
y1=-2*pi:0.01:2*pi;
g1=f1(x,y);
g2=f2(x,y);
contour(x1,y1,g1,[0.00, 0.00],'k')
hold on
contour(x1,y1,g2,[0.00, 0.00],'b')
grid on
zoom on
title('f1=sin(x.*y);f2=x.^2+2*cos(y)')
37
f 1( x, y )
f 2( x, y )
. Dados
E-12) Interceptar las superficies
x2
y2
y2
Resolucin
[x y]=meshgrid(-2:.1:2)
z=x.^2+y.^2;
mesh(x,y,z);
hold on;
z=4-y.^2;
mesh(x,y,z);
E-13) Interceptar las superficies
x2
sen( x)
y2
Resolucin
>>[x y]=meshgrid(-2:.1:2);
z=x.^2;
mesh(x,y,z);
hold on;
z=4-y.^2;
mesh(x,y,z);
hold on;
38
z=sin(x);
mesh(x,y,z);
xlabel('EJE X')
ylabel('EJE Y')
zlabel('EJE Z')
title('Interseccin de superficies')
E-14) Graficar :
x2
z
z
3(1 x) 2 e
y2
2
z=peaks
sen ( x)
x 2 ( y 1) 2
10( 5x
x3
y 5 )e
x2 y2
Resolucin
>> subplot(2,2,1)
[x y]=meshgrid(-2:.1:2);
z=x.^2;
mesh(x,y,z);
hold on;
z=4-y.^2;
mesh(x,y,z);
xlabel('EJE X')
ylabel('EJE Y')
zlabel('EJE Z')
title('GRFICA DE SUPERFICIES')
grid on
zoom on
subplot(2,2,2)
z= 3*(1-x).^2.*exp(-(x.^2) - (y+1).^2) ...
- 10*(x/5 - x.^3 - y.^5).*exp(-x.^2-y.^2) ...
- 1/3*exp(-(x+1).^2 - y.^2);
mesh(x,y,z);
xlabel('EJE X')
ylabel('EJE Y')
zlabel('EJE Z')
grid on
zoom on
subplot(2,2,3)
39
1
3
( x 1) 2 y 2
(funcin peaks)
z=x.^2;
surf(x,y,z);
hold on
z=4-y.^2;
mesh(x,y,z);
hold on
z=y.^2+sin(x);
mesh(x,y,z);
xlabel('EJE X')
ylabel('EJE Y')
zlabel('EJE Z')
grid on
zoom on
subplot(2,2,4)
[x y]=meshgrid(-3:.5:3,-3:.1:3);
z=peaks(x,y);
ribbon(y,z);
mesh(x,y,z);
xlabel('EJE X')
ylabel('EJE Y')
zlabel('EJE Z')
grid on
zoom on
40
E-15) Realizar el contorno de la superficie
3(1 x)2 e
x 2 ( y 1) 2
10( 5x
x3
y 5 )e
x2 y2
1
3
( x 1) 2 y 2
Resolucin
z=['3*(1-x).^2.*exp(-(x.^2) - (y+1).^2)', ...
'- 10*(x/5 - x.^3 - y.^5).*exp(-x.^2-y.^2)', ...
'- 1/3*exp(-(x+1).^2 - y.^2)'];
ezcontour(z,[-4,4],45)
grid on
xlabel('EJE X')
ylabel('EJEY')
title('Contornos')
E-16) Realizar el contorno rellenado de la superficie
3(1 x) 2 e
x 2 ( y 1) 2
10( 5x
x3
y 5 )e
x2 y2
1
3
Resolucin
z=[ ' 3*(1-x).^2.*exp(-(x.^2) - (y+1).^2)', ...
'- 10*(x/5 - x.^3 - y.^5).*exp(-x.^2-y.^2)', ...
'- 1/3*exp(-(x+1).^2 - y.^2) ' ];
ezcontourf(z,[-4,4],45)
grid on
xlabel('EJE X')
ylabel('EJEY')
title('Contorno rellenado')
41
( x 1) 2 y 2
E-17) Realizar la interseccin de los contornos de las superficies
f ( x; y )
g ( x; y )
sen x 2
x2
y2
y2
y
1 x
y2
Resolucin
f=sin(sqrt(x.^2+y.^2))./(sqrt(x.^2+y.^2));
contour(f,5,'r')
hold on
g=y./(1+x.^2+y.^2);
contour(g,6,'g')
grid on
xlabel('EJE X')
ylabel('EJEY')
title('Interseccin de contornos')
E-18) Grfico de la esfera

Resolucin
>>sphere
42
axis equal
xlabel('EJE X')
ylabel('EJE y')
zlabel('EJE Z')
title('La Esfera')
SESIN 09 - CLCULO SIMBLICO: ANLISIS MATEMTICO Y LGEBRA

CALCULO SIMBLICO CON MATLAB
MATLAB dispone del mdulo Symbolic Math toolbox, que permite manejar perfectamente el
clculo matemtico simblico, manipular con facilidad y rapidez las frmulas y expresiones
algebraicas y realizar la mayora de las operaciones con las mismas. Es posible expandir,
factorizar y simplificar polinomios y expresiones racionales y trigonomtricas; encontrar
soluciones algebraicas de ecuaciones polinmicas y sistemas de ecuaciones; evaluar
derivadas e integrales simblicamente y encontrar funciones solucin de ecuaciones
diferenciales; manipular series de potencias, lmites y muchas otras facetas de la matemtica
algebraica. Para realizar esta tarea, MATLAB requiere que todas las variables (o expresiones
algebraicas) sean declaradas como simblicas previamente con el comando syms (o con
sym).
El toolbox de matemtica simblica de MATLAB aporta varios comandos para definicin y
conversin de variables simblicas que se explican a continuacin:
FUCNCIN
DESCRIPCIN
syms x y z t
Convierte las variables x,y,z,,t en simblicas
syms x y z t real
Convierte las variables x,y,z,,t
en simblicas con valores
reales
syms x y z t unreal
43
Convierte las variables x, y, z ,, t en simblicas con valores
no reales
syms
Lista las variables simblicas en el espacio de trabajo
x=syms(x)
Convierte la variable x en simblica (equivale a syms x)
x=syms(x,real)
Convierte a x en una variable simblica real
x=syms(x,unreal)
Convierte a x en una variable simblica no real
S=syms(A)
Crea un objeto simblico a partir de A, donde A puede ser una

cadena, un escalar, una matriz, una expresin numrica, etc.
S=syms(A, opcion)
Convierte la matriz, escalar o expresin numrica A a smbolos

segn la opcin especificada. La opcin puede ser f para
punto flotante, r para racional, e para formato de error y d
para decimal.
numeric(x) o double(x)
Convierte la variable o expresin x a numrica de doble

precisin
syms2poly( poli )
Convierte el polinomio simblico poli en un vector cuyas

componentes son sus coeficientes.
poly2sym(vector)
Convierte
el
vector
en
un
polinomio
simblico
cuyos
coeficientes son las componentes del vector.

poly2sym(vector, v)
Convierte el vector en un polinomio simblico en la variable v

cuyos coeficientes son las componentes del vector.
char(X)
Convierte el array X de enteros ASCII a su valor como cadena
latex(S)
Convierte a codigo latex la expresin simblica S
ccode(S)
Convierte a cdigo C la expresin simblica S
pretty(expr)
Convierte la expresin simblica a escritura matemtica
vpa(expr)
Resultado numrico de la expresin con los dgitos decimales

de precisin situados en digits
vpa(expr, n)
Resultado numrico de la expresin con n dgitos decimales
vpa(expr, n)
Resultado numrico de la expresin con n dgitos decimales
findsym(S)
Devuelve todas las variables simblicas en la expresin

simblica o matriz simblica S
Isvarname(S)
Devuelve TRUE si S es una variable simblica vlida
vectorize(S)
Inserta un punto en la cadena S antes de cualquier smbolo ^ ,

* , / con la finalidad de operar vectorialmente de forma
correcta.
1. Calcular la funcin inversa de f ( x)
5 Log ( x 2)
Resolucin
>> syms x
44
f=5-log(x+2);
f1=finverse(f)
pretty(f1)
% f1 =exp(5-x)-2
2. Calcular la funcin inversa de g ( x)
x 5
Resolucin
>> syms x
g=4-sqrt(x+5);
g1=finverse(g)
pretty(g1)
% g1 =11-8*x+x^2
2
% 11 - 8 x + x
3. Calcular Lim
n
6n 3 4n 2 4n 5
2n 3 7 n 2
Resolucin
>> syms n
f=(6*n^3+4*n^2-4*n+5)/(2*n^3-7*n+2)
limit(f,inf)
%f =(6*n^3+4*n^2-4*n+5)/(2*n^3-7*n+2)
%ans =3
>> % en forma equivalente

limit((6*n^3+4*n^2-4*n+5)/(2*n^3-7*n+2),inf)
%ans =3
4. Calcular Lim
n
1 n
n
Resolucin
>> syms n
f=((1+n)/n)^(1/2);
limit(f,inf)
%ans =1
5. Calcular Lim
n
x3
3x 2 2 x 1
x3 2 x 5
x 1
Resolucin
45
>> syms x
f=((x.^3+3*x.^2+2*x-1)./(x.^3+2*x-5)).^(x+1);
limit(f,inf)
%ans =exp(3)
6. Calcular Lim
2x
2x
Resolucin
>> syms x
f=(sqrt(x+sqrt(2*x))-sqrt(x-sqrt(2*x)));
limit(f,inf)
%ans =2^(1/2)
7. Calcular Lim
x
x2
5x 6
x 2
Resolucin
>> syms x
f=(x^2-5*x+6)/(x-2);
limit(f,2)
%ans =-1
8. Calcular Lim
x
sen (7 x) sen (3 x)
x. cos( x)
Resolucin
>> syms x
f=(sin(7*x)-sin(3*x))./(x.*cos(x));
limit(f,0)
%ans =4
9. Calcular Lim
x
cos( x)
x2
cos( x)
Resolucin
>> syms x
f=(2-sqrt(cos(x))-cos(x))./(x.^2);
limit(f,0)
%ans =3/4
10. Calcular Lim
x
x 3
x 9
Resolucin
>> syms x
f=(sqrt(x)-3)/(x-9);
46
limit(f,9)
%ans =1/6
11. Calcular la derivada de f ( x)
x3
2e
2x
3sin( 3x)
Resolucin
>> syms x
f=x^3-2*exp(-2*x)+sin(3*x);
p=diff(f)
%p =3*x^2+4*exp(-2*x)+3*cos(3*x)
x3
12. Calcular la segunda derivada de f ( x)
2e
2x
3sin( 3x)
Resolucin
>> syms x
f=x^3-2*exp(-2*x)+sin(3*x);
p=diff(diff(f))
%p =6*x-8*exp(-2*x)-9*sin(3*x)
13. Calcular la tercera derivada de f ( x)
x3
2e
2x
3sin( 3x)
Resolucin
>> syms x
f=x^3-2*exp(-2*x)+sin(3*x);
p=diff(diff(diff(f)))
%p =6+16*exp(-2*x)-27*cos(3*x)
% En forma equivalente
>> syms x
f=x^3-2*exp(-2*x)+sin(3*x);
p=diff(f,3)
%calcula la tercera derivada
%p =6+16*exp(-2*x)-27*cos(3*x)
14. Sea f ( x; y )
cos( 2 x 3 y) 2 xy e xy . Calcular
f ( x; y )
x
Resolucin
>> syms x y
f=cos(2*x+3*y)-2*x*y-exp(x*y);
p=diff(f,x)
%derivada parcial respecto de x
pretty(p)
%p =-2*sin(2*x+3*y)-2*y-y*exp(x*y)
% -2 sin(2 x + 3 y) - 2 y - y exp(x y)
47
15. Sea f ( x; y )
f ( x; y )
y
Resolucin
>> syms x y
p=diff(f,y) %derivada parcial respecto de y
pretty(p)
%p =-3*sin(2*x+3*y)-2*x-x*exp(x*y)
% -3 sin(2 x + 3 y) - 2 x - x exp(x y)
16. Sea f ( x; y )
f ( x; y )
y x
Resolucin
>> syms x y
p=diff(diff(f,x),y)
%derivada parcial cruzada de df/dydx
pretty(p)
%p =-6*cos(2*x+3*y)-2-exp(x*y)-y*x*exp(x*y)
%
-6 cos(2 x + 3 y) - 2 - exp(x y) - y x exp(x y)
17 Calcular
( sen ( x) 3 x 2
2x
)dx
Resolucin
>> syms x
f=sin(x)-3*x^2+exp(-2*x);
int(f,x)
%ans =-cos(x)-x^3-1/2*exp(-2*x)
18. Calcular
x2 5x 9
dx
x2 5x 6
Resolucin
>> syms x
f=(x.^2-5*x+9)./(x.^2-5*x+6);
int(f,x)
%ans =x-3*log(x-2)+3*log(x-3)
19. Calcular
48
2 x 2 41x 91
dx
( x 1)( x 3)( x 4)
Resolucin
>> syms x
f=(2*x.^2+41*x-91)./((x-1).*(x+3).*(x-4));
int(f,x)
%ans =-7*log(x+3)+5*log(x-4)+4*log(x-1)
1
dx
8 4sen ( x) 7 cos( x)
20. Calcular
Resolucin
>> syms x
f=(1)./(8-4*sin(x)+7*cos(x));
int(f,x)
%ans =log(tan(1/2*x)-5)-log(tan(1/2*x)-3)
1 sen ( x) cos( x)
dx
1 sen ( x) cos( x)
21. Calcular
Resolucin
>> syms x
f=(1-sin(x)+cos(x))./(1+sin(x)-cos(x));
int(f,x)
%ans =-2*log(tan(1/2*x)+1)+2*log(tan(1/2*x))-x
(cos( x) 5 y 2
22. Calcular
2y
)dxdy
Resolucin
>> syms x y
f=cos(x)-5*y^2+exp(-2*y);
p=int(int(f,x),y)% integral doble de f respecto a x,luego a y
%p =sin(x)*y-5/3*y^3*x-1/2*exp(-2*y)*x
( sen ( x) 3 y 2
23. Calcular
2 z 3 )dxdydz
Resolucin
>> syms x y z
f=sin(x)-3*y^2+2*z^3;
p=int(int(int(f,x),y),z)
%p =-cos(x)*y*z-y^3*x*z+1/2*z^4*x*y
2
( x2
24. Calcular
2 x 3) dx
Resolucin
>> syms x
f=(x.^2-2*x+3);
49
p=int(f,x,1,2)
%p =7/3
4
25. Calcular
1
dx
x 3x 2
2
Resolucin
>> syms x
f=1./(x.^2-3*x+2);
p=int(f,x,3,4)
%p =2*log(2)-log(3)
4
( x2
26. Calcular
0
ysen( x)) dydx
Resolucin
>> syms x y
f=x^2+y*sin(x);
int(int(f,y,-2,5),x,0,4)
%ans =959/6-21/2*cos(4)
1 1
0 0
27. Calcular
1
dz dy dx
y z 1
Resolucin
>> syms x y z
f=1/sqrt(x+y+z+1);
p=int(int(int(f,z,0,1),y,0,1),x,1,0)
%p =-248/15+72/5*3^(1/2)-32/5*2^(1/2)
2
4 3
3 1
( x. y.sen( x z ) 6 y)dx dy dz
28. Calcular
Resolucin
>> syms x y z
f=x*y*(sin(x-z))+6*y;
p=int(int(int(f,x,1,3),y,-3,4),z,0,2)
%p =35/2*sin(1)+84-7/2*cos(3)-21/2*sin(3)+7/2*cos(1)
50

Ecuaciones Matlab

Cargado por

Copyright:

Formatos disponibles

Ecuaciones Matlab

Cargado por

Información del documento

Descripción original:

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

Ecuaciones Matlab

Cargado por

Copyright:

Formatos disponibles

TEORA DE ECUACIONES CON MATLAB

m.10 r entonces su equivalente es x

m.e.r as como por ejemplo:

Escalares, vectores o matrices

Escalares, vectores o matrices

a / b : Realiza la divisin de a entre b

a / b : Realiza la divisin de a entre b

a \ b : Realiza la divisin de b entre a

a .\ b : Realiza la divisin de b entre a

a ^ b : Eleva la base a al exponente b

a .^ b : Eleva la base a al exponente b

: Establece la condicin de a mayor que b

: Establece la condicin de a menor que b

a >= b : Establece la condicin de a mayor o igual que b

a <= b : Establece la condicin de a menor o igual que b

Dr. SORIA QUIJAITE JUAN JESS

: Establece la condicin de a igual a b

Negacin lgica (NOT) o complementario de A

Conjuncin lgica (AND) o interseccin de A y B

Disyuncin lgica inclusiva (OR) o unin de A y B

OR exclusivo (XOR) o diferencia simtrica de A y B

Dr. SORIA QUIJAITE JUAN JESS

psilon de la mquina =1.0000e-006

Mayor nmero (2 e)21023

Contador de operaciones de punto flotante

Nmero de argumentos de entrada de una funcin

Nmero de argumentos de salida de una funcin

7. COMANDOS DE LECTURA Y ESCRITURA

: indica que la impresin de la variable vol ser en la siguiente lnea

: formato de un nmero entero

: formato de un nmero real con 5 decimales.

Dr. SORIA QUIJAITE JUAN JESS

: Limpia la informacin de la ventana de comandos

8. FUNCIONES MATEMTICAS EN MATLAB

FUNCIONES EXPONENCIALES Y LOGARTMICAS

Funcin exponencial en base e (e^x)

Funcin logaritmo en base e de x

Funcin logaritmo en base 10 de x

Dr. SORIA QUIJAITE JUAN JESS

Funcin logaritmo en base 2 de x

Funcin potencia de base 2 de x

Funcin raz cuadrada de x

FUNCIONES ESPECFICAS DE VARIABLE NUMRICA

Valor absoluto del real x

El mayor valor entero o igual que el

El menor entero mayor o igual que el

El entero ms prximo al real x

Da el resto de la divisin entre los

Signo del real x (1 si x>0 ; -1 si x<0)

Dr. SORIA QUIJAITE JUAN JESS

str2mat('IEP "DE LA CRUZ"','MATEMATICA','LOS TRIANGULOS','PROPIEDADES','EJERCICIOS')

Transformar el numeral 3121(4) al

Dr. SORIA QUIJAITE JUAN JESS

Colocar V o F segn corresponda:

Transformar el numeral 101101(2) al

Si: abc (9) = 175

1. ENTORNO DE TRABAJO DEL MATLAB

Dr. SORIA QUIJAITE JUAN JESS