《2014秋2012级matlab程序设计实验报告三 - 副本》由会员分享，可在线阅读，更多相关《2014秋2012级matlab程序设计实验报告三 - 副本（10页珍藏版）》请在金锄头文库上搜索。
1、实验三 MATLAB 程序设计一、实验目的1、掌握建立和执行 M 文件的方法。2、掌握选择结构、多分支选择结构程序的编程方法。 3、掌握多种循环结构程序的编程方法。4、掌握定义函数文件和调用函数文件的方法。二、实验内容1、MATLAB M 文件的创建与使用。(1) 建立自己工作目录，如 D:MATLAB。在 File 菜单中，单击 Set Path 菜单项，在弹出对话框中建立目标文件夹，：单击 Add Folder 按钮，将自己的工作目录 D:MATLAB 添加到搜索目录之中，按 Save 按钮保存，如图 1。图 1 工作目录设置(2) 在 M 文件编辑器中编制命令行文件。在 File 菜单中
2、，单击 NewM-File，弹出 M 文件编辑器，如图 2：图 2 编辑器在编辑器中输入如下命令程序：%求小于 2000 且为 2 的整数次幂的正整数。f(1)=2;k=1;while f(k)3-x+6endend3、编写一个函数 M 文件，返回向量 x 中元素的和、平均值、最大值、最小值、均方根值。其中方均根值(Root Mean Square)的计算公式为： 21Nirms并用下面数据测试你写的函数：(1) x=sin(0:0.01:6*pi)(2) x=rand(1,200)，得到的 x 为 200 个(0,1)之间均匀分布的随机数。(3) function fun(x)(4) %UN
3、TITLED Summary of this function goes here(5) % Detailed explanation goes hereg(6) sum(x)(7) mean(x)(8) max(x)(9) min(x)(10) sqrt(meansqr(x)(11) end4、编写一个函数 M 文件，若给出一个向量 ，函数返回如12,nxxL下范德蒙矩阵。 122112nnnxxLL例如在命令窗口输入v =myvander(2 3 4 5)，得 ，同时1234596287v生成一些数据测函数。function fun4_vander(A)%UNTITLED Summary
4、of this function goes here% Detailed explanation goes hererot90(vander(A)end5、分别使用 while 和 for 编写命令文件，用循环语句编程计算 s 值:10322()ns同时对编写的文件设置断点，观察变量取值的变化情况。function fun5()%UNTITLED Summary of this function goes here% Detailed explanation goes heren=2;s=0;while n0&rem(x,1)=0disp(x)while x=1if rem(x,2)=1x=x
5、*3+1;disp(-) disp(x)endx=x/2;disp(-) disp(x)endendif rem(x,1)=0|x=190disp()endend8、Fibonacci 数组的元素满足 Fibonacci 规则：ak+2=ak+ak+1,（k=1 ，2，3，） ；且 a1= a2=1。请设计一段程序，求出该数组中第一个大于 10000 的元素要求显示的结果为：i=21a(i)=10946或者a(21)=10946function fun8()%UNTITLED4 Summary of this function goes here% Detailed explanation g
6、oes herea=1,1;i=2;while a(i)10100，并求出n!”function fun10(n)%UNTITLED2 Summary of this function goes here% Detailed explanation goes herem=n;while factorial(n)10100n=n-1;endnfor i=1:mif factorial(i)10100disp(i-1)breakendendend11、下列程序用来判断一个人的体温是否处于危险状态。调试程序是否正确，如果程序错误指出错误在哪里？并写出正确答案。temp=input(请输入体温 : t
7、emp = );if temp 36.5disp(体温正常);elseif temp 38.0disp(体温偏高！);elseif temp 39disp(体温高！);endfunction fun11()%UNTITLED3 Summary of this function goes here% Detailed explanation goes heretemp=input(: temp = );if temp 39disp();elseif temp 38.0disp();elseif temp 36.5disp();endend12、使用嵌套 for 循环命令创建下列矩阵：510015Afunction fun12()%UNTITLED4 Summary of this function goes here% Detailed explanation goes herea=5,1,0,0,0;toeplitz(a)b=zeros(5);for i=1:5b(i,i)=5;if(i+1)0b(i,i-1)=1;endenddisp(b)end
