涡格法计算气动力MATLAB

涡格法计算，适用范围0<=ma<=0.5

%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