NACA翼型Matlab问题

我在生成NACA翼型图时遇到了一个障碍。代码如下：

% NACA 6414%
clc
clear
format compact

m1=input('Please enter one digit number for maxi camber line:')
p1=input('Please enter one digit number for distance of m and chordline:')
t1=input('Please enter two digit number for thickness of airfoil:')

m=m1*0.01;
p=p1*0.1;
t=t1*0.01;

x1=[0:0.01:p];
x2=[p:0.01:1];
x=[x1,x2]


yc1=(m*(2*p*x1-x1.^2))/(p^2);
yc2=(m/(1-p)^2)*((1-2*p)+2*p*x2-x2.^2);
yt=(t/0.2)*(0.2969*sqrt(x)-0.1260*x-0.3516*x.^2+0.2843*x.^3-0.1015*x.^4)

yc=[yc1,yc2]

ang1=(atan((m/p^2)*(2*p-2*x1)))*180/pi
ang2=(atan((m/(1-p)^2)*(2*p-2*x2)))*180/pi
ang=[ang1,ang2]

xu=x-yt.*sin(ang)
yu=yc+yt.*cos(ang)
xl=x+yt.*sin(ang)
yl=yc-yt.*cos(ang)

xplot=[xu,xl]
yplot=[yu,yl]

plot(xplot,yplot)

它没有生成一个正常的翼型图，而是带有一个非常奇怪的图，图上有许多圆圈。有人能帮我解决这个问题吗

---

谢谢

我想问题出在这里：

 ang1=(atan((m/p^2)*(2*p-2*x1)))*180/pi
 ang2=(atan((m/(1-p)^2)*(2*p-2*x2)))*180/pi
 ang=[ang1,ang2]

 xu=x-yt.*sin(ang)
 yu=yc+yt.*cos(ang)
 xl=x+yt.*sin(ang)
 yl=yc-yt.*cos(ang)

函数sin和cos期望输入弧度。atan的输出也以弧度为单位，但你要将其转换为度，这意味着xu、yu等中的输出，你稍后绘制的将是错误的。尽量避免转换，除非您确实需要以度为单位的角度。或者，您可以使用这些函数的其他形式（sind、cosd、atand等），并始终保持所有函数的度数