twointerp.m

用matlab来编写的较复杂的程序

%二维插值函数，输入变量依次为：攻角数组，马赫数数组，待插参数数组，插值点攻角值，插值点马赫数值
%调用了一维插值函数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));