当前位置:文档之家› Matlab昆明理工大学习题

Matlab昆明理工大学习题

一、填空题:每空1分共32分1、M A T L A B是由英语Matrix(矩阵)和Laboratory(实验室)两个单词的前三个字母的组合而成的。

2.M A T L A B是以向量、矩阵为基本的数据单元,可以编写结构化、有面向对象特点的程序。

3.M A T L A B是以矩阵作为数据操作的基本单位,语言灵活,功能强大。

4.变量名可以由任意的字母、数字或下划线组成,但必须以字母打头;变量名区分字母大小写,变量名最多不超过19 个字符。

5.在MATLAB 中,主要有三种方式输入数据(1)直接生成数据;(2)使用M 文件将文件中的数据调入存;(3)用load 命令恢复工作空间中的变量。

6.在MATLAB 中数据的输出方式有三种(1)直接在终端上显示结果;(2)用diary 命令将控制窗口中的有关容以文本形式存入一个指定的文件;(3)用save 命令将工作空间中的变量保存到磁盘上。

7.MATLAB是以c语言为平台,属于第四代计算机语言。

优点是格式自由,语句简练,结构简单,学科性工具箱功能强大,可扩展性、可移植性强。

缺点是在命令窗口中解释一句,执行一句,所以运算速度慢。

8.在MATLAB中保存整个工作空间,保存个别变量文件的扩展名为.m a t。

9.数组是MATLAB的基础。

矩阵运算是MATLAB的核心。

10.MATLAB的构成:是由MATLAB的语言;MATLAB的工具箱MATLAB的工作环境构成的。

11.在创建多行多列数组时逗号或空格:用于分隔在同一行的元素,而分号用于分开不同的行。

12.多维数组的创建有将二维数组扩展为多维数组。

和用函数创建多维数组两种方法。

13.MATLAB 中三角函数只对弧度操作,所以自变量是角度的必须先转化为弧度,弧度= 角度量*pi/180。

14.在一维插值中只有三次样条插值可以外插,其他几种都只能插。

15.MATLAB中的语句序列由表达式语句语句、赋值语句语句、控制语句语句、调试语句和空语句组成。

二、选择题:每空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)plot3 (C)plot (D)surf4.MATLAB中 D 函数用于牛顿—科茨法计算函数的定积分。

(A)quad (B)trapz (C)leastq (D)quad85. MATLAB 中用 A 函数求二维离散傅立叶变换。

(A )fft2 (B )conv2 (C )del2 (D )interp26.MATLAB 中用C 函数实现一维多项式插值.(A )sun (B )prod (C )interp (D )diff7.MATLAB 中用 A 函数来完成线性拟合。

(A )p o l y f i t (B )sort (C )mean (D )trapz8.MATLAB 中用D 函数来数值计算函数微分。

(A )del2 (B )interp2 (C )p o l y f i t (D )diff9.MATLAB 中用D 函数计算卷积和多项式乘法。

(A )quad (B )angle (C )f i l t e r (D )conv三、应用及编程 1.求解方程组(10分)543210658320379530x x x x x x x x x x x x +++=⎧⎪+++=⎨⎪+++=⎩A=[5 4 3 2;6 5 8 3;3 7 9 5] B=[10 20 30]'X0=A\B %求解check=A*X0 %检验解是否正确norm_x0=norm(X0) %计算解x0的数X1=pinv(A)*B %用函数pinv求解线性方程组norm_x1=norm(X1)A =5 4 3 26 5 8 33 7 9 5B =102030X0 =2.55321.70210 check =10.000020.000030.0000 norm_x0 =3.2477 X1 =-0.79671.63311.3950norm_x1 =2.76802.用梯形法求积分30.50sin/6t e t dt ππ-+⎰(2) (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摄氏度乙醇溶液中,两种化合物反应的产物浓度与反应时间关系的原始数据为时间time 2.5 5.0 7.510.13.17.20.30.40.50.60.70.浓度res 0.290.560.771.051.361.522.02.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) %将图形窗口分成2×2个,激活第一个子图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 82.5000 5.0000 7.5000 10.0000 13.0000 17.0000 20.0000 30.0000Columns 9 through 1240.0000 50.0000 60.0000 70.0000res =Columns 1 through 80.2900 0.5600 0.7700 1.0500 1.36001.52002.0000 2.2700Columns 9 through 122.81003.0500 3.2500 3.5600 结果:4.在同一个窗口中的不同区域分别作652-3xy在[-10,10]间的图像和作=x+cos xy=在一个周期的图像。

(15分)(2sin x)(2y=、)(要求分别用红色‘+’和绿色的‘*’显示)程序:subplot(1,2,1) %将图形窗口分为1×2个,并激活第一个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*') %同时绘制两条曲title('y=sin(2*x)与y2=cos(2*x2)的图形') %给图形加标题subplot(1,2,2) %激活第二个图形窗口x=(-10:0.1:10);y=5*x.^2+3*x-6;plot(x,y,'r+') %在新建的图形窗口中绘y=5*x.^2+3*x-6 图title('y=5*x.^2+3*x-6的图形') %给图形加标题结果:图表。

相关主题