1.
Practica de laboratorio
1.1.-Parte práctica instrucciones básicas: desde la ventana de comandos digite las instrucciones y
pulse la tecla enter, en cada caso explicar el resultado:
>> sqrt(7) {raíz cuadrada de 7}
ans =2.6458
>> sqrt(7/5) {raíz cuadrada de 7/5}
ans =1.1832
>> a=2.1; sqrt(2*a) {raíz cuadrada de a}
ans =2.0494
>> exp (3) {épsilon elevada al cubo}
ans =20.0855
>> x=5*cos(pi/6) { 5 por coseno de pi entre seis}
x =4.3301
>> y=5*sin(pi/6) {5 por seno de pi entre seis)
y =2.5000
>> acos(x/5), asin(y/5) {arco coseno de ‘x’ entre cinco, arco seno de ‘y ‘entre cinco}
ans =0.5236
ans =0.5236
>> (3+4^2)/(2/(3^(1/5))-(1/(3.1-2))^(3/4)) {Matlab trabaja como una simple calculadora al escribir
expresiones aritméticas de acuerdo a las prioridades y formatos}
ans = 28.1702
>> (5+3^2)/((2*4^(1/3)-(1/3.1-2)^3))
ans = 1.7734
1.2.-Variables: Matlab no requieren ningún tipo de comandos para declarar variables, crea la
variable mediante la asignación de su valor, Procese las siguientes instrucciones:
>> fprintf('formato_e:%12.6e\n',12345,2)
formato_e:1.234500e+04
formato_e:2.000000e+00
>> v=4+6+7-4+12-7 {se crea la variable v}
v = 18
>> v^4 {la variable se eleva a la cuarta}
ans = 104976
>> b=pi {se asigna pi a la variable b}
b = 3.1416
>> clock {nos indica la fecha en valores decimales}
ans = 1.0e+03 *
2.0190 0.0040 0.0070 0.0130 0.0510 0.0213
>> fix(clock) {nos indica la fecha en valores reales}
ans = 2019 4 7 13 51 45
>> date {nos indica la fecha}
ans = '07-Apr-2019'
>> t_0=fix(clock)
t_0 = 2019 4 7 13 52 15
>> t_1=fix(clock)
t_1 = 2019 4 7 13 52 26
>> t_1-t_0
ans = 0 0 0 0 0 11
>> v=[3,4,5,6,13/7,2,6] {es un vector fila}
v = 3.0000 4.0000 5.0000 6.0000 1.8571 2.0000 6.0000
1.3.- Variables vectoriales: formato V=v1,v2,v3 ,....,vn o V= v1 v2 v3 .... vn , practique:
>> vector1=[5:20] {devuelve los elementos del vector situados del 5 al 20 y ambos incluidos}
vector1 = 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
>> vector2=[Link] {devuelve los elementos del vector situados del 5 al 20 y ambos incluidos
pero separados de 3 en 3 unidades}
vector2 = 2 5 8 11 14 17 20
>> vector3=linspace(10,30,6) {devuelve los elementos situados del 10 al 30 de la operación (30-
10)/(6-1) de 4 en 4 unidades }
vector3 = 10 14 18 22 26 30
>> vector4=[5:20]; (devuelve los elementos situados del 5 al 20 y ambos incluidos)
vector4 = 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
>> vector5=logspace(10,30,6)
vector5 = 1.0e+30 *
0.0000 0.0000 0.0000 0.0000 0.0001 1.0000
1.4.-Variables matriciales:
N [a11,a12,a13 ,a14 ,a21;a22,a23 ,a24;a31,a32,a33,a34 ]; M1=[1 2 3 4 5 6;3 4 5 6 7 9;11 3 5 7 15 8]
N=M1
Practique:
>> A=[1 2 3;4 5 6;7 8 9] {devuelve toda la matriz 3x3}
A=
1 2 3
4 5 6
7 8 9
>> x=[-1.3,sqrt(3),(1+2+3)*4/5]
x = -1.3000 1.7321 4.8000
>> x(5)=abs(x(1))
x = -1.3000 1.7321 4.8000 0 1.3000
>> r=[10 11 12] {creamos una matriz de una fila r}
r = 10 11 12
>> A=[A;r] {añadimos una fila a la matriz A la matriz que creamos r}
A=
1 2 3
4 5 6
7 8 9
10 11 12
>> x=[1 2 3];y=[4 5 6]; {realiza la multiplicación de dos matrices de una fila}
>> z=x*y
z= 4 10 18
>> x=1:5
x= 1 2 3 4 5
>> A=[1 2 3;4 5 6;7 8 9] {devuelve una matriz 3x3}
A=
1 2 3
4 5 6
7 8 9
>> A(3,3)=A(1,3)+A(3,1) {devuelve la suma de la primera columna tercera fila más la tercera
columna y primera fila en la tercera columna y tercera fila}
A=
1 2 3
4 5 6
7 8 10
1.5.- Investigue y explique qué resultados producirá las instrucciones, sabiendo que A es una
matriz de 10x10.
>> B=[1,2,3,4,5,6,7,8,9,10] {matriz 1x10}
B= 1 2 3 4 5 6 7 8 9 10
>> C=[1;2;3;4;5;6;7;8;9;10] {matriz 10x1}
C=
10
>> A=C*B {matriz 10x10}
A=
1 2 3 4 5 6 7 8 9 10
2 4 6 8 10 12 14 16 18 20
3 6 9 12 15 18 21 24 27 30
4 8 12 16 20 24 28 32 36 40
5 10 15 20 25 30 35 40 45 50
6 12 18 24 30 36 42 48 54 60
7 14 21 28 35 42 49 56 63 70
8 16 24 32 40 48 56 64 72 80
9 18 27 36 45 54 63 72 81 90
10 20 30 40 50 60 70 80 90 100
>> A(1:5,3) {indica los elementos de 1 al 5 en la columna 3}
ans =
12
15
>> A(1:5,2) {indica los elementos de 1 al 5 en la columna 2}
ans =
10
>> A(1:5,7:10) {indica los elementos de 1 al 5 en la columna de 7 al 10}
ans =
7 8 9 10
14 16 18 20
21 24 27 30
28 32 36 40
35 40 45 50
>> D=A
D=
1 2 3 4 5 6 7 8 9 10
2 4 6 8 10 12 14 16 18 20
3 6 9 12 15 18 21 24 27 30
4 8 12 16 20 24 28 32 36 40
5 10 15 20 25 30 35 40 45 50
6 12 18 24 30 36 42 48 54 60
7 14 21 28 35 42 49 56 63 70
8 16 24 32 40 48 56 64 72 80
9 18 27 36 45 54 63 72 81 90
10 20 30 40 50 60 70 80 90 100
>> D=A;
>> A(:,[3 5 10])=D(:,1:3) {las columnas 3,5 y 10 se cambian respectivamente por las columnas 1 al
3}
A=
1 2 1 4 2 6 7 8 9 3
2 4 2 8 4 12 14 16 18 6
3 6 3 12 6 18 21 24 27 9
4 8 4 16 8 24 28 32 36 12
5 10 5 20 10 30 35 40 45 15
6 12 6 24 12 36 42 48 54 18
7 14 7 28 14 42 49 56 63 21
8 16 8 32 16 48 56 64 72 24
9 18 9 36 18 54 63 72 81 27
10 20 10 40 20 60 70 80 90 30
1.6.- Dado B 1 2 3 4 5;6 7 8 9 10;2 4 5 6 7;3 6 1 1 3 obtenga:
>> B=[1 2 3 4 5;6 7 8 9 10;2 4 5 6 7; 3 6 1 1 3]
B=
1 2 3 4 5
6 7 8 9 10
2 4 5 6 7
3 6 1 1 3
>> C=[B,eye(4)] {C es la matriz aumentada de B}
C=
1 2 3 4 5 1 0 0 0
6 7 8 9 10 0 1 0 0
2 4 5 6 7 0 0 1 0
3 6 1 1 3 0 0 0 1
>> D=C(:,[Link]) {D matriz de elementos de C (fila:1 al 4,columna: 1 al 4)}
D=
1 2 3 4
6 7 8 9
2 4 5 6
3 6 1 1
>> F=C([1,2],3:5) {F matriz de elementos de C (columnas: 3 al 5, filas 1 y 2)}
F=
3 4 5
8 9 10
>> H=C([1,3],[2,3,5]) {H matriz de elementos de C (columnas: 2,3 y 5; filas:1 y 3)}
H=
2 3 5
4 5 7
>> [eye(5,4) zeros [5,4] ones(5,4)]
1.7.- comandos que manejas caracteres, practique:
>> setstr(32:126)
ans = ' !"#$%&'()*+,-
./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{|}~'
>> abs('{]}><#!¿?\')
ans = 123 93 125 62 60 35 33 191 63 92
[Link] para resolver en grupo.
2.1.- Evaluar las siguientes expresiones matemáticas en MATLAB.
a. log(2)
b. arcsen(-0.5)
c. 5+e2
d. tan(e)
>> log(2)
ans = 0.6931
>> asin(-0.5)
ans = -0.5236
>> sqrt(5)+exp(2)
ans = 9.6251
>> tan(exp(1))
ans = -0.4505
2.2.- Extraer las siguientes sub matrices de A=rand(20,10):
a. Las 5 primeras columnas
b. Las 15 últimas filas
c. Primera y quita fila
d. 10 últimos elementos de la tercera fila
e. La intersección de la fila: 2,3,17,19 con las columnas: 5,6,7,12
>> A=rand(20,10)
A=
0.6443 0.3111 0.0855 0.0377 0.0305 0.0596 0.1734 0.9516 0.0326 0.2518
0.3786 0.9234 0.2625 0.8852 0.7441 0.6820 0.3909 0.9203 0.5612 0.2904
0.8116 0.4302 0.8010 0.9133 0.5000 0.0424 0.8314 0.0527 0.8819 0.6171
0.5328 0.1848 0.0292 0.7962 0.4799 0.0714 0.8034 0.7379 0.6692 0.2653
0.3507 0.9049 0.9289 0.0987 0.9047 0.5216 0.0605 0.2691 0.1904 0.8244
0.9390 0.9797 0.7303 0.2619 0.6099 0.0967 0.3993 0.4228 0.3689 0.9827
0.8759 0.4389 0.4886 0.3354 0.6177 0.8181 0.5269 0.5479 0.4607 0.7302
0.5502 0.1111 0.5785 0.6797 0.8594 0.8175 0.4168 0.9427 0.9816 0.3439
0.6225 0.2581 0.2373 0.1366 0.8055 0.7224 0.6569 0.4177 0.1564 0.5841
0.5870 0.4087 0.4588 0.7212 0.5767 0.1499 0.6280 0.9831 0.8555 0.1078
0.2077 0.5949 0.9631 0.1068 0.1829 0.6596 0.2920 0.3015 0.6448 0.9063
0.3012 0.2622 0.5468 0.6538 0.2399 0.5186 0.4317 0.7011 0.3763 0.8797
0.4709 0.6028 0.5211 0.4942 0.8865 0.9730 0.0155 0.6663 0.1909 0.8178
0.2305 0.7112 0.2316 0.7791 0.0287 0.6490 0.9841 0.5391 0.4283 0.2607
0.8443 0.2217 0.4889 0.7150 0.4899 0.8003 0.1672 0.6981 0.4820 0.5944
0.1948 0.1174 0.6241 0.9037 0.1679 0.4538 0.1062 0.6665 0.1206 0.0225
0.2259 0.2967 0.6791 0.8909 0.9787 0.4324 0.3724 0.1781 0.5895 0.4253
0.1707 0.3188 0.3955 0.3342 0.7127 0.8253 0.1981 0.1280 0.2262 0.3127
0.2277 0.4242 0.3674 0.6987 0.5005 0.0835 0.4897 0.9991 0.3846 0.1615
0.4357 0.5079 0.9880 0.1978 0.4711 0.1332 0.3395 0.1711 0.5830 0.1788
>> B=A(:,1:5)
B=
0.6443 0.3111 0.0855 0.0377 0.0305
0.3786 0.9234 0.2625 0.8852 0.7441
0.8116 0.4302 0.8010 0.9133 0.5000
0.5328 0.1848 0.0292 0.7962 0.4799
0.3507 0.9049 0.9289 0.0987 0.9047
0.9390 0.9797 0.7303 0.2619 0.6099
0.8759 0.4389 0.4886 0.3354 0.6177
0.5502 0.1111 0.5785 0.6797 0.8594
0.6225 0.2581 0.2373 0.1366 0.8055
0.5870 0.4087 0.4588 0.7212 0.5767
0.2077 0.5949 0.9631 0.1068 0.1829
0.3012 0.2622 0.5468 0.6538 0.2399
0.4709 0.6028 0.5211 0.4942 0.8865
0.2305 0.7112 0.2316 0.7791 0.0287
0.8443 0.2217 0.4889 0.7150 0.4899
0.1948 0.1174 0.6241 0.9037 0.1679
0.2259 0.2967 0.6791 0.8909 0.9787
0.1707 0.3188 0.3955 0.3342 0.7127
0.2277 0.4242 0.3674 0.6987 0.5005
0.4357 0.5079 0.9880 0.1978 0.4711
>> C=A(6:end,:)
C=
0.9390 0.9797 0.7303 0.2619 0.6099 0.0967 0.3993 0.4228 0.3689 0.9827
0.8759 0.4389 0.4886 0.3354 0.6177 0.8181 0.5269 0.5479 0.4607 0.7302
0.5502 0.1111 0.5785 0.6797 0.8594 0.8175 0.4168 0.9427 0.9816 0.3439
0.6225 0.2581 0.2373 0.1366 0.8055 0.7224 0.6569 0.4177 0.1564 0.5841
0.5870 0.4087 0.4588 0.7212 0.5767 0.1499 0.6280 0.9831 0.8555 0.1078
0.2077 0.5949 0.9631 0.1068 0.1829 0.6596 0.2920 0.3015 0.6448 0.9063
0.3012 0.2622 0.5468 0.6538 0.2399 0.5186 0.4317 0.7011 0.3763 0.8797
0.4709 0.6028 0.5211 0.4942 0.8865 0.9730 0.0155 0.6663 0.1909 0.8178
0.2305 0.7112 0.2316 0.7791 0.0287 0.6490 0.9841 0.5391 0.4283 0.2607
0.8443 0.2217 0.4889 0.7150 0.4899 0.8003 0.1672 0.6981 0.4820 0.5944
0.1948 0.1174 0.6241 0.9037 0.1679 0.4538 0.1062 0.6665 0.1206 0.0225
0.2259 0.2967 0.6791 0.8909 0.9787 0.4324 0.3724 0.1781 0.5895 0.4253
0.1707 0.3188 0.3955 0.3342 0.7127 0.8253 0.1981 0.1280 0.2262 0.3127
0.2277 0.4242 0.3674 0.6987 0.5005 0.0835 0.4897 0.9991 0.3846 0.1615
0.4357 0.5079 0.9880 0.1978 0.4711 0.1332 0.3395 0.1711 0.5830 0.1788
>> D=A([1,5],:)
D=
0.6443 0.3111 0.0855 0.0377 0.0305 0.0596 0.1734 0.9516 0.0326 0.2518
0.3507 0.9049 0.9289 0.0987 0.9047 0.5216 0.0605 0.2691 0.1904 0.8244
>> E=A([3],1:end)
E = 0.8116 0.4302 0.8010 0.9133 0.5000 0.0424 0.8314 0.0527 0.8819 0.6171
>> F=A([2,3,17,19],[5,6,7,12])
Index exceeds matrix dimensions.
>> F=A([2,3,17,19],[5,6,7,9])
F=
0.7441 0.6820 0.3909 0.5612
0.5000 0.0424 0.8314 0.8819
0.9787 0.4324 0.3724 0.5895
0.5005 0.0835 0.4897 0.3846
2.3.- Crear los siguientes vectores:
x=[0 3 π e2 ]
Y=[0 0.1π 0.2π 0.3π 0.4π 0.5π 0.6π 0.7π 0.8π 0.9π π]
>> x=[0 sqrt(3) pi exp(2)]
x= 0 1.7321 3.1416 7.3891
>> Y=[0 0.1*pi 0.2*pi 0.3*pi 0.4*pi 0.5*pi 0.6*pi 0.7*pi 0.8*pi 0.9*pi pi]
Y = 0 0.3142 0.6283 0.9425 1.2566 1.5708 1.8850 2.1991 2.5133 2.8274 3.1416
2.4.- Crear un vector z de cuatro numeros complejos.
>> z=x+i*Y(1:4)
z = 0.0000 + 0.0000i 1.7321 + 0.3142i 3.1416 + 0.6283i 7.3891 + 0.9425i
2.5.- Listar el tercer element del vector z.
>> z(3)
ans = 3.1416 + 0.6283i
2.6.- Listar los 5 primeros elementos del vector Y.
>> Y(1:5)
ans = 0 0.3142 0.6283 0.9425 1.2566
2.7.- Listar los 5 ultimos elementos del vector Y.
>> size(Y)
ans = 1 11
>> Y((11-4):end)
ans = 1.8850 2.1991 2.5133 2.8274 3.1416
2.8.- Listar los elementos de posición impares del vector Y.
>> Y([1,3,5,7,9,11])
ans = 0 0.6283 1.2566 1.8850 2.5133 3.1416
2.9.- Listar los elementos de posiciones 2,4,5 y 7 del vector Y.
>> Y([2,4,5,7])
ans = 0.3142 0.9425 1.2566 1.8850
2.10.- Crear los vectores.
>> a=[1 2 3 4 5], b=[1 3 5 7 9]
a= 1 2 3 4 5
b= 1 3 5 7 9
2.11.- fusionar los vectores a y b en un vector c.
>> c=a+b
c= 2 5 8 11 14
2.12.- Obtener la transpuesta del vector c.
>> ct=c'
ct =
11
14
2.13.- Obtener la transpuesta del vector z.
>> zt=z'
zt =
0.0000 + 0.0000i
1.7321 - 0.3142i
3.1416 - 0.6283i
7.3891 - 0.9425i
2.14.- Crear las siguientes matrices.
>> g=[1 2 3 4;5 6 7 8],h=[1 1 1 1;2 2 2 2]
g=
1 2 3 4
5 6 7 8
h=
1 1 1 1
2 2 2 2
2.15.- Sumar las matrices g y h.
>> g+h
ans =
2 3 4 5
7 8 9 10
2.16.- Multiplicar las matrices g y h.
>> g*h
Error using *
Inner matrix dimensions must agree.
2.17.- multiplicar g con la transpuesta de h.
>> g*h'
ans =
10 20
26 52
2.18.- multiplicar g y h componente a componente.
>> g.*h
ans =
1 2 3 4
10 12 14 16
2.19.- Eleve 2 a cada elemento de g.
>> g.*g
ans =
1 4 9 16
25 36 49 64
2.20.- Obtener la inversa de cada elemento de g.
>> g./(g.*g)
ans =
1.0000 0.5000 0.3333 0.2500
0.2000 0.1667 0.1429 0.1250
2a 3b c 6
2.21.- Resolver el Sistema. 4a b 2c 7
6a b 7c 4
>> A=[2 3 1;4 1 2;6 1 7]
A=
2 3 1
4 1 2
6 1 7
>> B=[6;7;4]
B=
>> C=inv(A)*B
C=
2.2500
1.0000
-1.5000
2.22.- Resuelva el sistema anterior mediante \.
>> C=A\B
C=
2.2500
1.0000
-1.5000
2.23.- Utilizando MATLAB determine el valor de la expresión.
[(9.8-0.7)/(4-ln(2))]245/124
>> [(9.8-0.7)/(4-log(2))]^(245/124)
ans = 7.3895
2.24.- Crear la siguiente matriz.
>> T=[1 0 0 0 1 2 3 4;0 2 0 0 5 6 7 3;0 0 3 0 9 10 11 12;0 0 0 4 20 0 5 4;1 5 9 20 0 0 0 0;2 6 10 0 0 0 0
0;3 7 11 5 0 0 0 0;4 8 12 4 0 0 0 0]
T=
1 0 0 0 1 2 3 4
0 2 0 0 5 6 7 3
0 0 3 0 9 10 11 12
0 0 0 4 20 0 5 4
1 5 9 20 0 0 0 0
2 6 10 0 0 0 0 0
3 7 11 5 0 0 0 0
4 8 12 4 0 0 0 0
>> T1([3:6],[3:6])
ans =
3 0 9 10
0 4 20 0
9 20 0 0
10 0 0 0
2.25.- Escriba la matriz A.
>> A=[1 -2 0 0 0; -2 1 -2 0 0;0 0 -2 1 -2;0 0 0 -2 1]
A=
1 -2 0 0 0
-2 1 -2 0 0
0 0 -2 1 -2
0 0 0 -2 1
2.26.- Use la función diag para A.
>> diag(A)
ans =
-2
-2
2.27.-si b.
>> b=[1,2,3,4,5]'
b=
>> x=A.\b
Matrix dimensions must agree.
>> x=inv(A)*b
Error using inv
Matrix must be square.
3. Gráficos: construcción de gráficos
simples.
