第2章 MATALAB编程与作图VIP免费

E-mail:cxdzlh@163.comTel:49936260数学实验MATALAB编程与作图MATALAB控制流MATALAB作图M文件2.1M文件编写初步本质：就是将一些指令预先在M文件中编辑好，然后在需要时将M文件调出来执行即可。M文件的分类：一类为简单的脚本M文件；一类为函数M文件。♥脚本M文件就是就是一连串常用指令的组合，可在指令窗口中直接执行。♥函数M文件的开头一行需要有一个函数名称，然后用小括弧包括需要的输入参数；函数名称前可以有等号，在此等号的左方是它的输出参数MATLAB的内部函数是有限的，有时为了研究某一个函数的各种性态，需要为MATLAB定义新函数，为此必须编写函数文件.函数文件是文件名后缀为M的文件，这类文件的第一行必须是一特殊字符function开始，格式为：function因变量名=函数名（自变量名）函数值的获得必须通过具体的运算实现，并赋给因变量.1.函数M文件2121222121662,M1xxxxxxxxf文件，求解编写例▲单击MATALAB窗口里的［New-Mfile］图标▲输入代码：functionz=fun(x)z=x(1)*2+x(2)*2-2*x(1)*x(2)+6*x(1)-6*x(2)▲单击［File］－［Save］,把文件保存在［Work]文件夹下▲在指令窗口中输入如下指令即可看到运行结果。fun([2,3])例2定义函数f(x1,x2)=100(x2-x12)2+(1-x1)2functionf=fun(x)f=100*(x(2)-x(1)^2)^2+(1-x(1))^21.建立M文件：2.可以直接使用函数fun.m例如：计算f(1,2),只需在Matlab命令窗口键入命令：fun([12])2.2MATALAB程序设计初步MATALAB提供了四种控制程序流结构：◆for、while循环结构；◆if－else－end分支结构；◆switch－case结构；◆try－catch结构。1.for、while循环结构1.1for语句for循环：允许一组命令以固定的和预定的次数重复for变量=array（初始值：增量：终止值）运算指令（commands）end注1增量默认值为1，也可自己定义增量值，增量值可小于0。当增量值大于零时，程序将在变量大于终止值时终止；当增量小于零时，程序将在变量小于终止值时结束。的值求例n10sinx10,2,3,4.....,1n1nforn=1:10x(n)=sin(n*pi/10);endx单循环语句:的值求例-1nn12xx6,2,3,4.....n1,x2forn=2:6x(n)=2*x(n-1);endxx(1)=1;1000001)3121(S3iii求解级数求和问题例s=0;fori=1:100000s=s+1/2^i+1/3^i;ends多重嵌套循环语句:。阵对应的元素和构成的矩和求由例Hnm1,...5,n1,....5,m3form=1:5endHendH(m,n)=m+nforn=1:51.2while语句与for循环以固定次数求一组命令相反，while循环以不定的次数求一组语句的值.while表达式（expression）循环体{commands}end只要在表达式(expression)里的所有元素为真，就执行while和end语句之间的命令串{commands}，执行完成后再判断表达式是否为真，若不是则跳出循环体，向下继续执行。的值。的最小求出满足例m1000011miis=0;i=0;whiles<=10000iends=s+i;i=i+1;例2设银行年利率为11.25%。将10000元钱存入银行，问多长时间会连本带利翻一番？money=10000;years=0;whilemoney<20000years=years+1;money=money*(1+11.25/100);endyears1.3if－else－end分支结构1）有一个选择的一般形式是：if逻辑表达式（expression）执行语句{commands}end如果在表达式(expression)里的所有元素为真，就执行if和end语句之间的命令串{commands}.)1(),2(,1211)(12ffxxxxxf求设例先建立M文件f.m定义函数f（x），再在Matlab命令窗口输入f(2),f(-1)即可。functionf=f(x)Endifx<=1f=2*xendifx>1f=x^2+12)有多个选择的一般形式是：if逻辑表达式执行语句1elseif逻辑表达式执行语句2elseif逻辑表达式执行语句3elseif……else执行语句endend)1(),5.0(),2(010211232fff,xxxxxxf(x)求设例functionf=fun2(x)ifx>1f=x^2+1elseifx<=0f=x^3elsef=2*xendend1.3try－catch结构try语句用来检测程序代码是否会产生错误，一旦错误发生，MATALAB会立即跳入到相应的catch语句中去。trystatement1%命令块statement1被执行。若正确，则跳出此结构。catchstatement2%当命令块statement1出现执行错误则执行命令块则statement2。end输入矩阵A、B的值，执行上面的程序，当A的列数不等于B的行数时，程序...