Corts Franco Javier
Operaciones con
matrices
>
matrix(c(1,2,3,4,5,6,7,8,9),nrow
=3)
[,1] [,2] [,3]
[1,]
[2,]
[3,]
> matrix(1:8,ncol=2)
[2,]
[3,]
[4,]
[2,]
[3,]
> matriz1+2
[,1] [,2] [,3]
[1,]
[2,]
7 10
[3,]
8 11
> matriz2
[1,] 10 13 16
[2,] 11 14 17
[3,] 12 15 18
> matriz2*3
> vector10=(1:9)
> vector10=c(1:9)
> matrix(vector10,nrow=3)
[,1] [,2] [,3]
[1,]
[2,]
[3,]
> matriz1=matrix(1:9,nrow=3);
>
matriz2=matrix(10:18,nrow=3)
> matriz1
[1,]
[,1] [,2] [,3]
[,1] [,2]
[1,]
[,1] [,2] [,3]
[,1] [,2] [,3]
[1,] 30 39 48
[2,] 33 42 51
[3,] 36 45 54
> matriz1
[,1] [,2] [,3]
[1,]
[2,]
[3,]
> matriz2
[,1] [,2] [,3]
Corts Franco Javier
[1,] 10 13 16
[2,] 11 14 17
> matriz1+matriz2
[,1] [,2]
[3,] 12 15 18
[1,]
> matriz1+matriz2
[2,]
[,1] [,2] [,3]
[1,] 11 17 23
> matriz3-matriz1
[,1] [,2]
[2,] 13 19 25
[1,]
[3,] 15 21 27
[2,]
0 -8
............................
> 3*matriz2
> matrix(c(4,6,-1,9),nrow=2)
[,1] [,2]
[,1] [,2]
[1,]
9 -6
[1,]
4 -1
[2,]
[2,]
> 4*matriz2+2*matriz3
> matrix(c(0,3,3,-2),nrow=2)
[,1] [,2]
[1,]
[2,]
3 -2
> matrix(c(8,6,3,1),nrow=2)
[,1] [,2]
[1,]
[2,]
> matriz1=matrix(c(4,6,1,9),nrow=2)
> matriz2=matrix(c(0,3,3,2),nrow=2)
>
matriz3=matrix(c(8,6,3,1),nrow
=2)
[,1] [,2]
[1,] 16 18
[2,] 24 -6
Corts Franco Javier
matrix(c(1,2,-4,2,1,-2,1,1,1),nrow=3)
[,1] [,2] [,3]
Ejercicio con solve
1
matrix(c(2,1,3,3,1,-1,1,-1,2),
nrow=3)
[,1] [,2] [,3]
[1,]
2 -1
[2,]
[3,] -4 -2
> matrix(c(32,16,12),nrow=3)
[,1]
[1,] 32
[1,]
[2,] 16
[2,]
1 -1
[3,] 12
[3,]
3 -1
> matrix(c(12,4,6),nrow=3)
> matriz3=matrix(c(1,2,-4,2,1,2,-1,1,1),nrow=3)
[1,] 12
>
matriz4=matrix(c(32,16,12),nro
w=3)
[2,]
> solve(matriz3,matriz4)
[3,]
[,1]
[,1]
> matriz1=matrix(c(2,1,3,3,1,1,1,-1,2), nrow=3)
[1,] -14.66667
>
matriz2=matrix(c(12,4,6),nrow=
3)
[3,] 14.66667
> solve(matriz1,matriz2)
Hallar los valores y vectores
propios de la siguientes
matrices
[,1]
[2,] 30.66667
[1,] 2.3529412
[2,] 2.2352941
[3,] 0.5882353
2
> matrix(c(1,2,-4,2,1,-2,2,2,2),nrow=3)
[,1] [,2] [,3]
Corts Franco Javier
[1,]
2 -2
[2,]
[3,] -4 -2
> matrix(c(2,-2,1,3,-1,1,1,4,1),nrow=3)
[,1] [,2] [,3]
> a=matrix(c(1,2,-4,2,1,-2,2,2,2),nrow=3)
[1,]
[2,] -2 -1
> y=eigen(a)
[3,]
>y
> b=matrix(c(2,-2,1,3,-1,1,1,4,1),nrow=3)
$values
1 -1
> y=eigen(b)
[1] 3 3 -2
>y
$values
$vectors
[,1]
[,2]
[,3]
[1,] 0.4082483 -0.4082483
0.6666667
[2,] -0.4082483 0.4082483
-0.6666667
[3,] -0.8164966 0.8164966
0.3333333
> y$val
[1] -1.618034 1.000000 0.618034
$vectors
[,1]
[,2]
[,3]
[1,] 0.5257311 0.8890009
0.8506508
[2,] -0.7608452 -0.3810004
-0.4702282
[3,] 0.3804226 0.2540003
0.2351141
[1] 3 3 -2
> y$vec
[,1]
[,2]
[,3]
[1,] 0.4082483 -0.4082483
0.6666667
[2,] -0.4082483 0.4082483
-0.6666667
[3,] -0.8164966 0.8164966
0.3333333
> y$val
[1] -1.618034 1.000000 0.618034
> y$vec
[,1]
[,2]
[,3]
[1,] 0.5257311 0.8890009
0.8506508
Corts Franco Javier
[2,] -0.7608452 -0.3810004
-0.4702282
[2,] 0.3293570 -4.472136e-01
0.8834098
[3,] 0.3804226 0.2540003
0.2351141
[3,] -0.9297338 8.944272e-01
-0.1564734
[Link] la matriz inversa y los
valores y vectores propios de
las siguientes matrices
2
> matrix(c(1,-4,-1,3,1,1,2,2,1),nrow=3)
[,1] [,2] [,3]
[1,]
> matriz5=matrix(c(-2,2,2,2,1,2,1,2,-4),nrow=3)
[2,] -4 -1
> solve(matriz5)
[3,] -1
[,1]
[,2]
[,3]
[1,] 3.702550e-17 0.3333333
0.1666667
[2,] 6.666667e-01 0.3333333
0.3333333
[3,] -3.333333e-01 0.0000000
-0.3333333
> y=eigen(matriz5)
> matriz6=matrix(c(1,-4,-1,3,1,1,2,2,1),nrow=3)
> solve(matriz6)
[,1]
[,2]
[,3]
[1,] 0.4285714 0.1428571
-1.142857
[2,] -0.2857143 -0.4285714
1.428571
[3,] 0.7142857 0.5714286
-1.571429
>y
$values
> x=eigen(matriz6)
[1] -3.645751 -3.000000 1.645751
>x
$values
$vectors
[,1]
[,2]
[,3]
[1,] 0.1646785 -1.697156e-16
0.4417049
[1] 0.7933550+3.361774i
0.7933550-3.361774i
-0.5867099+0.000000i
$vectors
Corts Franco Javier
[,1]
[,2]
[,3]
Resolver los siguiente sistemas
de ecuaciones lineales.
[1,] 0.6977063+0.0000000i
0.6977063+0.0000000i
0.3940415+0i
En primera instancia sin utilizar
la funcin solve( ) y despus
usando este comando para
verificar las operaciones
realizadas
[2,] -0.1803899+0.6158444i
-0.1803899-0.6158444i
-0.6449483+0i
[3,] 0.1984960+0.2489988i
0.1984960-0.2489988i
0.6548076+0i
> matrix(c(1,4,7,2,5,8,3,6,9),nrow=3)
[,1] [,2] [,3]
> x$val
[1] 0.7933550+3.361774i
0.7933550-3.361774i
-0.5867099+0.000000i
[1,]
[2,]
[3,]
8 -9
> matrix(c(24,32,41),nrow=3)
> x$vec
[,1]
[,1]
[,3]
[1,] 0.6977063+0.0000000i
0.6977063+0.0000000i
0.3940415+0i
[2,] -0.1803899+0.6158444i
-0.1803899-0.6158444i
-0.6449483+0i
[3,] 0.1984960+0.2489988i
0.1984960-0.2489988i
0.6548076+0i
[,2]
[1,] 24
[2,] 32
[3,] 41
> d=matrix(c(1,4,7,2,5,8,3,6,9),nrow=3)
> e=matrix(c(24,32,41),nrow=3)
> solve(d,e)
[,1]
[1,] -18.72222222
[2,] 21.44444444
[3,] -0.05555556
