实验三 循环结构程序设计
一、实验目的
1. 掌握利用for语句实现循环结构的方法。 2. 掌握利用while语句实现循环结构的方法。 3. 熟悉利用向量运算来代替循环操作的方法。
二、实验环境
1．计算机
2．装有MATLAB6.0以上计算机一台
三 实验原理 1). 无条件循环
当需要无条件重复执行某些命令时，可以使用for循环： for 循环变量t=表达式1 : 达式2 : 表达式3 语句体 end
说明：表达式1为循环初值，表达式2为步长，表达式3为循环终值；当表达式2省略时则默认步长为1；for语句允许嵌套。
例1： 例2：矩阵输入程序 生成3×4阶的Hiltber矩阵。 m=input(‘矩阵行数：m=’); for i=1 : 3 n= input(‘矩阵列数：n=’); for j=1 : 4 for i=1:m H（i，j）=1/（i+j-1）； for j=1:n
end disp([‘输入第’,num2str(i),’行，第’,num2str(j),’ end 列元素’])
A(i, j) = input (‘ ’) end
end
2) while循环语句
while循环的一般使用形式为：
while 表达式 语句体 end 例3：
用迭代法计算多项式方程x3?2x?5= 0的一个实根。 解：
x0=3;r=1;n=0;
while r > 1e-5 x1=x0;
x0=3+5/x1.^2; r=abs(x0-x1); n=n+1; end,x0,n
运行结果为： x0 =
3.4260 n = 9
四 实验内容和步骤 （一）验证例题中的操作。
（二）完成以下各题，将答案写入实验报告，并在机器上运行结果验证。 1．编写程序，计算1+3+5+7+?+(2n+1)的值（用input语句输入n 值）。
2.根据
?26?1111，求?的近似值。当n分别取100、1 000、10 000时，?????122232n2结果是多少？
要求：分别用循环结构和向量运算（使用sum函数）来实现。 3.根据. y?1?111????，求： 352n?1（1）y<3时的最大n值。
（2）于（1）的n值对应的y值。
五 写出实验总结报告
实验四 函数与文件
一、实验目的
1. 掌握文件的建立、打开与关闭和文件的读写操作 2. 掌握定义和调用MATLAB函数的方法。 3. 掌握MATLAB文件的基本操作。
二、实验环境
1．计算机
2．装有MATLAB6.0以上计算机一台
三 实验原理 1. 建立M文件
将多个可执行的系统命令，用文本编辑器编辑后并存放在后缀为 .m 的文件．．．
中，若在MATLAB命令窗口中输入该m-文件的文件名（不跟后缀.m!），即可依次执行该文件中的多个命令。这个后缀为.m的文件，也称为Matlab的脚本文件(Script File)。
注意：文件存放路径必须在Matlab能搜索的范围内。 2. 建立函数文件
对于一些特殊用户函数，系统提供了一个用于创建用户函数的命令function，以备用户随时调用。
1．格式：
function [输出变量列表]=fun_name(输入变量列表)
用户自定义的函数体
2．函数文件名为：fun_name，注意：保存时文件名与函数名最好相同； ．．．．．．
3．存储路径：最好在系统的搜索路径上。 4. 调用方法：输出参量=fun_name (输入变量) 例：
计算s = n!，在文本编辑器中输入：
function s=pp(n); s=1;
for i=1:n s=s*i; end s;
在MATLAB命令窗口中输入：s=pp(5)
结果为: s = 120
四 实验内容和步骤
1. 定义一个函数文件，求给定复数的指数、对数、正弦和余弦，并在命令文件中调用该函数文件。
2、一个物理系统可用下列方程组来表示：
?m1cos??msin??1?0??0?m10m20?sin?cos??sin??cos?0??a1??0??a??mg?0???2???1? 0??N1??0??????1??N2??m2g?从键盘输入m1、m2和?的值，求a1、a2、N1、N2的值。其中g取 9.8,输入?时以角度为单位。
要求：定义一个求解线性方程组AX=B的函数文件，然后在命令文件中调用该函数文件。 3、设f?x??1?x?2?2?0.1?1?x?3?4?0.01，编写一个MATLAB函数文件fx.m,
使得调用f(x)时，x可用矩阵代入，得出的f(x)为同阶矩阵。
五 写出实验总结报告
实验五 MATLAB图形系统
一 实验目的
学会用MATLAB语言绘图 二 实验仪器和设备
装有MATLAB6.0以上计算机一台 三 实验原理 1、 plot
该命令在线性坐标系下绘制y对应于x的轨迹。若其中之一为矩阵则该命令将对应于矩阵的行或者列绘制一簇曲线。该命令也可以在同一坐标轴下绘制多条曲线。
例1 运行程序
t=[0:1:10];
x1=t; plot(x1) 可在图形窗口得到下面图形:
1098765432101234567891011
例2
x=[0:0.5:7]; y=sin(x); y1=2*x; y2=3*x;
plot(y,x,y1,x,y2,x)