R语言中向量和矩阵简单运算的实现

 更新时间:2021年03月11日 10:08:26   作者:init_bin  
这篇文章主要介绍了R语言中向量和矩阵简单运算的实现,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧

一、向量运算

向量是有相同基本类型的元素序列,一维数组,定义向量的最常用办法是使用函数c(),它把若干个数值或字符串组合为一个向量。

  1.R语言向量的产生方法

> x <- c(1,2,3)
> x
[1] 1 2 3

 2.向量加减乘除都是对其对应元素进行的,例如下面

> x <- c(1,2,3)
> y <- x*2
> y
[1] 2 4 6

(注:向量的整数除法是%/%,取余是%%。)

3.向量的内积,有两种方法。

第一种方法:%*%

> x <- c(1,2,3)
> y <- c(4,5,6)
> z <- x%*%y
> z
   [,1]
[1,]  32

第二种方法:crossprod(x,y).

> x <- c(1,2,3)
> y <- c(4,5,6)
> z <- crossprod(x,y)
> z
   [,1]
[1,]  32

4.向量的外积 ,有三种方法。

第一种方法:%o%

> x <- c(1,2,3)
> y <- c(4,5,6)
> x%o%y
   [,1] [,2] [,3]
[1,]  4  5  6
[2,]  8  10  12
[3,]  12  15  18

第二种方法:tcrossprod(x,y)

> x <- c(1,2,3)
> y <- c(4,5,6)
> tcrossprod(x,y)
   [,1] [,2] [,3]
[1,]  4  5  6
[2,]  8  10  12
[3,]  12  15  18

第三种方法:outer(x,y)

> x <- c(1,2,3)
> y <- c(4,5,6)
> outer(x,y)
   [,1] [,2] [,3]
[1,]  4  5  6
[2,]  8  10  12
[3,]  12  15  18

二、矩阵的运算

1. 矩阵的产生方式

> x <- matrix(1:9,3,3)
> x
   [,1] [,2] [,3]
[1,]  1  4  7
[2,]  2  5  8
[3,]  3  6  9

其中第一个3表示的是行数,第二个3表示的列数。 故产生一个3*3的矩阵。这里是将1到9按列排列,如果想按行排列,那么如下代码

> x <- matrix(1:9,3,3,byrow = TRUE)
> x
   [,1] [,2] [,3]
[1,]  1  2  3
[2,]  4  5  6
[3,]  7  8  9

2.矩阵对应元素的运算

> x <- matrix(1:9,3,3)
> y <- matrix(9:1,3,3)
> x*y
   [,1] [,2] [,3]
[1,]  9  24  21
[2,]  16  25  16
[3,]  21  24  9

3.矩阵的转置

> x <- matrix(1:9,3,3)
> x
   [,1] [,2] [,3]
[1,]  1  4  7
[2,]  2  5  8
[3,]  3  6  9
> t(x)
   [,1] [,2] [,3]
[1,]  1  2  3
[2,]  4  5  6
[3,]  7  8  9

4.矩阵乘法

> x <- matrix(1:9,3,3)
> y <- matrix(9:1,3,3)
> x%*%y
   [,1] [,2] [,3]
[1,]  90  54  18
[2,] 114  69  24
[3,] 138  84  30

5.矩阵 x乘y的转置,x的转置乘以y

> x <- matrix(1:9,3,3)
> y <- matrix(9:1,3,3)
> crossprod(x,y)
   [,1] [,2] [,3]
[1,]  46  28  10
[2,] 118  73  28
[3,] 190 118  46
#这个是x的转置乘以y
> tcrossprod(x,y)
   [,1] [,2] [,3]
[1,]  54  42  30
[2,]  72  57  42
[3,]  90  72  54
#这个是x乘以y的转置

6.求矩阵的行列式、对称矩阵的特征值、特征向量

> x <- matrix(1:9,3,3)
> x
   [,1] [,2] [,3]
[1,]  1  4  7
[2,]  2  5  8
[3,]  3  6  9
> det(x)
[1] 0
#这个是求特征值的
> x <- matrix(1:9,3,3)
> x
   [,1] [,2] [,3]
[1,]  1  4  7
[2,]  2  5  8
[3,]  3  6  9
> a <- crossprod(x,x)
> a1 <- eigen(a)  #这个是得到对称矩阵特征值、特征向量的主要函数
> a1
eigen() decomposition
$`values`  #这个是特征值
[1] 2.838586e+02 1.141413e+00 6.308738e-15
 
$vectors   #这个是特征向量
      [,1]    [,2]    [,3]
[1,] -0.2148372 0.8872307 0.4082483
[2,] -0.5205874 0.2496440 -0.8164966
[3,] -0.8263375 -0.3879428 0.4082483

到此这篇关于R语言中向量和矩阵简单运算的实现的文章就介绍到这了,更多相关R语言向量运算和矩阵运算内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家! 

相关文章

  • 详解R语言caret包trainControl函数

    详解R语言caret包trainControl函数

    这篇文章主要介绍了R语言caret包trainControl函数详解,本文通过源码分析给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友可以参考下
    2022-08-08
  • R语言绘制交互地图的实现示例

    R语言绘制交互地图的实现示例

    本文主要介绍了R语言绘制交互地图的实现示例,使用leaflet包绘制交互地图,具有一定的参考价值,感兴趣的可以了解一下
    2025-03-03
  • R语言基本运算的示例代码

    R语言基本运算的示例代码

    这篇文章主要介绍了R语言基本运算的示例代码,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧
    2021-03-03
  • R语言绘制数据可视化小提琴图Violin plot with dot画法

    R语言绘制数据可视化小提琴图Violin plot with dot画法

    这篇文章主要为大家介绍了R语言绘制数据可视化小提琴图Violin plot with dot画法的示例详解有需要的朋友可以借鉴参考下希望能够有所帮助,祝大家多多进步,早日升职加薪
    2022-02-02
  • R语言绘图基础教程(新手入门推荐!)

    R语言绘图基础教程(新手入门推荐!)

    数据作图是数据分析的重要方法之一,R提供了丰富的作图函,下面这篇文章主要给大家介绍了关于R语言绘图基础教程的相关资料,文中通过实例代码以及图文介绍的非常详细,需要的朋友可以参考下
    2022-11-11
  • R语言利用ggplot2绘制QQ图和箱线图详解

    R语言利用ggplot2绘制QQ图和箱线图详解

    这篇文章主要为大家介绍了R语言如何利用ggplot2绘制QQ图和箱线图,文中的示例代码讲解详细,对我们学习R语言有一定的帮助,需要的可以参考一下
    2022-06-06
  • 解决R语言报错:Error in y + 1:non-numeric argument to binary operator

    解决R语言报错:Error in y + 1:non-numeric argument to binary

    R语言编程中的常见错误有一些错误是R的初学者和经验丰富的R程序员都可能常犯的,下面这篇文章主要给大家介绍了关于解决R语言报错:Error in y + 1:non-numeric argument to binary operator的相关资料,需要的朋友可以参考下
    2022-11-11
  • 基于R语言xlsx安装遇到的问题及解决方案

    基于R语言xlsx安装遇到的问题及解决方案

    这篇文章主要介绍了基于R语言xlsx安装遇到的问题及解决方案,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
    2021-04-04
  • R语言对二进制文件操作详解

    R语言对二进制文件操作详解

    在本篇内容里小编给大家整理的是一篇关于R语言二进制文件操作实例内容,有兴趣的朋友们可以学习下。
    2021-05-05
  • R语言绘制小提琴图violin plot实现示例

    R语言绘制小提琴图violin plot实现示例

    这篇文章主要为大家介绍了R语言绘制小提琴图violin plot的实现方法示例,有需要的朋友可以借鉴参考下,希望能够有所帮助,祝大家多多进步
    2022-02-02

最新评论