第2讲_MATLAB入门（计算药物分析）

第2讲_MATLAB入门（计算药物分析）内容摘要：

第2讲_MATLAB入门（计算药物分析） + 0世纪七十年代后期：时任美国新墨西哥大学计算机科学系主任的 为学生设计了一组 “ 通俗易用 ” 的接口 ， 此即用 解决 工程 、 科学计算 和 数学学科 中许多问题。 量 、 数组 和 矩阵 的基础上，使用方便，人机界面直观，输出结果可视化。 矩阵 是 强大的 绘图、计算和统计 功能在欧美大学里 ， 诸如应用代数 、 数理统计 、 自动控制 、 数字信号处理 、 模拟与数字通信 、 时间序列分析 、 动态系统仿真等课程的教科书都把 这几乎成了九十年代教科书与旧版书籍的区别性标志。 在那里 ， 硕士生 、 博士生必须掌握的基本工具。 在国际学术界 ， 可靠的科学计算标准软件。 在许多国际一流学术刊物上 ， （ 尤其是信息科学刊物 ） ， 都可以看到 理工科大学生应该掌握的技术工具，它作为一种编程语言和可视化工具，可解决 工程 、科学计算 和 数学学科 中许多问题。 量 、 数组 和 矩阵 的基础上，使用方便，人机界面直观，输出结果可视化。 矩阵 是 1，提供了一种全新的编程语言，其语法规则简单，很贴近人的思维方式。 2， 3， 化学计量学程序调试和编制有很大好处。 4， 得特别简捷、高效和方便。 应用实例：1 简单的计算器使用法【 例 1求的算术运算结果。 （ 1）用键盘在 (12+2*(72 （ 2）在上述表达式输入完成后，按 【 ，该就指令被执行。 （ 3）在指令执行后， 2【 例 1简单矩阵的输入步骤。 （ 1）在键盘上输入下列内容A = 1,2,3; 4,5,6; 7,8,9 （ 2）按 【 键，指令被执行。 （ 3）在指令执行后， =1 2 34 5 67 8 9 【 例 1矩阵的分行输入A=1,2,34,5,67,8,9 （以下是显示结果）A =1 2 34 5 67 8 9 1） 矩阵 (算术 )运算符A表示矩阵 A+B ：表示矩阵 的和A *B：表示矩阵 的乘法A *B：表示矩阵 的对应元素相乘A.的元素除以矩阵 (2) 矩阵的创建a) 显示地输入一个元素序列；b) 用 c) 在 ： A=1 3 5 7; 12 3 5 3; 3 5 9 1 创建矩阵。 求解代数方程 问题。 教科书的基本叙述：当行数等于列数 (标量 )时 ， X=若行数大于列数满秩阵 ， X=(AA)；若列数大于行数时 ， 有无数解。 (，一般程序就必须按以上不同情况进行编程。 然而对 那只需一条简单指令： x=Ab。 计算的快速性 、 准确性和稳定性是普通程序所远不及的。 is an of % E= 0 0 0 = =)*A C=E 量 与 函 数二、数 组三、 矩 阵四、 量名规则 是：（ 1）变量名必须是不含空格的单个词；（ 2）变量名区分大小写；（ 3）变量名最多不超过 19个字符；（ 4）变量名必须以字母打头，之后可以是任意字母、数字或下划线，变量名中不允许使用标点符号 量 与 函 数特殊变量 取 值a n s 用 于 结 果 的 缺 省 变 量 名周率e p 机 的 最 小 数 ， 当 和 1 相 加 就 产 生 一 个 比1 大 的 数f l o p s 浮 点 运 算 数i n f 无 穷 大 ， 如 1 / 0N a N 不 定 量 ， 如 0 / 0i ， j i = j =1n a r g i n 所 用 函 数 的 输 入 变 量 数 目n a r g o u t 所 用 函 数 的 输 出 变 量 数 目r e a l m i n 最 小 可 用 正 实 数r e a l m a x 最 大 可 用 正 实 数特殊变量表2、数学运算符号及标点符号+ 加 法 运 算 ， 适 用 于 两 个 数 或 两 个 同 阶 矩 阵 相 加 . 减法运算* 乘法运算.* 点乘运算/ 除法运算./ 点除运算 乘幂运算. 点 乘 幂 运 算 反 斜 杠 表 示 左 除 .（ 1） 为 逗号 或 无标点 符号，则显示命令的结果；若命令后为 分号 ，则禁止显示结果 . （ 2）“ %” 后面所有文字为注释 . （ 3） “ .”表示续行 数 名 称 函 数 名 称s i n ( x ) 正弦函数 a s i n ( x ) 反 正 弦 函 数c o s ( x ) 余弦函数 a c o s ( x ) 反 余 弦 函 数t a n ( x ) 正切函数 a t a n ( x ) 反 正 切 函 数a b s ( x ) 绝 对 值 m a x ( x ) 最 大 值m i n ( x ) 最 小 值 s u m ( x ) 元 素 的 总 和s q r t ( x ) 开 平 方 e x p ( x ) 以 e 为 底 的 指 数l o g ( x ) 自然对数)(lo 0 为 底 的 对 数s i g n ( x ) 符号函数 f i x ( x ) 取整3、数学函数时为了研究某一个函数的各种性态，需要为 此必须编写函数文件 . 函数文件是文件名后缀为 类文件的 第一行必须是一特殊字符格式为： 因变量名 =函数名（自变量名）函数值的获得必须通过具体的运算实现，并赋给因变量 1. 在 :在编辑窗口中输入程序内容3. 点： 盘， 文件保存。 例：定义函数 f(x1,100(+(1 f=x)f=100*(x(2)2)2+(1)文件： . 可以直接使用函数 算 f(1,2), 只需在 x=1 2x)返回x=n）创建从开始，到结束，有 1、创建简单的数组二、数 组 x=a b c d e f 创建包含指定元素的行向量x= 1计数，到 n）创建从 数组元素的访问 3） 直接使用元素编址序号 . x(a b c d) 表示提取数组 a、 b、 c、 x(a) x(b) x(c) x(d). （ 2） 访问一块元素 ： x(a ： b ： c)表示访问数组 步长为 不超过 c）， . （ 1） 访问一个元素 ： x(i)表示访问数组 3、数组的方向前面例子中的数组都是一行数列，是行方向分布的 . 称之为 行向量 . 数组也可以是 列向量 ，它的数组操作和运算与行向量是一样的，唯一的区别是结果以列形式显示 接产生 例 c=1； 2； 3； 4转置产生 例 b=1 2 3 4; c=b 说明：以 空格 或 逗号 分隔的元素指定的是不同列的元素，而以 分号 分隔的元素指定了不同行的元素 . 4、数组的运算（ 1）标量 、乘、除、乘方是数组的每个元素对该标量施加相应的加、减、乘、除、乘方运算 . 设： a=a1,a n, c=标量则： a+c=a1+c,a2+c,a n+ca.*c=a1*c,a2*c,a n*ca./c= a1/c,a2/c,a n/c(右除）a.c= c/a1,c/,c/a n (左除）a.c= a1c,a2c,a ncc.a= ca1,c,ca n 2） 数组 、减、乘、除、幂运算可按 元素对元素 方式进行的，不同大小或维数的数组是不能进行运算的 . 设： a=a1,a n, b=b1,b n则： a+b= a1+b1,a2+,a n+bna.*b= a1*b1,a2*,a n*bna./b= a1/b1,a2/,a n/bna.b=b1/a1,b2/,b n/ana.b=a1b1,a2,a n回2. 控制指令： 例 2用 查 在指令窗中运行以下指令 ， 就可看到内存变量。 Y x y 例 2键入 获得驻留。