本文整理汇总了C++中corr函数的典型用法代码示例。如果您正苦于以下问题:C++ corr函数的具体用法?C++ corr怎么用?C++ corr使用的例子?那么恭喜您, 这里精选的函数代码示例或许可以为您提供帮助。
在下文中一共展示了corr函数的20个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于我们的系统推荐出更棒的C++代码示例。
示例1: r
int HrPwWindow::fit(QVector<double> &Xi,QVector<double> &Yi,int n)
{
QVector<double> r(5),lnXi(100),lnYi(100),logXi(100),logYi(100),invXi(100);
//Linear correlation
r[0]=val_abs(corr(Xi,Yi,n));
//Exponential correlation
lnarray(Yi,lnYi,n);
r[1]=val_abs(corr(Xi,lnYi,n));
//Power correlation
logarray(Xi,logXi,n);
logarray(Yi,logYi,n);
r[2]=val_abs(corr(logXi,logYi,n));
//Inverse correlation
invarray(Xi,invXi,n);
r[3]=val_abs(corr(invXi,Yi,n));
//Logarithmic correlation
lnarray(Xi,lnXi,n);
r[4]=val_abs(corr(lnXi,Yi,n));
//Best fit test
return rmax(r);
}
开发者ID:GoldenCheetah,项目名称:GoldenCheetah,代码行数:33,代码来源:HrPwWindow.cpp
示例2: r
int
Statistic::ajustement(QVector<double> &Xi,QVector<double> &Yi,int n)
{
QVector<double> r(5),lnXi(100),lnYi(100),logXi(100),logYi(100),invXi(100);
//corrélation pour linéaire
r[0]=val_abs(corr(Xi,Yi,n));
//corrélation pour exponetielle
lntab(Yi,lnYi,n);
r[1]=val_abs(corr(Xi,lnYi,n));
//corrélation pour puissance
logtab(Xi,logXi,n);
logtab(Yi,logYi,n);
r[2]=val_abs(corr(logXi,logYi,n));
//corrélation pour inverse
invtab(Xi,invXi,n);
r[3]=val_abs(corr(invXi,Yi,n));
//corrélation pour logarithmique
lntab(Xi,lnXi,n);
r[4]=val_abs(corr(lnXi,Yi,n));
//Test du meilleur ajustement
return rmax(r);
}
开发者ID:AartGoossens,项目名称:GoldenCheetah,代码行数:34,代码来源:Statistic.cpp
示例3: NCC
double NCC(const Image<float>& I1,const Image<float>& meanI1,Point m1,const Image<float>& I2,const Image<float>& meanI2,Point m2,int n) {
if (m1.x<n || m1.x>=I1.width()-n || m1.y<n || m1.y>=I1.height()-n) return -1;
if (m2.x<n || m2.x>=I2.width()-n || m2.y<n || m2.y>=I2.height()-n) return -1;
double c1=corr(I1,meanI1,m1,I1,meanI1,m1,n);
if (c1==0) return -1;
double c2=corr(I2,meanI2,m2,I2,meanI2,m2,n);
if (c2==0) return -1;
return corr(I1,meanI1,m1,I2,meanI2,m2,n)/sqrt(c1*c2);
}
开发者ID:momolenghun,项目名称:disparity,代码行数:9,代码来源:image.cpp
示例4: NCC
double NCC(const Image<float>& I1,const Image<float>& meanI1,const Image<float>& corrI1,Point m1,const Image<float>& I2,const Image<float>& meanI2,const Image<float>& corrI2,Point m2,int n) {
if (m1.x<n || m1.x>=I1.width()-n || m1.y<n || m1.y>=I1.height()-n) return -1;
if (m2.x<n || m2.x>=I2.width()-n || m2.y<n || m2.y>=I2.height()-n) return -1;
if (corrI1(m1)==0) return -1;
if (corrI2(m2)==0) return -1;
return corr(I1,m1,I2,m2,n)/sqrt(corrI1(m1)*corrI2(m2));
}
开发者ID:KevinBienvenu,项目名称:ImageProcessingProjects,代码行数:7,代码来源:image.cpp
示例5: delayHr
int
HrPwWindow::findDelay(QVector<double> &wattsArray, QVector<double> &hrArray, int rideTimeSecs)
{
int delay = 0;
double maxr = 0;
if (rideTimeSecs>= 60) {
for (int a = 10; a <=60; ++a) {
QVector<double> delayHr(rideTimeSecs);
for (int j = a; j<rideTimeSecs; ++j) {
delayHr[j-a] = hrArray[j];
}
for (int j = rideTimeSecs-a; j<rideTimeSecs; ++j) {
delayHr[j] = 0.0;
}
double r = corr(wattsArray, delayHr, rideTimeSecs-a);
//fprintf(stderr, "findDelay %d: %.2f \n", a, r);
if (r>maxr) {
maxr = r;
delay = a;
}
}
}
delayEdit->setText(QString("%1").arg(delay));
rDelayEdit->setText(QString("%1").arg(delay));
delaySlider->setValue(delay);
rDelaySlider->setValue(delay);
return delay;
}
开发者ID:Gerryboy1964,项目名称:GoldenCheetah,代码行数:35,代码来源:HrPwWindow.cpp
示例6: bmfm_fancy
void bmfm_fancy(float *disp, // output disparities image (dx, dy)
float *errc, // output error image
float *a, // input image A
float *b, // input image B
int w, // width
int h, // height
int pd, // pixel dimension
double fm[9], // fundamental matrix
float *disp_init, // initialization (optional)
float *search_radius // optional, w.r.t. initialization
)
{
int maxpoints = 2 * (w+h), (*p)[2] = xmalloc(maxpoints*sizeof*p);
for (int j = 0; j < h; j++)
for (int i = 0; i < w; i++) {
float rad = NAN, ini[2] = {0, 0};
if (search_radius) rad = search_radius[j*w+i];
if (disp_init) ini[0] = disp_init[2*(j*w+i)+0];
if (disp_init) ini[1] = disp_init[2*(j*w+i)+1];
int np = plot_epipolar_fancy(p, fm, w, h, i, j, ini, rad);
float mincorr = INFINITY;
int minidx = 0;
for (int k = 0; k < np; k++) {
float c = corr(a, b, w, h, pd, i, j, p[k][0], p[k][1]);
if (c < mincorr) {
mincorr = c;
minidx = k;
}
}
if (disp) disp[2*(j*w+i) + 0] = p[minidx][0] - i;
if (disp) disp[2*(j*w+i) + 1] = p[minidx][1] - j;
if (errc) errc[j*w+i] = mincorr;
}
free(p);
}
开发者ID:Fahdben,项目名称:imscript,代码行数:35,代码来源:bmfm_fancier.c
示例7: neigh_indices
void shot_detector::compare(pcl::PointCloud<DescriptorType>::Ptr model_descriptions, pcl::PointCloud<DescriptorType>::Ptr scene_descriptions)
{
model_scene_corrs->clear();
pcl::KdTreeFLANN<DescriptorType> match_search;
match_search.setInputCloud (model_descriptions);
std::cerr << scene_descriptions->size() << " and " << model_descriptions->size() << std::endl;
model_good_keypoints_indices.clear();
scene_good_keypoints_indices.clear();
// For each scene keypoint descriptor, find nearest neighbor into the model keypoints descriptor cloud and add it to the correspondences vector.
for (size_t i = 0; i < scene_descriptions->size (); ++i) {
std::vector<int> neigh_indices (1);
std::vector<float> neigh_sqr_dists (1);
if (!pcl_isfinite (scene_descriptors->at (i).descriptor[0])) { //skipping NaNs
continue;
}
int found_neighs = match_search.nearestKSearch (scene_descriptors->at (i), 1, neigh_indices, neigh_sqr_dists);
if(found_neighs == 1 && neigh_sqr_dists[0] < corr_dist_) { // add match only if the squared descriptor distance is less than 0.25 (SHOT descriptor distances are between 0 and 1 by design)
pcl::Correspondence corr (neigh_indices[0], static_cast<int> (i), neigh_sqr_dists[0]);
model_scene_corrs->push_back (corr);
}
}
pcl::copyPointCloud (*model_keypoints, model_good_keypoints_indices, *model_good_kp);
pcl::copyPointCloud (*scene_keypoints, scene_good_keypoints_indices, *scene_good_kp);
std::cerr << "Correspondences found: " << model_scene_corrs->size () << std::endl;
}
开发者ID:ehuang3,项目名称:apc_ros,代码行数:27,代码来源:shot_detector.cpp
示例8: operator
double operator()(configuration const& c) const {
auto delta = (rhs - kern(c)) / error_bars;
int M = first_dim(delta);
double kappa = 0;
for(int i = 1; i < M; ++i) kappa += corr(delta(i), delta(i-1));
kappa /= M - 1;
return kappa;
}
开发者ID:krivenko,项目名称:som,代码行数:8,代码来源:fit_quality.hpp
示例9: corrImage
// Renvoit la matrice avec toutes les auto-correlations
Image<float> corrImage(const Image<float>& I, const Image<float>& meanI, int n){
Image<float> corrI(I.width(), I.height(),CV_32F);
for(int i=n+1; i<I.height()-n-1; i++){
for(int j=n+1; j<I.width()-n-1; j++){
//cout<< i << " " << j << endl;
corrI.at<float>(i,j) = (float)(corr(I, meanI, Point(j,i), I, meanI, Point(j,i),n));
}
}
return corrI;
}
开发者ID:KevinBienvenu,项目名称:ImageProcessingProjects,代码行数:11,代码来源:image.cpp
示例10: InitialGuess
int InitialGuess (double *s, double *p, int nphase, int nchn, int *chn)
{
int index;
double frac_off = 0.05; // set to be 0.05
double ptemp[nphase];
double temp[nchn];
int i, h;
for (i = 0; i < nchn; i++)
{
for (h = 0; h < nphase; h++)
{
ptemp[h] = p[i*nphase+h];
}
int x;
x = def_off_pulse (nphase, ptemp, frac_off);
double ptemp_out[nphase];
pre_diff (ptemp, nphase, x, frac_off, ptemp_out);
temp[i] = find_peak_value(nphase,ptemp_out);
}
int peak;
find_peak (0, nchn, temp, &peak);
//printf ("%d\n",peak);
(*chn) = peak;
double p_use[nphase];
double s_use[nphase];
for (h = 0; h < nphase; h++)
{
p_use[h] = p[peak*nphase+h];
s_use[h] = s[peak*nphase+h];
}
// remove the baseline of template
index = def_off_pulse (nphase, s_use, frac_off);
double s_out[nphase];
pre_diff (s_use, nphase, index, frac_off, s_out);
// remove the baseline of profile
index = def_off_pulse (nphase, p_use, frac_off);
double p_out[nphase];
pre_diff (p_use, nphase, index, frac_off, p_out);
// Guess the phase shift
int d;
d = corr (s_out, p_out, nphase);
return d;
}
开发者ID:shidai,项目名称:ptimeD,代码行数:54,代码来源:fdt.c
示例11: main
int main()
{
int a, i, j;
gets(string);
a = 1;
i = 0;
for(j = 0; string[j] != '\0'; j++);
j--;
while(j - i >= 1)
{
while(corr(string[i])) i++;
while(corr(string[j]) && (j >= 1)) j--;
if ((reg(string[i]) != reg(string[j])) && (j - i >= 1)) a = 0;
i++;
j--;
}
if (a) printf("Yes");
else printf("No");
return 0;
}
开发者ID:Lozov-Petr,项目名称:programming-semester-1,代码行数:21,代码来源:Palindrom.c
示例12: roguejoin
static void roguejoin(struct level *lev, int x1, int y1, int x2, int y2, int horiz)
{
int x,y,middle;
#ifndef MAX
#define MAX(a,b) (((a) > (b)) ? (a) : (b))
#endif
#ifndef MIN
#define MIN(a,b) (((a) < (b)) ? (a) : (b))
#endif
if (horiz) {
middle = x1 + rn2(x2-x1+1);
for (x=MIN(x1,middle); x<=MAX(x1,middle); x++)
corr(lev, x, y1);
for (y=MIN(y1,y2); y<=MAX(y1,y2); y++)
corr(lev, middle,y);
for (x=MIN(middle,x2); x<=MAX(middle,x2); x++)
corr(lev, x, y2);
} else {
middle = y1 + rn2(y2-y1+1);
for (y=MIN(y1,middle); y<=MAX(y1,middle); y++)
corr(lev, x1, y);
for (x=MIN(x1,x2); x<=MAX(x1,x2); x++)
corr(lev, x, middle);
for (y=MIN(middle,y2); y<=MAX(middle,y2); y++)
corr(lev, x2,y);
}
}
开发者ID:mbi,项目名称:NitroHack,代码行数:27,代码来源:extralev.c
示例13: fit
int fit() {
cout << "start of fit procedure" << endl;
readsource("run369/CRS/CRS_QDC_01!qdc01.dat",7214);
readbackgr("run371/CRS/CRS_QDC_01!qdc01.dat",56948);
corr();
TCanvas* c1 = new TCanvas("c1","c1",800,600);
int maxbin=-1;
Double_t maxvalue=-1;
Double_t tempbin=0;
TString text="";
hcorr->GetXaxis()->SetRange(200,4000);
maxbin=hcorr->GetMaximumBin();
maxvalue=hcorr->GetBinContent(maxbin);
hcorr->GetXaxis()->SetRange(0,4095);
TF1* f1 = new TF1("f1","[0]*exp(-((x-[1])^2)/[2])");
f1->SetParameter(0,maxvalue);
f1->SetParameter(1,maxbin);
f1->SetParameter(2,1);
f1->SetLineColor(2);
hcorr->SetMinimum(0);
hcorr->SetMaximum(1.2*maxvalue);
hcorr->Fit("f1","","",maxbin-maxbin/6,4095);
tempbin=sqrt(-log(maxvalue*2/3/(f1->GetParameter(0)))*(f1->GetParameter(2)))+f1->GetParameter(1);
cout << "2/3 of maximum height is in bin " << tempbin << "->" <<int(tempbin+0.5) << endl;
cout << endl << "max bin = " << maxbin << " with content = " << maxvalue << endl;
}
开发者ID:phystory,项目名称:newStandalone,代码行数:48,代码来源:fit.C
示例14: spatial_mat
static void
spatial_mat(double *par, double *dist, longint *n, longint *nug,
double (*corr)(double ), double *mat)
{
longint i, j, np1 = *n + 1;
double aux, *sdist, ratio = 1.0;
sdist = dist;
if (*nug) ratio = par[1];
for(i = 0; i < *n; i++) {
mat[i * np1] = 1.0;
for(j = i + 1; j < *n; j++, sdist++) {
aux = *sdist / *par;
*(mat + i + j * (*n)) = *(mat + j + i * (*n)) = ratio * corr(aux);
}
}
}
开发者ID:csilles,项目名称:cxxr,代码行数:16,代码来源:corStruct.c
示例15: corr
std::vector<Joystick::CalibrationData>
Joystick::get_calibration()
{
std::vector<struct js_corr> corr(get_axis_count());
if (ioctl(fd, JSIOCGCORR, &*corr.begin()) < 0)
{
std::ostringstream str;
str << filename << ": " << strerror(errno);
throw std::runtime_error(str.str());
}
else
{
std::vector<CalibrationData> data;
std::transform(corr.begin(), corr.end(), std::back_inserter(data), corr2cal);
return data;
}
}
开发者ID:dino0815,项目名称:jstest-gtk,代码行数:18,代码来源:joystick.cpp
示例16: insert_cmd
void AddEntryDlg::OnButtonEditInsertPressed()
{
//QString insert_cmd("UPDATE 'musicdb'.'main' SET 'titel' = '");
QString insert_cmd("UPDATE 'main' SET 'titel' = '");
insert_cmd.append(corr(le_title.text()));
insert_cmd.append("', 'kuenstler' = '");
insert_cmd.append(corr(le_artist.text()));
insert_cmd.append("', 'album' = '");
insert_cmd.append(corr(le_album.text()));
insert_cmd.append("', 'tag' = '");
insert_cmd.append(corr(le_tag.text()));
insert_cmd.append("', 'genre' = '");
insert_cmd.append(corr(le_genre.text()));
insert_cmd.append("', 'jahr' = '");
insert_cmd.append(le_year.text());
insert_cmd.append("', 'others' = '");
if(cb_interest_others.isChecked())
insert_cmd.append("1");
else
insert_cmd.append("0");
insert_cmd.append("', 'yours' = '");
if(cb_interest_yours.isChecked())
insert_cmd.append("1");
else
insert_cmd.append("0");
insert_cmd.append("', 'dateityp' = '");
insert_cmd.append(le_filetype.text());
insert_cmd.append("', 'qualitaet' = '");
insert_cmd.append(le_quality.text());
insert_cmd.append("', 'bew_yours' = '");
insert_cmd.append(QString::number(sb_vote_yours.value()));
insert_cmd.append("', 'bew_others' = '");
insert_cmd.append(QString::number(sb_vote_others.value()));
insert_cmd.append("', 'pfad' = '");
insert_cmd.append(le_path.text());
insert_cmd.append("', 'url' = '");
insert_cmd.append(corr(le_source.text()));
insert_cmd.append("' WHERE 'main'.'id' =");
insert_cmd.append(QString::number(editnum));
printf("Query command: %s\n", insert_cmd.toAscii().data());
sqlhelper.exec(insert_cmd);
accept();
}
开发者ID:JohannesLorenz,项目名称:LoQuality,代码行数:51,代码来源:AddEntries.cpp
示例17: main
int main(){
int i,j,st;
int sp = n();
double data[500000];
double theta[5] = {0.3,0.2,0.3,0,0};
double phi[5] = {1.,0.3,0.7,0,0};
double nofphot[5] = {1000.,1000.,1000,0,0};
st = corr(data,3,theta,phi,0.25,50,5,nofphot,0.,50);
FILE* f = fopen("T.txt","w+");
for(i=0;i<256;i++){
for(j=0;j<256;j++){
fprintf(f,"%f ",data[sp*256*256+j*256+i]);
}
fprintf(f,"\n");
}
fclose(f);
return 0;
}
开发者ID:pssncp142,项目名称:dtu-space-rmc,代码行数:19,代码来源:corr.c
示例18: corr
void DSPEngine::dcOffset(SampleVector::iterator begin, SampleVector::iterator end)
{
double count;
int io = 0;
int qo = 0;
Sample corr((qint16)m_iOffset, (qint16)m_qOffset);
// sum and correct in one pass
for(SampleVector::iterator it = begin; it < end; it++) {
io += it->real();
qo += it->imag();
*it -= corr;
}
// moving average
count = end - begin;
m_iOffset = (15.0 * m_iOffset + (double)io / count) / 16.0;
m_qOffset = (15.0 * m_qOffset + (double)qo / count) / 16.0;
}
开发者ID:zoutepopcorn,项目名称:rtl-sdrangelove,代码行数:19,代码来源:dspengine.cpp
示例19: density
doublereal WaterPropsIAPWS::psat(doublereal temperature, int waterState) {
doublereal densLiq = -1.0, densGas = -1.0, delGRT = 0.0;
doublereal dp, pcorr;
if (temperature >= T_c) {
densGas = density(temperature, P_c, WATER_SUPERCRIT);
setState_TR(temperature, densGas);
return P_c;
}
doublereal p = psat_est(temperature);
bool conv = false;
for (int i = 0; i < 30; i++) {
if (method == 1) {
corr(temperature, p, densLiq, densGas, delGRT);
doublereal delV = M_water * (1.0/densLiq - 1.0/densGas);
dp = - delGRT * Rgas * temperature / delV;
} else {
corr1(temperature, p, densLiq, densGas, pcorr);
dp = pcorr - p;
}
p += dp;
if ((method == 1) && delGRT < 1.0E-8) {
conv = true;
break;
} else {
if (fabs(dp/p) < 1.0E-9) {
conv = true;
break;
}
}
}
// Put the fluid in the desired end condition
if (waterState == WATER_LIQUID) {
setState_TR(temperature, densLiq);
} else if (waterState == WATER_GAS) {
setState_TR(temperature, densGas);
} else {
throw Cantera::CanteraError("WaterPropsIAPWS::psat",
"unknown water state input: " + Cantera::int2str(waterState));
}
return p;
}
开发者ID:hkmoffat,项目名称:cantera,代码行数:42,代码来源:WaterPropsIAPWS.cpp
示例20: go_appendix_backgr
int go_appendix_backgr(){
readsource("run370/CRS/CRS_QDC_00!qdc00.dat",7215);
readbackgr("run371/CRS/CRS_QDC_00!qdc00.dat",56948);
corr();
TPaveText *pt = new TPaveText(0.65,0.7,0.85,0.85, "NDC");
pt->AddText("run371");
pt->AddText("background");
pt->AddText("box I (QDC ch00)");
hbackgr->GetXaxis()->SetTitle("QDC energy [bin]");
hbackgr->SetTitle("Background - box I");
TCanvas* c1=new TCanvas("c1","c1",800,300);
hbackgr->Draw();
pt->Draw("same");
c1->SaveAs("appendix/backgr_boxI.jpg");
}
开发者ID:phystory,项目名称:newStandalone,代码行数:20,代码来源:run370-Cs137-boxI.C
注:本文中的corr函数示例由纯净天空整理自Github/MSDocs等源码及文档管理平台,相关代码片段筛选自各路编程大神贡献的开源项目,源码版权归原作者所有,传播和使用请参考对应项目的License;未经允许,请勿转载。 |
请发表评论