大佬教程收集整理的这篇文章主要介绍了Kronecker 积函数定义,大佬教程大佬觉得挺不错的,现在分享给大家,也给大家做个参考。
我有这段代码来定义一个函数,该函数可以做一个 kronecker 产品,但之后不起作用
kron <- function(A,B){
nfila = nrow(A)*nrow(B)
n@R_262_8620@ = ncol(A)*ncol(B)
m = matrix(NA,nrow = nfila,ncol = n@R_262_8620@)
for (i in 1:nfila){
for (j in 1:ncol(A)){
m[i,j]=(A[i,j]*B)
}
}
return(m)}
我使用了这两个矩阵
A <- matrix(c(1,2,3,1),nrow = 2,ncol = 2,byrow = TRUE)
B <- matrix(c(0,byrow = TRUE)
要检查 de process 是否正常,运行的 kron(A,B) 必须等于 kronecker(A,B)
如果 c<-kronecker(A,B)
和 @H_286_5@m<-kron(A,B)。然后,C==m
。
在 R 中你可以这样做:
library(purrr)
kron <- function(A,B){
map(seq_len(nrow(A)),~do.call("cbind",map(A[.x,],~.x * B))) %>%
do.call("rbind",.)
}
内循环(映射):将矩阵 A 行中的每个元素与矩阵 B 相乘。最后绑定列
外循环(map):通过矩阵A的行迭代内循环。最后绑定行
您的代码的问题在于 A[i,j]*B
是与 B
具有相同维度的矩阵,并且您想将其放入单单元格 (i,j)
中。
您可以将代码修改为如下所示:
kron <- function(A,B){
m <- matrix(0,nrow = nrow(A)*nrow(B),ncol = ncol(A)*ncol(B))
for (i in seq_len(nrow(A))){
cur_rows <- seq_len(nrow(B)) + (i-1)*nrow(B)
for (j in seq_len(ncol(A)))
m[cur_rows,seq_len(ncol(B)) + (j-1)*ncol(B)] <- A[i,j]*B
}
m
}
以上是大佬教程为你收集整理的Kronecker 积函数定义全部内容,希望文章能够帮你解决Kronecker 积函数定义所遇到的程序开发问题。
如果觉得大佬教程网站内容还不错,欢迎将大佬教程推荐给程序员好友。
本图文内容来源于网友网络收集整理提供,作为学习参考使用,版权属于原作者。
如您有任何意见或建议可联系处理。小编QQ:384754419,请注明来意。