Coeficiente De Sustenta O

%DETERMINAÇÃO DA LINHA DE CÂMBER
%Os vetores e1, e2, i1, e i2 foram importados do arquivo NACA63212.txt
%correspondendo, respectivamente, a A[1:26],B[1:26], A[27:52] e B[27:52],
%representando as linhas de extradorso e intradorso do perfil camber=mean([e2 i2]')'; z=polyfit(e1,camber,6); der=polyder(z);
%MÉTODO ANALÍTICO
%transformação da derivada de z em uma função simbólica x=sym('x'); f(x)=der(1)*x^5+der(2)*x^4+der(3)*x^3+der(4)*x^2+der(5)*x+der(6);
%transformação da variável x em theta theta=sym('theta'); g(theta)=(1-cos(theta))/2; h(theta)=compose(f,g); %cálculo dos coeficientes A0, A1 e A2
%cálculo de Cl para gerar a curva Cl x alfa analítica
A1=2*int(h(theta)*cos(theta),0,pi)/pi;
A2=2*int(h(theta)*cos(2*theta),0,pi)/pi; alfa(1)=degtorad(-8); A0=zeros(17,1);
Cl=zeros(17,1);
for i=1:17 if i~=1 alfa(i)=alfa(i-1)+degtorad(1); end
A0(i)=alfa(i)-int(h(theta),0,pi)/pi;
Cl(i)=pi*(2*A0(i)+A1); end alfa_g=radtodeg(alfa); plot(alfa_g,Cl,'*') %plota a curva Cl x alfa analítica hold on
%segura o gráfico para acrescentar as outras curvas posteriormente
%cálculo do Cm a 1/4 da corda
Cm_c4a=pi*(A2-A1)/4
%cálculo do Cl e CP para um ângulo de 3 graus alfa=degtorad(3); A0=alfa-int(h(theta),0,pi)/pi;
Cl_3a=pi*(2*A0+A1)
CP_3a=(1+pi*(A1-A2)/Cl_3a)/4
%cálculo do Cl e CP para um ângulo de 4 graus alfa=degtorad(4); A0=alfa-int(h(theta),0,pi)/pi;
Cl_4a=pi*(2*A0+A1)
CP_4a=(1+pi*(A1-A2)/Cl_4a)/4
%MÉTODO NUMÉRICO (VORTEX LATTICE)
%cálculo do número N de elementos para o qual o método converge
%avaliado para alfa igual a 5 graus alfa=degtorad(5); %definição arbitrária de Cl para iniciar o comando while, com N inicial igual a 2
Cl_anterior=1;
Cl=100;
N=2;
%crítério de convergência: diferença percentual entre Cl e Cl_anterior menor que 0.2% while abs(Cl-Cl_anterior)/Cl_anterior>0.002
%definição dos vetores C e V
C=zeros(N,1);
V=zeros(N,1);
C(1)=3/(4*N);
V(1)=1/(4*N); for i=2:N
C(i)=C(i-1)+1/N;
V(i)=V(i-1)+1/N; end %definição da matriz A
A=zeros(N,N);
for k=1:N for j=1:N