twointerp.m

来自「采用Matlab编写的弹道分析算法」· M 代码 · 共 22 行

%二维插值函数，输入变量依次为：攻角数组，马赫数数组，待插参数数组，插值点攻角值，插值点马赫数值
%调用了一维插值函数oneinterp2

function F = twointerp(arg1,arg2,arg3,arg4,arg5)
[m n]=size(arg3);
%确定要用到哪二行马赫数
if arg5<=arg2(1)%马赫数小于或者等于最小给定马赫数 
    i=1;
end    

if arg5>arg2(m)%马赫数大于最大给定马赫数
    i=m-1;
else
     i=1;
     while arg5>arg2(i+1)
          i=i+1;
     end    
end     

fy1=oneinterp2(arg1,arg3,arg4,i);
fy2=oneinterp2(arg1,arg3,arg4,i+1);
F=fy1+(fy1-fy2)*(arg5-arg2(i))/(arg2(i)-arg2(i+1));