《MATLAB及用》实验指导书-mbook2008a.doc

《MATLAB及应用》实验指导书 班 级： 姓 名： 学 号： 总评成绩： 汽车工程系 电测与汽车数字应用中心 目录 实MATLAB语言基础 1 实MATLAB科学计算及绘图 3 实MATLAB综合实例编程 5  实MATLAB语言基础 操作成绩 报告成绩 实验目的 熟悉MATLAB的 掌握MATLAB的矩阵和数组的运算 掌握MATLAB符号表达式的创建熟悉符号方程的求解 randn产生一个（10×10）的矩阵A 计算B=（A+A’）/2。请注意，B一定是一个对称矩阵 计算矩阵B的特征向量e1 e2 … e10（使用函数eig） 验证在i不等于j的情况下，ei和ej的内积必定为0 下列 MATLAB 语句用于画出函数 在 [0, 10] 的值。 x = 0:0.1:10; y = 2*exp(-0.2*x); plot(x, y); 利用 MATLAB的M-文件编辑器创建一个新的 M-文件，输入上述代码，保存为文件test1.m。然后在命令窗口中输入 test1 执行这个文件。看得到什么结果？ 熟悉MATALB帮助的使用： 通过以下两种方式得到关于 exp函数的帮助 在命令窗口中输入 help exp 命令； 运用帮助空间窗口。 使用lookfor 命令查找出以10为底的对数函数，并计算1、10、1000的对数值。 假设u=1和v=3，用 MATLAB 符号计算功能计算下列表达式的值： a. b. c. d. a = 3, A = 4, b = a2, B = b2-1, c = a+A-2B, C = a+B+2c，求C。（使用符号函数solve） 求解常微分方程x’’=-x’+x+1，x’(0)=1，x(0)=0（使用符号函数dsolve） 应用 MATLAB 帮助窗口查找显示当前目录的命令。 MATLAB 启动时的当前目录是什么？ 创建一个新的目录mynewdir，并将其改变为当前目录。然后打开M-文件编辑窗口，增加以下语句： % create an input array from -2*pi to 2*pi t = -2*pi:pi/10:2*pi; % calculate|sin(t)| x=abs(sin(t)); %plot result plot(t, x); 把此文件以 test2.m 为文件名保存，然后在命令窗口中输入 test2 执行此文件。 关闭绘图窗口，返回到原始目录，然后在命令窗口中输入 test2。看有何情况发生，为什么。 思考题 MATLAB中，数组与矩阵在表示与应用上有哪些区别。 实MATLAB科学计算及绘图 操作成绩 报告成绩 实验目的 熟悉MATLAB所提供的常用数值计算的函数(方程（组）的求解、插值、拟合)； 掌握MATLAB二维图形绘制命令及其图形控制（plot、loglog、contour、polar等）； 熟悉MATLAB三维图形绘制命令及其图形控制(mesh、surf等)。 实验内容 求f(x) = 4x6-x+x3-95在（0，100）上最大值、最小值，根。（使用函数fminbnd、roots） 请用梯形法、辛普森法分别计算积分值 （trapz、quad） ，用两种方法求函数的根，并求其极值与零点。 计算二重积分 （使用函数dblquad） 求解常微分方程x’’=-x’+x+1，x’(0)=1，x(0)=0（使用函数ode45，需建立M函数） 矩阵M = [1, 2, 6; 4, 2, 7; 8, 9, 3]，求M的LU分解，QR分解，对角阵，特征值分解。（使用函数lu、qr、svd、eig） 建立一个5×5矩阵，分别用polyval函数和polyvalm函数将矩阵代入函数 f(x) = 1.35+0.0668x+0.436x2+0.695x3 ，计算出结果并比较。 已知t = an2 + bn，测得对应数据如下：（多项式插值interp1，多项式拟合polyfit） t = [0, 20, 40, 60, 80, 100, 120, 140, 160, 183.5]; n = [0, 1153, 2045, 2800, 3466, 4068, 4621, 5135, 5619, 6152]; 通过插值求解若 n = 0:5:6152时相应的 t值。 通过多项式拟合试求a和b的值。