matlab插值法原理,五孔探针的原理和插值法数据处理源代码

clear

xydata=load('Calib-1.dat'); % 载入标定数据

xydata1=load('Meas-1.dat'); % 载入试验数据

fid=fopen('Pout-1.dat','w'); % 输出文件

ma=xydata(:,1)';

beta=xydata(:,2)';

alpha=xydata(:,3)';

pt=xydata(:,4)';

p1=xydata(:,6)';

p2=xydata(:,7)';

p3=xydata(:,8)';

p4=xydata(:,9)';

p5=xydata(:,10)';

BB=p1-(p2 p3 p4 p5)/4;

C_alpha=(p2-p3)./BB;

C_beta=(p4-p5)./BB;

C_ma=(p1-BB)./p1;

C_pt=(pt-p1)./BB;

p1_t=xydata1(:,1)';

p2_t=xydata1(:,2)';

p3_t=xydata1(:,3)';

p4_t=xydata1(:,4)';

p5_t=xydata1(:,5)';

BB_t=p1_t-(p2_t p3_t p4_t p5_t)/4;

C_alpha_t=(p2_t-p3_t)./BB_t;

C_beta_t=(p4_t-p5_t)./BB_t;

C_ma_t=(p1_t-BB_t)./p1_t;

% 插值

alpha_out=griddata3(C_alpha,C_beta,C_ma,alpha,C_alpha_t,C_beta_t,C_ma_t);

beta_out=griddata3(C_alpha,C_beta,C_ma,beta,C_alpha_t,C_beta_t,C_ma_t);

ma_out=griddata3(C_alpha,C_beta,C_ma,ma,C_alpha_t,C_beta_t,C_ma_t);

C_pt_out=griddata3(C_alpha,C_beta,C_ma,C_pt,C_alpha_t,C_beta_t,C_ma_t);

pt_out=C_pt_out.*BB_t p1_t;

fprintf(fid,' alpha beta Ma Pt\n');

for i=1:length(ma_out)

fprintf(fid,'%8.2f %8.2f %8.3f %9.3f\n',alpha_out(i), beta_out(i),ma_out(i),pt_out(i));

end

fclose(fid);