add: phase shift corrected with average
This commit is contained in:
parent
9e64a106af
commit
1ab07d335f
9
main.cpp
9
main.cpp
|
@ -203,12 +203,15 @@ int main() {
|
|||
}
|
||||
|
||||
cout << "[INFO] Starting the calculate the phase shift" << endl;
|
||||
double delta = 0;
|
||||
double deltas;
|
||||
int ndeltas = 0;
|
||||
for(int i = 0; i < channelsZeros[0].size(); i++) {
|
||||
if(channelsZeros[1].size() <= i) continue;
|
||||
cout << "[INFO] Phase shift for zero " << i << " is " << ((channelsXData[0][channelsZeros[0][i]] - channelsXData[1][channelsZeros[1][i]]) / channelOneT) * 2 * M_PI << endl;
|
||||
delta += channelsZeros[1][i] - channelsZeros[0][i];
|
||||
// cout << "[INFO] Phase shift for zero " << i << " is " << ((channelsXData[0][channelsZeros[0][i]] - channelsXData[1][channelsZeros[1][i]]) / channelOneT) * 360 << endl;
|
||||
deltas += channelsXData[1][channelsZeros[1][i]] - channelsXData[0][channelsZeros[0][i]];
|
||||
ndeltas++;
|
||||
}
|
||||
cout << "[INFO] Average phase shift (CHA0 to CHA1) " << ((deltas/ndeltas) / channelOneT) * 360 << endl;
|
||||
|
||||
plt::plot(channelsXData[0], channelsYData[0], channelsXData[1], channelsYData[1]);
|
||||
plt::xlabel("X Label");
|
||||
|
|
Reference in New Issue