MATLAB批量处理生成profili生成的翼型数据

话不多说，从profili光滑处理后导出的翼型数据存在NACA4415.txt里。

load NACA4415.txt
YX(136:261,:)=flipud(YX(136:261,:));  //手动删除文字行，便于三维软件绘制，将后半部分数据倒置处理
YX=[YX zeros(261,1)]//加入z轴坐标

YXZ=YX*8.5; //调整为自己需要的弦长
//批量处理..
chord0=max(YXZ(:,1))-min(YXZ(:,1));
chord=chord0;
for i=1:9
    y = chord0-(100-sqrt(100^2-(10*i)^2))-(200-200*sqrt(1-i^2/27^2));
    chord = [chord y];
end
for i=10:23
    y=578.83-(10*i-95.76)*(578.83/174.24)+200*sqrt(1-i^2/27^2);
    chord = [chord y];
end
//调整z坐标
for i = 1:24
    %str1 = ['NACA4415S',num2str(i+1),'.txt'];
    YXS(:,:,i) = YXZ * chord(i) / chord0;
    %NACA4415S(:,1,i) = NACA4415S(:,1,i) + x(i);
    YXS(:,3,i) = 10 * (i-1) * ones(261,1);
    %q = NACA4415S(:,:,i);
    %save(str1,'q','-ascii');
end
//批量处理
for i=1:10
    YXS(:,1,i) = YXS(:,1,i) + 100 - sqrt(100^2 - (10*(i-1))^2);
end
for i=11:24
    YXS(:,1,i) = YXS(:,1,i) + 71.17+(10*(i-1)-95.76)*(578.83/174.24);
end
//批量保存
for i = 1:24
    str1 = ['NACA4415S',num2str(i),'.txt'];
%     NACA4415S(:,:,i) = NACA4415 * chord(i) / chord0;
%     NACA4415S(:,1,i) = NACA4415S(:,1,i) + x(i);
%     NACA4415S(:,3,i) = 10 * (i-1) * ones(261,1);
    q = YXS(:,:,i);
    save(str1,'q','-ascii');
end