matlab在地质中的应用讲义.ppt

（1）if 逻辑表达式 执行语句 end 2.1 二维三维图形的绘制 为了获取这一地区的点的高度，在原图上，以左下角为坐标原点，图的下边界为X轴，图的左边界为Y轴，分别在两轴上每隔1cm取点，根据图的比例尺1:20000，图上的1cm在实际中表示200米，根据取的点作方格，如图2，再根据地质图上原有的等高线的信息，就可以得到一系列的（X,Y,H）。其中（X,Y）表示在这张地形地质图上某点的位置坐标，H表示该点的高程（通过图上的等高线可以得出来）。下表就是得到的数据： Matlab在构造地质学中的应用 2.1 二维三维图形的绘制 在编程中为了导入的方便，将相同的数据改用txt格式进行导入(h.txt)。 2.1 二维三维图形的绘制 用原始数据作该地区的地貌图。 x=0:200:3000; y=0:200:2400; load h.txt; meshz(x,y,h); xlabel(X),ylabel(Y),zlabel(H) 运行结果如图3。 明显山区表面很粗糙，效果不好。 分号在matlab语言中所起的作用是作为分隔行的标志，这里包括在矩阵定义中对行的区分以及在程序代码中对行的划分。在矩阵定义中，分号之间的数据被认为是矩阵的同一行。对于以分号结尾的行语句，matlab语言不会把其运算结果显示在命令窗口中。 矩阵的定义需要使用中括号。此外，中括号可以表示空矩阵，并通过该空矩阵实现对矩阵的部分删除。 1.3 Matlab的简单操作 1.3 Matlab的简单操作 a=[1,2,3,4;5,6,7,8;4,3,2,1;8,7,6,5] a = 1 2 3 4 5 6 7 8 4 3 2 1 8 7 6 5 a(:,1:2)=[] a = 3 4 7 8 2 1 6 5 b=a([1 3],:) b = 3 4 2 1 1.3 Matlab的简单操作 编程语句结构 Matlab语言最基本的赋值语句结构为： 变量名列表=表达式 注1：整个赋值语句以；结束，则不在屏幕上返回 结果，否则立即返回结果。 注2：多个语句可在同一行，用逗号分开。 注3：表达是太长可以用续行符号… 1.3 Matlab的简单操作 编程语句结构 循环语句 for, while 条件转移 if end, if elseif else end 开关语句 switch case 注释语句 % 中断语句 break 暂停语句 pause 回显语句 echo on/off 1.3 Matlab的简单操作 for循环语句 for循环语句其功能是重复执行一条或一组语句的固定次数(也可以是零次)。 表达式实际上是一个矩阵，因为在Matlab中基本元素就是矩阵。当矩阵的每一列元素依次赋给变量时，循环体语句就执行一次。 1.3 Matlab的简单操作 for循环语句 基本格式： for 循环变量＝起始值：步长：终止值 循环体 end 步长缺省值为1，可以在正实数或负实数范围内任意指定。对于正数，循环变量的值大于终止值时，循环结束；对于负数，循环变量的值小于终止值时，循环结束。循环结构可以嵌套使用。 1.3 Matlab的简单操作 for循环语句 for i=1:5 for j=1:5 a(i,j)=1/(i+j-1); end end a= 1.0000 0.5000 0.3333 0.2500 0.2000 0.5000 0.3333 0.2500 0.2000 0.1667 0.3333 0.2500 0.2000 0.1667 0.1429 0.2500 0.2000 0.1667 0.1429 0.1250 0.2000 0.1667 0.1429 0.1250 0.1111 1.3 Matlab的简单操作 while循环语句 基本格式： while 表达式 循环体 end 若表达式为真，则执行循环体的内容，执行后再判断表达式是否为真，若不为真，则跳出循环体，向下继续执行。