逆发动机模型

niuju=xlsread('\throttle.xls','B2:L18');
zhuansu=xlsread('\throttle.xls','A2:A18');
jieqimen=xlsread('\throttle.xls','B1:L1');
%zhuansu
speed_engine=[];
speed_row=[];
    for i=1:length(zhuansu)
        speed_engine(i,:)=ones(length(jieqimen),1)*zhuansu(i);
    end
    for i=1:length(zhuansu)
        for j=1:length(jieqimen)
          speed_row=[speed_row,speed_engine(i,j)];
        end
    end
 speed_row=speed_row';
  %jieqimen
  jieqimen_table=[]; jieqimen_row=[];
  for i=1:length(zhuansu)
       jieqimen_table(i,:)=jieqimen;
  end
   for i=1:length(zhuansu)
        for j=1:length(jieqimen)
          jieqimen_row=[jieqimen_row,jieqimen_table(i,j)];
        end
   end
     jieqimen_row= jieqimen_row';
     %zhuanju
     zhuanju_row=[]
      for i=1:length(zhuansu)
        for j=1:length(jieqimen)
          zhuanju_row=[zhuanju_row,niuju(i,j)];
        end
      end
   zhuanju_row=zhuanju_row';
   F=scatteredInterpolant(zhuanju_row,speed_row,jieqimen_row);
   su=min(speed_row):100:max(speed_row);
   T=min(zhuanju_row):10:max(zhuanju_row);
   
   for i=1:length(su)
       for j=1:length(T)
           table(i,j)=F(T(j),su(i));
           if table(i,j)>=1
               table(i,j)=1;
           elseif table(i,j)<=0
               table(i,j)=0;
           else
               table(i,j)=table(i,j);
           end
           end
                          
       end

   [TU,SU]=meshgrid(T,su);
   mesh(SU,TU,table)