FDTD中及MATlAB编程基础.pptVIP

* FDTD中的MATlAB编程基础 一、系统预定义变量 pi ： 圆周率 ? ，其值为 imag(log(-1)) inf，Inf ：无穷大 nan，NaN ：Not-a-Number，一个不定值，如 0/0 eps ：浮点运算相对精度 i，j ：虚部单位，即 应尽量避免给系统预定义变量重新赋值！ Matlab 中的数默认是双精度实数，表示方法同 C 语言 3, -9, 0.4, 1.603e-12, 3.23e+20 浮点运算的相对误差为 eps 浮点数表示范围为：10-308 ～ 10308 输出格式 Matlab 以双精度执行所有的运算，运算结果可以在屏幕上输出，同时赋给指定变量；若无指定变量，则系统会自动将结果赋给变量 “ans” Matlab 中数的输出格式可以通过 format 命令指定 format 只改变变量的输出格式， 但不会影响变量的值！ 各种 format 格式 3.141592653589793e+000 自由格式 format loose 3.1416 短格式（缺省显示格式），同short format 例 解释 格式 format + / format bank / format rat / format hex (详情查看联机帮助) 压缩格式 format compact 3.14159265358979 长格式g方式 format long g 3.1416 短格式g方式 format short g 长格式e方式 format long e 3.1416e+000 短格式e方式（科学计数格式） format short e 3.14159265358979 长格式，双精度数15位，单精度数7位 format long 3.1416 短格式（缺省显示格式），只显示5位 format short 二、变量的读取 清除当前工作空间中的变量 clear 清除当前工作空间中的所有变量 clear A x 清除指定的变量 将数据文件中的变量载入当前工作空间 load mydata 载入数据文件中的所有变量 load mydata A x 从数据文件中提取指定变量 三、建立矩阵的函数 常用函数有： eye(size(A)) 产生与A矩阵同阶的单位矩阵 zeros(m,n) 产生0矩阵 ones(m,n) 产生幺矩阵 rand (m,n) 产生随机元素的矩阵 Size(a) 返回包含两个元素的向量。 Length(a) 返回向量的长度。 常见矩阵生成函数 其它特殊矩阵生成函数：magic、hilb、pascal 产生 0～1 间均匀分布的随机矩阵 m=n 时简写为 rand(n) rand(m,n) 产生均值为0，方差为1的标准正态分布随机矩阵m=n 时简写为 randn(n) randn(m,n) 提取一个矩阵的上三角部分 triu(A) 提取一个矩阵的下三角部分 tril(A) 若 X 是矩阵，则 diag(X) 为 X 的主对角线向量 若 X 是向量，diag(X) 产生以 X 为主对角线的对角矩阵 diag(X) 生成一个主对角线全为 1 的 m 行 n 列矩阵, m=n 时可简写为 eye(n)，即为 n 维单位矩阵 eye(m,n) 生成一个 m 行 n 列的元素全为 1 的矩阵,  m=n 时可写为 ones(n) ones(m,n) 生成一个 m 行 n 列的零矩阵，m=n 时可简写为 zeros(n) zeros(m,n) 可利用冒号提取矩阵 的整行或整列。 例： A(1, :) A(:, 1:3) A(:, :) 四、常用命令 1. input函数：用于向计算机输入一个参数。 调用格式： A=input(提示信息，选项)； 注：‘s’选项，则允许用户输入一个字符串。 例如想输入一个人的姓名，可采用命令 xm=input(Whats your name:,s) 【例2】 求一元二次方程a2 +bx+c=0的根。 a=input(a=?); b=input(b=?); c=input(c=?); d=b*b-4*a*c; x=[(-b+sqrt(d))/(2*a),(-b-sqrt(d))/(2*a)] 将该程序以aa.m文件存盘，然后运行aa.m文件。 2． pause函数：暂停程序的执行。 调用格式： pause(延迟秒数) 注：如果省略延迟时间，直接使用pause，则将暂停程序，直到用户按任一键后程序继续执行。 3.Drawnow函数： 将还未处理完的图像实时的显示出