导弹总体设计《防空导弹主级燃料相对质量因数计算》大作业.

N*s/kgР#define g 9.801 //重力加速度g,单位:m/s^2Р#define P 2.2 //推重比PР#define p0 5880.0 //翼载p0,单位:N/m^2Р#define vt 420.0 //飞机速度vt,单位:m/sР#define yt0 15000.0 //飞机初始纵坐标yt,单位:mР#define dt0 34200.0 //弹目斜距Dt0,单位:mР#define xt0 31569.39345 //飞机初始橫坐标xt0,单位:m xt0=sqrt(dt0^2-(yt0-y0)^2)+x0Р#define cotq0 2.10462623//初始目标对制导站的航向角的余切,cotq0=xt0/yt0;Рdouble inter_linear(double x0,double x1,double y0,double y1,double x){//线性插值函数Р double a0,a1,y;Р a0=(x-x1)/(x0-x1);Р a1=(x-x0)/(x1-x0);Р y=a0*y0+a1*y1;Р return y;Р}Рdouble interp1(double x[],double y[],int n,double a){//一维线性插值函数Р double z;Р int i=0;Р for(i=0;i<(n-1);i++){Р if((a>=x[i])&&(a<=x[i+1])) break;//确定a的位置Р }Р z = inter_linear(x[i],x[i+1],y[i],y[i+1],a);Р return z;Р}Рdouble interp2(double x[],double y[],double z[6][6],int m,int n,double a,double b){//双线性插值函数