化工微机应用---第三讲.pptVIP

化工微机应用(第三讲)Applications of Computer in Chemical Engineering 一、数值数组及其运算 数组：由实数或复数排成的长方阵列 数组运算：任何运算均是对数组中的每个元素平等的操作 Example 1 绘制函数 在 时的函数曲线 程序: x=0:0.1:1 y=x.*exp(-x) plot(x,y),xlabel(x),ylabel(y),title(y=x*exp(-x)) 一维数组的创建 逐个元素输入法 冒号生成法 指令： x=a: inc: b a是数组的第一个元素；inc是采样点之间的间隔，即步长。 若(b-a)是inc的整数倍，则生成数组的最后一个元素等于b；否则小于b。 a,inc,b之间使用冒号”:”(英文)隔开。 inc可以省略。省略时，默认值为1 Inc可正可负，取负值时ab 定数线性采样法 指令：x=linspace(a,b,n) a,b分别是生成数组的第一个和最后一个元素。n是采样点总数。 该指令生成(1*n)数组 该指令等价于x=a:(b-a)/n-1:b Example 2 一维数组的寻访与赋值 子数组的赋值 二维数组的创建 直接输入法 a、数组元素置于方括号 “[ ]”中 b、数组的行与行之间必须用分号“;”，或者回车键“Enter”分隔 c、每一行内的元素必须由逗号“,”或者空格分隔 d、数组元素可以既有实数也有复数 Example 3 复数数组的输入 利用M文件输入和保存数组 操作步骤： a、打开脚本文件编辑器 b、输入数组，如： AM=[ 2, 3, 5, 7,11,13;… 17,19,23,29,31,37 ] c、保存文件，如 MN.m d、在命令窗中运行MN，可以将数组AM调入MATLAB内存 二维数组元素的标识 “全下标”标识： 数组元素采用第几行、第几列标识 A(m,n)?第m行、第n列 “单下标”标识： 对二维数组采用一维编号，编号方法：将二维数组的所有列，按照从左至右的顺序首尾相接排成“一维长列”，然后自上而下对元素的位置编号 如：(m*n)二维数组A的第r行，c列元素的单下标为l=(c-1)*m+r “逻辑1”标识 便于寻找数组中大于某值的元素 Example 4 数组运算的常用函数 函数数组的运算规则 对于(m*n)数组 X=[xij]m*n , 函数 f ( . )的数组运算规则是: f (X)=[f (xij)]m*n 常用函数举例2 指数函数 常用函数举例3 复数函数 Example 5　 数组运算和矩阵运算 数组和矩阵的联系与区别 a、外形--- 一致 b、运算规则--- 不同 在编制程序时要分清数组操作与矩阵操作 数组与矩阵运算异同 Example 6 标准生成数组 Example 6 数组操作函数 Example 7 数组构作技法综合 为了生成比较复杂的数组，或者为了更改和扩展已有数组，MATLAB提供了翻转、插入提取、收缩和重组等操作，构作技法方便灵活 数组的赋值扩展 多次寻访扩展法 合成扩展法 提取子数组；合成新数组 单下标寻访和reshape指令演示 “对列（或行）同加一个数”三种的操作方法 逻辑函数的运用 按指定次序，对矩阵各维的次序进行重组 permute 按照指定维上的数目铺放模块数组 repmat 以数组的“水平中线”为对称轴，上下交换 flipud 以数组的“垂直中线”为对称轴，左右交换 fliplr 以指定维的对称轴进行翻转 flipdim 提取对角线元素,或者生成对角阵 diag 含义 指令 提取数组上三角部分，生成上三角阵 triu 提取数组下三角部分，生成下三角阵 tril 逆时针旋转二维矩阵900 rot90 总元素不变的前提下，改变各维的大小 reshape 含义 指令 * * 常用函数举例1 三角函数 双曲余弦 cosh 反正弦 asin 双曲正弦 sinh 正切 tan 反正切 atan 余弦 cos 反余弦 acos 正弦 sin 含义 名称 含义 名称 平方根 sqrt 2的幂 pow2 以2为底的对数 log2 常