MATLAB实验指导书.docVIP

《MATLAB语言及应用》 实验指导书 目 录 实验一 矩阵的运算 1 实验二 选择结构和循环结构程序设计 3 实验三 函数与文件 5 实验四 线性代数中的数值计算问题 7 实验五 数据处理和多项式计算 9 实验六 数值微积分与方程数值求解 10 实验七 绘图功能 11 实验八 图形用户界面设计 12 实验九 Simulink的应用 14 实验十 综合实验 15  实验一 矩阵的运算 一、实验目的 熟悉启动和退出MATLAB的方法。 熟悉MATLAB的桌面环境。 掌握建立矩阵的方法。 掌握MATLAB各种表达式的书写规则以及常用函数的使用。 二、实验内容 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。 已知 求下列表达式的值： （1）和（其中I为单位矩阵） （2）和 （3）A^3和A.^3 （4）A/B和B\A （5）[A，B]和[A([1,3]，：）；B^2] （1）求它们的乘积C。 （2）将矩阵C的右下角3x2子矩阵赋给D。 （3）查看MATLAB工作空间的使用情况。 完成下列操作： （1）求[100，999]之间能被21整除的数的个数。 （2）将101~125这25个数填入一个5行5列的表格中，使其每行每列及对角线的和均为565。 三、实验报告 实验程序（命令）清单 运行结果 实验收获与体会 四、实验学时 3学时 实验二 选择结构和循环结构程序设计 一、实验目的 掌握建立和执行M文件的方法。 掌握利用if语句或switch语句实现选择结构的方法。 掌握利用for语句或while语句实现循环结构的方法。 熟悉利用向量运算来代替循环操作的方法。 二、实验内容 假定某地区电话收费标准为：通话时间在3分钟以下，收费0.50元；3分钟以上，则每超过1分钟加收0.15元；在7:00~22:00之间通话者，按上述收费标准全价收费，在其他时间通话者，按上述收费标准半价收费。计算某人在t1时间通话至t2时间，应缴多少电话费。 提示： （1）t1、t2t2-t1,相减时可以将t1、t2（2） 求?1~?100中： （1）最大值、最小值、各数之和。 （2）正数、零、负数的个数。 提示：可以考虑使用MATLAB有关函数来实现。 三、实验报告 1. 实验程序（命令）清单 2. 运行结果 3. 实验收获与体会 四、实验学时 3学时 实验三 函数与文件 一、实验目的 掌握定义和调用MATLAB函数的方法。 掌握MATLAB文件的基本操作。 二、实验内容 定义一个求线性方程组AX=B的解的函数文件，然后在命令文件中调用该函数文件，求如下的方程组的解。 （1）（2） （3）提示：求AX=B的特解利用x=A\B。 rank(A) 函数用来求矩阵A的秩。 norm(V) 函数用来求向量V的2-范数，即向量的模。 null(A,’r’) 函数可求得满足AX=0的解空间，实际上是求出解空间的一组基（基础解系）。 x=0:0.1:1时,y=exp(x)的值,把结果[x;y]按指定的格式即浮点型,每行两个数据，写入到文件myfile.txt中，并输出文件的大小。 三、实验报告 1. 实验程序（命令）清单 2. 运行结果 3. 实验收获与体会 四、实验学时 3学时 实验四 线性代数中的数值计算问题 一、实验目的 掌握生成特殊矩阵的方法。 掌握矩阵分析的方法。 掌握线性方程组的求解方法。 实验内容 产生5阶希尔伯特矩阵H和5阶帕斯卡矩阵P，且求其行列式的值Hh和Hp以及它们的条件数Th和Tp，判断哪个矩阵性能更好，为什么？ 已知 求A的特征值及特征向量，并分析其数学意义。 分别用3种不同的数值方法解线性方程组。 不用rot90函数，实现方阵左旋90°或右旋90°的功能。例如，原矩阵为A，A左旋后得到B，右旋后得到C 要求：定义一个函数文件，然后在命令文件中调用该函数文件。函数文件的调用格式为myrot(A,t)，其中A在调用前为原矩阵，调用后为经旋转后的矩阵，t确定是左旋还是右旋，当t=1时为左旋，t=-1时为右旋，t为其他值时显示出错信息。 提示：先将A转置，再作上下翻转，则完成左旋90°；若将A转置后作左右翻转，则完成右旋90°。 实验五 数据处理和多项式计算 一、实验目的 掌握数据插值与曲线拟合的方法及其应用。 掌握多项式的常用运算。 6：00~18：00之间每隔2h的室内外温度（℃）如表1所示。 表1室内外温度观测结果（℃） 时间h 6 8 10 12 14 16 18 室内温度t1 18.0 20.0 22.0 25.0 30.0 28.0 24.0 室外温度t2 15.0 19.0 24.0 28.