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)