matlab POD处理

admin

% 本文件为主函数，不执行具体的计算，仅调用各子函数完成需要的功能
clear all;
clc;
row = 25096;                                              % 每个文件共10772
col = 134;                                                % 共104 组模拟数据
Velocity_U = zeros(row,col);                           
Velocity_V = zeros(row,col);   
Velocity_W = zeros(row,col);   
path = '.\Fluent\dawo_';   % POD模态路径

filename=sprintf('%s%d.txt',path,1);
CFD_data = importdata(filename);                                  % 加载脉动数据
CFD_X= CFD_data.data(:,1:4);                                      % 存储 坐标值
for i=1:col
    filename=sprintf('%s%d.txt',path,i);
    CFD_data= importdata(filename);                             % 加载脉动数据
    Velocity_U(:,i)=CFD_data.data(:,6);                         % 提取速度 U
    Velocity_V(:,i)=CFD_data.data(:,5);                         % 提取速度 V
    Velocity_W(:,i)=CFD_data.data(:,4);                         % 提取速度 W
end

[U_M1, U_Phi_POD, U_lam_POD, U_xmean, U_Coeff_POD]=POD2(Velocity_U);
[V_M1, V_Phi_POD, V_lam_POD, V_xmean, V_Coeff_POD]=POD2(Velocity_V);
[W_M1, W_Phi_POD, W_lam_POD, W_xmean, W_Coeff_POD]=POD2(Velocity_W);
POD_Data_chuli2(CFD_X,U_Phi_POD,U_Coeff_POD,V_Phi_POD,V_Coeff_POD,W_Phi_POD,W_Coeff_POD,1,134);  % 输出四阶模态速度矢量
POD_Data_chuli2(CFD_X,U_Phi_POD,U_Coeff_POD,V_Phi_POD,V_Coeff_POD,W_Phi_POD,W_Coeff_POD,2,134);  % 输出二阶模态速度矢量
POD_Data_chuli2(CFD_X,U_Phi_POD,U_Coeff_POD,V_Phi_POD,V_Coeff_POD,W_Phi_POD,W_Coeff_POD,3,134);  % 输出三阶模态速度矢量
POD_Data_chuli2(CFD_X,U_Phi_POD,U_Coeff_POD,V_Phi_POD,V_Coeff_POD,W_Phi_POD,W_Coeff_POD,4,134);  % 输出四阶模态速度矢量

POD_Data_chuli3(CFD_X,U_Phi_POD,U_Coeff_POD,V_Phi_POD,V_Coeff_POD,W_Phi_POD,W_Coeff_POD,4,134);

fprintf('全部结束 !!!\n');