%Cy-alpha%
clear
clc
Cy=[];
mz=[];
XXX=[];
l0=input（‘請輸入翼根弦長‘）
ll=input（‘請輸入翼梢弦長‘）
L=input（‘請輸入展長‘）
ka=input（‘請輸入前緣后掠角（以弧度計量）‘）
tic
ka=ka*100000;
ka=round（ka）;
ka=ka/100000;
for?ii=0:1:10
alpha=ii;
S=（ll+l0）*0.5*L;
bA=2/S*（L/6*（l0-ll）^2+L/2*ll*（l0-ll）+L/2*ll^2）;
u=20;
[QXZX1X2Z1Z2W]=fenkuai（l0llLkau）;
%列方程組%
[Tqi]=qiujie（XZX1X2Z1Z2u）;
%以下部分求解方程組%
????v=100*sin（alpha/180*pi）;
????V（1:2*u*u）=v*4*pi;
????H=inv（T）*V‘;
????P（u2*u）=0;
????for?i=1:1:u
????????for?j=1:1:2*u
????????????P（ij）=H（2*u*（i-1）+j）;
????????end
????end
????%以下部分計算答案%
????zi（1:u1:2*u）=L/2/u;
????F=1.112*100*sum（sum（P.*zi））;
????M=1.112*100*sum（sum（P.*zi.*W））;
????Mz=1.112*100*sum（sum（P.*zi.*W））/（0.5*1.112*S*bA*100^2）;
????Cy=[Cy?F/（0.5*1.112*S*100^2）];
????mz=[mz?M/（0.5*1.112*S*bA*100^2）];
????XXX=[XXX?Mz/（F/（0.5*1.112*S*100^2））];
????ii
end
K=[];
for?i=1:1:10
????K（i）=（Cy（i+1）-Cy（i））;
end
K=K/pi*180;
i=1:1:10;
plot（iK）
grid?on
xlabel（‘迎角‘）
ylabel（‘升力線斜率‘）
toc

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1073??2014-10-24?10:20??渦格法程序包\Cy_alpha.m

?????文件???????2267??2014-10-22?23:51??渦格法程序包\fenkuai.m

?????文件????????974??2014-10-24?09:43??渦格法程序包\main.m

?????文件???????1352??2014-10-24?09:02??渦格法程序包\qiujie.m

?????文件????????261??2014-10-24?10:29??渦格法程序包\程序說明.txt

?????目錄??????????0??2014-11-22?16:20??渦格法程序包

-----------?---------??----------?-----??----

?????????????????5927????????????????????6