%%     Code written by WBX 2022/08/16
% 对比仅用FAST计算所得的浮式和固定式风机的响应
%
clear all
clc
% Step1:读取.out文件，注意只读取"-1"文件(Span1to9叶片段1-9)即可
% 1)指定文件名 Filename
OpenFASTout_Floating='Test24_1.out';
OpenFASTout_Fixed   ='Test24_1_fixed.out';

% 2-1) 将浮式风机叶片段1-9所在的FAST结果文件数据存入矩阵   Floating
file_Floating = fopen(OpenFASTout_Floating,'r');
Output_Floating=textscan(file_Floating, '%f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f', 'headerlines',8);fclose(fileID1);
Floating=cell2mat(Output_Floating);               %使用cell2mat函数将元胞数组Output_Floating转化为矩阵，便于后续运算
% 2-2) 将固定式风机叶片段1-9所在的FAST结果文件数据存入矩阵 Fixed
file_Fixed = fopen(OpenFASTout_Fixed,'r');
Output_Fixed=textscan(file_Fixed, '%f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f', 'headerlines',8);fclose(fileID1);
Fixed=cell2mat(Output_Fixed);               %使用cell2mat函数将元胞数组Output_Fixed转化为矩阵，便于后续运算

% 3) 绘图对比结果 *
% 3-1) 6-DOF 运动 

