資源簡介
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);
!xfoil clc
[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);
!xfoil clc
[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
- 上一篇:LTE系統中頻率調度算法程序
- 下一篇:PSO_dg_prog
評論
共有 條評論