1、第一,启动MATLAB,新建脚本(Ctrl+N),在脚本编辑区输入如下代码:close a盟敢势袂ll;clear all;clc稆糨孝汶;dbstop if error;[x,y,z,v]=flow;sx=1:9;sy=[];sz=0;cvals=linspace(-8,2,10)figurecontourslice(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
2、第二,保存和运行上述脚本,得到如下流体流动的轮廓切片图(Contour slice of fluid flow)。
3、第三,关于脚本中的flow,linspace,contourslice等命令的介绍可以查看MATLAB帮助文档。在命令行窗口输入doc flow,查看flow命令的介绍。[x,y,z,v]=flow是返回一个25*50*25的矩阵v及其相应坐标x,y,z值。
4、第四,在命令行窗口输入doc linspace,查看linspace的介绍。cvals=l坡纠课柩inspace(-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
5、第五,在命令行窗口输入doc contourslice,查看contourslice的介绍。contourslice(x,y,z,v,sx,sy,sz)是沿着sx,sy,sz定义的二维平面,根据v值画轮廓切片图。contourslice(...,cvals)是根据length(cvals)的值在每个切片上画轮廓线。关于daspect,campos等命令,同样通过doc XXX查看其介绍,不再赘述。