数学实验_第二章

云南大学数学与统计学实验教学中心
实 验 报 告 课程名称：数学实验
学期：大二下 成绩： 指导教师：
学生姓名： 学生学号： 实验名称：M 文件编写
实验要求： 实验学时： 实验编号：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 ：用迭代法求a x =的值。

求平方根的迭代公式为： 112n n n a x x x +⎛⎫=+ ⎪⎝⎭
迭代的终止条件为前后两次求出的x 的差的绝对值小于10-5。

四、程序源代码：
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 k^3+j^3+i^3==100*k+10*j+i
w=100*k+10*j+i;
fprintf('w=%d\n', w); % fprintf 为matlab 中的输出函数同c 中的printf ，而input （'***'）也为输入命令，类似c 中get （）和scanf 。

end % if-end
end % for-end
end % 中for-end
end % 外for-end
实验结果：
w=153
w=370
w=371
w=407
法（2）%M-file,qiushuixianhuashu.m
for m=100:999
m1=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==m1^3+m2^3+m3^3 % 因为是数字，所以用^（幂运算符）
disp(m) % DISP(X) displays the array, without printing the array name.
end % if-end
end % for-end
实验结果：
153
370
371
407
2.编写函数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-end
x=xn;
i % 输出循环次数和所求结果
实验结果：
输入a=10;x=SQRT(a) 后回车
运行结果：
i =
3
x =
3.1623
五、实验过程
1.命令M文件的编写及调用：
1）在Matlab主菜单上选择菜单命令File→New→M-File；
2）在弹出的编辑框中编写源代码如上；
3）保存：编辑框中的菜单命令File→Save as ，文件名必须与源代码中函数名相同。

4）调用命令M文件：在command 对话框中直接输入源代码中函数名，然后回车即可得到答案。

2.函数M文件的编写及调用：
1）在Matlab主菜单上选择菜单命令File→New→M-File；
2）在弹出的编辑框中编写源代码如上；
3）保存：编辑框中的菜单命令File→Save ，或者直接点击保存按钮即可
六、实验总结
初步掌握函数M文件的应用，但还不能熟练应用。