《matlab编程初步.ppt》由会员分享，可在线阅读，更多相关《matlab编程初步.ppt（43页珍藏版）》请在三一办公上搜索。
1、MATLAB编程初步,2023年11月8日,2,Matlab 编程入门,Matlab 编程简介,Matlab 作为一种广泛应用于科学计算的工具软件，不仅具有强大的数值计算、符号计算、矩阵运算能力和丰富的绘图功能，同时也具有和 C、FORTRAN 等高级语言一样进行程序设计。,在 Matlab 程序设计中，要充分利用 Matlab 数据结构的特点，提高编程效率。,利用 Matlab 的程序控制功能，可以将有关 Matlab 命令编成程序存储在一个文件中（M 文件），然后在命令窗口中运行该文件，Matlab 就会自动依次执行文件中的命令，直到全部命令执行完毕。,2023年11月8日,3,用 Mat
2、lab 语言编写的程序称为 M 文件。,M 文件是由若干 Matlab 命令组合在一起构成的，它可以完成某些操作，也可以实现某种算法。事实上，Matlab 提供的内部函数以及各种工具箱，都是利用 Matlab 语言开发的 M 文件。用户也可以结合自己的工作需要，开发自己的程序或工具箱。,M 文件,M 文件以.m 为扩展名。,M脚本文件,M函数文件,脚本文件与函数文件的区别,2023年11月8日,7,M 文件是一个文本文件，可以用任何文本编辑器来建立和编辑，通常使用 Matlab 自带的 M 文件编辑器。,M 文件的建立与打开,M文件路径的加载,通过修改current directory选项将M
3、文件所在底层文件夹加载进来才可以使用。点击后面的符号即可修改路径。,2023年11月8日,9,在命令窗口中输入 f2cs，即可执行该 M 文件。,例：编写一个脚本文件将华氏温度转化为摄氏温度,脚本文件举例,脚本文件运行后，文件中所产生的所有变量都驻留在 Matlab 当前工作空间，可以用 whos 查看。,2023年11月8日,10,关系运算,比较大小，如果结论是 真 则返回 1，否则返回 0,注意=与=的区别,关系操作符能用来比较两个同样大小的数组，或用来比较一个数组和一个标量，在后一种情况，标量和数组中的每一个元素相比较，比较结果与数组大小一样。,关系运算符,2023年11月8日,11,例
4、：,关系运算举例,2+2=4 ans=1,23 ans=0,A=1 3 5;2 0 6;A=1 3 5 B=3 1 0 ans=0 1 1 B=3 1 0;2 4 6;2 0 6 2 4 6 1 0 1 A=B,x=5 0;4 2;x=5 0 ans=0 1 x4 4 2 0 1,2023年11月8日,12,逻辑运算,逻辑运算符,A&B 等价于 and(A,B)A|B 等价于 or(A,B)A 等价于 not(A),2023年11月8日,13,逻辑运算表,2023年11月8日,14,all(x)如果向量 X 中所有元素都非零，则返回 1，否则返回 0,若 x 为矩阵，则 any 和 all 按
5、列运算，返回一个 0-1 向量,逻辑运算,相关函数,any(x)如果向量 X 中存在非零元素，则返回 1，否则返回 0A=0 2;3 4;5 6 any(A)any按列运算A=0 2 ans=1 1 3 4 all(A)5 6 ans=0 1 all也按列运算,2023年11月8日,15,程序控制结构有三种：顺序结构、选择结构和循环结构。任何复杂的程序都由这三种基本结构组成。,顺序结构,M文件控制流,2023年11月8日,16,数据的输入：input,数据输入 input,2023年11月8日,17,disp(X),输出变量 X 的值，可以是数值矩阵或字符串,B=1 2 3;4 5 6;7 8
6、 9;disp(B),数据输出 disp,数据的输出：disp,2023年11月8日,18,数据输出 fprintf,format 中还可以使用一些特殊格式，如：,n(换行)t(制表符)b(退格)(反斜杆)%(百分号),2023年11月8日,19,数据输出 fprintf,a=Hello;b=2.4;c=100*pi;fprintf(a=%s,b=%f,c=%en,a,b,c)a=Hello,b=2.400000,c=3.141593e+002,例：,format 中的输出格式要与输出变量一一对应,可以没有输出变量,2023年11月8日,20,pause,程序的暂停：pause,pause 或
7、 pause(n),若想强行终止程序的运行，可以使用 Ctrl+c,其中 n 是延迟时间，以秒为单位；若缺省，则将暂停程序，直到用户按任意键后继续,pause off 屏蔽程序中所有 pause 的作用 pause on 打开 pause 的作用,2023年11月8日,21,选择结构 是根据给定的条件成立或不成立，分别执行不同的语句。Matlab 用于实现选择结构的语句有 if 语句和 switch 语句。,选择结构,2023年11月8日,22,单分支结构,if expression（条件）statements（语句组）end,if 条件语句,2023年11月8日,23,if expressi
8、on1（条件1）statements1（语句组1）elseif expression2（条件2）statements2（语句组2）.elseif expressionm（条件m）statementsm（语句组m）else statements（语句组）end,多分支结构,if 条件语句,2023年11月8日,24,if 语句举例,数论中的一个有趣问题：任取一个正整数，如果是偶数，用 2 除，如果是奇数，用 3 乘再加 1，反复这个过程，直到所得到的数为 1。问：是否存在使该过程永不中止的整数？(exp01.m),2023年11月8日,25,while 1 n=input(Please ente
9、r n(nonpositive quit):);m=n;if n1 if rem(n,2)=0 n=n/2;else n=3*n+1;end fprintf(n n=%d,n);end fprintf(n%d is not we need!continue.n,m);end,if 语句举例,2023年11月8日,26,根据表达式的不同取值，分别执行不同的语句。,switch expression（表达式）case value1（表达式1）statement1（语句组1）case value2（表达式2）statement2（语句组2）.case valuem（表达式m）statementm（语
10、句组m）otherwise statement（语句组）end,switch 语句,2023年11月8日,27,Matlab 首先计算 expression 的值，然后将它依次与各个 case 指令后的检测值进行比较，当比较结果为真时，就执行相应的语句组，然后跳出 switch 结构。,switch 后面的表达式 expression 的值可以是一个标量或字符串。,switch 语句,otherwise 指令可以不出现。,如果所有的比较结果都为假，则执行 otherwise 后面的语句组，然后跳出 switch 结构。,2023年11月8日,28,循环结构 是按照给定的条件，重复执行指定的语句
11、。Matlab 用于实现循环结构的语句有 for 语句和 while 语句,循环结构,2023年11月8日,29,for variable=expression statement（循环体）end,循环变量,表达式 expression 可以是行向量，也可以是矩阵,for 循环,clear;y=0;n=100;for k=1:n y=y+1/(2*k-1);end,2023年11月8日,30,循环语句可以嵌套使用,while 循环,while expression（条件）statement（循环体）end,不能在 for 循环体内改变循环变量的值,为了提高代码的运行效率，应尽可能提高代码的向量
12、化程度，避免 for 循环的使用,如果预先就知道循环的次数，则可以采用 for 循环；否则，如果预先无法确定循环的次数，则可以使用 while 循环。,2023年11月8日,31,这个例子给出了计算 eps 的一种方法。这里我们用大写 EPS，因此系统中 eps 的值不会被覆盖。当 EPS=eps 时，条件仍成立；而当 EPS=eps/2 时，条件为假，退出循环，所以最后应有 EPS=eps/2。num 用来记录循环次数。,计算 Matlab 中 eps 值(exp02.m),while 循环举例,num=0;EPS=1;while(1+EPS)1 EPS=EPS/2;num=num+1;en
13、d eps为浮点运算的相对精度,2023年11月8日,32,break 和 continue,break 语句用于终止循环的执行，即跳出最内层循环 continue 语句用于结束本次循环，进行下一次循环 break 和 continue 一般与 if 语句配合使用,其它流控制语句,2023年11月8日,33,首先由计算机随机产生一个 1,100 之间的一个整数，然后由用户猜测所产生的这个数。根据用户猜测的情况给出不同的提示，如果猜测的数大于产生的数，则显示“High”，小于则显示“Low”，等于则显示“You win！”，同时退出游戏。用户最多有 7 次机会。(exp03.m),程序示例,猜数
14、游戏,2023年11月8日,34,%猜数游戏clear;x=round(100*rand);%由计算机随机产生一个0,100的整数n=7;%有7次机会test=1;disp(欢迎参加猜数游戏！你共有,num2str(n),次机会);disp(请输入一个0到100之间的一个整数);for k=1:n guess=input(Enter your guess:);if guess x disp(high);else disp(You win!);test=0;break;end n=n-1;disp(你还有,num2str(n),次机会！);endif test=1 disp(You lost!)
15、end,程序示例,2023年11月8日,35,function 输出形参列表=函数名(输入形参列表)%注释说明部分(可选)函数体语句(必须),第一行为引导行，表示该 M文件是函数文件 函数名的命名规则与变量名相同(必须以字母开头)当输出形参多于一个时，用方括号括起来 函数必须是一个单独的 M文件 函数文件名必须与函数名一致 以百分号开始的语句为注释语句,函数文件,函数文件由 function 语句引导,2023年11月8日,36,function y=myfliplr(x)%returns X with row preserved and columns%flipped in the left
16、/right direction%X=1 2 3 becomes 3 2 1%4 5 6 6 5 4%num_col=size(x,2);%取 x 的列数y=x(:,num_col:-1:1);先输入矩阵x，而后写上myfliplr(x),函数文件举例,myfliplr.m,2023年11月8日,37,function printyh(n)%打印杨辉三角形,本函数没有输出参数yh=1;disp(yh);if n=1,return;endyh=1,1;disp(yh);for k=3:n yh_old=yh;k2=ceil(k/2);for l=2:k2 yh(l)=yh_old(l-1)+yh
17、_old(l);end yh(k2+1:k)=yh(k-k2:-1:1);disp(yh);end（在matlab中可输入printyh(5)）,函数文件举例,printyh.m,2023年11月8日,38,函数调用的一般格式,输出实参列表=函数名(输入实参列表),函数调用时，实参的顺序应与函数定义时的形参的顺序一致。函数可以嵌套调用，即一个函数可以被其它函数调用，甚至可以被它自身调用，此时称为递归调用。函数所传递的参数具有可调性，Matlab 用两个永久变量 nargin 和 nargout 分别记录调用该函数时的输入实参和输出实参的个数。,函数调用,2023年11月8日,39,%函数文件
18、factor.mfunction f=factor(n)if(n=1)f=1;else f=n*factor(n-1);end,%main.m%计算 s=1!+2!+3!+4!+5!%clear;s=0;for i=1:5 s=s+factor(i);endfprintf(s=%g n,s),例：利用函数的递归调用计算 n!,递归函数举例,2023年11月8日,40,%ex4argin.mfunction out=ex4argin(a,b)if(nargin=1)out=a;elseif(nargin=2)out=a*b;end,参数的可调性举例,%ex4argout.mfunction p,q=ex4argout(a,b)if(nargout=1)p=a+b;elseif(nargout=2)p=a+b;q=a-b;end,2023年11月8日,41,上机任务,2.建立M函数文件作以下计算:1)自然数n的阶乘;2)n中取m的组合。,（1）建立M脚本文件求1到100的和。（2）建立M函数文件求 并求n=50,100,150的值。,2023年11月8日,42,上机任务,3.,4.,5,
