0

我对h这个实验中的值有点困惑。在cpp中,

int h,J=3,n=200,p=3,h_m=(n+p+1)/2;
float rt=(float)h_m/n;

for(int j=0,j<J,j++){
   h=floor((j+1)/J*rt*(n-p-1)+p+1);
   std::cout<<h<<" "<<j<<" "<<rt<<" "<<n-p-1<<" "<<h_i<<" "<<J<<std::endl;
}

给出:

4 0 0.51 196 7  3
4 1 0.51 196 7  3
103 2 0.51 196 7  3

我试图得到(这是在 R 中):

n<-200
p<-3
h_m<-as.integer((n+p+1)/2)
J<-3
j<-0:(J-1)
rt<-h_m/n
floor((j+1)/J*rt*(n-p-1)+p+1)
[1]  37  70 103

cpp 构造有什么问题?

4

1 回答 1

5

(j+1)/J

是整数除法,因此被截断。明确地将操作数之一转换为浮点数:

(j+1)/(float)J
于 2012-08-02T11:40:47.813 回答