Matlab

PARTE EXPERIMENTAL
Considere la siguiente matriz:

1.- Indique el resultado de las siguientes operaciones y que significa cada una de ellas.

A=[11 12 13; 21 22 23; 31 32 33]

a) A(:,1)

Muestre todos los elementos de la primera columna.

ans =
11
21
31
b) A(2,:)

Muestretodos los elementos de la segunda fila.

ans =
21 22 23
c) A(:,2:3)

Muestre todos los elementos desde la segunda hasta la tercera columna.

ans =
12 13
22 23
32 33
d) A(:,1:2:3)

Muestre todos los elementos desde la primera hasta la terceracolumna, pero cada 2 columnas.

ans =
11 13
21 23
31 33
e) B=[A,[ones(1,2),eye(2)]]

La matriz B está conformada por los elementos de la matriz A mas otra columna compuesta por una matriz de unos (de una fila por dos columnas), como primera fila; y una matriz identidad de 2*2, en la siguientefila.

B =
11 12 13 1 1
21 22 23 1 0
31 32 33 0 1
f) A(:,:)

Mostrar los elementos de todas las filas y todas las columnas.

ans =
11 12 13
21 22 23
31 32 33g) A(2:4)=[]
Muestra los elementos de A, a excepción de la segunda a cuarta entrada como columna.
A =
11 22 32 13 23 33
h) Ones(3,3)
Genera una matriz de unos de 3*3
ans =
1 1 1
1 1 1
1 1 1
i) Eye(3)Genera una matriz identidad de 3*3

ans =
1 0 0
0 1 0
0 0 1

j) Diag(A)
Muestra los elementos de la diagonal principal de A.
ans =
11
22
33
k) Zeros(size(A))

Genera una matriz deceros del tamaño de A.

ans =
0 0 0
0 0 0
0 0 0
l) Rand(size(A))
Genera una matriz de valores aleatorios del tamaño de A.

ans =
0.8147 0.9134 0.2785
0.9058 0.6324 0.54690.1270 0.0975 0.9575

m) magic(length(A))

Genera una matriz mágica del tamaño de A. O sea q la suma de los elementos de cada fila o columna es el mismo valor siempre.

ans =
8 1 6
3 5 7
4 9 2

2.- Grafique lassiguientes funciones sobre un mismo gráfico. Los vectores deberían tener al menos 200 puntos para su graficación.
[pic]
Código:
t = -5 : 0.007 : 5;
y = (1-sin(t).^2)./t;
z = (t+1)./(t.^2+4);
plot(t,y,t,z);

Grafica:

[pic]

3.- Utilizando el comando subplot, grafique las funciones del numeral anterior.

Código:
t = -5 : 0.007 :5;
y = (1-sin(t).^2)./t;
z = (t+1)./(t.^2+4);
subplot(2,1,1);
plot(t,y);
subplot(2,1,2);
plot(t,z);
xlabel ('t');
ylabel ('z(t)');
grid on;
hold;

Gráfica:

[pic]

4.- Grafique la órbita de Mercurio (x vs. y)...