欢迎来到淘文阁 - 分享文档赚钱的网站! | 帮助中心 好文档才是您的得力助手!
淘文阁 - 分享文档赚钱的网站
全部分类
  • 研究报告>
  • 管理文献>
  • 标准材料>
  • 技术资料>
  • 教育专区>
  • 应用文书>
  • 生活休闲>
  • 考试试题>
  • pptx模板>
  • 工商注册>
  • 期刊短文>
  • 图片设计>
  • ImageVerifierCode 换一换

    第3章MATLAB程序设计课件.ppt

    • 资源ID:92175583       资源大小:327.50KB        全文页数:73页
    • 资源格式: PPT        下载积分:20金币
    快捷下载 游客一键下载
    会员登录下载
    微信登录下载
    三方登录下载: 微信开放平台登录   QQ登录  
    二维码
    微信扫一扫登录
    下载资源需要20金币
    邮箱/手机:
    温馨提示:
    快捷下载时,用户名和密码都是您填写的邮箱或者手机号,方便查询和重复下载(系统自动生成)。
    如填写123,账号就是123,密码也是123。
    支付方式: 支付宝    微信支付   
    验证码:   换一换

     
    账号:
    密码:
    验证码:   换一换
      忘记密码?
        
    友情提示
    2、PDF文件下载后,可能会被浏览器默认打开,此种情况可以点击浏览器菜单,保存网页到桌面,就可以正常下载了。
    3、本站不支持迅雷下载,请使用电脑自带的IE浏览器,或者360浏览器、谷歌浏览器下载即可。
    4、本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰。
    5、试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓。

    第3章MATLAB程序设计课件.ppt

    第第3章章 MATLAB程序设计程序设计 3.1 M文件文件 3.2 程序控制结构程序控制结构 3.3 函数文件函数文件 3.4 程序举例程序举例 3.5 程序调试程序调试MATLAB是一种科学计算语言,但具有与是一种科学计算语言,但具有与BASIC、FORTRAN、C等传统高级语言类似的等传统高级语言类似的语言特征,可以方便地实现程序设计。语言特征,可以方便地实现程序设计。利用利用MATLAB的程序控制功能,将有关的的程序控制功能,将有关的MATLAB命令编成程序存储在一个文件中命令编成程序存储在一个文件中(称为称为M文件文件),然后运行该程序,然后运行该程序,MATLAB就会自动就会自动按照程序的要求,执行程序中的命令。以后需按照程序的要求,执行程序中的命令。以后需要这些命令时,只需再次运行该程序。这就是要这些命令时,只需再次运行该程序。这就是MATLAB程序设计。程序设计。充分利用充分利用MATLAB数据结构的特点,可使数据结构的特点,可使程序结构简单,提高编程效率。程序结构简单,提高编程效率。3.1 M文件文件3.1.1 M文件概述文件概述 M文件文件:用用MATLAB语言编写的程序。语言编写的程序。M文件根据调用方式的不同分为两类:文件根据调用方式的不同分为两类:1.命令命令(脚本脚本)文件文件(Script File)2.函数文件函数文件(Function File)命令文件与函数文件的主要区别:命令文件与函数文件的主要区别:1)命令文件没有输入参数,也不返回输出参数;命令文件没有输入参数,也不返回输出参数;而函数文件可以带输入参数,也可以返回输出而函数文件可以带输入参数,也可以返回输出参数。参数。2)命令文件以工作空间中的变量进行操作,文件命令文件以工作空间中的变量进行操作,文件中所有命令的执行结果也完全返回到工作空间中所有命令的执行结果也完全返回到工作空间中;而函数文件中定义的变量为局部变量,当中;而函数文件中定义的变量为局部变量,当函数文件执行完毕时,这些变量被清除。函数文件执行完毕时,这些变量被清除。3)命令文件可以直接运行;而函数文件不能直接命令文件可以直接运行;而函数文件不能直接运行,要以函数调用的方式运行。运行,要以函数调用的方式运行。例例3-1 分别建立命令文件和函数文件,将华氏温分别建立命令文件和函数文件,将华氏温度度f转换为摄氏温度转换为摄氏温度c。程序程序1:首先建立命令文件并以文件名首先建立命令文件并以文件名f2c1.m存盘。存盘。clear;f=input(Input Fahrenheit temperature:);c=5*(f-32)/9 然后在然后在MATLAB的命令窗口中输入的命令窗口中输入f2c1,将会执行该命令文件,执行情况为:将会执行该命令文件,执行情况为:Input Fahrenheit temperature:73c=22.7778程序程序2:首先建立函数文件首先建立函数文件f2c.m。function c=f2c(f)c=5*(f-32)/9 然后在然后在MATLAB的命令窗口调用该函数文件。的命令窗口调用该函数文件。clear;y=input(Input Fahrenheit temperature:);x=f2c(y)输出情况为:输出情况为:Input Fahrenheit temperature:70c=21.1111x=21.1111 3.1.2 M文件的建立与打开文件的建立与打开 M文件是一个文本文件,它可以用任何编辑程文件是一个文本文件,它可以用任何编辑程序来建立和编辑,而一般常用且最为方便的是使用序来建立和编辑,而一般常用且最为方便的是使用MATLAB提供的提供的文本编辑器文本编辑器。1建立新的建立新的M文件文件 为建立新的为建立新的M文件,启动文件,启动MATLAB文本编辑器有文本编辑器有3种方法:种方法:(1)菜单操作菜单操作。从。从MATLAB主窗口的主窗口的File菜单中选择菜单中选择New菜单项,再选择菜单项,再选择M-file命令,屏幕上将出现命令,屏幕上将出现MATLAB 文本编辑器窗口。文本编辑器窗口。(2)命令操作命令操作。在。在MATLAB命令窗口输入命令命令窗口输入命令edit,启动启动MATLAB文本编辑器后,输入文本编辑器后,输入M文件的内容文件的内容并存盘。并存盘。(3)命令按钮操作命令按钮操作。单击。单击MATLAB主窗口工具栏上主窗口工具栏上的的New M-File命令按钮,启动命令按钮,启动MATLAB文本编辑文本编辑器后,输入器后,输入M文件的内容并存盘。文件的内容并存盘。2打开已有的打开已有的M文件文件 打开已有的打开已有的M文件,也有文件,也有3种方法:种方法:(1)菜单操作菜单操作。从。从MATLAB主窗口的主窗口的File菜菜单中选择单中选择Open命令,则屏幕出现命令,则屏幕出现Open对话对话框,在框,在Open对话框中选中所需打开的对话框中选中所需打开的M文文件。在文档窗口可以对打开的件。在文档窗口可以对打开的M文件进行文件进行编辑修改,编辑完成后,将编辑修改,编辑完成后,将M文件存盘。文件存盘。(2)命令操作命令操作。在。在MATLAB命令窗口输入命命令窗口输入命令:令:edit 文件名文件名,则打开指定的,则打开指定的M文件。文件。(3)命令按钮操作命令按钮操作。单击。单击MATLAB主窗口工主窗口工具栏上的具栏上的Open File命令按钮,再从弹出的命令按钮,再从弹出的对话框中选择所需打开的对话框中选择所需打开的M文件。文件。3.2 程序控制结构程序控制结构程序的控制结构有程序的控制结构有3 3种:顺序结构、选种:顺序结构、选择结构和循环结构。择结构和循环结构。任何复杂的程序都可以由这任何复杂的程序都可以由这3 3种基本结种基本结构构成。构构成。MATLABMATLAB提供了实现控制结构的语提供了实现控制结构的语句,利用这些语句可以编写解决实际问题句,利用这些语句可以编写解决实际问题的程序。的程序。3.2.1 3.2.1 顺序结构顺序结构 按照程序中语句的排列顺序从上到下依次执行,直按照程序中语句的排列顺序从上到下依次执行,直到程序的最后一个语句。到程序的最后一个语句。1数据的输入数据的输入 从键盘输入数据,则可以使用从键盘输入数据,则可以使用input函数函数来进行,该函数的调用格式为:来进行,该函数的调用格式为:A=input(提示信息,选项提示信息,选项);其中提示信息为一个字符串,用于提示用其中提示信息为一个字符串,用于提示用户输入什么样的数据。户输入什么样的数据。如果在如果在input函数调用时采用函数调用时采用s选项,则允选项,则允许用户输入一个字符串。例如,想输入一个许用户输入一个字符串。例如,想输入一个人的姓名,可采用命令:人的姓名,可采用命令:xm=input(Whats your name?,s);2数据的输出数据的输出 MATLAB提供的命令窗口输出函数主要提供的命令窗口输出函数主要有有disp函数,其调用格式为函数,其调用格式为:disp(输出项输出项)其中输出项既可以为其中输出项既可以为字符串字符串,也可以为,也可以为矩阵矩阵。例如:。例如:A=Hello,Tom;disp(A)又如:又如:A=1 2 3;4 5 6;7 8 9;disp(A)例例3-2 输入输入x,y的值,并将它们的值互换后的值,并将它们的值互换后输出。输出。程序如下程序如下(example3_2):nclc,clearnx=input(Input x please.);ny=input(Input y please.);nz=x;nx=y;ny=z;ndisp(x);ndisp(y);例例3-3 求一元二次方程求一元二次方程ax2+bx+c=0的根的根。程序如下程序如下(example3_3):nclc,clearna=input(a=?);nb=input(b=?);nc=input(c=?);nd=b*b-4*a*c;nx=(-b+sqrt(d)/(2*a),(-b-sqrt(d)/(2*a);ndisp(x1=,num2str(x(1),x2=,num2str(x(2);3程序的暂停程序的暂停 暂停程序的执行可以使用暂停程序的执行可以使用pause函数,其函数,其调用格式为:调用格式为:pause(延迟秒数延迟秒数)如果省略延迟时间,直接使用如果省略延迟时间,直接使用pause,则,则将暂停程序,直到用户按任一键后程序继续将暂停程序,直到用户按任一键后程序继续执行。执行。若要强行中止程序的运行可使用若要强行中止程序的运行可使用:Ctrl+C3.2.2 选择结构选择结构 根据给定的条件成立或不成立,分别执行不根据给定的条件成立或不成立,分别执行不同的语句。有同的语句。有if、switch和和try三种语句。三种语句。1if语句语句 在在MATLAB中,中,if语句有语句有3种格式。种格式。(1)单分支单分支if语句:语句:if 条件条件 语句组语句组end 当条件成立时,则执行语句组,执行完之后继续当条件成立时,则执行语句组,执行完之后继续执行执行if语句的后继语句,若条件不成立,则直接语句的后继语句,若条件不成立,则直接执行执行if语句的后继语句。语句的后继语句。(2)双分支双分支if语句:语句:if 条件条件 语句组语句组1else 语句组语句组2end 当条件成立时,执行语句组当条件成立时,执行语句组1,否则执行语,否则执行语句组句组2,语句组,语句组1或语句组或语句组2执行后,再执行执行后,再执行if语句的后继语句。语句的后继语句。例例3-4 计算分段函数的值。计算分段函数的值。程序如下程序如下(example3_4):nclc,clearnx=input(请输入请输入x的值的值:);nif x=0n y=(x+sqrt(pi)/exp(2);nelsen y=log(x+sqrt(1+x*x)/2;nendny也可以用单分支也可以用单分支if语句实现语句实现:nclc,clearnx=input(请输入请输入x的值的值:);nif x0n y=log(x+sqrt(1+x*x)/2;nendny或者:或者:nclc,clearnx=input(请输入请输入x的值的值:);ny=(x+sqrt(pi)/exp(2);nif x0n y=log(x+sqrt(1+x*x)/2;nendny(3)多分支多分支if语句:语句:if 条件条件1 语句组语句组1 elseif 条件条件2 语句组语句组2 elseif 条件条件m 语句组语句组m else 语句组语句组n end语句用于实现多分支选择结构。语句用于实现多分支选择结构。例例3-5 输入一个字符,若为大写字母,则输出其对应输入一个字符,若为大写字母,则输出其对应的小写字母;若为小写字母,则输出其对应的大写的小写字母;若为小写字母,则输出其对应的大写字母;若为数字字符则输出其对应的数值,若为其字母;若为数字字符则输出其对应的数值,若为其他字符则原样输出。他字符则原样输出。nclc,clearnc=input(请输入一个字符请输入一个字符,s);nif c=A&c=a&c=0&c=9n disp(abs(c)-abs(0);nelsen disp(c);nend2switch语句语句 switch语句根据表达式的取值不同,分别执行不语句根据表达式的取值不同,分别执行不同的语句,其语句格式为:同的语句,其语句格式为:switch 表达式表达式 case 表达式表达式1 语句组语句组1 case 表达式表达式2 语句组语句组2 case 表达式表达式m 语句组语句组m otherwise 语句组语句组n end 当表达式的值等于表达式当表达式的值等于表达式1的值时,执行语的值时,执行语句组句组1;当表达式的值等于表达式当表达式的值等于表达式2的值时,执行语的值时,执行语句组句组2,;当表达式的值等于表达式当表达式的值等于表达式m的值时,执行的值时,执行语句组语句组m;当表达式的值不等于当表达式的值不等于case所列的表达式的所列的表达式的值时,执行语句组值时,执行语句组n;当任意一个分支的语句执行完后,直接执当任意一个分支的语句执行完后,直接执行行switch语句的下一句。语句的下一句。例例3-6 某商场对顾客所购买的商品实行打某商场对顾客所购买的商品实行打折销售,标准如下折销售,标准如下(商品价格用商品价格用price来表示来表示):price200 没有折扣没有折扣 200price500 3%折扣折扣 500price1000 5%折扣折扣 1000price2500 8%折扣折扣 2500price 0)n sumn mean=sum/cntnend3break语句和语句和continue语句语句 与循环结构相关的语句还有与循环结构相关的语句还有break语句和语句和continue语句。它们一般与语句。它们一般与if语句配合使用。语句配合使用。break语句用于终止循环的执行语句用于终止循环的执行。当在循环体内执。当在循环体内执行到该语句时,程序将跳出循环,继续执行循环语行到该语句时,程序将跳出循环,继续执行循环语句的下一语句。句的下一语句。continue语句控制跳过循环体中的某些语句语句控制跳过循环体中的某些语句。当在。当在循环体内执行到该语句时,程序将跳过循环体中所循环体内执行到该语句时,程序将跳过循环体中所有剩下的语句,继续下一次循环。有剩下的语句,继续下一次循环。例例3-12 求求100,200之间第一个能被之间第一个能被21整除的整数。整除的整数。程序如下:程序如下:nclc,clearnfor n=100:200n if rem(n,21)=0n continuen endnbreaknendnn4循环的嵌套循环的嵌套 如果一个循环结构的循环体又包括一个循环结构,就称如果一个循环结构的循环体又包括一个循环结构,就称为循环的嵌套,或称为多重循环结构。为循环的嵌套,或称为多重循环结构。例例3-13 若一个数等于它的各个真因子之和,则称该数为完若一个数等于它的各个真因子之和,则称该数为完数,如数,如6=1+2+3,所以,所以6是完数。求是完数。求1,500之间的全部完之间的全部完数。数。nclc,clearnfor m=1:500ns=0;nfor k=1:m/2n if rem(m,k)=0n s=s+k;n endnendnif m=sn disp(m);nendnend 3.3 函数文件函数文件 3.3.1 函数文件的基本结构函数文件的基本结构 函数文件由函数文件由function语句引导,其基本结构为:语句引导,其基本结构为:function 输出形参表输出形参表=函数名函数名(输入形参表输入形参表)注释说明部分注释说明部分 函数体语句函数体语句 其中以其中以function开头的一行为引导行,表示开头的一行为引导行,表示该该M文件是一个函数文件。函数名的命名规则与文件是一个函数文件。函数名的命名规则与变量名相同。输入形参为函数的输入参数,输出变量名相同。输入形参为函数的输入参数,输出形参为函数的输出参数。当输出形参多于一个,形参为函数的输出参数。当输出形参多于一个,则应该用方括号括起来。则应该用方括号括起来。例例3-14 编写函数文件求半径为编写函数文件求半径为r的圆的面积和周长。的圆的面积和周长。函数文件如下:函数文件如下:nfunction s,p=fcircle(r)n%CIRCLE calculate the area and perimeter of a circle of radii rn%r 圆的半径圆的半径n%s 圆的面积圆的面积n%p 圆的周长圆的周长n%2010年年4月月1日编写,浙江师范大学工学院日编写,浙江师范大学工学院.ns=pi*r*r;np=2*pi*r;在在MATLABMATLAB命令窗口调用该函数:命令窗口调用该函数:s,p=fcircle(10)s,p=fcircle(10)输出如下结果:输出如下结果:s=s=314.1593 314.1593p=p=62.8319 62.8319利用利用helphelp命令可以查询命令可以查询fcirclefcircle函数的注释函数的注释说明:说明:help fcircle help fcircle屏幕显示:CIRCLE calculate the area and perimeter of a circle of radii r r 圆的半径 s 圆的面积 p 圆的周长再用再用lookfor命令在第一注释行查询指定的关命令在第一注释行查询指定的关键词:键词:lookfor perimeter屏幕显示:屏幕显示:fcircle -CIRCLE calculate the area and perimeter of a circle of radii r 3.3.2 函数调用函数调用 函数调用的一般格式是:函数调用的一般格式是:输出实参表输出实参表=函数名函数名(输入实参表输入实参表)要注意的是,函数调用时各实参出现的要注意的是,函数调用时各实参出现的顺序、个数,应与函数定义时形参的顺序、顺序、个数,应与函数定义时形参的顺序、个数一致,否则会出错。个数一致,否则会出错。函数调用时,先将函数调用时,先将实参传递给相应的形参,从而实现参数传递,实参传递给相应的形参,从而实现参数传递,然后再执行函数的功能。然后再执行函数的功能。例例3-15 利用函数文件,实现直角坐标利用函数文件,实现直角坐标(x,y)与极坐标与极坐标(,)之间的转换。之间的转换。函数文件函数文件tran.m:nfunction rho,theta=tran(x,y)nrho=sqrt(x*x+y*y);ntheta=atan(y/x);调用调用tran.m的命令文件的命令文件main1.m:nclc,clearnx=input(Please input x=:);ny=input(Please input y=:);nrho,the=tran(x,y);nrhonthe 在在MATLAB中,函数可以嵌套中,函数可以嵌套调用,即一个函数可以调用别的函调用,即一个函数可以调用别的函数,甚至调用它自身。数,甚至调用它自身。一个函数调用它自身称为函数一个函数调用它自身称为函数的的递归调用。递归调用。例例3-16 利用函数的递归调用,求利用函数的递归调用,求n!。!。n!本身就是以递归的形式定义的。求本身就是以递归的形式定义的。求n!需需要求要求(n-1)!,这时可采用递归调用。,这时可采用递归调用。递归调用函数文件递归调用函数文件factor.m如下:如下:nfunction f=factor(n)nif n=1n f=1;nelsen f=factor(n-1)*n;%递归调用求递归调用求(n-1)!nend命令文件命令文件main2.m调用函数文件调用函数文件factor.m求:求:s=1!+2!+3!+4!+5!程序如下:程序如下:nclc,clearns=0;nfor i=1:5n s=s+factor(i);%调用调用factor.mnendns 3.3.3 函数参数的可调性函数参数的可调性 在调用函数时,在调用函数时,MATLAB用两个永久用两个永久变量变量nargin和和nargout分别记录调用该函分别记录调用该函数时的输入实参和输出实参的个数。只数时的输入实参和输出实参的个数。只要在函数文件中包含这两个变量,就可要在函数文件中包含这两个变量,就可以准确地知道该函数文件被调用时的输以准确地知道该函数文件被调用时的输入输出参数个数,从而决定函数如何进入输出参数个数,从而决定函数如何进行处理。行处理。例例3-17 nargin用法示例。用法示例。函数文件函数文件examp.m:nfunction fout=charray(a,b,c)nif nargin=1n fout=a;nelseif nargin=2n fout=a+b;nelseif nargin=3n fout=(a*b*c)/2;nend命令文件命令文件mydemo.m:nclc,clearnx=1:3;ny=1;2;3;nexamp(x)nexamp(x,y)nexamp(x,y,3)ans=1 2 3ans=2 4 6ans=21 3.3.4 全局变量与局部变量全局变量与局部变量在在Matlab中,函数文件中的变量是局部中,函数文件中的变量是局部的,与其他函数文件及的,与其他函数文件及Matlab工作空间相互工作空间相互隔离。隔离。如果在若干函数中,都把某一变量定义如果在若干函数中,都把某一变量定义为全局变量,那么这些函数将公用这一变量。为全局变量,那么这些函数将公用这一变量。全局变量用全局变量用global命令定义,格式为:命令定义,格式为:global 变量名变量名例例3-18 全局变量应用示例。全局变量应用示例。先建立函数文件先建立函数文件wadd.m,该函数将输入的参数,该函数将输入的参数加权相加。加权相加。nfunction f=wadd(x,y)nglobal ALPHA BETAnf=ALPHA*x+BETA*y;在命令窗口中输入:在命令窗口中输入:nclc,clearnglobal ALPHA BETAnALPHA=1;nBETA=2;ns=wadd(1,2)3.4 程序举例程序举例例例3-19 猜数游戏。首先由计算机产猜数游戏。首先由计算机产生生1,100之间的随机整数,然后由之间的随机整数,然后由用户猜测所产生的随机数。根据用用户猜测所产生的随机数。根据用户猜测的情况给出不同提示,如猜户猜测的情况给出不同提示,如猜测的数大于产生的数,则显示测的数大于产生的数,则显示“High”,小于则显示,小于则显示“Low”,等于,等于则显示则显示“You won”,同时退出游戏。,同时退出游戏。用户最多可以猜用户最多可以猜7次。次。n%Play the game of guessing the numbernx=fix(100*rand);%a random number calculated by the computernn=7;ntest=1;nfor k=1:7n numb=int2str(n);n disp(You have a right to,numb,guesses)n disp(A guess is a number between 0 and 100)n guess=input(Enter your guess:);n if guessxn disp(High)n elsen disp(You won)n test=0;n break;n endn n=n-1;nendnif test=1n disp(You lost)nend例例3-20 用筛选法求某自然数范围内的用筛选法求某自然数范围内的全部素数。全部素数。素数是大于素数是大于1,且除了,且除了1和它本身以外,不能被其他任何整和它本身以外,不能被其他任何整数所整除的整数。用筛选法求素数数所整除的整数。用筛选法求素数的基本思想是:要找出的基本思想是:要找出2m之间的之间的全部素数,首先在全部素数,首先在2m中划去中划去2的倍的倍数数(不包括不包括2),然后划去,然后划去3的倍数的倍数(不不包括包括3),由于,由于4已被划去,再找已被划去,再找5的的倍数倍数(不包括不包括5),直到再划去不,直到再划去不超过的数的倍数,剩下的数都超过的数的倍数,剩下的数都是素数。是素数。nclc,clearnm=input(m=);np=1:m;p(1)=0;nfor i=2:sqrt(m)n for j=2*i:i:mn p(j)=0;n endnendnn=find(p=0);np(n)从从p中划去中划去i的倍数(不包括的倍数(不包括i),可利用矩阵),可利用矩阵运算一步完成,得到更为简洁的程序:运算一步完成,得到更为简洁的程序:nclc,clearnm=input(m=);np=2:m;nfor i=2:sqrt(m)n n=find(rem(p,i)=0&p=i);n p(n)=;nendnp例例3-21 设,求设,求 求函数求函数f(x)在在a,b上的定积分,其几何意义就是上的定积分,其几何意义就是求曲线求曲线y=f(x)与直线与直线x=a,x=b,y=0所围成的曲边所围成的曲边梯形的面积。为了求得曲边梯形面积,先将积分区梯形的面积。为了求得曲边梯形面积,先将积分区间间a,b分成分成n等分,每个区间的宽度为等分,每个区间的宽度为h=(b-a)/n,对应地将曲边梯形分成对应地将曲边梯形分成n等分,每个小部分即是一等分,每个小部分即是一个小曲边梯形。近似求出每个小曲边梯形面积,然个小曲边梯形。近似求出每个小曲边梯形面积,然后将后将n个小曲边梯形的面积加起来,就得到总面积,个小曲边梯形的面积加起来,就得到总面积,即定积分的近似值。近似地求每个小曲边梯形的面即定积分的近似值。近似地求每个小曲边梯形的面积,常用的方法有:积,常用的方法有:矩形法矩形法、梯形法梯形法以及以及辛普生辛普生法法等。等。nclc,clearna=0;b=3*pi;nn=1000;h=(b-a)/n;nx=a;s=0;nf0=exp(-0.5*x)*sin(x+pi/6);nfor i=1:nn x=x+h;n f1=exp(-0.5*x)*sin(x+pi/6);n s=s+(f0+f1)*h/2;n f0=f1;nendns利用向量运算,可使程序更加简洁,更赋有利用向量运算,可使程序更加简洁,更赋有Matlab的特点。的特点。nclc,clearna=0;b=3*pi;nn=1000;h=(b-a)/n;nx=a:h:b;nf=exp(-0.5*x).*sin(x+pi/6);nfor i=1:n n s(i)=(f(i)+f(i+1)*h/2;nendns=sum(s)事实上,事实上,Matlab提供了有关数值积分的标准提供了有关数值积分的标准函数,实际应用中可以直接调用这些函数。函数,实际应用中可以直接调用这些函数。nclc,clearns=quad(exp(-0.5*x).*sin(x+pi/6),0,3*pi)例例3-22 Fibonacci数列定义如下:数列定义如下:f1=1f2=1fn=fn-1+fn-2 (n2)求求Fibonacci数列的第数列的第20项。项。建立函数文件建立函数文件ffib.mnfunction f=ffib(n)n%用于求用于求Fibonacci数列的函数文件数列的函数文件n%f=ffib(n)n%2004年年9月月30日编日编nif n2n f=ffib(n-1)+ffib(n-2);nelsen f=1;nend也可由也可由Matlab命令实现,程序如下:命令实现,程序如下:nclc,clearnf(1)=1;f(2)=1;nn=3;nwhile n0n E=E+F;n F=F*X/n;n n=n+1;nendnEnexpm(X)%调用调用MATLAB矩阵指数函数求矩阵指数矩阵指数函数求矩阵指数 3.5 程序调试程序调试3.5.1 程序调试概述程序调试概述 一般来说,应用程序的错误有两类一般来说,应用程序的错误有两类:一类是语法错误,另一类是运行时的错误一类是语法错误,另一类是运行时的错误。语法错误包括词法或文法的错误,例如函语法错误包括词法或文法的错误,例如函数名的拼写错、表达式书写错等。数名的拼写错、表达式书写错等。程序运行时的错误是指程序的运行结果有程序运行时的错误是指程序的运行结果有错误,这类错误也称为程序逻辑错误。错误,这类错误也称为程序逻辑错误。运行程序运行程序err1.mnclc,clearnA=87;nB=9.3;nC=A+*B;错误信息:错误信息:n?Error:File:F:课程教学课程教学-Matlab程序设程序设计教程计教程Matlab程序设计教程机械程序设计教程机械061第第3章程序章程序err1.m Line:5 Column:5nMissing variable or function.3.5.2 调试菜单调试菜单1Debug菜单项菜单项 该菜单项用于程序调试,需要与该菜单项用于程序调试,需要与Breakpoints菜菜单项配合使用。单项配合使用。2Breakpoints菜单项菜单项 该菜单项共有该菜单项共有6个菜单命令,前两个是用于在程个菜单命令,前两个是用于在程序中设置和清除断点的,后序中设置和清除断点的,后4个是设置停止条件的,个是设置停止条件的,用于临时停止用于临时停止M文件的执行,并给用户一个检查局文件的执行,并给用户一个检查局部变量的机会,相当于在部变量的机会,相当于在M文件指定的行号前加入文件指定的行号前加入了一个了一个keyboard命令。命令。例例3-24在在MATLAB编辑窗口,有一个判断编辑窗口,有一个判断m是否为完数的程序是否为完数的程序perfect.m,试设置断点,试设置断点来控制程序的执行。来控制程序的执行。m=input(请输入一个数:请输入一个数:);s=0;for k=1:m/2 if rem(m,k)=0 s=s+k;endend if m=s disp(m);endperfect.m程序:程序:3.5.3 调试命令调试命令 除了采用调试器调试程序外除了采用调试器调试程序外MATLAB还提供了一些命令用于程还提供了一些命令用于程序调试。命令的功能和调试器菜单序调试。命令的功能和调试器菜单命令类似,命令类似,表表3-1对这些命令进行了总结。对这些命令进行了总结。本章作业本章作业实验内容:实验内容:134思考练习:思考练习:245

    注意事项

    本文(第3章MATLAB程序设计课件.ppt)为本站会员(飞****2)主动上传,淘文阁 - 分享文档赚钱的网站仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知淘文阁 - 分享文档赚钱的网站(点击联系客服),我们立即给予删除!

    温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载不扣分。




    关于淘文阁 - 版权申诉 - 用户使用规则 - 积分规则 - 联系我们

    本站为文档C TO C交易模式,本站只提供存储空间、用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。本站仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知淘文阁网,我们立即给予删除!客服QQ:136780468 微信:18945177775 电话:18904686070

    工信部备案号:黑ICP备15003705号 © 2020-2023 www.taowenge.com 淘文阁 

    收起
    展开