r - 在 R 中编写适应函数 - 输出计算中的错误
问题描述
我有两个用 MATLAB 编写的函数(distitpf
以及distchpf
来自语音箱工具箱),我正在尝试将它们转换为 R。这些函数计算 Itakura 距离和 COSH 距离。下面是这些函数的代码以及我迄今为止在 R 中编写的代码。我的问题是当我测试我的 R 函数时出现错误。任何想法是什么错误以及如何解决它们?
MATLAB 中的函数:
function d=distitpf(pf1,pf2)
[nf1,p2]=size(pf1);
p1=p2-1;
nf2=size(pf2,1);
nx=min(nf1,nf2);
r=pf1(1:nx,:)./pf2(1:nx,:);
q=log(r);
d=log((sum(r(:,2:p1),2)+0.5*(r(:,1)+r(:,p2)))/p1)-(sum(q(:,2:p1),2)+0.5*(q(:,1)+q(:,p2)))/p1;
end
function d=distchpf(pf1,pf2)
[nf1,p2]=size(pf1);
p1=p2-1;
nf2=size(pf2,1);
nx=min(nf1,nf2);
r=pf1(1:nx,:)./pf2(1:nx,:);
q=r+r.^(-1);
d=(2*sum(q(:,2:p1),2)+q(:,1)+q(:,p2))/(4*p1)-1;
end
R中的功能:
distitpf <- function(pf1,pf2){
nf1=nrow(pf1)
p2=length(pf1);
p1=p2-1;
nf2=nrow(pf2);
nx=min(nf1,nf2);
r=pf1[1:nx,]/pf2[1:nx,];
qx=log(r);
d=log((sum(r[,2:p1],2)+0.5*(r[,1]+r[,p2]))/p1)-(sum(qx[,2:p1],2)+0.5*(qx[,1]+qx[,p2]))/p1;
return(d)}
distchpf <- function(pf1,pf2) {
nf1=nrow(pf1)
p2=length(pf1);
p1=p2-1;
nf2=nrow(pf2);
nx=min(nf1,nf2);
r=pf1[1:nx,]/pf2[1:nx,];
qx=r+r^(-1);
d=(2*sum(qx[,2:p1],2)+qx[,1]+qx[,p2])/(4*p1)-1;
return(d)}
测试 MATLAB:
pf1=[2,4,5,6];
pf2=[7,8,9,1];
>> distitpf(pf1,pf2)
ans =
0.6732
>> distchpf(pf1,pf2)
ans =
0.6386
测试 R:
pf1=c(2,4,5,6);
pf2=c(7,8,9,1);
> distchpf(pf1,pf2)
Error in 1:nx : result would be too long a vector
In addition: Warning message:
In min(nf1, nf2) : no non-missing arguments to min; returning Inf
> distitpf(pf1,pf2)
Error in 1:nx : result would be too long a vector
In addition: Warning message:
In min(nf1, nf2) : no non-missing arguments to min; returning Inf
(编辑)我已经解决了这个问题。我在下面发布了更正的 R 代码。
distitpf <- function(pf1,pf2){
pf1<-as.matrix(pf1)
pf2<-as.matrix(pf2)
nf1=nrow(pf1)
p2=length(pf1);
p1=p2-1;
nf2=nrow(pf2);
nx=min(nf1,nf2);
r=pf1/pf2
qx=log(r);
d=log((sum(r[2:p1])+0.5*(r[1]+r[p2]))/p1)-(sum(qx[2:p1])+0.5*(qx[1]+qx[p2]))/p1;
return(d)}
distchpf <- function(pf1,pf2) {
pf1<-as.matrix(pf1)
pf2<-as.matrix(pf2)
nf1=nrow(pf1)
p2=length(pf1);
p1=p2-1;
nf2=nrow(pf2);
r=pf1/pf2;
qx=r+r^(-1);
d=(2*sum(qx[2:p1])+qx[1]+qx[p2])/(4*p1)-1;
return(d)}
解决方案
答案是菲利帕。但是我会在seeway
和proxy
库中使用已经实现的函数:
distitpf <- function(pf1, pf2) {
pf1 <- as.matrix(pf1)
pf2 <- as.matrix(pf2)
nf1 <- nrow(pf1)
p2 <- length(pf1)
p1 <- p2 - 1
nf2 <- nrow(pf2)
nx <- min(nf1, nf2)
r <- pf1 / pf2
qx <- log(r)
log((sum(r[2:p1]) + 0.5 * (r[1] + r[p2])) / p1) - (sum(qx[2:p1]) + 0.5 * (qx[1] + qx[p2])) / p1
}
distchpf <- function(pf1, pf2) {
pf1 <- as.matrix(pf1)
pf2 <- as.matrix(pf2)
nf1 <- nrow(pf1)
p2 <- length(pf1)
p1 <- p2 - 1
nf2 <- nrow(pf2)
r <- pf1 / pf2
qx <- r + 1 / r
(2 * sum(qx[2:p1]) + qx[1] + qx[p2]) / (4 * p1) - 1
}
推荐阅读
- python - 如何在 python 3.6.9 中将 dbf 转换为 csv?
- javascript - jQuery - 检查目标是否是特定的 div
- android - 在将 Android 应用程序集成到 Firebase 远程配置时,如何使用 2 个 google-services.json(一个用于 PROODUCTION,另一个用于 QA 环境)?
- python-3.x - 使用 Python 和 PyDrive 在谷歌驱动器上成功上传文件,但文件已损坏
- string - 如何连接 Rust Vector 的两个 `&str` 类型元素并添加到零索引
- smt - smt2 中的类型不匹配
- c# - 当集合为空时,UWP XAML 引发未处理的异常
- php - 在 MySql 中使用 JSON+PHP 插入 \ 字符
- node.js - 如何在 node.js 中使用 jsconfig 设置路径别名
- java - 如何从 Firebase 获取特定键的特定子值