matlab实验五报告

《matlab实验五报告》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验五MATLAB程序设计实验者：祝松年级：机设092学号：09405701002一、实验目的掌握MATLAB程序设计的主要方法，熟练编写MATLAB函数。二、实验内容(1)M文件的编辑。(2)程序流程控制结构。(3)子函数调用和参数传递。(4)局部变量和全局变量。三、实验步骤1.M文件的编辑选择MATLAB的菜单File

2、New

3、M-file，打开新的M文件进行编辑，然后输入以下内容，并保存文件名为exp1.m。%M脚本文件%功能：计算自然数列1～100的数列和s=0;forn=1:100s=s+n;ends保存好

4、文件后，在命令窗口输入exp1即可运行该脚本文件，注意观察变量空间。接着创建M函数文件，然后输入以下内容，并保存文件名为exp2.m。%这是M函数文件%功能：计算自然数列1～x的数列和functions=exp2(x)s=0;forn=1:xs=s+n;end保存好文件后，在命令窗口输入>>clear>>s=exp2(100)以open命令可以打开M文件进行修改>>openconv%打开conv函数2.程序流程控制结构1)for循环结构forn=1:10nend另一种形式的for循环：n=10:-1:5fori=n%

5、循环的次数为向量n的列数iend2)while循环结构在命令窗口输入：>>clear,clc;x=1;while1x=x*2end将会看到MATLAB进入死循环，因为while判断的值恒为真，这时须按下Ctrl＋C键来中断运行，并且可看到x的值为无穷大。练习：(1)请把exp2.m函数文件用while循环改写。(2)用π/4≈1-1/3+1/5-1/7+…公式求π的近似值，直到最后一项的绝对值小于为止，试编写其M脚本文件。3)if-else-end分支结构if-else-end分支有如下3种形式。(a)if表达式语句

6、组1end(b)if表达式语句组1else语句组2end(c)if表达式A语句组1elseif表达式B语句组2elseif语句组3……else语句组nend4)switch-case结构创建M脚本文件exp3.m，输入以下内容并在命令窗口中运行。%功能：判断键盘输入的数是奇数还是偶数n=input('n=');ifisempty(n)error('pleaseinputn')endswitchmod(n,2)case1A='奇数'case0A='偶数'end3.子函数和参数传递有一个函数g(x)=(x＝1,2,3…)

7、，试编写实现该函数的函数文件。functiong=exp4(x)%主函数g=0;forn=1:xg=g+fact(n);%调用子函数endfunctiony=fact(k)%子函数y=1;fori=1:ky=y*i;end输入参数可以由函数nargin计算，下面的例子sinplot2()，当只输入一个参数w时，sinplot2()函数会给p赋予默认值0。functiony=sinplot(w,p)ifnargin>2erro(‘toomanyinput’)endifnargin==1p=0;endx=linspace

8、(0,2*pi,500);z=sin(x.*w+p);练习：(1)编写求矩形面积函数rect，当没有输入参数时，显示提示信息；当只输入一个参数时，则以该参数作为正方形的边长计算其面积；当有两个参数时，则以这两个参数为长和宽计算其面积。(2)编写一个字符串加密函数nch=my_code(ch,x)，其中ch是字符串参数，x为整数；加密方法是：把ch的每一个字符的ASCII码值加上x，得到的即为加密后的新的字符串nch。由于可显示ASCII码值是有范围的(32，126)，因此当得到的ASCII码值大于126时，需要减去9

9、3。同理，再编写一个解码函数nch=my_dcode(ch,x)。提示：char(32:126)可获得ASCII码值为32～126的字符。4.局部变量和全局变量自程序执行开始到退出MATLAB，始终存放在工作空间，可被任何命令文件和数据文件存取或修改的变量即是全局变量，全局变量可用于函数之间传递参数，全局变量用关键字global声明。编写一个求和的函数文件，其名为summ.m。程序如下：functions=summglobalBEGENDk=BEG:END;s=sum(k);再编写调用M脚本文件use.m来调用sum

10、m.m函数文件，它们之间通过全局变量传递参数。程序如下：globalBEGENDBEG=1;END=10;s1=summ;BEG=1;END=20;s2=summ;