matlab调用xfoil程序计算翼型升阻力力矩

matlab调用xfoil程序，计算翼型升阻力力矩，将结果生成文件，程序对三种翼型进行了计算，用户可根据需要删减或复制相同部分即可

delete（‘input-E186.dat.dat‘）;
delete（‘output-E186.dat‘）;
delete（‘input-E334.dat.dat‘）;
delete（‘output-E334.dat‘）;
delete（‘input-E637.dat.dat‘）;
delete（‘output-E637.dat‘）;

fid=fopen（‘input-E186.dat‘‘wt‘）;
fprintf（fid‘LOAD \n‘）;
fprintf（fid‘E186.DAT \n‘）;
fprintf（fid‘PANE \n‘）;
fprintf（fid‘OPER \n‘）;
fprintf（fid‘VISC %12d \n‘2500000）;
fprintf（fid‘M %12.3f \n‘0.11）;
fprintf（fid‘ITER %12.3f \n‘50）;
fprintf（fid‘PACC \n‘）;
fprintf（fid‘output-E186.dat \n‘）;
fprintf（fid‘ \n‘）;
fprintf（fid‘ASEQ %12.3f %12.3f %12.3f \n‘ [1 18 1]）;
%fprintf（fid‘A %f \n‘2）;
fprintf（fid‘PACC \n‘）;
fprintf（fid‘CPWR \n‘）;
fprintf（fid‘E186cp.txt \n‘）;
fprintf（fid‘PACC \n‘）;
fprintf（fid‘QUIT \n‘）;
%fprintf（fid‘polar.txt \n‘）;
fclose（fid）;
!xfoilclc
[ABCDEFG]=textread（‘output-E186.dat‘‘%12.5f %12.5f   %12.5f  %12.5f %12.5f %12.5f %12.5f‘‘delimiter‘‘‘‘headerlines‘12）;
display（‘    迎角        CL         CD     CL/CD      CDp       CM      TOP_Xtr     BOT_Xtr‘）
H1=[ABCB./CDEFG];
J1=unique（H1‘rows‘）  %返回不同行组成的矩阵

fid=fopen（‘input-E334.dat‘‘wt‘）;
fprintf（fid‘LOAD \n‘）;
fprintf（fid‘E334.dat \n‘）;
fprintf（fid‘PANE \n‘）;
fprintf（fid‘OPER \n‘）;
fprintf（fid‘VISC %12d \n‘2500000）;
fprintf（fid‘M %12.3f \n‘0.11）;
fprintf（fid‘ITER %12.3f \n‘50）;
fprintf（fid‘PACC \n‘）;
fprintf（fid‘output-E334.dat \n‘）;
fprintf（fid‘ \n‘）;
fprintf（fid‘ASEQ %12.3f %12.3f %12.3f \n‘ [1 18 1]）;
%fprintf（fid‘A %f \n‘2）;
fprintf（fid‘PACC \n‘）;
fprintf（fid‘CPWR \n‘）;
fprintf（fid‘E334cp.txt \n‘）;
fprintf（fid‘PACC \n‘）;
fprintf（fid‘QUIT \n‘）;
%fprintf（fid‘polar.txt \n‘）;
fclose（fid）;
!xfoilclc
[ABCDEFG]=textread（‘output-E334.dat‘‘%12.5f %12.5f   %12.5f  %12.5f %12.5f %12.5f %12.5f‘‘delimiter‘‘‘‘headerlines‘12）;
display（‘    迎角        CL         CD     CL/CD      CDp       CM      TOP_Xtr     BOT_Xtr‘）
H2=[ABCB./CDEFG];
J2=unique（H2‘rows‘）  %返回不同行组成的矩阵


%figure（1）
subplot（221）
plot（AJ1（:2）‘r‘AJ2（:2）‘g‘）
legend（‘E186‘‘E334‘）
title（‘CL-a‘）;

%figure（2）
subplot（222）
plot（AJ1（:3）‘r‘AJ2（:3）‘g‘）
legend（‘E186‘‘E334‘）
title（‘CD-a‘）;

%figure（3）
subplot（223）
plot（AJ1（:6）‘r‘AJ2（:6）‘g‘）
legend（‘E186‘‘E334‘）
title（‘CM-a‘）;

subplot（224）
plot（AJ1（:4）‘r‘AJ2（:4）‘g‘）
legend（‘E186‘‘E334‘）
title（‘L/D-a‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1412  2012-05-28 15:16  matlab-xfoil\E186.DAT

     文件       3380  2012-05-28 15:23  matlab-xfoil\E186cp.txt

     文件       1735  2012-05-28 15:17  matlab-xfoil\E334.dat

     文件       3380  2012-05-28 15:23  matlab-xfoil\E334cp.txt

     文件       1351  2012-05-14 22:44  matlab-xfoil\E637.dat

     文件        204  2012-05-28 15:23  matlab-xfoil\input-E186.dat

     文件        204  2012-05-28 15:23  matlab-xfoil\input-E334.dat

     文件       1619  2012-05-28 15:23  matlab-xfoil\output-E186.dat

     文件       1619  2012-05-28 15:23  matlab-xfoil\output-E334.dat

     文件        431  2012-05-28 15:23  matlab-xfoil\QUIT

     文件       2432  2012-05-28 14:13  matlab-xfoil\xfoil.asv

     文件    1339392  2006-05-02 22:10  matlab-xfoil\xfoil.exe

     文件       2434  2012-05-28 15:23  matlab-xfoil\xfoil.m

     目录          0  2014-08-20 22:02  matlab-xfoil

----------- ---------  ---------- -----  ----

              1359593                    14