2013 MATLAB实验六.docVIP

实验六 Matlab程序设计 一、实验目的： 1、熟练掌握M文件的结构。 2、熟练掌握MATLAB的程序流程控制结构。 3、学会编写基本程序。 二、实验内容与步骤： 用户如想灵活应用matlab去解决实际问题，充分调用matlab的科学技术资源，就需要编辑m文件 包含matlab语言代码的文件称为m文件，其扩展名为m。 编辑m文件可使用各种文本编辑器。 1.M文件的创建方法： 在matlab命令窗口点击file 菜单new ( m-file，可以打开空白的文件编辑器，也可以通过打开已有的M文件来打开文件编辑器。如图所示，为打开已建的M文件。 2. M文件的编写 2.1 M脚本文件 （1）在MATLAB桌面打开M文件编辑器； （2）将以下命令全部写入编辑器，为了能标志该文件名称，在第一行写入包含文件名的注释。 %exam01二阶系统时域响应曲线 x=0:0.1:20; y1=1-1/sqrt(1-0.3^2)*exp(-0.3*x).*sin(sqrt(1-0.3^2)*x+acos(0.3)) plot(x,y1,r) hold on y2=1-1/sqrt(1-0.707^2)*exp(-0.707*x).*sin(sqrt(1-0.707^2)*x+acos(0.707)) plot(x,y2,g) y3=1-exp(-x).*(1+x) plot(x,y3,b) 编写好之后，将文件存在当前路径下，文件名exam01.m，在MATLAB命令窗口执行命令： exam01 2.2 M函数文件 MATLAB函数文件（function file） function 返回变量=函数名（输入变量） 注释说明语句段 程序语句段 特定规则： (1)函数m文件第一行必须以单词function作为引导词，必须遵循如下形式： function 因变量=函数名(自变量) (2)m文件的文件名必须是函数名 .m。 (3)程序中的变量均为局部变量，不保存在工作空间中。其变量只在函数运行期间有效。 创建M函数文件并调用的步骤如下： （1）编写函数代码 function y=exam02(zeta) %exam02 x=0:0.1:20; if(zeta0)(zeta1) y=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt(1-zeta^2)*x+acos(zeta)); elseif zeta==1 y=1-exp(-x).*(1+x); end plot(x,y) （2）在M文件编辑窗口编写上述程序，保存为exam02.m文件 （3在MATLAB命令窗口输入： zeta=0.707； exam02（zeta） 3.MATALB程序设计 matlab语言的程序结构与其它高级语言是一致的，分为顺序结构，循环结构，分支结构。 (1).顺序结构 —— 依次顺序执行程序的各条语句 (2).循环结构 —— 被重复执行的一组语句，循环是计算机解决问题的主要手段。 循环语句主要有： for — end 语法为： for i = 表达式 （表达式为一个向量） 可执行语句 end 利用for循环求1！+2！+3！+ …… +5！的值( sum=0; for i=1:5 pdr=1; for k=1:i pdr=pdr*k; end sum=sum+pdr; end 在M文件编辑窗口编写上述程序，保存为exam03.m文件；在MATLAB命令窗口执行命令： exam03 sum for循环语句计算sum=的程序 sum=0; for n=1:10 sum = n^n + sum; end 在M文件编辑窗口编写上述程序，保存为exam04.m文件；在MATLAB命令窗口执行命令： exam04 sum while — end 循环 while循环将循环体中的语句循环执行不定次数。 语法为：while 表达式 循环体语句 end 表达式一般是由逻辑运算和关系运算以及一般运算组成的，以判断循环的进行和停止；只要表达式的值 非0，继续循环；直到表达式值为0，循环停止。 用while循环求1~100间整数的和 sum=0; i=1; while i=100 s