From 1e157e25e191295493d375a7fbfe7f46341b974b Mon Sep 17 00:00:00 2001
From: Thomas Lambert <t.lambert@uliege.be>
Date: Fri, 9 Dec 2022 00:03:15 +0100
Subject: [PATCH] feat(analysis): add comare frequencies function
---
.gitignore | 1 +
utils/comparefrequencies.m | 44 ++++++++++++++++++++++++++++++++++++++
2 files changed, 45 insertions(+)
create mode 100644 utils/comparefrequencies.m
diff --git a/.gitignore b/.gitignore
index 89bb900..1f51cd2 100644
--- a/.gitignore
+++ b/.gitignore
@@ -4,6 +4,7 @@
*.eps
# Folders to ignore
+figures/
data/wtFiles/*/*
data/arduinoFiles/*
data/ORIGINAL_RAW_DATA
diff --git a/utils/comparefrequencies.m b/utils/comparefrequencies.m
new file mode 100644
index 0000000..f672197
--- /dev/null
+++ b/utils/comparefrequencies.m
@@ -0,0 +1,44 @@
+function [allRelDiff, meanDiff, stdDiff] = comparefrequencies(ResData)
+ % COMAPREFREQUENCIES Compare the frequencies of the two sets (should be identical)
+ % Todo
+
+ % ----------------------------------------------------------------------------------------------
+ % (c) Copyright 2022 University of Liege
+ % Author: Thomas Lambert <t.lambert@uliege.be>
+ % ULiege - Aeroelasticity and Experimental Aerodynamics
+ % MIT License
+ %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+ allRelDiff = zeros(size(ResData));
+
+ for i = 1:size(ResData, 1)
+ for j = 1:size(ResData, 2)
+
+ allFreqF(i, j) = ResData(i, j).Front.trueFreq;
+ allFreqA(i, j) = ResData(i, j).Aft.trueFreq;
+
+ allRelDiff(i, j) = abs(ResData(i, j).Front.trueFreq - ResData(i, j).Aft.trueFreq) / ...
+ ResData(i, j).Front.trueFreq;
+
+ end
+ end
+
+ allRelDiff = abs(allFreqF - allFreqA) ./ allFreqF;
+
+ meanDiff = mean(mean(allRelDiff));
+ stdDiff = mean(std(allRelDiff));
+
+ figure('Name', 'Frequency differences');
+ hold on;
+ plot(allFreqF, allFreqA, 'ok', 'MarkerFaceColor', 'k');
+ xlim([1.5, 3.5]);
+ ylim(xlim);
+ plot(xlim, xlim, '--');
+
+ hold off;
+ grid on;
+
+ xlabel('f_{aft} [Hz]');
+ ylabel('f_{front} [Hz]');
+
+end
--
GitLab