《MATLAB语言与应用》实验测验高成。.doc

MATLAB语言及其应用实验报告.doc一、实验目的1、熟悉MATLAB语言及其基本操作；2、掌握利用MATLAB进行数据分析和可视化的方法；3、掌握MATLAB应用于科学计算的基本方法。

二、实验环境MATLAB R2018a三、实验内容及方法3.1 实验内容1、打开MATLAB环境；2、读入数据文件，并对数据进行分析和处理；3、通过MATLAB绘制图表，对数据进行可视化。

四、实验细节及流程4.1 读入数据文件本实验使用的数据文件为sales.xlsx，其中包括2018年各个季度的销售数据。

首先，通过以下命令导入数据：data=xlsread('sales.xlsx');4.2 数据分析和处理1、计算各季度销售总额通过以下代码计算每个季度的销售额并求和，得到每年的销售总额：Q1=data(:,2);TotalSales=sum([Q1 Q2 Q3 Q4],2);2、计算增长率根据每年的销售总额，计算出每年的增长率。

具体代码如下：GrowthRate(1)=0;for i=2:length(TotalSales)GrowthRate(i)=((TotalSales(i)-TotalSales(i-1))/TotalSales(i-1))*100; endSalesGrowth=[TotalSales GrowthRate];3、计算每个季度的均值和标准差meanQ1=mean(Q1);stdQ1=std(Q1);4、计算出每年第一个季度的销售额所占比例首先，我们将第一个季度的销售额单独提出来，具体代码如下：4.3 数据可视化1、柱形图对于销售总额，使用柱形图进行可视化，具体代码如下：结果如下图所示：图1 销售总额2、线性图3、箱形图boxplot([Q1 Q2 Q3 Q4],{'Q1','Q2','Q3','Q4'});图3 每个季度的销售额4、饼图pie(FirstQSalesRatio(:,1));五、结论本实验通过对销售数据的分析和可视化，得出以下结论：1、2018年销售总额呈逐年上升趋势，其中2017年到2018年的增长率最高；2、每年第一个季度的销售额所占比例在40%至45%之间，与其他季度相比，显著高于其他季度；3、2018年第二季度的销售额经过调整后，表现出了相对较高的波动。

《MATLAB 语言与应用》实验课程任务书一、 实验教学目标与基本要求上机实验是本课程重要的实践教学环节；实验的目的不仅仅是验证理论知识，更重要的是通过上机实验，加强学生的实验手段与实践技能，掌握应用MATLAB 语言求解问题的方法，培养学生分析问题、解决问题、应用知识的能力和创新精神，全面提高学生的综合素质。

上机实验共8学时。

主要实验内容是基于理论课所学知识对课后典型习题进行MATLAB 求解，基本掌握常见数学问题的求解方法与命令调用，更深入地认识和了解MATLAB 语言强大的计算功能。

上机实验最终以书面报告的形式提交，并作为期末成绩考核内容的一部分。

二、 实验内容（8学时）第一部分MATLAB 语言编程、科学绘图与基本数学问题求解（4学时）主要内容：掌握MATLAB 语言编程基础、科学绘图方法、微积分问题、线性代数问题等基本数学问题的求解与应用。

练习题：1、安装MATLAB 软件，应用demo 命令了解主要功能，熟悉基本功能，会用help 命令。

2、用MATLAB 语句输入矩阵A 和B⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=1423143212344321A ， ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡++++++++++++++++=4j 11j43j 22j34j 11j 42j 33j 24j 13j 22j 31j 41j 42j 33j 24j 1B 前面给出的是44⨯矩阵，如果给出5)6,5(=A 命令将得出什么结果？ 代码：A=[1 2 3 4;4 3 2 1;2 3 4 1;3 2 4 1]B=[1+4j 2+3j 3+2j 4+1j;4+1j 3+2j 2+3j 1+4j;2+3j 3+2j 4+1j 1+4j;3+2j 2+3j 4+1j 1+4j] A(6,5)=53、假设已知矩阵A，试给出相应的MATLAB命令，将其全部偶数行提取出来，赋给B矩阵，用magic(8)A 命令生成A矩阵，用上述命令检验一下结果是不是正确。

第二章 MATLAB 语言及应用实验项目实验一 MATLAB 数值计算三、实验内容与步骤1.创建矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a（1（2）用（3）用（42.矩阵的运算（1）利用矩阵除法解线性方程组。

⎪⎪⎩⎪⎪⎨⎧=+++=-+-=+++=+-12224732258232432143214321421x x x x x x x x x x x x x x x 将方程表示为AX=B ，计算X=A\B 。

（2）利用矩阵的基本运算求解矩阵方程。

已知矩阵A 和B 满足关系式A -1BA=6A+BA ，计算矩阵B 。

其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=7/10004/10003/1A ，Ps: format rata=[1/3 0 0;0 1/4 0;0 0 1/7];b=inv(a)*inv(inv(a)-eye(3))*6*a（3）计算矩阵的特征值和特征向量。

已知矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=1104152021X ，计算其特征值和特征向量。

（4）Page:322利用数学函数进行矩阵运算。

已知传递函数G(s)=1/(2s+1),计算幅频特性Lw=-20lg(1)2(2w )和相频特性Fw=-arctan(2w),w 的范围为[0.01，10]，按对数均匀分布。

3.多项式的运算（1）多项式的运算。

已知表达式G(x)=(x-4)(x+5)(x 2-6x+9),展开多项式形式，并计算当x 在[0，20]内变化时G(x)的值，计算出G(x)=0的根。

Page 324（2）多项式的拟合与插值。

将多项式G(x)=x 4-5x 3-17x 2+129x-180，当x 在[0，20]多项式的值上下加上随机数的偏差构成y1，对y1进行拟合。

对G(x)和y1分别进行插值，计算在5.5处的值。

Page 325 四、思考练习题1.使用logspace 函数创建0～4π的行向量，有20个元素，查看其元素分布情况。

Ps: logspace(log10(0),log10(4*pi),20) (2) sort(c,2) %顺序排列 3.1多项式1）f(x)=2x 2+3x+5x+8用向量表示该多项式，并计算f(10)值. 2）根据多项式的根[-0.5 -3+4i -3-4i]创建多项式。

一．填空题：(每空1分,共20分)1. MATLAB的通用命令中列出内存中的变量目录可用whos 命令；若对函数功能和用法不熟悉，可用help 命令；显示变量、阵列或文本时，可用disp 命令；清除工作空间中的变量时，可用clear 命令；保存和恢复工作空间，或者工作空间中的特定变量时，可用save 和load 命令。

2．在MATLAB图形系统中，如果需要给已经绘制好的图形加上栅格线，可以使用grid on 命令，如果需要在指定的坐标位置标注文本时使用的函数为text ，需要绘制多条曲线，保持当前图形窗口可使用hold on 命令。

3．用MATLAB绘制条形图时可使用bar 函数，绘制柱状图，用以表示数据值的分布情况时可使用hist 函数。

4．MATLAB中用于控制不定次重复循环语句为while ，若在循环执行过程中需要终止该循环时采用的语句为break 。

5．对于MATLAB中的多维阵列，在保持所有元素个数和内容不变的前提下可以使用reshape 函数改变其尺寸和维数；利用cell 函数可以定义单元阵列。

6. 在MATLAB中，可用内部函数与工具箱函数产生矩阵，其中如果要产生单位阵可用eye 函数，产生对角阵可用diag 函数，产生均匀分布的随机阵可用rand 函数。

7. 在MATLAB中，将含有小数的数值转换成字符串可用num2str 命令，在字符串中查找子字符串可用findstr 命令。

二．选择题（每题2分，共20分，在每题的四个备选答案中选择一个正确的答案。

）1. 设A=[2 4 3; 5 3 1; 3 6 7]，则sum(A)，length(A)和size(A)的结果( D )A. [10 13 11] 9 [3 3]B. [9 9 16] 3 [3 3]C. [9 9 16] 9 [3 3]D. [10 13 11] 3 [3 3]2. 下列关于脚本文件和函数文件的描述中不正确的是( B )A. 去掉函数文件第一行的定义行可转变成脚本文件;B. 函数文件可以在命令窗口直接运行;C. 脚本文件可以调用函数文件;D. 函数文件中的第一行必须以function开始;3. 在Command Window窗口中分别输入下列命令，对应输出结果错误的是( C )A. x=[-3:2] x=[-3 -2 1 0 1 2]B. x=zeros(1,2);x>0 ans=[0 0]C. y=diag(eye(3),2).’ y=[0 0]D. 3-2*rand(1,2) ans=[1.0997 2.5377]4. 对于矩阵B，统计其中大于A的元素个数，可以使用的语句是( B )A. length(B) - length(find(B<=A))B. sum(sum(B>A))C. length(sum(B>A))D. sum(length(B>A))5. 已知str1=’find’, str2=’fund’, str3=’I like you’,有：k1=sum(str1==str2), k2=sum(strr ep(str1,’i’,’u’)==str2), k3=findstr(str3,’ke’), 则k1，k2，k3的结果分别为：( B )A. 3, 3, 5B. 3, 4, 5C. 4, 3, 5D. 3, 4, 66. 工作空间workspace中有变量a,b,c，按如下步骤进行操作，(1)将变量a存成数据文件adata.mat，(2)清除工作区间所有变量，(3)再将数据文件调入工作区间，下面( A )组命令是正确的：A save adata a; clear; load adata;B save adata a -ascii; clear; load adata;C save adata; clear all; load adata;D save adata; clear all; load –ascii adata;7. 在Command Window窗口中分别输入下列命令，对应输出结果正确的是( C )A. j=sqrt(-1) j=-1B. x=ones(1,2);x>0 ans=1C. y=diag(eye(2)).’ y=[1 1]D. rand(1,3) ans=[0.9501 1.2311 0.6068]8. 对应MATLAB提供的绘制平面曲线、离散序列、三维曲线、三维网格曲线、等高线图的函数为( B )A. plot, stem, mesh, plot3, surfB. plot, stem, plot3, mesh, contourC. plot, hist, mesh, plot3, contourD. plot, hist, plot3, meshc, contour9. 下列对条件语句、分支语句、固定循环次数语句、不定循环次数语句关键字排序正确的是( B )A. if else end; for end; switch case end; while end;B. if else end; switch case end; for end; while end;C. for end; if else end; while end; switch case end;D. while end; for end; switch case end; if else end;10. 在MATLAB的Editor中输入下列语句：x=[1 3;2 4];y=norm(x(:,2));z=det(x);u=diag(x);v=polyval(u,1);运行后，在Command Window 中对应y,z,u,v的输出结果为( A )A.5 -2 [1;4] 5B.5 2 [1;2] 3C.4 2 [1;2] 3D.4 -2 [1;4] 5三．简答题(每题10分，共60分)1．下面三种表示方法有什么不同的含义？（1）f=3*x^2+5*x+2（2）f='3*x^2+5*x+2'（3）x=sym('x')f=3*x^2+5*x+2答：（1）f=3*x^2+5*x+2表示在给定x时，将3*x^2+5*x+2的数值运算结果赋值给变量f，如果没有给定x则指示错误信息。

《MATLAB及应用》实验指导书班级：姓名：学号：总评成绩：汽车工程系电测与汽车数字应用中心目录实验04051001 MATLAB语言基础 (3)实验04051002 MATLAB科学计算及绘图 (12)实验04051001 MATLAB语言基础实验目的1)熟悉MATLAB的运行环境2)掌握MATLAB的矩阵和数组的运算3)掌握MATLAB符号表达式的创建4)熟悉符号方程的求解实验内容（任选6题）1.利用rand等函数产生下列矩阵：产生一个均匀分布在（-5，5）之间的随机阵（50×2），要求显示精度为精确到小数点后一位（精度控制指令为format）。

format banka=-5; b=5;r = a + (b-a).* rand(50,2)r =3.15 -2.244.06 1.80-3.73 1.554.13 -3.371.32 -3.81-4.02 -0.02-2.22 4.600.47 -1.604.58 0.854.65 -2.76-3.42 2.514.71 -2.454.57 0.06-0.15 1.993.00 3.91-3.58 4.59-0.78 0.474.16 -3.612.92 -3.514.59 -2.421.56 3.41-4.64 -2.463.49 3.144.34 -2.562.58 -1.502.43 -3.03-1.08 -2.491.55 1.16-3.29 -0.272.06 -1.48-4.68 3.31-2.23 0.85-4.54 0.50-4.03 4.173.23 -2.141.952.57-1.83 2.544.50 -1.20-4.66 0.68-0.61 -4.24-1.18 -4.462.66 0.312.95 2.79-3.13 4.34-0.10 -3.70-0.54 0.691.46 -0.312.09 -4.882.55 -1.632.在一个已知的测量矩阵T（100×100）中，删除整行数据全为0的行，删除整列数据全为0的列（判断某列元素是否为0方法：检查T(: , i) .* (T(: , i))是否为0）。

国家开放大学《Matlab语言及其应用》实验报告
——Matlab基本程序设计
姓名：学号：
1．M文件
M文件是一个文本文件，可以用任何文本编辑器来建立和编辑，通常使用 Matlab 自带的 M 文件编辑器。

例：编写一个脚本文件将华氏温度转化为摄氏温度。

c=5
9
(f−32)
●新建一个 M 文件 f2cs.m，内容如下：
clear; % 清除当前工作空间中的变量
f=input('Please input Fahrenheit temperature:');
c=5*(f-32)/9;
fprintf('The centigrade temperature is %g\n',c);
●在命令窗口中输入 f2cs，即可执行该 M 文件。

脚本文件运行后，文件中所产生的所有变量都驻留在 Matl ab 当前工作空间，可以用 whos 查看。

2.程序控制结构
例：计算函数y=cos(x+1)+√x2+1 x≥10d的值。

clear
x=input('请输入x的值:');
if x>=10
y=cos(x+1)+sqrt(x*x+1);
end
y
请输入x的值:12
y =12.9490
3.程序调试
例：拼写错误
clear
t=0:0.1:10;
y=son(t);
例：函数错误
调试器的使用：MatLab的调试器和M文件编辑器集成一体。

实验结果及分析实验1：程序如下x=1:10y=2*x;plot(x,y)仿真结果：实验结果分析：仿真结果是条很规则的直线，X轴和Y轴一一对应，清楚明了，而序又特别简单。

所以用Maltab 软件很方便地画出规则的直线，方便研究。

实验结果及分析1、A=2、A=1A=实验结果及分析实验三 Matlab在信号与系统中的应用实验名称实验1、掌握信号与系统课程中基本知识的Matlab编程、仿真方法目的实验原理实验1程序：b=[1];a=[1 1];p=;t=0:p:5;x=exp(-3*t);subplot(1,2,1);impulse(b,a,0:p:5);title('冲激响应');subplot(1,2,2);step(b,a,0:p:5);title('阶跃响应');实验内容<设计性实验>1、用MATLAB在时域中编程求解y′(t)+y(t)=f(t), f(t)= exp(-3t)ε(t)的冲激响应、阶跃响应。

在simulink仿真环境下，设计系统框图，分析系统的冲激响应、阶跃响应。

<设计性实验>（选做）2、用MATLAB在时域中编程求解y′(t)+y(t)=f(t), f(t)=(1+exp(-3t))ε(t)的冲激响应、阶跃响应，要求用conv编程实现系统响应。

在simulink仿真环境下，设计系统框图，分析系统的冲激响应、阶跃响应。

实验结果及分析实验1仿真结果：simulink仿真环境下冲激响应阶跃响应实验名称实验四 Matlab在数字信号处理中的应用实验结果及分析实验1仿真结果:6khz12kHZ。

《Matlab语言与应用》课程实验报告*名：**班级：电信114学号：************指导老师：***二〇一三年十一月二十一日Matlab实验报告实验一一、实验问题求[12 + sin(2)×( 22 −4)]÷3^2的算术运算结果。

二、问题的分析该题主要熟悉Matlab环境下的对数学运算的熟悉，如何输入数据、建立函数输出结果。

三、上机实验结果如图四、实验的总结与体会通过本次实验，我学会了用Matlab来计算数学运算中的复杂技术。

我们也可以自己编写一个可以调用的函数，首先我们要了解Matlab语言函数的基本结构，结构如下：Function [返回变量列表]=函数名（输入变量列表）注释说明语句段，由%引导输入、返回变量格式的检测函数体语句例如：输入变量为k，返回的变量为m和s,其中s为前m项的和Function [m,s]=findsum(k)s=0;m=0;while(s<=k),m=m+1;s=s+m;end编写了函数，就可以将其存为findsum.m文件，这样就可以在Matlab环境中对不同的k值调用该函数了。

这样就可以灵活的实现我们想要的数据。

实验二一、实验问题二、问题分析输入矩阵时，空格或逗号表示间隔，分号表示换行，比如上面的矩阵A应写为A=[1,2,3;4,5,6;7,8,9]三、上机实验结果如图四、实验的总结与体会通过对本次上机实验了解到在Matlab中对矩阵的代数运算矩阵转置、矩阵的加减法运算、矩阵乘法、矩阵的左、右除、矩阵翻转、矩阵乘方运算、点运算等。

实验三一、实验要求画图，理解plot函数用法二、代码如下clear; clf;t=0:pi/20:2*pi;R=5;x=R*sin(t); y=R*cos(t);plot(x,y,'b:'), gridhold onrrr=[x;y;x+y];plot(rrr(1,:),'.','MarkerSize',10,'Color','r')plot(rrr(2,:),rrr(3,:),'o','MarkerSize',15,'Color','b'); axis([-8,20,-8,8]), % axis squarehold off三、Matlab运行结果如图实验四一、实验要求二、代码如下t=0:pi/50:4*pi;y0=exp(-t/3);y=exp(-t/3).*sin(3*t);plot(t,y,'-r',t,y0,':b',t,-y0,':b') Grid三、Matlab运行结果如图实验五一、实验要求傅里叶频谱分析二、代码及分析(1)首先生成数据，包含50Hz和120Hz频率的正弦波x >>t = 0:.001:.25;>>x = sin(2*pi*50*t) + sin(2*pi*120*t);(2)再生成噪音信号yy = x + 2*randn(size(t));plot(y(1:50))title('Noisy time domain signal')(3)对y进行快速傅里叶变换Y = fft(y,256);(4)计算功率谱Pyy = Y.*conj(Y)/256;f = 1000/256*(0:127);plot(f,Pyy(1:128))title('Power spectral density')xlabel('Frequency (Hz)')(5)只查看200Hz以下频率段plot(f(1:50),Pyy(1:50))title('Power spectral density')xlabel('Frequency (Hz)')三、上机结果如下图实验六一、实验要求FIR数字滤波器设计代码如下clear;close allf=[0,0.6,0.6,1]; m=[0,0,1,1]; % 设定预期幅频响应b=fir2(30,f,m); n=0:30; % 设计FIR 数字滤波器系数subplot(3,2,1),stem(n,b,'.')xlabel('n'); ylabel('h(n)');axis([0,30,-0.4,0.5]),line([0,30],[0,0])[h,w]=freqz(b,1,256);subplot(3,2,2),plot(w/pi,20*log10(abs(h)));gridaxis([0,1,-80,0]),xlabel('w/pi'); ylabel('幅度(dB)');二、上级结果如图实验七二、实验要求用guide实验一个简单的加减乘除计算器二、实验步骤在命令行输入guide命令，进入guide界面新建一个空白guide文件在空白文件中设置好功能模块如图模块建立好后，就要把编写好的加减乘除代码加入到各自的回调函数中，见下图两个被加数代码如下图加模块代码如下图减模块代码如下图乘模块代码如下图除模块代码如下图各模块的回调函数加完后就可以运行了，运行结果如下图总结：Matlab一个高级的距阵/阵列语言，它包含控制语句、函数、数据结构、输入和输出和面向对象编程特点。

实验一：MATLAB 运算基础学 院班级 姓 名 学 号一、 实验目的1.熟悉启动和退出MATLAB 方法。

2.熟悉MATLAB 命令窗口的组成。

3.掌握建立矩阵的方法。

4.掌握MATLAB 各种表达式的书写规则以及常用函数的使用。

二、 实验内容1.先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量；（1）122sin 851z e︒=+ （2）(21ln 2z x =+，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦（3）()0.30.330.3sin 0.3ln , 3.0, 2.9, 2.8,,2.8,2.9,3.022a a e e az a a --+=++=--- 提示：利用冒号表达式生成a 向量，求各点的函数值时用点乘运算。

（4）2242,011,12,21,23t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩其中0:0.5:2.5t =2.已知：1234413134787,2033657327A B --⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦求下列表达式：（1）A +6*B 和A -B +I (其中I 为单位矩阵) （2）A *B 和A .*B （3）A ^3和A .^3 （4）A /B 和B \A（5）[A ,B ]和[A ([1,3],:);B ^2]完成时间：得 分：3.设有矩阵A 和B1234530166789101769,111213141502341617181910970212223242541311A B ⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦（1）求它们的乘积C 。

（2）将矩阵C 的右下角3x2子矩阵赋给D 。

（3）查看MATLAB 工作空间的使用情况。

4.完成下列操作：（1）求[100,999]之间能被21整除的数的个数。

提示：先利用冒号表到式，在利用find 和length 函数。

实验二：MATLAB 矩阵分析与处理学 院班级 姓 名 学 号一、 实验目的1.掌握生成特殊矩阵的方法。

2.掌握矩阵分析的方法。

3.用矩阵求逆法解线性方程组。

二、 实验内容习题三1.写出完成下列操作的命令。

（1）建立3阶单元矩阵A 。

（2）建立5X6随机矩阵A ,其元素为[100,200]范围内的随机整数。

（3）产生均值为1，方差为0.2的500个正态分布的随机数。

（4）产生和A 同样大小的幺矩阵。

（5）将矩阵A 对角线的元素加30.（6）从矩阵A 提取主对角线元素，并以这些元素构成对角阵B 。

2.使用函数，实现方阵左旋90。

或右旋90的功能。

例如，原矩阵为A ,A 左旋后得到B ，右旋后得到C 。

1011123211471078965425811,,45698736912123121110A B C ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥===⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦3.建立一个方阵A ,求A 的逆矩阵和A 的行列式的值，并验证A 和A -1是互逆。

4.求下列线性方程组的解。

1231231242232101238x x x x x x x x +-=⎧⎪-+=⎨⎪+=⎩5.求下列矩阵的主对角线元素、上三角阵、下三角阵、秩、范数、条件数和迹。

（1）11235142305211159A -⎡⎤⎢⎥-⎢⎥=⎢⎥⎢⎥⎣⎦（2）0.434328.9421B ⎡⎤=⎢⎥-⎣⎦6.求矩阵A 的特征值和相应的特征向量。

110.5110.250.50.252A ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦完成时间：得 分：实验二内容1.设有分块矩阵33322322ER A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦,其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和零矩阵，试通过数值计算验证22E R RS A O S +⎡⎤=⎢⎥⎣⎦。

2.产生5阶希尔伯特矩阵H 和5阶帕斯卡矩阵P ，且求其行列式的值Hh 和Hp以及它们的条件数Th 和Tp ，判断哪个矩阵性能更好。

试卷总分:100 得分:100
判断题（共17小题，1至15小题每小题6分，16至17小题每小题5分，共100分）1.matlab中有两种描述系统的方法：传递函数模型和零极点增益模型。

答案:正确
2.matlab中两种表示模型的方法之间可以相互转换。

答案:正确
3.matlab中支持的系统互联方式有三种：串联、并联和反馈。

答案:正确
4.方波信号的产生函数是square()。

答案:正确
5.冲激信号是特异信号。

答案:正确
6.系统反馈的函数为feedback。

答案:正确
7.step函数用于计算系统的冲激响应。

答案:错误
8.impulse函数用于计算系统的阶跃响应。

答案:错误
9.lism函数可用于求解系统对任意输入的响应。

答案:正确
10.一个线性系统对正弦信号的响应中，其频率、幅值和相角均有可能发生改变。

答案:错误
11.Matlab 中freqs函数用于绘制系统的幅频特性。

答案:正确
12.bode图可用于判断系统是否稳定。

答案:正确
13.M-函数是利用Matlab语言编写的特定功能函数。

答案:正确
14.simulink的S函数，可以是用C语言编写的。

答案:正确。

国家开放大学《Matlab语言及其应用》实验报告——熟悉Matlab 操作环境姓名：学号：1. Matlab桌面2.指令窗口历史指令窗记录着用户在命令指令窗中所输入过的所有指令行，且所有这些被记录的指令行都能被复制，并送到指令窗中再运行。

3.当前目录流量器（Current Directory Browser）当前目录窗口中显示了MATLAB当前工作目录下的所有文件夹与文件，以便用户对当前目录下的文件进行管理。

工作空间窗口以列表形式显示了MATLAB工作区中当前所有变量的名称及其属性，属性中包括了变量的类型、长度及其占用空间大小。

工作空间由运行MATLAB的程序或命令所生成的所有变量和MATLAB 本身提供的常量构成。

MATLAB 6.x后才提供。

早期版本是在命令窗口中键入who查看变量名，或者键入whos命令查看变量名、变量大小和数据类型等内容。

4.文件编辑/调试器（Editor/debugger）进入MATLAB文件编辑器的方法如下：①点击MATLAB桌面上的图标，打开空白的M文件编辑器。

②点击MATLAB桌面上的图标，填写所选文件名后，再点[打开]，即可展示相应 M文件编辑器。

③用鼠标左键双击当前目录窗中的所需M文件，可直接引出相应的M文件编辑器。

MATLAB通过自带的M文件编辑/调试器来创建和编辑M文件。

M文件(带.m扩展名的文件)类似于其它高级语言的源程序。

M文件编辑器可用来对M文件进行编辑和交互调试，也可阅读和编辑其它ASCⅡ码文件。

M文件编辑/调试器窗口由菜单栏、工具条和文本编辑区组成，是标准的windows 风格。

5.使用MATLAB帮助两个重要的帮助命令:1．help命令使用方法：help 帮助主题“帮助主题”可以是命令名（函数名）、目录名或者部分路径名。

回车后，命令窗口中就会出现与帮助主题的“help”有关的内容。

注意：命令窗口显示的MATLAB帮助信息中，是用大写字母来突出函数名的，但在使用这些函数时，应该用小写字母。

《MATLAB 语言及其用》实验指导书目录实验一Matlab 使用方法和程序设计........................实验二控制系统的模型及其转换.............................实验三控制系统的时域、频域和根轨迹分析...........实验四动态仿真集成环境-Simulink.........................实验一Matlab使用方法和程序设计一、实验目的1、掌握Matlab软件使用的基本方法；2、熟悉Matlab的数据表示、基本运算和程序控制语句3、熟悉Matlab绘图命令及基本绘图控制4、熟悉Matlab程序设计的基本方法二、实验内容：1、帮助命令使用help命令，查找 sqrt（开方）函数的使用方法；在 CommandWindowL里输入help，接在在search里输入sqr即可。

sqrtSquare rootSyntaxB = sqrt(X)DescriptionB = sqrt(X) returns the square root of each element of the array X. For the elements of X that are negative or complex, sqrt(X) produces complex results.TipsSee sqrtm for the matrix square root. Examplessqrt((-2:2)')ans =0 + 1.4142i0 + 1.0000i1.00001.4142See Alsonthroot | realsqrt | sqrtm2、矩阵运算（1）矩阵的乘法已知A=[1 2;3 4]; B=[5 5;7 8];求A^2*BA=[1 2;3 4];B=[5 5;7 8];C=A^2*B>> format compactC =105 115229 251（2）矩阵除法已知 A=[1 2 3;4 5 6;7 8 9];B=[1 0 0;0 2 0;0 0 3];A\B,A/BA=[1 2 3;4 5 6;7 8 9];B=[1 0 0;0 2 0;0 0 3];C=A\B,D=A/BC =1.0e+016 *0.3152 -1.2609 0.9457-0.6304 2.5218 -1.89130.3152 -1.2609 0.9457D =1.0000 1.0000 1.00004.0000 2.5000 2.00007.0000 4.0000 3.0000（3）矩阵的转置及共轭转置已知A=[5+i,2-i,1;6*i,4,9-i];求A.', A'A=[5+i,2-i,1;6*i,4,9-i];B=A.', C=A'B =5.0000 + 1.0000i 0 +6.0000i2.0000 - 1.0000i 4.00001.0000 9.0000 - 1.0000iC =5.0000 - 1.0000i 0 -6.0000i2.0000 + 1.0000i 4.00001.0000 9.0000 + 1.0000i（4）使用冒号表达式选出指定元素已知： A=[1 2 3;4 5 6;7 8 9];求A中第3列前2个元素；A中所有列第2，3行的元素；A=[1 2 3;4 5 6;7 8 9];B1=A([1,2],[3])B2=A([2,3],:)B1 =36B2 =4 5 67 8 9方括号[]用magic函数生成一个4阶魔术矩阵，删除该矩阵的第四列A=magic(4)B=A(:,[1,2,3])或A=magic(4)A(:,4)=[]A =16 2 3 135 11 10 89 7 6 124 14 15 1B =16 2 35 11 109 7 64 14 153、多项式（1）求多项式4xxp的根=x)2(3--Y=[1 0 -2 -4];S=roots(Y)S =2.0000-1.0000 + 1.0000i-1.0000 - 1.0000i（2）已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ，求矩阵A的特征多项式;把矩阵A作为未知数代入到多项式中；A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4]P=poly(A)polyval(P,A)A =1.2000 3.0000 5.0000 0.90005.0000 1.7000 5.00006.00003.0000 9.0000 0 1.00001.00002.00003.00004.0000P =1.0000 -6.9000 -77.2600 -86.1300 604.5500ans =1.0e+003 *0.3801 -0.4545 -1.9951 0.4601-1.9951 0.2093 -1.9951 -2.8880-0.4545 -4.8978 0.6046 0.43530.4353 0.0840 -0.4545 -1.16174、基本绘图命令（1）绘制余弦曲线 y=cos(t)，t∈[0，2π]（2）在同一坐标系中绘制余弦曲线y=cos(t-0.25)和正弦曲线y=sin(t-0.5)，t∈[0，2π]（1）t=[0:0.05:2*pi];y=cos(t);plot(t,y)-1-0.8-0.6-0.4-0.20.20.40.60.81（2）t=[0:0.05:2*pi];y1=cos(t-0.25);y2=sin(t-0.5);plot(t,y1)hold onplot(t,y2)5、基本绘图控制绘制[0，4π]区间上的x1=10sint曲线，并要求：（1）线形为点划线、颜色为红色、数据点标记为加号；（2）坐标轴控制：显示范围、刻度线、比例、网络线（3）标注控制：坐标轴名称、标题、相应文本；t=[0:0.1:4*pi];x1=10*sin(t);plot(t,x1,'r-.+'); %画图，显示红色、点划线、标记加号；axis([0,15,-10,10]); %定义显示范围，横轴为[0,15]，纵轴为[-10,10]；title('曲线x1=10sint'); %显示标题；xlabel('T轴');ylabel('X1轴'); %显示坐标轴名称；set(gca,'xminortick','on');set(gca,'yminortick','on'); %显示刻度线；grid on %显示网络线T 轴X 1轴6、基本程序设计（1）编写命令文件：计算1+2+…+n<2000 时的最大n 值； （2）编写函数文件：分别用for 和while 循环结构编写程序，求2的0到n 次幂的和。

《MATLAB 语言与应用》实验课程任务书一、 实验教学目标与基本要求上机实验是本课程重要的实践教学环节；实验的目的不仅仅是验证理论知识，更重要的是通过上机实验，加强学生的实验手段与实践技能，掌握应用MATLAB 语言求解问题的方法，培养学生分析问题、解决问题、应用知识的能力和创新精神，全面提高学生的综合素质。

上机实验共8学时。

主要实验内容是基于理论课所学知识对课后典型习题进行MATLAB 求解，基本掌握常见数学问题的求解方法与命令调用，更深入地认识和了解MATLAB 语言强大的计算功能。

上机实验最终以书面报告的形式提交，并作为期末成绩考核内容的一部分。

二、 实验内容（8学时）第一部分MATLAB 语言编程、科学绘图与基本数学问题求解（4学时）主要内容：掌握MATLAB 语言编程基础、科学绘图方法、微积分问题、线性代数问题等基本数学问题的求解与应用。

练习题：1、安装MATLAB 软件，应用demo 命令了解主要功能，熟悉基本功能，会用help 命令。

2、用MATLAB 语句输入矩阵A 和B⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=1423143212344321A ， ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡++++++++++++++++=4j 11j 43j 22j34j 11j 42j 33j 24j 13j 22j 31j 41j 42j 33j 24j 1B 前面给出的是44⨯矩阵，如果给出5)6,5(=A 命令将得出什么结果？ InputA=[1,2,3,4;4,3,2,1;2,3,4,1;3,2,4,1];B=[1+4j,2+3j,3+2j,4+1j;4+1j,3+2j,2+3j,1+4j;2+3j,3+2j,4+1j,1+4j;3+2j,2+3j,4+1j,1+4j]; A(5,6)=5 Answer= A =1 2 3 4 0 04 3 2 1 0 02 3 4 1 0 03 24 1 0 00 0 0 0 0 53、假设已知矩阵A，试给出相应的MATLAB命令，将其全部偶数行提取出来，赋给B矩阵，用magic(8)A=命令生成A矩阵，用上述命令检验一下结果是不是正确。

国开《Matlab语言及其应用》形考任务123答案题目如为随机抽题请用CTRL+F来搜索试题一、单项选择题（每小题3分，共36分）假设数组A=[123456789]，则mod(A，2)为(B)。

正确答案是: B．A=[101010101]M文件根据调用方式的不同可以分为(C)。

正确答案是: C．脚本文件和函数文件Matlab中的程序控制结构不包含(D)。

正确答案是: D．过程结构符号运算得到的是(B) 正确答案是: B．表达式Matlab中的变量名的第一个字母是(A)。

正确答案是: A．字母plot绘制图形时若想画出分格线，需进行的操作是(C)。

正确答案是: C．gridon图形控制中，legend函数的作用是(D)。

正确答案是: D．在图右上角建立图例饼图的绘制函数为(D) 正确答案是: D．piex=linspace(5，1，5)得到的数组是(C)。

正确答案是: C．54321X=5*ones(1，6)得到的数组是(B)。

正确答案是: B．555555一般情况下，绘图命令每执行一次就刷新当前图形窗口，图形窗口原有图形将不复存在。

若希望在已绘图形上叠加新的图形，可以使用(A)命令。

正确答案是: A．holdon下列哪条指令是求矩阵的行列式的值？(C) 正确答案是: C．det三阶幻方又称为九宫图，提取三阶幻方矩阵对角元并构造对角阵用(C) 正确答案是: (C)diag(diag(magic(3)))data=rand(1000,2);x=data(:,1);y=data(:,2);II=find(yx.^2);的功能是(B) 正确答案是: (B)统计1000个随机点落入特殊区域的点的索引值；MATLAB计算二项分布随机变量分布律的方法是(C) 正确答案是:(C)binopdf(x,n,p);。

MATLAB命令symse2;f=sqrt(1-e2*cos(t)^2);S=int(f,t,0,pi/2)功能是(D) 正确答案是: (D)计算f(t)定积分的符号结果。

《MATLAB及应用》实验报告3实验报告实验项目名称 MATLAB符号计算所属课程名称 MATLAB及应用实验类型上机实验实验日期12月日指导教师谢芹班级学号姓名成绩一、实验名称MATLAB符号计算及应用二、实验目的（1）掌握定义符号对象的方法（2）掌握符号表达式的运算法则以及符号矩阵运算（3）掌握求符号函数极限及导数的方法（4）掌握求符号函数定积分和不定积分的方法三、实验原理（略）四、实验内容1. 一无阻力抛射体的飞行，给定初速v，试计算物体在真空中飞行的时间和距离，并绘出其运行轨迹。

2y??x?115，在x?0到x?10之间所围面积，并讨论步长和积分方法对精度2. 用数值积分法求的影响。

3. 计算二重积分22(x?y)dxdy，积分区域D为由x?1，y?x及y?0所围成的闭合区域，并??D作出积分区域D的示意图。

4. 计算三重积分???xy2z2dxdydz，积分区域V为由x?1，y?x，z?xy及z?0所围成的闭合V区域，并作出积分区域V的示意图。

五、实验过程及结果（含源代码）1.h=1000*rand(); m=1000*rand(); g=9.8t=sqrt(2*h/g) v=100*rand(); x=v*tt=0:0.1:sqrt(2*h/g); x=v*t;h=h-(1/2).*g.*t.^2; comet(x,h) g =9.8000 t =12.8946 x =163.74412.syms x; y='-x.^2+115'; i=50; for n=1:ix0=linspace(0,10,n); y0=-(x0).^2+115; s1(n)=trapz(x0,y0);s2(n)=quad(y,0,10,1/n); end plot(s1)hold on plot(s2) hold off3.syms x y;int(int(x^2+y^2,x,0,1),y,0,1) x=0:0.001:1; y=x;ezplot('0',[0 2]); hold on;ezplot('x',[0 2]); plot([1,1],[0,1],'r');area(x,y) %蓝色区域为填充区域 ans = 2/34.syms x y zs='x*y^2*z^2';int(int(int(s,z,0,x*y),y,0,x),x,0,1)ans = 1/198z=-1:0.01:1; [Y,Z]=meshgrid(z); z=ones(201); X=ones(201); mesh(X,Y,Z) z=-1:0.01:1;[X1,Z1]=meshgrid(z); Y1=X1; hold on mesh(X1,Y1,Z1) z=-1:0.01:1;[X2,Y2]=meshgrid(z); Z2=X2.*Y2; hold on mesh(X2,Y2,Z2) z=-1:0.01:1;[X3,Y3]=meshgrid(z); Z3=zeros(201); hold onmesh(X3,Y3,Z3)六、实验总结通过本次matlab的上机实验课程，使我可以运用matlab求解函数的定积分与不定积分，并画出相应的图形，根据题目所给的条件列出函数并求解函数，对符号对象的定义使用的更加熟练，积分与不懂积分的方法掌握的更加熟练，较为轻松的解决该类问题。

《Matlab语言》实验报告《Matlab语言》实验（报告）题库1、TDOA信号分析类（1）已给出一段事先采集的信号，该信号为进行TDOA定位使用的基本信号，其格式为GPS+IQ + GPS+IQ … ，即每包数据由GPS头文件和IQ信号构成，GPS头文件共58B，其数据格式为$HT,20130114,084556,N3606.82273,E10343.59311,M1538.7, 11,0*，每包IQ数据共8192B，其数据格式为I0,Q0,I1,Q1,I2,Q2…,I2047,Q2047，即I数据2048点、Q数据2048点交叉出现。

换言之，每包数据实际内容为：$HT,20130114,084556,N3606.82273,E10343.59311,M1538.7,11,0 * I0 Q0 I1 Q1 I2 Q2 … I2047 Q2047，程序前期已实现读取IQ数据文件并进行关键信息读取分解，请根据程序提醒，完成相关功能（数据及程序见“1-实际IQ信号实验”文件夹）。

2、TDOA时差估计仿真类（2）在TDOA定位技术中，时差估计是一个非常重要的环节。

自行仿真2个具有一定时差的信号，用广义互相关法（GCC）计算该2个信号的时差，并与设定时差进行对比分析（需给出详细过程及适当的仿真图）。

（3）在TDOA定位技术中，时差估计是一个非常重要的环节。

自行仿真2个具有一定时差的信号，用互模糊函数法计算该2个信号的时差，并与设定时差进行对比分析（需给出详细过程及适当的仿真图）。

（4）在TDOA定位技术中，时差估计是一个非常重要的环节。

自行仿真2个具有一定时差的信号，用广义互相关（GCC）结合多项式拟合方法计算该2个信号的时差，并比较广义互相关法估计时差和广义互相关结合多项式拟合方法估计时差的结果，进行分析（需给出详细过程及适当的仿真图）。

（5）在TDOA定位技术中，时差估计是一个非常重要的环节。

自行仿真2个具有一定时差的信号，自选方法计算该2个信号的时差，并与设定时差进行对比分析（需给出详细过程及适当的仿真图）。