气动参数的拟合

气动参数一般有两种使用方法，插值和拟合，拟合计算较快，写程序简单。matlab有自带函数regress可以进行多元参数最小二乘拟合，可以拟合出想要的形式，得到每一项的系数。

clc;clear;close;
alpha1=xlsread('SR72数据.xlsx','sheet1','C4:C19');
CD1=xlsread('SR72数据.xlsx','sheet1','D4:D19');
plot(alpha1,CD1)
alpha2=xlsread('SR72数据.xlsx','sheet1','E4:E17');
CD2=xlsread('SR72数据.xlsx','sheet1','F4:F17');
hold on
plot(alpha2,CD2)
alpha3=xlsread('SR72数据.xlsx','sheet1','G4:G17');
CD3=xlsread('SR72数据.xlsx','sheet1','H4:H17');
hold on
plot(alpha3,CD3)
hold off
%% regress函数调用
Ma=[0.3*ones(1,length(alpha1)) 0.4*ones(1,length(alpha2)) 0.5*ones(1,length(alpha3))]';
alpha=[alpha1;alpha2;alpha3];
CD=[CD1;CD2;CD3];
X=[ones(size(CD)) Ma.^2 alpha.^2 Ma alpha Ma.*alpha];
[b,bint,r,rint,stats]=regress(CD,X);
% CD=b(2)Ma^2+b(3)alpha^2+b(4)Ma+b(5)alpha+b(6)Ma*alpha+b(1);
Ma=0:0.1:0.5;
alpha=0:30;
[Ma alpha]=meshgrid(Ma,alpha);
CD=b(2)*Ma.^2+b(3)*alpha.^2+b(4)*Ma+b(5)*alpha+b(6)*Ma.*alpha+b(1);
mesh(Ma,alpha,CD)
hold on
plot3(0.3*ones(1,length(alpha1)),alpha1,CD1,'o')
plot3(0.4*ones(1,length(alpha2)),alpha2,CD2,'*')
plot3(0.5*ones(1,length(alpha3)),alpha3,CD3,'s')
xlabel('Ma')
ylabel('alpha')
zlabel('CD')