Fayorg
/
gymnyo-phy-tp8
Archived
1
0
Fork 0
Code

add: phase shift

This commit is contained in:
Elie Baier 2024-02-25 21:25:13 +01:00
parent 3951adadfb
commit 9e64a106af
1 changed files with 13 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
main.cpp
View File

@ -110,6 +110,7 @@ int main() {
vector<vector<double>> channelsXData, channelsYData, channelsZeros; vector<vector<double>> channelsXData, channelsYData, channelsZeros;
vector<double> channelsP2P; vector<double> channelsP2P;
double channelOneT = 0;
for(int channel = 0; channel < channels.size(); channel++) { for(int channel = 0; channel < channels.size(); channel++) {
string line; string line;
@ -170,6 +171,10 @@ int main() {
avgFrequency *= 2; avgFrequency *= 2;
cout << "[INFO] Average peak time is " << avgFrequency << "s" << endl; cout << "[INFO] Average peak time is " << avgFrequency << "s" << endl;
cout << "[INFO] Average frequency is " << 1/avgFrequency << "Hz" << endl; cout << "[INFO] Average frequency is " << 1/avgFrequency << "Hz" << endl;
if(channel == 0) {
channelOneT = avgFrequency;
cout << "Set chan 1 freq: " << channelOneT << endl;
}
vector<double> mins, maxs; vector<double> mins, maxs;
for(int i = 1; i < channelsZeros[channel].size(); i++) { for(int i = 1; i < channelsZeros[channel].size(); i++) {
@ -197,6 +202,14 @@ int main() {
channelsP2P.push_back(abs(avgMin) + abs(avgMax)); channelsP2P.push_back(abs(avgMin) + abs(avgMax));
} }
cout << "[INFO] Starting the calculate the phase shift" << endl;
double delta = 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];
}
plt::plot(channelsXData[0], channelsYData[0], channelsXData[1], channelsYData[1]); plt::plot(channelsXData[0], channelsYData[0], channelsXData[1], channelsYData[1]);
plt::xlabel("X Label"); plt::xlabel("X Label");
plt::ylabel("Y Label"); plt::ylabel("Y Label");