example5_8.m

MATLAB7辅助控制系统设计与仿真 第5章

i=0;
a=-0.5:0.02:0.5;
for x=a
  i=i+1; j=0;
  for y=a
    j=j+1;
    z(i,j)=corana([x y]);
  end
end
clf
%x平面内的图形
plot(z(:,1)) 
xlabel('x')  
%x的范围为[250.0-250.25]
clf
%y平面内的图形
plot(z(1,:))
xlabel('y') 
%y的范围为[0-250]
%三维图形
mesh(a,a,z);
view(30,60);
grid;