《MATLAB及imulink应用》实验指导书+答案.doc

《MATLAB及应用》实验指导书 班 级： 姓 名： 学 号： 总评成绩： 汽车工程系 电测与汽车数字应用中心 目录 实MATLAB语言基础 1 实MATLAB科学计算及绘图 3 实MATLAB综合实例编程 5  实MATLAB语言基础 操作成绩 报告成绩 1实验目的 熟悉MATLAB的 掌握MATLAB的矩阵和数组的运算 掌握MATLAB符号表达式的创建熟悉符号方程的求解 2、 a.给出以下程序的伪代码描述：计算当米时的V和A，并同时绘制V相对于A的图形。 b.编写和运行部分a中所描述的程序。 解：clc,clear; r=0:0.1:3; V=4/3*pi*r.^3; A=4*pi*r.^2; plot(A,V);grid on; xlabel(Surface Area); ylabel(Volume); 在一个已知的测量矩阵T（100×100）中，删除整行全为0的行，删除整列内容全为0的整列（判断某列元素是否为0方法：检查T(: , i) .* (T(: , j))是否为0）。 假设汽车系在下列各年度的人口统计如表所示 年份 类 别 大一新生 学士毕业生 2002 102 82 2003 120 100 2004 98 94 2005 105 97 2006 121 110 试用一个二维矩阵STU表示上述数据，并请计算出下列各数值： 汽车系在2002~2006年之间的每年平均新生、毕业生人数。 5年来汽车系共有多少毕业生？ 在哪几年，新生数目比毕业生多？ 5年来每年的毕业生对新生的比例平均值为何？ 完成下列矩阵运算： 使用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为底的对数函数。 假设u=1和v=3,用 MATLAB 语句编写下列语句 a. b. c. d. 应用 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。看有何情况发生，为什么。 3思考题 MATLAB中，数组与矩阵在表示与应用上有哪些区别。 实MATLAB科学计算及绘图 操作成绩 报告成绩 1实验目的 熟悉MATLAB所提供的常用数值计算的函数(方程（组）的求解、插值、拟合)； 掌握MATLAB二维图形绘制命令及其图形控制（plot、loglog、contour、polar等）； 熟悉MATLAB三维图形绘制命令及其图形控制(mesh、surf等)。 2实验内容 说明：从下面的题目、第5~8章课后习题中任选10题，学号为奇数的选择奇数号题，学号为偶数的选择偶数号题，并在答案之前标注题目来源，例如“课本第7章题2”、“指导书题3”等。可以采用notebook格式输出结果（关于notebook请查阅其他参考资料，或上网搜索） 求f(x) = 4x6-x+x3-95在（0，100）上最大值、最小值，根。（使用函数fminbnd、roots） 求解常微分方程x’’=-x’+x+1，x’(