matlab编程必备.doc

《matlab编程必备.doc》由会员分享，可在线阅读，更多相关《matlab编程必备.doc（13页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、一编程原则 1. 正确：能准确实现原仿真目的； 2. 高效：循环向量化，少用或不用循环，尽量调用MATLAB 自带函数； 3. 清晰：养成良好的编程习惯，程序具有良好的可读性； 4. 通用：程序具高度的可移植性，同时易于扩展，避免同学之间重复工作。 二. 编程规则 1. 定义变量，以英文单词小写缩写开头表示类别名，再接具体变量的英文名称英文单词， 具体变量首字母大写。 教研室常用类别缩写：最大值 max，最小值min 例如：定义变量存贮临时数组TempArray 的最大值maxTempArray 依照工程大小确定变量名长短，小范围应用的变量应该用短的变量名。定义务必清晰，避免 混淆。 2. 循

2、环变量使用常用变量i,j,k；程序中使用复数时，采用i,j以外的循环变量以避免和虚数单位冲突，同时要在注释部分说明变量意义。 3. 程序应高内聚、低耦合、模块函数化，便于移植、重复使用。 4. 使用if 语句判断变量是否等于某一常数时，将常变量数写在等号之前，常数写在等号之后。 例如判断变量a 是否等于100 写作if a =100。 5. 用常数代替数字，少用或不用数字。 例如上一条：写作if a = =100 就不标准。应先定义meanConst=100；为期望常量。同时在注释中说明，然后在程序部分写作：if a = =const；如果要修改期望常量的话，只要在程序定义部分修改就可以，不必

3、再逐行修改，此外还易于在该程序的基础上增加其他功能。 三. 注释规则 1. 注释要用汉语，且用一个“%”连接文字作为注释开始标志。程序的起始要注明作者、时间、程序功能、复杂的算法需要加上流程说明。对于比较复杂的程序，由几个人协作完成,要在每一个模块的起始分别注明作者、时间、程序功能。注释采用matlab 的注释符号。同时第一行增加关键字，便于将来程序文件较多时，使用help 和lookfor 查找。例如： %关键字% %此处为函数名，如果不是函数文件就取消此行 % %创建人： % 日期： %修改人： % 日期： %函数变量及功能的简单描述 % 2. 对于所有有物理含义的变量、常量、数据结构声明

4、，在声明时都必须加以注释，说明其 物理含义，单位。 3. 源程序有效注释量必须在 20%以上，边写代码边注释，修改代码同时修改相应的注释， 以保证注释与代码的一致性，不再有用的注释要删除。注释语言必须准确、易懂、简洁，尽 量不要用缩写。如果接手前人的程序，在维护过程中需要对变量以及源程序作部分修改及补 充，在修改变量及程序注释的同时需要在修改的位置标明修改人的名字及修改日期。 对原程序修改要及时修改注释。否则有注释比没有注释还要糟糕。 4. 程序内部不同功能模块要空一行，同时附注释说明模块的作用、功能。这样程序结构清 晰易懂，便于同学间学习交流。 例如： 模块一代码 此处为空行 %注释说明下一

5、模块的功能% 模块二代码 5. 一行代码不宜过长，应限制在80 列之内。一旦有过长代码须用MATLAB 转接符号“.” 连接两行代码，以便于同学之间打印程序交流。因为大多是打印机的行长度为80。 例 ： 第一行代码 nVariable=1+1+1+1+1+1+1+1+1+11+1+1+1+1+1+1+1+1+1+1+1+1. 第二行代码 +1+1+1; 编程入门 一前言 Matlab 作为一种广泛应用于科学计算的工具软件，不仅具有强大的数值计算、符号计算、 矩阵运算能力和丰富的绘图功能，同时也具有和 C、FORTRAN等高级语言一样进行程序设 计。利用 Matlab 的程序控制功能，可以将有关

6、 Matlab 命令编成程序存储在一个文件中（M 文 件），然后在命令窗口中运行该文件，Matlab就会自动依次执行文件中的命令，直到全部命令执行完毕.在 Matlab 程序设计中，要充分利用 Matlab 数据结构的特点，提高编程效率。 二M 文件 M 文件以 .m为扩展名。M 文件是由若干 Matlab 命令组合在一起构成的，它可以完 成某些操作，也可以实现某种算法。事实上，Matlab 提供的内部函数以及各种工具箱，都 是利用 Matlab 语言开发的 M 文件。用户也可以结合自己的工作需要，开发自己的程序或工具箱。 M 文件根据调用方式的不同可以分为两类：Script：脚本文件/命令文

7、件 ，Function：函数文件。它们均是普通的ASCII码构成的文件。M 脚本文件中包含一族由 MATLAB 语言所支持的语句，它类似于 DOS 下的批处理文件，它的执行方式很简单，用户只需在 MATLAB的提示符 下键入该 M 文件的文件名，这样 MATLAB 就会自动执行该 M 文件中的各条语句，并将结果直接返回到 MATLAB 的工作空间。M 函数格式是MATLAB程序设计的主流，一般情况下，不建议您使用 M 脚本文件格式编程。 MATLAB 的M 函数是由function语句引导的，其基本格式如下： function 返回变量列表 = 函数名(输入变量列表) 注释说明语句段, 由%引

8、导 输入、返回变量格式的检测 函数体语句 建议：M 文件名与文件内主函数名相同（matlab是以文件名做区分的，当然函数名和文件名也可以不同）。 三程序实例 1.脚本文件举例：编写一个脚本文件将华氏温度转化为摄氏温度。 5 c = - ( f -32) 9 新建一个 M 文件f2cs.m，内容如下： clear; % 清除当前工作空间中的变量 f=input(Please input Fahrenheit temperature:); c=5*(f-32)/9; fprintf(The centigrade temperature is %gn,c); 在命令窗口中输入f2cs，即可执行该M

9、文件。 2.function 文件举例：计算两个数字之和。 新建一个m 文件，命名为MySum.m，内如如下： function rt = MySum(x,y) rt=x+y; end 在命令窗口输入MySum(10,20)即可输出计算结果。 程序控制结构 一顺序结构 按排列顺序依次执行各条语句，直到程序的最后，这是最简单的一种程序结构，一般涉及数 据的输入输出、数据的计算或处理等。 二选择结构 选择结构 是根据给定的条件成立或不成立，分别执行不同的语句。Matlab用于实现选择结 构的语句有if语句和switch语句。 If 条件语句 单分支结构 if expression （条件） sta

10、tements （语句组） end 双分支结构 if expression （条件） statements1 （语句组1） else statements2 （语句组2） end 多分支结构 if expression1 （条件1） statements1 （语句组1） elseif expression2 （条件2） statements2 （语句组2） . . elseif expressionm （条件m） statementsm （语句组m） else statements （语句组） end if 语句举例 数论中的一个有趣问题：任取一个正整数，如果是偶数，用2除，如果是奇数，用3乘

11、再加 1，反复这个过程，直到所得到的数为1。 问：是否存在使该过程永不中止的整数？ while 1 n=input(Please enter n(nonpositive quit):); if n1 if rem(n,2)=0 n=n/2; else n=3*n+1; end fprintf(n n=%d,n); end fprintf(n n=%d is not we need! continue . n, nt); end switch 语句 根据表达式的不同取值，分别执行不同的语句。 switch expression （表达式） case value1 （表达式1） statement

12、1 （语句组1） case value2 （表达式2） statement2 （语句组2） . . case valuem （表达式m） statementm （语句组m） otherwise statement （语句组） end l Matlab 首先计算 expression 的值，然后将它依次与各个case指令后的检测值进行比较，当比较结果为真时，就执行相应的语句组，然后跳出 switch结构。 l 如果所有的比较结果都为假，则执行otherwise后面的语句组，然后跳出switch结构。 l otherwise 指令可以不出现。 l switch后面的表达式 expression的值

13、可以是一个标量或字符串。 三循环结构 循环结构 是按照给定的条件，重复执行指定的语句。Matlab 用于实现循环结构的语句有for 语句和while语句 for 循环 for variable=expression statement （循环体） end 表达式 expression 可以是行向量，也可以是矩阵 while 循环 while expression （条件） statement （循环体） end l 循环语句可以嵌套使用 l 不能在for循环体内改变循环变量的值 l 为了提高代码的运行效率，应尽可能提高代码的向量化程度，避免 for循环的使用 l 如果预先就知道循环的次数，则可

14、以采用for 循环；否则，如果预先无法确定循环 的次数，则可以使用while 循环。 四其它流控制语句 break 和 continue l break 语句用于终止循环的执行，即跳出最内层循环 l continue 语句用于结束本次循环，进行下一次循环 l break 和continue 一般与if 语句配合使用 return l return 语句用于退出正在运行的脚本或函数，通常用在函数文件中。 Matlab 常用函数及函数实例 一常用函数列表 = 小整理：MATLAB 常用的基本数学函数 abs(x)：纯量的绝对值或向量的长度 angle(z)：复数z 的相角(Phase angle)

15、 sqrt(x)：开平方 real(z)：复数z 的实部 imag(z)：复数z 的虚部 conj(z)：复数z 的共轭复数round(x)：四舍五入至最近整数 fix(x)：无论正负，舍去小数至最近整数 floor(x)：地板函数，即舍去正小数至最近整数 ceil(x)：天花板函数，即加入正小数至最近整数 rat(x)：将实数x 化为分数表示 rats(x)：将实数x 化为多项分数展开 sign(x)：符号函数(Signum function)。 当x0 时，sign(x)=1。 rem(x,y)：求x 除以y 的馀数 gcd(x,y)：整数x 和y 的最大公因数 lcm(x,y)：整数x

16、和y 的最小公倍数 exp(x)：自然指数 pow2(x)：2 的指数 log(x)：以e 为底的对数，即自然对数或 log2(x)：以2 为底的对数 log10(x)：以10 为底的对数 = 小整理：MATLAB 常用的三角函数 sin(x)：正弦函数 cos(x)：馀弦函数 tan(x)：正切函数 asin(x)：反正弦函数 acos(x)：反馀弦函数 atan(x)：反正切函数 atan2(x,y)：四象限的反正切函数 sinh(x)：超越正弦函数 cosh(x)：超越馀弦函数 tanh(x)：超越正切函数 asinh(x)：反超越正弦函数 acosh(x)：反超越馀弦函数 atanh(

17、x)：反超越正切函数 = 小整理：适用於向量的常用函数有： min(x): 向量x 的元素的最小值 max(x): 向量x 的元素的最大值 mean(x): 向量x 的元素的平均值 median(x): 向量x 的元素的中位数 std(x): 向量x 的元素的标准差 diff(x): 向量x 的相邻元素的差 sort(x): 对向量x 的元素进行排序（Sorting ） length(x): 向量x 的元素个数 norm(x): 向量x 的欧氏（Euclidean）长度 sum(x): 向量x 的元素总和 prod(x): 向量x 的元素总乘积 cumsum(x): 向量x 的累计元素总和 c

18、umprod(x): 向量x 的累计元素总乘积 dot(x, y): 向量x 和y 的内积 cross(x, y): 向量x 和y 的外积 （大部份的向量函数也可适用於矩阵，详见下述。） = 小整理：MATLAB 的查询命令 help：用来查询已知命令的用法。例如已知inv 是用来计算反矩阵，键入 help inv 即可得知有关inv 命令的用法。（键入help help 则显示help 的用 法，请试看看！） lookfor：用来寻找未知的命令。例如要寻找计算反矩阵的命令，可键入 lookfor inverse，MATLAB 即会列出所有和关键字inverse 相关的指令。找到所需的 命令後

19、，即可用help 进一步找出其用法。（lookfor 事实上是对所有在搜寻 路径下的M 档案进行关键字对第一注解行的比对，详见後叙 = 下表即为MATLAB 常用到的永久常数。 小整理：MATLAB 的永久常数 i 或j ：基本虚数单位（即） eps：系统的浮点（Floating-point）精确度 inf：无限大，例如1/0 nan 或NaN：非数值（Not a number），例如0/0 pi：圆周率p （= 3.1415926.） realmax：系统所能表示的最大数值 realmin：系统所能表示的最小数值 nargin: 函数的输入引数个数 nargin: 函数的输出引数个数 = 1

20、-7、结束MATLAB 有三种方法可以结束MATLAB： 1.键入exit 2.键入quit 3.直接关闭MATLAB 的命令视窗（Command window） = 小整理：MATLAB 基本绘图函数 plot: x 轴和y 轴均为线性刻度（Linear scale） loglog: x 轴和y 轴均为对数刻度（Logarithmic scale） semilogx: x 轴为对数刻度，y 轴为线性刻度 semilogy: x 轴为线性刻度，y 轴为对数刻度 = 若要画出多条曲线，只需将座标对依次放入plot 函数即可： plot(x, sin(x), x, cos(x); 若要改变颜色，在

21、座标对後面加上相关字串即可： plot(x, sin(x), c, x, cos(x), g); 若要同时改变颜色及图线型态（Line style），也是在座标对後面加上相 关字串即可： plot(x, sin(x), co, x, cos(x), g*); = 小整理：plot 绘图函数的叁数 字元 颜色字元 图线型态 y 黄色 . 点 k 黑色o 圆 w 白色x x b 蓝色+ + g 绿色 * * r 红色- 实线 c 亮青色 : 点线 m 锰紫色-. 点虚线 - 虚线 = 图形完成後，我们可用axis(xmin,xmax,ymin,ymax)函数来调整图轴的范围： axis(0, 6,

22、 -1.2, 1.2); 此外，MATLAB 也可对图形加上各种注解与处理： xlabel(Input Value); % x 轴注解 ylabel(Function Value); % y 轴注解 title(Two Trigonometric Functions); % 图形标题 legend(y = sin(x),y = cos(x); % 图形注解 grid on; % 显示格线 我们可用subplot 来同时画出数个小图形於同一个视窗之中： subplot(2,2,1); plot(x, sin(x); subplot(2,2,2); plot(x, cos(x); subplot(

23、2,2,3); plot(x, sinh(x); subplot(2,2,4); plot(x, cosh(x); MATLAB 还有其他各种二维绘图函数，以适合不同的应用，详见下表。 = 小整理：其他各种二维绘图函数 bar 长条图 errorbar 图形加上误差范围 fplot 较精确的函数图形 polar 极座标图 hist 累计图 rose 极座标累计图 stairs 阶梯图 stem 针状图 fill 实心图 feather 羽毛图 compass 罗盘图 quiver 向量场图 二常用函数举例 以下我们针对每个函数举例。 当资料点数量不多时，长条图是很适合的表示方式： close

24、all; % 关闭所有的图形视窗 x=1:10; y=rand(size(x); bar(x,y); 如果已知资料的误差量，就可用errorbar 来表示。下例以单位标准差来做 资料的误差量： x = linspace(0,2*pi,30); y = sin(x); e = std(y)*ones(size(x); errorbar(x,y,e) 对於变化剧烈的函数，可用fplot 来进行较精确的绘图，会对剧烈变化处进 行较密集的取样，如下例： fplot(sin(1/x), 0.02 0.2); % 0.02 0.2是绘图范围 若要产生极座标图形，可用polar： theta=linspac

25、e(0, 2*pi); r=cos(4*theta); polar(theta, r); 对於大量的资料，我们可用hist 来显示资料的分 情况和统计特性。下面 几个命令可用来验证randn 产生的高斯乱数分 ： x=randn(5000, 1); % 产生5000 个 ?=0，?=1 的高斯乱数 hist(x,20); % 20 代表长条的个数 rose 和hist 很接近，只不过是将资料大小视为角度，资料个数视为距离，?眉昊嬷票硎荆? x=randn(1000, 1); rose(x); stairs 可画出阶梯图： x=linspace(0,10,50); y=sin(x).*exp(-

26、x/3); stairs(x,y); stems 可产生针状图，常被用来绘制数位讯号： x=linspace(0,10,50); y=sin(x).*exp(-x/3); stem(x,y); stairs 将资料点视为多边行顶点，并将此多边行涂上颜色： x=linspace(0,10,50); y=sin(x).*exp(-x/3); fill(x,y,b); % b为蓝色 feather 将每一个资料点视复数，并以箭号画出： theta=linspace(0, 2*pi, 20); z = cos(theta)+i*sin(theta); feather(z); compass 和feat

27、her 很接近，只是每个箭号的起点都在圆点： theta=linspace(0, 2*pi, 20); z = cos(theta)+i*sin(theta); compass(z); - 3.基本XYZ 立体绘图命令 在科学目视表示（Scientific visualization）中，三度空间的立体图是 一个非常重要的技巧。本章将介绍MATLAB 基本XYZ 三度空间的各项绘图命 令。 mesh 和plot 是三度空间立体绘图的基本命令，mesh 可画出立体网状图， plot 则可画出立体曲面图，两者产生的图形都会依高度而有不同颜色。下 列命令可画出由函数形成的立体网状图: x=linsp

28、ace(-2, 2, 25); % 在x 轴上取25 点 y=linspace(-2, 2, 25); % 在y 轴上取25 点 xx,yy=meshgrid(x, y); % xx 和yy 都是21x21 的矩阵 zz=xx.*exp(-xx.2-yy.2); % 计算函数值，zz 也是21x21 的矩阵 mesh(xx, yy, zz); % 画出立体网状图 surf 和mesh 的用法类似： x=linspace(-2, 2, 25); % 在x 轴上取25 点 y=linspace(-2, 2, 25); % 在y 轴上取25 点 xx,yy=meshgrid(x, y); % xx

29、和yy 都是21x21 的矩阵 zz=xx.*exp(-xx.2-yy.2); % 计算函数值，zz 也是21x21 的矩阵 surf(xx, yy, zz); % 画出立体曲面图 为了方便测试立体绘图，MATLAB 提供了一个peaks 函数，可产生一个凹凸有 致的曲面，包含了三个局部极大点及三个局部极小点，其方程式为： 要画出此函数的最快方法即是直接键入peaks： peaks z = 3*(1-x).2.*exp(-(x.2) - (y+1).2) . - 10*(x/5 - x.3 - y.5).*exp(-x.2-y.2) . - 1/3*exp(-(x+1).2 - y.2) 我们

30、亦可对peaks 函数取点，再以各种不同方法进行绘图。meshz 可将曲面加上围裙： x,y,z=peaks; meshz(x,y,z); axis(-inf inf -inf inf -inf inf); waterfall 可在x 方向或y 方向产生水流效果： x,y,z=peaks; waterfall(x,y,z); axis(-inf inf -inf inf -inf inf); 下列命令产生在y 方向的水流效果： x,y,z=peaks; waterfall(x,y,z); axis(-inf inf -inf inf -inf inf); meshc 同时画出网状图与等高线：

31、x,y,z=peaks; meshc(x,y,z); axis(-inf inf -inf inf -inf inf); surfc 同时画出曲面图与等高线： x,y,z=peaks; surfc(x,y,z); axis(-inf inf -inf inf -inf inf); contour3 画出曲面在三度空间中的等高线： contour3(peaks, 20); axis(-inf inf -inf inf -inf inf); contour 画出曲面等高线在XY 平面的投影： contour(peaks, 20); plot3 可画出三度空间中的曲线： t=linspace(0,2

32、0*pi, 501); plot3(t.*sin(t), t.*cos(t), t); 亦可同时画出两条三度空间中的曲线： t=linspace(0, 10*pi, 501); plot3(t.*sin(t), t.*cos(t), t, t.*sin(t), t.*cos(t), -t); y(2:4)-1 % 取出y 的第二至第四个元素来做运算 同样的方法可用於产生公差为1 的等差数列：x = 7:16 x = 7:3:16 % 公差为3 的等差数列 x = linspace(4, 10, 6) % 等差数列：首项为4,末项为10,项数为6 若要重新安排矩阵的形状，可用reshape 命令

33、：B = reshape(A, 4, 2) % 4 是新矩阵的列数，2是新矩阵的行数 举例来说，下列命令会产生一个长度为6 的调和数列（Harmonic sequence）： x = zeros(1,6); % x 是一个16 的零矩阵 for i = 1:6, x(i) = 1/i; end for圈可以是多层的，下例产生一个16 的Hilbert矩阵h，其中为於第i列、第j行的元素为： h = zeros(6); for i = 1:6, for j = 1:6, h(i,j) = 1/(i+j-1); end end format rat % 使用分数来表示数值 disp(x) 1 1/

34、2 1/3 1/4 1/5 1/6 function output = fact(n) % FACT Calculate factorial of a given positive integer. output = 1; for i = 1:n, output = output*i; end 其中fact 是函数名，n 是输入引数，output 是输出引数，而i 则是此函数用到的暂时变数。要使用此函数，直接键入函数名及适当输入引数值即可： MATLAB 的函数也可以是递式的（Recursive），也就是说，一个函数可以呼叫它本身。举例来说，n! =n*(n-1)!，因此前面的阶乘函数可以改成

35、递式的写法： function output = fact(n)% FACT Calculate factorial of a given positive integer recursively. if n = 1, % Terminating condition output = 1; return; end output = n*fact(n-1); 在写一个递函数时，一定要包含结束条件（Terminating condition），否则此函数将会一再呼叫自己，永远不会停止，直到电脑的记忆体被耗尽为止。以上例而言，n=1 即满足结束条件，此时我们直接将output 设为1，而不再呼叫此函数本身。