1.龙格库塔法1）建立old_y数组储存最开始的。2）不考虑t,直接加步长后进行右函数迭代计算。3）设定下有多少需要积分的。2.插值法（1维，2维，3维）1）1维找到位置后，进行一维插值计算2）2维中的某一维找到位置，进行插值计算，再调用一维插值。3）三维找位置计算，调2维。

参考《导弹飞行力学》对部分参数的解释:dx/dt=f(t,x):   之前一直看不懂f（t,x）到底指的哪个式子，其实在开头就提出来了，f是dy  K2=△t*f(tk+△t/2,xk+1/2*K1)：t在导弹飞行力学，是y(0)(M中是1)，所以可忽略。然后xk+1/2*K1即dx/dt的变量进行变换，在飞行力学没有直接加，需要通过y来改变dy而y的变换则是借助y+h*dy去改变

问：1.Matlab读数据  （C读数据）         2.超出边界如何处理                自：按旁边斜率处理    众：？见： 3.减少重复程序interp,单独拉出来        4.减少输入量，矩阵长度自己计算测试案例：clc;clear;%%%一维插值测试A1=[0;1;2;3;4;5];A2=[1.2;2.4;3.6;4.8;

1.JB变量2.导引控制q_old_y = atan((pos_Target[1] - y[8]) / (pos_Target[0] - y[7])); q_old_z = atan((pos_Target[2] - y[9]) / (pos_Target[0] - y[7]));if (flag4 == 1 && y[0] > T5) { q_new_

1.Wz主要由攻角、舵角和滚动阻尼项有关（Wy则是侧滑角）。Wx则主要是舵和滚转阻尼项。2.注意：计算的时候，是弧度制；插值的时候，有可能是角度制。3.对于弹簧刀，速度不影响气动参数。