M文件编写.doc

云南大学数学与统计学实验教学中心实验报告云南大学数学与统计学实验教学中心实 验 报 告课程名称：数学实验学期：大二下成绩：指导教师：朱娟萍学生姓名：学生学号：实验名称：M文件编写实验要求：用MATLAB 实验学时：实验编号：2实验日期： 完成日期：学院： 数学与统计学院专业 ： 年级： 一、实验目的熟悉MATLAB基本界面，熟悉M文件的概念和组成，学会一些简单的matlab语言编写，for循环，while循环，if-else-end结构和switch-case-end结构，达到对M文件的熟练掌握！二、实验环境Matlab7.0版本三、实验内容 P31 1.7.3 实验三：编写M文件1. 建立一个命令M文件：求所有的“水仙花数”，所谓“水仙花数”是指一个三位数，其各位数字的立方和等于该数本身。例如，153是一个水仙花数，因为153=13+53+33。2.编写函数M-文件SQRT.m：用迭代法求的值。求平方根的迭代公式为：迭代的终止条件为前后两次求出的x的差的绝对值小于10-5。 4、 程序源代码：1.命令M文件源程序：法（1） %M-file,qiushuixianhuashu.m % 函数名必须与文件名相同for k=1:1:9 % 循环变量k从1开始，到9结束，增量为1。k表示百位数 for j=0:1:9 for i=0:1:9 if k3+j3+i3=100*k+10*j+iw=100*k+10*j+i;fprintf(w=%dn, w); % fprintf为matlab中的输出函数同c中的printf，而input（*）也为输入命令，类似c中get（）和scanf。 end % if-end end % for-end end % 中for-endend % 外for-end法（2） %M-file,qiushuixianhuashu.m for m=100:999m1=fix(m/100); %求m的百位数字，命令fix为取整。m2=rem(fix(m/10),10); %求m的十位数字，命令rem为求余，具体格式为rem（a，b）表示a除以b求余m3=rem(m,10); %求m的个位数字if m=m13+m23+m33 % 因为是数字，所以用（幂运算符）disp(m) % DISP(X) displays the array, without printing the array name.end % if-endend % for-end2. 编写函数M文件源程序：function x=SQRT(a) % M文件必须包含function，指定函数名为SQRT，输入一个参数为a的值xn=1;xn1=2; % 初始值i=0; % 记录迭代次数while abs(xn1-xn)1e-5 % 循环条件xn1=(xn+a/xn)/2;xn=(xn1+a/xn1)/2;i=i+1;end % while-endx=xn;i,x % 输出循环次数和所求结果 五、实验过程1.命令M文件的编写及调用： 1）在Matlab主菜单上选择菜单命令FileNewM-File； 2）在弹出的编辑框中编写源代码如上； 3）保存：编辑框中的菜单命令FileSave as ，文件名必须与源代码中函数名相同。 4）调用命令M文件：在command 对话框中直接输入源代码中函数名，然后回车即可得到答案。2. 函数M文件的编写及调用： 1）在Matlab主菜单上选择菜单命令FileNewM-File； 2）在弹出的编辑框中编写源代码如