所以我正在尝试用C ++中的MATLAB做我的一个项目,但是我一直困扰着。
这是我想在MATLAB中转换为C ++的代码部分。它适用于MATLAB但不适用于C ++
RelRough = [0, 1E-6, 5E-6, 1E-5, 5E-5, 0.0001, 0.0002, 0.0004, 0.0006, 0.0008, 0.001];
ReT = [4000, 5000, 6000, 7000, 8000, 9000, 10000, 20000, 30000, 40000, 50000, 60000, 70000, 80000, 90000, 100000, 200000, 300000, 400000, 500000];
for i = 1:length(ReT)
for j = 1:length(RelRough)
FCT_guess = 1;
tolerance = 1;
while tolerance > 1e-14
FCT_cal = 1/(-2*log10((RelRough(j)/3.7) + (2.51/(ReT(i)*sqrt(FCT_guess)))))^2;
tolerance = abs(FCT_cal-FCT_guess);
FCT_guess = FCT_cal;
FCT(i,j) = FCT_cal*1000;
end
end
end
这是我的C ++版本,我不断收到变量g的错误,例如“表达必须有整数或未整合的枚举类型”
double RelRough[] = { 0, 1E-6, 5E-6, 1E-5, 5E-5, 0.0001, 0.0002, 0.0004, 0.0006, 0.0008, 0.001 };
const int lengthRelRough = sizeof(RelRough) / sizeof(RelRough[0]);
double ReT[] = { 4000, 5000, 6000, 7000, 8000, 9000, 10000, 20000, 30000, 40000, 50000, 60000, 70000, 80000, 90000, 100000, 200000, 300000, 400000, 500000 };
const int lengthReT = sizeof(ReT) / sizeof(ReT[0]);
double fct[lengthReT][lengthRelRough] = { 0 };
double fct_guess = 1;
double tolerance = 1;
double fct_cal = 0;
for (int ii = 0; ii < lengthReT; ++ii) {
for (int jj = 0; jj < lengthRelRough; ++jj) {
while (tolerance > 1e-14) {
double h = (RelRough[jj] / 3.7), w = (2.51 / (ReT[ii] * sqrt(fct_guess)));
double g = (-2*log10(h+w))^2;
fct_cal = 1/g;
tolerance = abs(fct_cal - fct_guess);
fct_guess = fct_cal;
fct[ii][jj] = fct_cal;
std::cout << fct[ii][jj] << "\t";
}
}
}
return 0;
}
有没有人帮助看出什么是错的。提前致谢!
改变这个:
double g = (-2*log10(h+w))^2;
成:
double g = pow(-2*log10(h+w),2.0);
正如@Eljay在他的评论中指出的那样,运营商^
在XOR
中执行C++
,而不是取幂。欲获得更多信息: