该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

总有问题，现附上信息：

N=100;

%线段份数

c=1;

%假设的弦长

x=0:c/N:c;

%x的取值范围

m=2/100;p=3/10;t=10/100;

%NACA2310的所包含的翼型数据信息

yc=zeros(1,N+1);

angle_of_attack=0;

U=10;

ux=U*cos(angle_of_attack);

uy=U*sin(angle_of_attack);

%这是攻角angle_of_attack和来流速度U

%ux,uy是水平来流u∞和竖直来流v∞，根据攻角求得

for i=1:N+1

if x(i)<=p*c

yc(i)=m*x(i)/p^2*(2*p-x(i)/c);

else

yc(i)=m*(c-x(i))/(1-p)^2*(1-2*p+x(i)/c);

end

end

%计算出翼型的中心线

%plot(x,yc);

%axis equal;

%xlabel('x');

%ylabel('y');

%上面是检验中心线是否正确

yt=zeros(1,N+1);

yt_dill=zeros(1,N+1);

th=zeros(1,N+1);

for i=1:N+1

yt(i) =t*c*(0.2969*(x(i)/c)^0.5 -0.126*(x(i)/c)^1 -0.3516*(x(i)/c)^2 +0.2843*(x(i)/c)^3 -0.101036*(x(i)/c)^4)/0.2;

yt_dill(i)=t*c*(0.2969*0.5*(x(i)/c)^(-0.5)/c-0.126*(x(i)/c)^0/c-0.3516*2*(x(i)/c)^1/c+0.2843*3*(x(i)/c)^2/c-0.101036*4*(x(i)/c)^3)/0.2/c;

th(i)=atan(yt_dill(i));

end

%计算翼厚yt，yt_dill是翼厚导数

%plot(yt,x); %检验yt的信息

figure;

plot(x,th);

for i=1:N+1

xu(i)=x(i)-yt(i)*sin(th(i));

xl(i)=x(i)+yt(i)*sin(th(i));

yu(i)=yc(i)+yt(i)*cos(th(i));

yl(i)=yc(i)-yt(i)*cos(th(i));

end

%上下翼型的函数

for I=1:N

XU(I)=(xu(I)+xu(I+1))/2;

YU(I)=(yu(I)+yu(I+1))/2;

%控制点为上边界时的中心点

XL(I)=(xl(I)+xl(I+1))/2;

YL(I)=(yl(I)+yl(I+1))/2;

%控制点为下边界时的中心的 %因为下面的积分的控制点是段的重点，所以控制点和段数的数量是相同的 %N个上控制点，N个下控制点

end

%定义上下边界的中点YU,YL

%————————————————上面是基础信息，不会错的

% plot(x,yc,'r',xu,yu,'b',xl,yl,'b');

% axis equal;

% xlabel('x');

% ylabel('y'); %画出上下边界线

% plot(x,yc,'r',XU,YU,'b',XL,YL,'b');

% axis equal;

% xlabel('x');

% ylabel('y'); %检验上下重点控制点的正确性

shuzu_xishu=zeros(2*N+1,2*N+1); %定义了系数数组，共2N+1，2n+1项

for i=1:N;

feiqicixiang(i)=-ux*YU(i)+uy*XU(i);

feiqicixiang(N+i)=-ux*YL(i)+uy*XL(i);

end

feiqicixiang(2*N+1)=0; %非齐次项的系数凑齐了，2N+1个

n=99;

%n是积分过程中用的总次数，就是把xu(i)到xl(i+1)中间分成99分，共100个点

for I=1:N %对上边界控制点的积分系数

for i=1:N %对控制点积分的点的移动

%上边界对上控制点的积分

ang=atan((yu(i+1)-yu(i)) / (xu(i+1)-xu(i)));

kuandu=(xu(i+1)-xu(i))/n;

x=xu(i) : kuandu : xu(i+1);

jifen_u=0;

for k=1:n+1

y(k)=log(((x(k)-XU(I))^2+((yu(i)+(x(k)-xu(i))*tan(ang))-YU(I))^2)^(1/2))/cos(ang); %积分的内容,yu(i)=yu(i)+(x-xu(i))*tan(ang)

jifen_u=y(k)*kuandu+jifen_u;

end

shuzu_xishu(I,i)=jifen_u/(2*pi);

%下边界对上控制点的积分

ang=atan((yl(i+1)-yl(i)) / (xl(i+1)-xl(i)));

kuandu=(xl(i+1)-xl(i))/n;

x=xl(i) : kuandu : xl(i+1);

jifen_u=0;

for k=1:n+1

y(k)=log(((x(k)-XU(I))^2+((yl(i)+(x(k)-xl(i))*tan(ang))-YU(I))^2)^(1/2))/cos(ang); %积分的内容,yu(i)=yu(i)+(x-xu(i))*tan(ang)

jifen_u=y(k)*kuandu+jifen_u;

end

shuzu_xishu(I,i+N)=1/(2*pi)*jifen_u; %对上控制点N的积分已经好了，现在赋值到系数数组里面

end

end

for I=1:N %对于下边界控制点的积分,I是下控制点的移动

for i=1:N %i是对下控制点积分的点的移动

%上边界对下控制点的积分

ang=atan((yu(i+1)-yu(i)) / (xu(i+1)-xu(i)));

kuandu=(xu(i+1)-xu(i))/n;

x=xu(i) : kuandu : xu(i+1);

jifen_l =0;

for k=1:n+1

y(k)=log(((x(k)-XL(I))^2+((yu(i)+(x(k)-xu(i))*tan(ang))-YL(I))^2)^(1/2))/cos(ang); %积分的内容,yu(i)=yu(i)+(x-xu(i))*tan(ang)

jifen_l=y(k)*kuandu+jifen_l;

end

shuzu_xishu(I+N,i)=1/(2*pi)*jifen_l;

%下边界对下控制点的积分

ang=atan((yl(i+1)-yl(i)) / (xl(i+1)-xl(i)));

kuandu=(xl(i+1)-xl(i))/n;

x=xl(i) : kuandu : xl(i+1);

jifen_l =0;

for k=1:n+1

y(k)=log(((x(k)-XL(I))^2+((yl(i)+(x(k)-xl(i))*tan(ang))-YL(I))^2)^(1/2))/cos(ang); %积分的内容,yu(i)=yu(i)+(x-xu(i))*tan(ang)

jifen_l=y(k)*kuandu+jifen_l;

end

shuzu_xishu(I+N,i+N)=1/(2*pi)*jifen_l; %对下控制点N的积分已经好了，现在赋值到系数数组里面

end

end