%% This read me is best read using Notepad++ %%

Dear user,

This repository contains the data used for the writing of the manuscript titled "Wake Mixing Control for Floating Wind Farms" which is submitted for a special edition of the Control Systems Magazine. In the manuscript the interaction between the Helix wake mixing method and IEA 15MW turbine mounted on the VolturnUS platform is analyzed. The investigation is conducted using two different methods. The first is a linear frequency analysis. By performing and identification experiment linear input-output relations can be identified. The data for these experiments is contained in the folder 'Frequency_Identification'. The second method is an analysis through full non-linear simulations using QBlade as simulation suite. With this method the downstream wind speed is analyzed. The data for these experiments is contained in the folder 'Simulation Results'. A third folder named 'Functions' is also included and contains matlab functions required for the system identification. 

%% Frequency_Identification %%
Within the folder 'Frequency_Identification' a subfolder 'Identification_Data' can be found as well as a matlab script and function. When run, the Matlab script will import the data contained in the subfolder (line 34-68 in the code) and run identification functions that process the data. This happens in line 79 tot 87, for example the line:

[G_Surge,sys_Surge,~] = id_funcs(detrend(squeeze(PitchYaw)),detrend(squeeze(Turbine1_data(:,12))),dt,4,6) ;

uses the matlab function id_funcs to calculate the input-output spectrum between the blade pitch signals and the floating turbine surge motion. This spectrum is contained in the sys_Surge variable. G_surge contains a state space model that is fitted to this spectrum. The last three inputs are the sampling time, resampling and model order. The last two are only required for model fitting. The contents of the id_funcs are taken from: https://www.dcsc.tudelft.nl/~jwvanwingerden/pbsid/pbsidtoolbox_product_page.html. After identification has been done (this can take some time), the same figures as those shown in the manuscript are created.

%% Simulation Data %%
The data from the simulations carried out in QBlade can be found in 'Simulation Results' which contains two sub folders ('PostProcessing_2_degrees' and 'PostProcessing_4_degrees'), imported data as .mat files and two matlab scripts. The two subfolder contain the raw data, divided in turbine data and the wind speed data. The Matlab script 'ImportData.m' imports the raw data and saves it to the .mat files also found in the same folder. The script 'ProcessData' opens the .mat files and plots the figures as seen in the manuscript. The function crameri.m together with the CrameriColourMaps7.0.mat is required to get the same colourmaps. These colourmaps are designed to best represent the data, for more information see: https://www.fabiocrameri.ch/colourmaps/.

All of the scripts were tested in MATLAB R2022b.