面元法计算机身的气动参数的matlab程序代码.docVIP

面元法计算超音速机身的气动参数的编程实验 机身参数如下图所示： 面元法计算机身的气动参数的matlab程序代码（仅供参考） function [] = surface(a,b,D,n) P=0; Q=0; R=0; L=D*4.2; A=16.5*pi/180;B=20*pi/180; S=pi*(L/4*tan(A)+D/2)^2 for i=1:11 X(i)=-a*(i-1)/10; r1=sqrt(1-(X(i)+a)*(X(i)+a)/(a*a))*b; for j=1:(n+1) Y(i,j)=r1*cos(2*pi*(j-