MATLAB学习与使用：流体流动的轮廓切片图

第一，启动MATLAB，新建脚本（Ctrl+N），在脚本编辑区输入如下代码：

close all;clear all;clc;dbstop if error;

[x,y,z,v]=flow;

sx=1:9;sy=[];sz=0;

cvals=linspace(-8,2,10)

figure

contourslice(x,y,z,v,sx,sy,sz,cvals);

axis([0,10,-3,3,-3,3]);

daspect([1,1,1]);

campos([0,-20,7]);

xlabel('x');ylabel('y');zlabel('z');

box on

第二，保存和运行上述脚本，得到如下流体流动的轮廓切片图（Contour slice of fluid flow）。

第三，关于脚本中的flow，linspace，contourslice等命令的介绍可以查看MATLAB帮助文档。

在命令行窗口输入doc flow，查看flow命令的介绍。[x,y,z,v]=flow是返回一个25*50*25的矩阵v及其相应坐标x,y,z值。

第四，在命令行窗口输入doc linspace，查看linspace的介绍。cvals=linspace(-8,2,10)主要是得到-8到2之间的10个线性值。

即cvals =

  Columns 1 through 6

  -8.0000   -6.8889   -5.7778   -4.6667   -3.5556   -2.4444

  Columns 7 through 10

   -1.3333   -0.2222    0.8889    2.0000

第五，在命令行窗口输入doc contourslice，查看contourslice的介绍。contourslice(x,y,z,v,sx,sy,sz)是沿着sx,sy,sz定义的二维平面，根据v值画轮廓切片图。contourslice（...，cvals）是根据length（cvals）的值在每个切片上画轮廓线。

关于daspect，campos等命令，同样通过doc XXX查看其介绍，不再赘述。