Matlab昆明理工大学习题.doc

一、 填空题：每空1分 共32分1、MATLAB是由英语Matrix（矩阵）和Laboratory（实 验室）两个单词的前三个字母的组合而成的。2MATLAB是以 向量 、 矩阵 为基本的数据单元，可以编写结构化、有面向对象特点的程序。3MATLAB是以 矩阵 作为数据操作的基本单位，语言灵活，功能强大。4变量名可以由任意的字母、数字或下划线组成，但必须以 字母 打头；变量名区分字母大小写，变量名最多不超过 19 个字符。5在MATLAB 中，主要有三种方式输入数据（1）直接生成数据 ； （2）使用 M 文件将文件中的数据调入内存 ；（3）用load 命令恢复工作空间中的变量 。6在MATLAB 中数据的输出方式有三种 （1）直接在终端上显示结果 ；（2）用diary 命令将控制窗口中的有关内容以文本形式存入一个指定的文件 ；（3）用 save 命令将工作空间中的变量保存到磁盘上 。7MATLAB是以 c 语言为平台，属于第四代计算机语言。优点是格式自由，语句简练，结构简单，学科性 工具箱 功能强大，可扩展性、可移植性强。缺点是 在命令窗口中解释一句，执行一句，所以运算速度慢。 8在MATLAB中保存整个工作空间，保存个别变量，文件的扩展名为 .mat 。9 数组 是MATLAB的基础。 矩阵运算 是MATLAB的核心。10MATLAB的构成：是由MATLAB的 语言 ；MATLAB的工具箱MATLAB的 工作环境 构成的。11在创建多行多列数组时 逗号 或 空格 ：用于分隔在同一行的元素,而 分号 用于分开不同的行。12多维数组的创建有 将二维数组扩展为多维数组。 和用函数创建多维数组 两种方法。13MATLAB 中三角函数只对弧度操作，所以自变量是角度的必须先转化为弧度，弧度= 角度量*pi/180 。14在一维插值中只有 三次样条 插值可以外插，其他几种都只能内插。15MATLAB中的语句序列由 表达式语句 语句、 赋值语句 语句、 控制语句 语句、 调试语句 和空语句组成。二、选择题：每空2分共18分1清除命令窗口中的文字（但不改变工作空间）用 B 命令。（A）home （B）clc （C）esc （D）del2在M文件的编写与管理中，显示或修改MATLAB搜索路径用 D 命令。（A）cd （B）cd path （C）path （D）p=cd3在MATLAB中，用 c 命令绘制简单的二维图形。（A）fplot （B）plot （C）plot3 （C）plot （D）surf4MATLAB中 D 函数用于牛顿科茨法计算函数的定积分。 （A）quad （B）trapz （C）leastq （D）quad85 MATLAB中用 A 函数求二维离散傅立叶变换。 （A）fft2 （B）conv2 （C）del2 （D）interp26MATLAB中用 C 函数实现一维多项式插值.（A）sun （B）prod （C）interp （D）diff7MATLAB中用 A 函数来完成线性拟合。 （A）polyfit （B）sort （C）mean （D）trapz8MATLAB中用 D 函数来数值计算函数微分。 （A）del2 （B）interp2 （C）polyfit （D）diff9MATLAB中用 D 函数计算卷积和多项式乘法。（A）quad （B）angle （C）filter （D）conv三、应用及编程1求解方程组（10分）A=5 4 3 2;6 5 8 3;3 7 9 5B=10 20 30X0=AB 求解check=A*X0 %检验解是否正确norm_x0=norm(X0) %计算解x0的范数X1=pinv(A)*B %用函数pinv求解线性方程组norm_x1=norm(X1)A = 5 4 3 2 6 5 8 3 3 7 9 5B = 10 20 30X0 = -1.0638 2.5532 1.7021 0check = 10.0000 20.0000 30.0000norm_x0 = 3.2477X1 = -0.7967 1.6331 1.5537 1.3950norm_x1 = 2.76802用梯形法求积分 （10分）f=(exp(-0.5*t).*sin(2*t+pi/6); %定义函数(exp(-0.5*t).*sin(2*t+pi/6)fplot(f,0,4*pi) %绘制曲线t1=0:pi/1000:3*pi; %定义积分范围0，3*pi和步长pi/100y1=(exp(-0.5*t1).*sin(2*t1+pi/6); %定义函数(exp(-0.5*t)*sin(2*t+pi/6)re=trapz(t1,y1) %求积分pause%积分范围不变，但改变步长t2=0:pi/2000:3*pi; %定义积分范围0，3*pi和步长3pi/40y2=(exp(-0.5*t2).*sin(2*t2+pi/6); %定义函数(exp(-0.5*t).*sin(2*t+pi/6)re=trapz(t2,y2)re = 0.4622re = 0.46223 一次化学动力学试验中，在95。15摄氏度乙醇溶液中，两种化合物反应的产物浓度与反应时间关系的原始数据为时间time2.55.07.510.013.017.020.030.040.050.060.070.0浓 度res0.290.560.771.051.361.522.002.272.813.053.253.56对其进行各种插值。（15分）程序：time=2.5 5.0 7.5 10.0 13.0 17.0 20.0 30.0 40.0 50.0 60.0 70.0%时间res=0.29 0.56 0.77 1.05 1.36 1.52 2.00 2.27 2.81 3.05 3.25 3.56%反应浓度timei=2.5:0.5:70; %时间插值范围及插值点resi1=interp1(time,res,timei,nearest); %用最近点插值subplot(2,2,1) %将图形窗口分成22个，激活第一个子图plot(time,res,*,timei,resi1) %绘制最近点插值曲线title(最近点插值)resi2=interp1(time,res,timei,linear); %用线性插值subplot(2,2,2) %激活第二个子图plot(time,res,+,timei,resi2) %绘制线性插值曲线title(线性插值)resi3=interp1(time,res,timei,spline); %用三次样条插值subplot(2,2,3) %激活第三个子图plot(time,res,timei,resi3) %绘制三次样条插值曲线title(三次样条插值)resi4=interp1(time,res,timei,cubic); %用三次多项式插值subplot(2,2,4) %激活第四个子图plot(time,res,x,timei,resi4) %绘制三次多项式插值曲线title(三次多项式插值)time = Columns 1 through 8 2.5000 5.0000 7.5000 10.0000 13.0000 17.0000 20.0000 30.0000 Columns 9 through 12 40.0000 50.0000 60.0000 70.0000res = Columns 1 through 8 0.2900 0.5600 0.7700 1.0500 1.3600 1.5200 2.0000 2.2700 Columns 9 through 12 2.8100 3.0500 3.2500 3.5600 结果：4在同一个窗口中的不同区域分别作在-10，10间的图像和作、在一个周期内的图像。（15分）（要求分别用红色+和绿色的*显示）程序：subplot(1,2,1) %将图形窗口分为12个，并激活第一个x1=(0:pi/1000:pi); %定义数组x1y1=sin(2*x1); %定义函数y1=sin(2*x1)x2=(0:pi/1000:pi); %定义数组x2y2=cos(2*x2); %定义函数y2=cos(2*x)plot(x1,y1,y*,x2,y2,g*)