matab画图助手.docVIP

clear;clc a=load(at1.txt); x=a(:,2); y=a(:,3); h=a(:,4); gnq=a(:,5); xi=linspace(min(x),max(x),319); yi=linspace(min(x),max(x),319); zi=griddata(x,y,gnq,xi,yi,cubic); %插值 hold on; contourf(xi,yi,zi); %填充图 contour(xi,yi,zi); %画等高线 title(城市功能区分布图); %标题插入 colormap(1-gray); %灰度设置 axis([0 3.0*10^4 0 18000]); %x，y轴设置 hc=colorbar; hc=colorbar(ytick,[1:5]) %划分为五个范围 newstr={1生活区,2工业区,3山区,4交通区,5生活区}; set(hc,Yticklabel,newstr) gtext(城市功能类别图例) %自己添加命令 hold off surf(x,y,,z) 三维图 hold on contour(x,y,z,30) 等高线 30指的是花多少条等高线 修改函数可以把等高线投影到自己想要的平面 海拔值改成整数 for n=p set(n,string,sprintf(%.0f,get(n,userdata))) End 函数contour 描述：矩阵的等高线图 调用语法： contour(Z)contour(Z,n)contour(Z,v)contour(X,Y,Z)contour(X,Y,Z,n)contour(X,Y,Z,v)contour(...,LineSpec)[C,h] = contour(...) [C,h] = contour(v6,...) 详细解释：等高线图显示了矩阵Z的等值线，标记等高线可以用函数clabel. contour(Z)画出了矩阵Z的等值线图，图中Z以高度的形式在x-y平面上表示出来。Z必须至少是2X2的矩阵。等高线划分的数目和等高线的数值根据Z中最小和最大值自动确定。x和y轴的范围分别是[1:n]和[1：m]，其中[m,n]=size(Z). contour(Z,n) 画矩阵Z有n条线的等值图。 contour(Z,v) 画矩阵Z由指定向量v确定的数值的等值图。等高线的数目等于length(v).画单独的等高线阶数为i的用法为contour(Z,[i i]). contour(X,Y,Z),contour(X,Y,Z,n),和contour(X,Y,Z,v)画Z的等高线图。X和Y明确指定了x和y轴的范围。当X和Y是矩阵时，必须大小和Z相同。例如在确定平面或者由surf函数定义平面时。 如果X或者Y是不规则间隔时，contour计算等高线是利用规则间隔的等高线网格，然后转换到数据X或者Y。 contour(...,lineSpec) 依据LineSpec使用指定的线型和颜色来画等高线图。等高线忽略标记符号 [C,h]=contour(...)返回等高线矩阵C和等高线组对象的句柄。clabel函数利用等高线矩阵C来创建标签。 [C,h]=contour(v6...)涉及patch图形对象和Matlab版本问题，不做介绍。 具体实例： 润滑油膜厚度等高线图显示，源数据f1.dat f2.dat f3.dat f1.dat其中数据为x轴上的取值：-1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 大小为1X21 double f2.dat其中数据为y轴上的数值，大小为1X151 double f3.da为z轴方向的数据大小为21X151 double 要求：等高线数目显示10条，颜色黑，有数值显示. 前提：将f1 f2 f3数据通过Load data file图标导入Matlab 程序修改过程： 1.初始程序 [X,Y]=meshgrid(f2,f1); [C,h]=contour(X,Y,f3); clabe(C,h,manual); 出来的结果如图 分析：根据f3.dat中数据，主要是想把其中下图等高线图表示出来 结果出来的显然不符合预期要求。因此要修改等高线间隔大小。 2 修改后程序 [xx,yy]=meshgrid(f2,f1);max_z=max(max(H_3D));min_z=min(min(H_3D));logb=log10(m