R語言中兩個數組(或向量)的外積怎樣計算

所謂數組(或向量)a和b的外積,指的是a的每個元素和b的每個元素搭配在一起相乘得到的新元素.當然運算規則也可自己定義.外積運算符為 %o%(注意:百分號中間的字母是小寫的字母o).比如:

> a <- 1:2
> b <- 3:5
> d <- a %o% b
> d
     [,1] [,2] [,3]
[1,]    3    4    5
[2,]    6    8   10

注意維數公式為:

dim(d) = c( dim(a) , dim(b) )

實際上R語言提供了一個更為一般化得外積函數outer(X,Y,Fun),比較簡單,看個演示樣例吧.

以下的樣例是為了繪制 z = x**2 +50cos(y)的3維表面圖:

> x <- seq(from = -5, to = 5, by = 0.25)
> y <- seq(from = -2, to = 2, by = 0.1)
> f <- function(x,y){x**2+50*cos(y)}
> z <- outer(x,y,f)
> persp(x,y,z)

繪制圖像例如以下:

其他看幫助文檔吧.