Matlab编程必备手册(适合初学者).pdf

《Matlab编程必备手册(适合初学者).pdf》由会员分享，可在线阅读，更多相关《Matlab编程必备手册(适合初学者).pdf（14页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、本人精通MATLAB等编程语言，可以提供以下方向的帮助 本人精通MATLAB等编程语言，可以提供以下方向的帮助 1.1.MATLAB/GUI/SIMULINK/C+/VC+编程问题；2.MATLAB/GUI/SIMULINK/C+/VC+编程问题；2.线性与非线性控制、智能控制、模糊控制；3.线性与非线性控制、智能控制、模糊控制；3.数值计算问题、小波分析算法、有限元问题；4.数值计算问题、小波分析算法、有限元问题；4.电机控制、电力系统、机器人路径优化、机器人控制；5.电机控制、电力系统、机器人路径优化、机器人控制；5.粒子群算法、神经网络、模拟退火算法等智能优化算法；6.粒子群算法、神经网

2、络、模拟退火算法等智能优化算法；6.图像处理、信号处理、语音信号处理、电子通信等方向；图像处理、信号处理、语音信号处理、电子通信等方向；有问题的朋友，有问题的朋友，可以将问题直接发到我的邮箱可以将问题直接发到我的邮箱，24小时内给您答复！非常欢迎大家加我为QQ好友，欢迎访问我的空间！，24小时内给您答复！非常欢迎大家加我为QQ好友，欢迎访问我的空间！联系方式：联系方式：QQ：626815632QQ：626815632 邮箱： 邮箱： QQ空间：http:/ 编程必备手册编程必备手册 北京索为高科系统技术有限公司 雒海涛 2011-02-11 编程格式规范编程格式规范 这儿把编程格式规范放在最前

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

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

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

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

7、为空行%注释说明下一模块的功能%模块二代码 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 的程

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

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

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

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

12、f expression（条件）（条件）statements（语句组）语句组）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 语句举例语句举例 数论中的一个有趣问题：任取一个

13、正整数，如果是偶数，用 2 除，如果是奇数，用 3 乘再加 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 val

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

15、都为假，则执行如果所有的比较结果都为假，则执行 otherwise 后面的语句组，然后跳出后面的语句组，然后跳出 switch 结结构。构。otherwise 指令可以不出现。指令可以不出现。switch 后面的表达式后面的表达式 expression 的值可以是一个标的值可以是一个标量或字符串。量或字符串。三循环结构三循环结构 循环结构循环结构 是按照给定的条件，重复执行指定的语句。Matlab 用于实现循环结构的语句有 for 语句语句和 while 语句语句 for 循环循环 for variable=expression statement（循环体）end 表达式 expression

16、 可以是行向量，也可以是矩阵 while 循环循环 while expression（条件）statement（循环体）end 循环语句可以嵌套使用循环语句可以嵌套使用 不能在不能在 for 循环体循环体内改变循环变量的值内改变循环变量的值 为了提高代码的运行效为了提高代码的运行效率，应尽可能提高代码的向量化程度，避免率，应尽可能提高代码的向量化程度，避免 for 循环循环的使用的使用 如果预先就知道循环的次数，则可以采用如果预先就知道循环的次数，则可以采用 for 循环循环；否则，如果预先无法确定循环；否则，如果预先无法确定循环的次数，则可以使用的次数，则可以使用 while 循环循环。四四

17、其它流控制语句其它流控制语句 break 和和 continue break 语句用于终止循环的执行，即跳出最内层循环 continue 语句用于结束本次循环，进行下一次循环 break 和 continue 一般与 if 语句配合使用 return return 语句用于退出正在运行的脚本或函数，通常用在函数文件中。Matlab 常用常用函数及函数实例函数及函数实例 一一常用函数列表常用函数列表=小整理：MATLAB 常用的基本数学函数 abs(x)：纯量的绝对值或向量的长度 angle(z)：复数 z 的相角(Phase angle)sqrt(x)：开平方 real(z)：复数 z 的实部

18、 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 和 y 的最小公倍数 exp(x)：自

19、然指数 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(x)：反超越正切函数=小整理：适用

20、於向量的常用函数有：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 的累计元素总和 cumprod(x):向量 x 的累计元素总乘

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

22、ookfor 事实上是对所有在搜寻 路径下的 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-7、结束 MATLAB 有三种方法可以结束 MA

23、TLAB：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);若要改变颜色，在座标对後面加上相关字串即可：plot(x,sin(x),c,x,c

24、os(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,-1.2,1.2);此外，MATLAB 也可对图形加上各种注解与处理：xlabel(Input Value);%x 轴注解

25、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(2,2,3);plot(x,sinh(x);subplot(2,2,4);plot(x,cosh(x);MATLAB 还有其他各种二维绘图函数，以适合不同的应用，详见下表。

26、=小整理：其他各种二维绘图函数 bar 长条图 errorbar 图形加上误差范围 fplot 较精确的函数图形 polar 极座标图 hist 累计图 rose 极座标累计图 stairs 阶梯图 stem 针状图 fill 实心图 feather 羽毛图 compass 罗盘图 quiver 向量场图 二常用函数举例二常用函数举例 以下我们针对每个函数举例。当资料点数量不多时，长条图是很适合的表示方式：close all;%关闭所有的图形视窗 x=1:10;y=rand(size(x);bar(x,y);如果已知资料的误差量，就可用 errorbar 来表示。下例以单位标准差来做 资料的误

27、差量：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=linspace(0,2*pi);r=cos(4*theta);polar(theta,r);对於大量的资料，我们可用 hist 来显示资料的分 情况和统计特性。下面 几个命令可用来验证 randn 产生的高斯乱数分：x=ran

28、dn(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(-x/3);stairs(x,y);stems 可产生针状图，常被用来绘制数位讯号：x=linspace(0,10,50);y=sin(x).*exp(-x/3);stem(x,y);stairs 将资料点视为多边行顶点，并将此多边行涂上颜色：x

29、=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 和 feather 很接近，只是每个箭号的起点都在圆点：theta=linspace(0,2*pi,20);z=cos(theta)+i*sin(theta);compass(z);-3.基本 XYZ 立体绘图命令 在科学目视表示（Scientific visualization）中，三度空

30、间的立体图是 一个非常重要的技巧。本章将介绍 MATLAB 基本 XYZ 三度空间的各项绘图命 令。mesh 和 plot 是三度空间立体绘图的基本命令，mesh 可画出立体网状图，plot 则可画出立体曲面图，两者产生的图形都会依高度而有不同颜色。下 列命令可画出由函数 形成的立体网状图:x=linspace(-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 的矩

31、阵 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 和 yy 都是 21x21 的矩阵 zz=xx.*exp(-xx.2-yy.2);%计算函数值，zz 也是 21x21 的矩阵 surf(xx,yy,zz);%画出立体曲面图 为了方便测试立体绘图，MATLAB 提供了一个 peaks 函数，可产生一个凹凸有 致的曲面，包含了三个局部极大点及三个局部极小点，其方程式为：要画出此

32、函数的最快方法即是直接键入 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)我们亦可对 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

33、 inf);下列命令产生在 y 方向的水流效果：x,y,z=peaks;waterfall(x,y,z);axis(-inf inf-inf inf-inf inf);meshc 同时画出网状图与等高线：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

34、 inf);contour 画出曲面等高线在 XY 平面的投影：contour(peaks,20);plot3 可画出三度空间中的曲线：t=linspace(0,20*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=linspa

35、ce(4,10,6)%等差数列：首项为 4,末项为 10,项数为 6 若要重新安排矩阵的形状，可用 reshape 命令：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);en

36、d end format rat%使用分数来表示数值 disp(x)1 1/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），也就是说，一个函数可以 呼叫它本身。举例来说

37、，n!=n*(n-1)!，因此前面的阶乘函数可以改成递式的写法：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，而不再呼叫此函数本身。