求螺旋桨推力系数(KT2.m)

【核心代码】

function kt = KT2(J, PD, AEAO)
	conn = database("Propdb","aaa","123456");
    cursor1 = exec(conn, "select * from MAU4_40KT");
    result = fetch(cursor1);
    temp = result.Data;
    mau4_40_pd_value =cell2mat(temp(:,5));
    mau4_40_pd_value=mau4_40_pd_value(end:-1:1);
    MAU4_40KT =cell2mat(temp(:,7:2:end));
    MAU4_40KT(find(isnan(MAU4_40KT)==1))= 0;
    close(cursor1);
    cursor2 = exec(conn, "select * from MAU4_55KT");
    result = fetch(cursor2);
    temp = result.Data;
    mau4_55_pd_value =cell2mat(temp(:,5));
    mau4_55_pd_value=mau4_55_pd_value(end:-1:1);
    MAU4_55KT =cell2mat(temp(:,7:2:end));
    MAU4_55KT(find(isnan(MAU4_55KT)==1))= 0;
    close(cursor2);   
    cursor3 = exec(conn, "select * from MAU4_70KT");
    result = fetch(cursor3);
    temp = result.Data;
    mau4_70_pd_value =cell2mat(temp(:,5));
    mau4_70_pd_value=mau4_70_pd_value(end:-1:1);
    MAU4_70KT =cell2mat(temp(:,7:2:end));
    MAU4_70KT(find(isnan(MAU4_70KT)==1))= 0;
    close(cursor3);  
    close(conn);		
    [c40,r40] = size(MAU4_40KT);
    z1 = [];
    z2 = [];
    z3 = [];
    for i=1:c40
        z1(c40 1-i,:) = MAU4_40KT(i,:);
    end
    [c55,r55] = size(MAU4_55KT);
    for i=1:c55
        
        z2(c55 1-i,:) = MAU4_55KT(i,:);
    end
    [c70,r70] = size(MAU4_70KT);
    for i=1:c70
        
        z3(c55 1-i,:) = MAU4_70KT(i,:);
    end
    
    x_kt40 = 0:0.1:0.1*(r40-1);
    Kt40=[];
    for i=1:c40
        Kt40(i) = interp1(x_kt40,z1(i,:),J);
    end
    KT40_res = interp1(mau4_40_pd_value,Kt40,PD);
    x_kt55 = 0:0.1:0.1*(r55-1);
    Kt55=[];
    for i=1:c55
        Kt55(i) = interp1(x_kt55,z2(i,:),J);
    end
    KT55_res = interp1(mau4_55_pd_value,Kt55,PD);
    x_kt70 = 0:0.1:0.1*(r70-1);
    Kt70=[];
    for i=1:c70
        Kt70(i) = interp1(x_kt70,z3(i,:),J);
    end
    KT70_res = interp1(mau4_70_pd_value,Kt70,PD);
    AE = [0.4,0.55,0.7];
    tempres = interp1(AE,[KT40_res,KT55_res,KT70_res],AEAO);

    kt = tempres;
 
    
end