2014《MATLAB及应用》实验报告4 -

第四章N=21时：M4-3程序如下：N=input('N=');W=blackman(N);wc=0.4*pi;alpha=(N-1)/2;k=0:N-1;hd=(wc/pi)*sinc((wc/pi)*(k-alpha));h=hd.*W';omega=linspace(0,pi,512);mag=freqz(h,[1],omega);magdb=20*log10(abs(mag));plot(omega/pi,magdb);grid; N=51时：axis([0 1 -100 5]);set(gca,'ytick',[-80 -75 -60 -40 -20 0]);结论：滤波器的长度增加，即N增加时，通阻带的震荡波纹数量随之增加，波纹宽度变窄，过渡带变窄。

但最小阻带衰减不变，即最小阻带衰减与N无关。

m4-6（b）程序如下：wp=input('wp=');ws=input('ws=');N=61;alpha=(N-1)/2;dw=ws-wp;wc=(ws+wp)/2;h1=sinc((dw/pi).*((0:29)-alpha)/2)*(wc/pi).*sinc((wc/pi).*((0:29)-alpha));h2=wc/pi;h3=sinc((dw/pi).*((31:60)-alpha)/2)*(wc/pi).*sinc((wc/pi).*((31:60)-alpha));h=[h1 h2 h3];omega=linspace(0,pi,512); （1）ΔΩ=0.1*pi时，As= -36db: mag=freqz(h,[1],omega);magdb=20*log10(abs(mag));plot(omega/pi,magdb);grid;axis([0 1 -100 5]);（2）ΔΩ=0.3*pi时，As= -48db: （3）ΔΩ=0.5*pi时，As= -52db:可见：随着过渡带的增宽，阻带衰减值增加。

信号与系统MATLAB第一次实验报告一、实验目的1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。

2.学会运用MATLAB表示常用连续时间信号的方法3.观察并熟悉一些信号的波形和特性。

4.学会运用MATLAB进行连续信号时移、反折和尺度变换。

5.学会运用MATLAB进行连续时间微分、积分运算。

6.学会运用MATLAB进行连续信号相加、相乘运算。

7.学会运用MATLAB进行连续信号的奇偶分解。

二、实验任务将实验书中的例题和解析看懂，并在MATLAB软件中练习例题，最终将作业完成。

三、实验内容1.MATLAB软件基本运算入门。

1). MATLAB软件的数值计算：算数运算向量运算：1.向量元素要用”[ ]”括起来，元素之间可用空格、逗号分隔生成行向量，用分号分隔生成列向量。

2.x=x0:step:xn.其中x0位初始值，step表示步长或者增量，xn为结束值。

矩阵运算：1.矩阵”[ ]”括起来；矩阵每一行的各个元素必须用”,”或者空格分开；矩阵的不同行之间必须用分号”;”或者ENTER分开。

2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。

3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。

举例：计算一个函数并绘制出在对应区间上对应的值。

2).MATLAB软件的符号运算：定义符号变量的语句格式为”syms 变量名”2.MATLAB软件简单二维图形绘制1).函数y=f(x)关于变量x的曲线绘制用语：>>plot(x,y)2).输出多个图像表顺序：例如m和n表示在一个窗口中显示m行n列个图像，p表示第p个区域，表达为subplot(mnp)或者subplot(m,n,p)3).表示输出表格横轴纵轴表达范围：axis([xmax,xmin,ymax,ymin])4).标上横轴纵轴的字母：xlabel(‘x’),ylabel(‘y’)5).命名图像就在subplot写在同一行或者在下一个subplot前：title(‘……’)6).输出：grid on举例1：举例2：3.matlab程序流程控制1).for循环：for循环变量=初值:增量:终值循环体End2).while循环结构：while 逻辑表达式循环体End3).If分支：(单分支表达式)if 逻辑表达式程序模块End(多分支结构的语法格式)if 逻辑表达式1程序模块1Else if 逻辑表达式2程序模块2…else 程序模块nEnd4).switch分支结构Switch 表达式Case 常量1程序模块1Case 常量2程序模块2……Otherwise 程序模块nEnd4.典型信号的MATLAB表示1).实指数信号：y=k*exp(a*t)举例：2).正弦信号：y=k*sin(w*t+phi)3).复指数信号：举例：4).抽样信号5).矩形脉冲信号:y=square(t,DUTY) (width默认为1)6).三角波脉冲信号：y=tripuls(t,width,skew)(skew的取值在-1~+1之间，若skew取值为0则对称)周期三角波信号或锯齿波:Y=sawtooth(t,width)5.单位阶跃信号的MATLAB表示6.信号的时移、反折和尺度变换:Xl=fliplr(x)实现信号的反折7.连续时间信号的微分和积分运算1).连续时间信号的微分运算：语句格式：d iff(function,’variable’,n)Function:需要进行求导运算的函数，variable:求导运算的独立变量，n：求导阶数2).连续时间信号的积分运算：语句格式：int(function,’variable’,a,b)Function:被积函数variable:积分变量a:积分下限b:积分上限(a&b默认是不定积分)8.信号的相加与相乘运算9.信号的奇偶分解四、小结这一次实验让我能够教熟悉的使用这个软件，并且能够输入简单的语句并输出相应的结果和波形图，也在一定程度上巩固了c语言的一些语法。

4)MATLAB实验报告MATLAB实验报告一、引言MATLAB是一种强大的数学软件，广泛应用于科学研究、工程设计和数据分析等领域。

本实验报告旨在介绍使用MATLAB进行数据处理和可视化的基本方法，并通过实例展示其应用。

二、数据处理1. 数据导入与读取在MATLAB中，可以使用load函数或importdata函数导入外部数据文件。

load函数适用于导入MATLAB格式的数据文件，而importdata函数可以导入多种格式的数据文件，如文本文件、Excel文件等。

2. 数据清洗与转换数据清洗是指对原始数据进行去除异常值、缺失值填充等处理，以保证数据的准确性和完整性。

MATLAB提供了丰富的函数和工具箱，如统计工具箱、优化工具箱等，可以方便地进行数据清洗和转换操作。

3. 数据分析与建模MATLAB具有强大的数学和统计分析功能，可以进行数据的描述性统计、回归分析、时间序列分析等。

通过使用相关函数和工具箱，可以对数据进行深入分析，并建立相应的数学模型。

三、数据可视化1. 统计图表MATLAB提供了丰富的绘图函数和工具箱，可以绘制各种统计图表，如直方图、散点图、箱线图等。

这些图表可以直观地展示数据的分布、关系和趋势，有助于更好地理解和解释数据。

2. 三维可视化除了二维图表外，MATLAB还支持三维数据的可视化。

通过使用plot3函数、mesh函数等，可以将三维数据以曲线、曲面等形式展示出来，进一步提供对数据的深入理解。

3. 动态可视化MATLAB还支持动态可视化，即通过动画或交互式图形来展示数据的变化过程。

通过使用animate函数、interactiveplot函数等，可以将数据的变化以动态的方式展示出来，增加数据分析和呈现的趣味性。

四、实例应用以某电商平台销售数据为例，展示如何使用MATLAB进行数据处理和可视化。

首先，导入销售数据文件，清洗数据，去除异常值和缺失值。

然后，通过统计分析，计算销售额、销量、平均价格等指标，并绘制相应的统计图表。

MATLAB及应用实验指导书实验一MATLAB使用基础一、实验目的1、熟悉MATLAB基本的操作界面。

2、掌握MATLAB中变量、数组、向量等对象的生成方法。

3、掌握符号矩阵的生成和基本运算方法。

4、掌握MATLAB中平面绘图方法。

5、掌握MATLAB中的常用帮助命令使用方法。

二、实验设备计算机，MATLAB语言环境三、实验基础理论1、常见数学函数2、系统的在线帮助（1）help 命令：1）.当不知系统有何帮助内容时，可直接输入help以寻求帮助:>> help（回车）2）.当想了解某一主题的内容时，如输入：>> help syntax （了解Matlab 的语法规定）3）.当想了解某一具体的函数或命令的帮助信息时，如输入：>> help sqrt （了解函数sqrt 的相关信息）（2）lookfor 命令现需要完成某一具体操作，不知有何命令或函数可以完成，如输入：>> lookfor line （查找与直线、线性问题有关的函数）3、常量与变量系统的变量命名规则：变量名区分字母大小写；变量名必须以字母打头，其后可以是任意字母，数字，或下划线的组合。

此外，系统内部预先定义了几个有特殊意义和用途的变量，见下表：（1）数值型向量（矩阵）的输入1）．任何矩阵（向量），可以直接按行方式输入每个元素：同一行中的元素用逗号（，）或者用空格符来分隔；行与行之间用分号（；）分隔。

所有元素处于一方括号（[ ]内；例1：>> Time = [11 12 1 2 3 4 5 6 7 8 9 10]>> X_Data = [2.32 3.43；4.37 5.98]2）．系统中提供了多个命令用于输入特殊的矩阵：4、数组（矩阵）的点运算运算符：+（加）、-（减）、./（右除）、.\（左除）、.^（乘方），例2：>> g = [1 2 3 4]；h = [4 3 2 1]；>> s1 = g + h, s2 = g.*h, s3 = g.^h, s4 = g.^2, s5 = 2.^h5、矩阵的运算运算符：+（加）、-（减）、*（乘）、/（右除）、\（左除）、^（乘方）、’（转置）等；常用函数：det（行列式）、inv（逆矩阵）、rank（秩）、eig（特征值、特征向量）、rref （化矩阵为行最简形）例3：>> A=[2 0 –1;1 3 2]; B=[1 7 –1;4 2 3;2 0 1];>> M = A*B % 矩阵A与B 按矩阵运算相乘>> det_B = det(B) % 矩阵A的行列式>> rank_A = rank（A）% 矩阵A的秩>> X = A/B % 矩阵右除，A/B是B左乘A的逆>> Y = B\A % 矩阵左除，B\A其计算结果是矩阵B的逆和矩阵A相乘6、平面绘图MA TLAB 能绘制表现非常丰富的平面图形，以下是一些常用的MATLAB 绘图函数，以绘制不同平面图元和实现不同的平面图形绘制功能：可以使用grid on 绘图网格线和hold on 同时绘制两个图形，还可以给图形加上各种注解和处理：例4：>>x=-2*pi:pi/10:2*pi;>>y=sin(x);>>plot(x,y,’r*:’)>>grid on % 绘制网格线>>hold on>>y2=2*cos(x);>>plot(x,y2,’b+-‘)>>xlabel('x 轴'); %x 轴注解>>ylabel('y 轴'); %y 轴注解>>title('正弦和余弦函数');%图形标题>>legend({'y = sin(x)','y2 = cos(x)'});%图形注解（2）．MA TLAB 还可以在同一图形窗口完成多个子图的功能，如：>>x=-2*pi:pi/10:2*pi;>>y=sin(x);>>y2=cos(x);>>y3=2*cos(x);>>y4=sin(x)+cos(x);>>subplot(2,2,1);>>plot(x,y,‘r*:’)>>subplot(2,2,2);>>plot(x,y2,‘b+-’)>>subplot(2,2,3);>>plot(x,y3,‘gp:’)>>subplot(2,2,4)>>plot((x,y4,‘kx’)四、实验内容与步骤：1．练习数据和符号的输入方式，将前面的命令在命令窗口中执行通过；2．输入A=[7 1 5；2 5 6；3 1 5]，B=[1 1 1; 2 2 2; 3 3 3]，在命令窗口中执行下列表达式，掌握其含义：A(2, 3) A*B A.*B A^2 A.^2 B/A B./A 3．输入C=1:2:20，则C（i）表示什么？其中i=1,2,3, (10)4. 建立全1矩阵（ones)；建立对角线为1的矩阵（eye)zeros(3) ones(2,3)5．. 练习把y = sin( x ). / x ，z = cos(x ) ，u = 2* sin( x ) ，v = sin( x ) ./ cos( x ) 在[0 , 2 π] 区间内的四个子图分别用不同的颜色、点型和线型绘制在同一个窗口中，并给其中加上纵坐标,标题、图例和网格线。

电气工程学院实验报告实验项目名称MATLAB数值计算所属课程名称MATLAB及应用实验类型上机实验实验日期2014-03- 5指导教师盛义发班级XXXXXX学号2012XXXXXXXXX 姓名XXXXXX成绩一、实验名称MATLAB数值计算二、实验目的（1）掌握MATLAB变量的使用（2）掌握MATLAB数组的创建（3）掌握MATLAB数组和矩阵的运算（4）熟悉MATLAB多项式的运用三、实验原理1. 矩阵分析矩阵转置：单引号（’）矩阵的旋转：rot90(A,k)，功能是将矩阵A旋转90度的k倍，缺省值是1 矩阵的左右翻转：fliplr(A)矩阵的上下翻转：flipud(A)矩阵的逆：inv(A)，与A^(-1)等价矩阵的行列式：det(A)矩阵的秩：rank(A)矩阵的迹：trace(A)将矩阵化为最简式：rref(A)矩阵的特征值与特征向量：(1) E = eig(A)；矩阵A的所有特征值构成向量E；(2) [V,D]=eig(A)；A的所有特征值构成对角阵D，A的特征向量构成V的列向量；2. 多项式多项式的建立：若多的项的全部根构成的向量为X，则以X为根的多项式为poly(X)多项式的根：roots(p)计算以向量p为系数的多项式的根，包括重根，复根多项式求值：polyval(p,x)，p是多项式的系数，x可以是一个数也可以是一个矩阵多项式求拟合次数：polyfit(x,y,n)，x可以是一个数也可以是一个矩阵，y 是x对应的数或矩阵多项式的四则运算：(1)P1+P2；(2)P1-P2；(3)conv(P1,P2), (4)deconv(P1,P2)四、实验内容1. 已知矩阵11 12 13 1421 22 23 24A= 31 32 33 344142 43 44(1) A(:,1) (2) A(2,:) (3) A(:,2:3)(4) A(2:3,2:3) (5) A(:,1:2:3) (6) A(2:3)(7) A(:) (8) A(:,:) (9) ones(2,2)(10) eye(2) (11) [A,[ones(2,2);eye(2)]] (12) diag(A)(13) diag(A,1) (14) diag(A,-1) (15) diag(A,2)2. 已知2354127354217623A-⎡⎤⎢⎥--⎢⎥=⎢⎥-⎢⎥--⎣⎦.求（1）A的逆；（2）A的行列式；（3）A的迹；（4）A的所有特征向量和特征值。

信号与系统MATLAB第一次实验报告一、实验目的1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。

2.学会运用MATLAB表示常用连续时间信号的方法3.观察并熟悉一些信号的波形和特性。

4.学会运用MATLAB进行连续信号时移、反折和尺度变换。

5.学会运用MATLAB进行连续时间微分、积分运算。

6.学会运用MATLAB进行连续信号相加、相乘运算。

7.学会运用MATLAB进行连续信号的奇偶分解。

二、实验任务将实验书中的例题和解析看懂，并在MATLAB软件中练习例题，最终将作业完成。

三、实验内容1.MATLAB软件基本运算入门。

1). MATLAB软件的数值计算：算数运算向量运算：1.向量元素要用”[ ]”括起来，元素之间可用空格、逗号分隔生成行向量，用分号分隔生成列向量。

2.x=x0:step:xn.其中x0位初始值，step表示步长或者增量，xn为结束值。

矩阵运算：1.矩阵”[ ]”括起来；矩阵每一行的各个元素必须用”,”或者空格分开；矩阵的不同行之间必须用分号”;”或者ENTER分开。

2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。

3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。

举例：计算一个函数并绘制出在对应区间上对应的值。

2).MATLAB软件的符号运算：定义符号变量的语句格式为”syms 变量名”2.MATLAB软件简单二维图形绘制1).函数y=f(x)关于变量x的曲线绘制用语：>>plot(x,y)2).输出多个图像表顺序：例如m和n表示在一个窗口中显示m行n列个图像，p表示第p个区域，表达为subplot(mnp)或者subplot(m,n,p)3).表示输出表格横轴纵轴表达范围：axis([xmax,xmin,ymax,ymin])4).标上横轴纵轴的字母：xlabel(‘x’),ylabel(‘y’)5).命名图像就在subplot写在同一行或者在下一个subplot前：title(‘……’)6).输出：grid on举例1：举例2：3.matlab程序流程控制1).for循环：for循环变量=初值:增量:终值循环体End2).while循环结构：while 逻辑表达式循环体End3).If分支：(单分支表达式)if 逻辑表达式程序模块End(多分支结构的语法格式)if 逻辑表达式1程序模块1Else if 逻辑表达式2程序模块2…else 程序模块nEnd4).switch分支结构Switch 表达式Case 常量1程序模块1Case 常量2程序模块2……Otherwise 程序模块nEnd4.典型信号的MATLAB表示1).实指数信号：y=k*exp(a*t)举例：2).正弦信号：y=k*sin(w*t+phi)3).复指数信号：举例：4).抽样信号5).矩形脉冲信号:y=square(t,DUTY) (width默认为1)6).三角波脉冲信号：y=tripuls(t,width,skew)(skew的取值在-1~+1之间，若skew取值为0则对称)周期三角波信号或锯齿波:Y=sawtooth(t,width)5.单位阶跃信号的MATLAB表示6.信号的时移、反折和尺度变换:Xl=fliplr(x)实现信号的反折7.连续时间信号的微分和积分运算1).连续时间信号的微分运算：语句格式：d iff(function,’variable’,n)Function:需要进行求导运算的函数，variable:求导运算的独立变量，n：求导阶数2).连续时间信号的积分运算：语句格式：int(function,’variable’,a,b)Function:被积函数variable:积分变量a:积分下限b:积分上限(a&b默认是不定积分)8.信号的相加与相乘运算9.信号的奇偶分解四、小结这一次实验让我能够教熟悉的使用这个软件，并且能够输入简单的语句并输出相应的结果和波形图，也在一定程度上巩固了c语言的一些语法。

《MATLAB程序设计与应用》实验报告专业：机械设计制造及其自动化班级：机械三班学号：17428320姓名：刘羽新2019年上期实验四：MATLAB 绘图基础一、实验目的通过本实验掌握MATLAB 二维、三维作图的基本操作，了解绘制图形的辅助操作。

二、实验要求1. 根据实验内容，编写相应的MATLAB 程序，并将程序及结果放置于相应位置。

2. 实验内容采用小四号宋体、1.5倍行距进行排版。

3. 实验报告以word 文档的.doc 格式保存；文件文件名为：学号_姓名，如“15312345_张三”。

三、实验内容1. 绘制下列曲线 (1) 21001y x =+, 55x -≤≤ （用绿色的虚线绘制，添加图形说明“二维曲线绘图”） x=-5:0.001:5;y=100./(1+x.*x);plot(x,y,'b:');title('二维曲线绘图');(2)235x ty t⎧=⎪⎨=⎪⎩, 1010x-≤≤（用红色的双划线绘制，添加图例“参数曲线绘图”）t=-10:0.001:10;x=t.*t;y=5*t.*t.*t;plot(x,y,'r--');legend('参数曲线绘图');(3)cossinx ty tz t=⎧⎪=⎨⎪=⎩, 02xπ≤≤（用黑色实线绘制，添加坐标轴名称）t = 0:0.001:2*pi; x=cos(t);y=sin(t);plot3(x,y,z,'k-'); xlabel('X');ylabel('Y');zlabel('Z');2. 绘制下列分段函数的曲线，并在各段曲线附近分别添加文本说明“21y x x =++”与“21sin y x π=+”：221,013sin ,12x x x y x x π⎧++≤<⎪=⎨+≤≤⎪⎩ t_1=0:0.001:1;t_2=1:0.001:2;y_1=x_1.*x_1+x_1+1;y_2=3+(sin(pi*x_2)).^2;plot(x_1,y_1,x_2,y_2);gtext('y=x^2+x+1');gtext('y=1+(sin(pi*x))^2');3. 已知241(11)1xy xx+=-≤≤+，sin3cos(0)sin3sinx t tty t tπ=⎧≤≤⎨=⎩，以子图形式绘制两条曲线。

国家开放大学《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文件编辑器集成一体。

M A T L A B应用实验指导书结果公司内部档案编码：[OPPTR-OPPT28-OPPTL98-OPPNN08]???????????????????????? MATLAB语言实验指导书????????????????????????中国矿业大学信息与电气工程学院2014年3月实验一 MATLAB 工作环境熟悉及基本运算一、实验目的：熟悉MATLAB 的工作环境，学会使用MATLAB 进行一些简单的运算。

掌握基本的矩阵运算及常用的函数。

二、实验内容：MATLAB 的启动和退出，熟悉MATLAB 的桌面（Desktop ），包括菜单（Menu ）、工具条 （Toolbar ）、命令窗口(Command Window)、历史命令窗口、工作空间(Workspace)等；完成一些基本的矩阵操作；学习使用在线帮助系统。

三、实验步骤：1、启动MATLAB ，熟悉MATLAB 的桌面。

2、在命令窗口执行命令完成以下运算，观察workspace 的变化，记录运算结果。

（1）（365-52?2-70）?3 =（2）area=pi*^2 =（3）已知x=3，y=4，在MATLAB 中求z ：()232y x y x z -== 576 （4）将下面的矩阵赋值给变量m1,在workspace 中察看m1在内存中占用的字节数。

m1=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡11514412679810115133216 执行以下命令>>m1( 2 , 3 )=10>>m1( 11 )=6>>m1( : , 3 )= 3 10 6 15>>m1( 2 : 3 , 1 : 3 )=[ 5 11 10;9 7 6]>>m1( 1 ,4 ) + m1( 2 ,3 ) + m1( 3 ,2 ) + m1( 4 ,1)=34（5）执行命令>>help abs查看函数abs 的用法及用途，计算abs( 3 + 4i )=5（6）执行命令>>x=0::6*pi;>>y=5*sin(x);>>plot(x,y)（7）运行MATLAB 的演示程序，>>demo ，以便对MATLAB 有一个总体了解。

实验四文件操作一、实验目的掌握MATLAB文件操作的各个函数。

二、实验内容1. 验证本章所举的例子例 5.1>> Fid=fopen('std.dat','w');>> Fid=fopen('std.dat','r')Fid =3>> A=fread(Fid);>> Sat=fclose(Fid)Sat =例 5.2>> magic=[1 2 3 4;5 6 7 8]magic =123 45678>> Fid=fopen('magic5.bin','w') Fid =3>> fwrite(Fid,magic,'int32') ans =8>> Sat=fclose(Fid)Sat =>> A=[-0.6515 -0.2727 -0.4354 -0.3190 -0.9047;-0.7534 -0.4567 -0.3212 -0.4132 -0.3583;-0.9264 -0.8173 -0.7823 -0.3265 -0.0631;-0.1735 -0.7373 -0.0972 -0.3267 -0.6298;-0.4768 -0.6773 -0.6574 -0.1923 -0.4389]A =-0.6515-0.2727-0.4354-0.3190-0.9047-0.7534-0.4567-0.3212-0.4132-0.3583-0.9264-0.8173-0.7823-0.3265-0.0631-0.1735-0.7373-0.0972-0.3267-0.6298-0.4768-0.6773-0.6574-0.1923-0.4389>> Fid=fopen('test.dat','w')Fid =3>> cnt=fwrite(Fid,A,'float')25>> fclose(Fid)ans =程序段将矩阵A的数据以二进制浮点数格式写入文件test.dat 中。

《MATLAB 语言及其应用》实验报告实验指导老师：学院：电气与信息工程学院班级：自动化姓名：学号：学号：《MATLAB 语言及其应用》实验指导目录实验一Matlab 使用方法和程序设计........................实验二控制系统的模型及其转换.............................实验三控制系统的时域、频域和根轨迹分析...........实验四动态仿真集成环境-Simulink.........................实验五直流电机自动调速系统控制器设计实验一Matlab使用方法和程序设计一、实验目的1、掌握Matlab软件使用的基本方法；2、熟悉Matlab的数据表示、基本运算和程序控制语句3、熟悉Matlab绘图命令及基本绘图控制4、熟悉Matlab程序设计的基本方法二、实验内容：1、帮助命令使用help命令，查找 sqrt（开方）函数的使用方法；2、矩阵运算（1）矩阵的乘法已知A=[1 2;3 4]; B=[5 5;7 8];求A^2*B>> A=[1 2;3 4];>> B=[5 5;7 8];>> C=A^2*BC =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/B>> A=[1 2 3;4 5 6;7 8 0];>> B=[1 0 0;0 2 0;0 0 3];>> A\B,A/Bans =-1.7778 1.7778 -0.33331.5556 -1.5556 0.6667-0.1111 0.4444 -0.3333ans =1.0000 1.0000 1.00004.0000 2.5000 2.00007.0000 4.0000 0（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];>> A.',A'ans =5.0000 + 1.0000i 0 +6.0000i2.0000 - 1.0000i 4.00001.0000 9.0000 - 1.0000ians =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行的元素；方括号[]用magic函数生成一个4阶魔术矩阵，删除该矩阵的第四列>> A=[1 2 3;4 5 6;7 8 9];>> B=A(1:2,3)B =6 >> B1=A(2:3,:) B1 =4 5 6 7 8 9 >> B1=A(2:3,:) B1 =4 5 6 7 8 9B(:,4)=[ ] B =16 2 3 5 11 10 9 7 6 4 14 153、多项式（1）求多项式 42)(3--=x x x p 的根 >> p=[1 0 -2 -4]; >> r=roots(p) r =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) p =1.0000 -6.9000 -77.2600 -86.1300 604.5500 >> y=sym(A)[ 6/5, 3, 5, 9/10][ 5, 17/10, 5, 6][ 3, 9, 0, 1][ 1, 2, 3, 4]4、基本绘图命令（1）绘制余弦曲线 y=cos(t)，t∈[0，2π]（2）在同一坐标系中绘制余弦曲线y=cos(t-0.25)和正弦曲线y=sin(t-0.5)，t∈[0，2π]>> t=0:0.01:2*pi;>> y=cos(t);>> plot(t,y);>> title('余弦函数y=cos(t)');>> x=0:0.01:2*pi;y=cos(x-0.25);plot(x,y)hold ony=sin(x-0.5);plot(x,y)>> title('y=cos(t-0.25),y=sin(t-0.5)')5、基本绘图控制绘制[0，4π]区间上的x1=10sint曲线，并要求：（1）线形为点划线、颜色为红色、数据点标记为加号；（2）坐标轴控制：显示范围、刻度线、比例、网络线（3）标注控制：坐标轴名称、标题、相应文本；>> t=0:0.01:4*pi;>> x1=10*sin(t);>> plot(t,x1,'-.r+');>> xlabel('X');>> ylabel('Y');>> grid on>> title('x1=10sin(t)')>> legend('y',4)6、基本程序设计（1）编写命令文件：计算1+2+…+n<2000 时的最大n值；>> s=0;n=0;>> while s<2000;s=s+n;n=n+1;end>> mn=n-1mn =63（2）编写函数文件：分别用for和while循环结构编写程序，求2的0到n次幂的和。

MATLAB应用实验指导书1234-结果************************ MATLAB语言实验指导书************************中国矿业大学信息与电气工程学院2014年3月实验一 MATLAB 工作环境熟悉及基本运算一、实验目的：熟悉MATLAB 的工作环境，学会使用MATLAB 进行一些简单的运算。

掌握基本的矩阵运算及常用的函数。

二、实验内容：MATLAB 的启动和退出，熟悉MATLAB 的桌面（Desktop ），包括菜单（Menu ）、工具条（Toolbar ）、命令窗口(Command Window)、历史命令窗口、工作空间(Workspace)等；完成一些基本的矩阵操作；学习使用在线帮助系统。

三、实验步骤：1、启动MATLAB ，熟悉MATLAB 的桌面。

2、在命令窗口执行命令完成以下运算，观察workspace 的变化，记录运算结果。

（1）（365-522-70） 3 = （2）area=pi*^2 = （3）已知x=3，y=4，在MATLAB 中求z ：()232y x y x z -== 576 （4）将下面的矩阵赋值给变量m1,在workspace 中察看m1在内存中占用的字节数。

m1=11514412679810115133216 执行以下命令>>m1( 2 , 3 )=10 >>m1( 11 )=6>>m1( : , 3 )= 3 10 6 15>>m1( 2 : 3 , 1 : 3 )=[ 5 11 10;9 7 6]>>m1( 1 ,4 ) + m1( 2 ,3 ) + m1( 3 ,2 ) + m1( 4 ,1)=34（5）执行命令>>help abs查看函数abs 的用法及用途，计算abs( 3 + 4i )=5 （6）执行命令>>x=0::6*pi; >>y=5*sin(x); >>plot(x,y)（7）运行MATLAB 的演示程序，>>demo ，以便对MATLAB 有一个总体了解。

matlab实验报告4实验四M⽂件的编写⼀、实验环境计算机MA TLAB软件⼆、实验⽬的1.学习MA TLAB中的关系运算和逻辑运算，掌握它们的表达形式和⽤法。

2.掌握MA TLAB中的选择结构和循环结构。

3.学会⽤MATLAB进⾏M⽂件的编写和调⽤。

三、预备知识1.关系和逻辑运算关系运算符⽤来完成关系运算，在控制程序流程⽅⾯有着极为重要的作⽤。

MATLAB 常⽤的关系符有：<、>、<=⼩于或等于、>=⼤于或等于、==等于、~=不等于。

关系运算符可以⽤来⽐较两个数值，若所描述的关系成⽴，则结果为1，表⽰逻辑真，反之，若所描述的关系不成⽴，结果为0，表⽰逻辑假。

MATLAB中的逻辑运算符有&与、|或、~⾮、逻辑运算法则A B A&B A|B xor(A,B) ~A0 0 0 0 0 10 1 0 1 1 11 0 0 1 1 01 1 1 1 0 02.选择结构if语句和switch语句if语句的⼀般形式如下：if A1 %表达式1B1 %命令2else if A2 %表达式2B2 %命令2else B3 %命令3endswitch语句的⼀般结构如下：switch a %读⼊⼀个语句case A1 %情况1B1 %命令1case A2case …………other case %其余情况Bn %最后⼀个命令3.循环结构for语句⼀半⽤于循环次数已知的情况，⽽while语句⼀般⽤于循环次数未知的情况。

for语句的格式为：for 变量=表达式命令1命令2…… endwhile 语句的格式为：while 表达式命令 end 四、实验内容和步骤1.创建⼀个矩阵，⽤函数all 和any 作⽤于该矩阵，⽐较结果。

创建⼀个矩阵a=[1,2,0,3；2,0,1,4];b=all(a),c=any(a),d=all(b),e=any(b) ⽐较结果b=1 0 0 1 c=1 1 1 1 d=0 e=1结果分情况讨论：当函数为all 时，在所选的数组中，只要有数等于0，那么，返回值为0；当不为0时，返回值为1。

《MATLAB 原理及应用》实验报告三.课后练习题答案1.为 ⎪⎩⎪⎨⎧<-=>+=01001x x x xx x y 编写赋值程序。

程序如下：①建立如下的M 文件：x=input('x=');%让用户通过键盘输入数值、字符串或表达式if x>0y=x+1;elseif x==0y=x;else x<0y=x-1;e nd程序执行结果如下>> kh1 %在当前工作目录下，文件名为“kh.1.m ” x=1>> yy =22.使用for ... end循环的array向量编程求出1+3+5...+100 的值程序如下：sum=0;>> for k=1:2:100sum=sum+k;end>> sumsum =25003.计算1+3+5...+100 的值，当和大于1000时终止计算。

程序如下：sum=0;for m=1:2:100; %建立1 3 5….100的向量if sum<=1000 %如果sum小于1000则可以继续加sum=sum+m; %累加elsebreak; %若sum的结果不符合条件就跳出整个循环endend结果为：sum =1024k =653.1计算从1开始多少个自然数之和超过100。

程序如下：>> sum=0;n=0;>> while sum<=100n=n+1;sum=sum+n;end结果为：n =14sum =1054.求1！+2！+3！+……+8！的值程序如下：n=1;sum=1;for m=2:8; %循环7次使得得到各次阶乘n=n*m;sum=sum+n; %累加end结果为:sum =462335.写程序，判断一年是否为闰年，符合下面两条件之一：（1990~2014）A、能被4整除，不能被100整除B、能被400整除程序如下count=0;for y=1990:2014;if((rem(y,4)==0&rem(y,100)~=0)|(rem(y,4)==0&rem(y,400)~=0));count=count+1;endend结果为:count =5。

学生实验报告开课学院及实验室: 机电学院2012年12月21日学院机电学院年级、专业、班姓名学号实验课程名称MATLAB程序设计成绩实验项目名称实验4: 数据和函数的可视化指导老师一、实验目的1、掌握MATLAB绘图的基本步骤和相关指令调用的先后顺序。

2、掌握MATLAB绘图指令的调用方法。

二、实验内容数学函数从形式上可以分为离散函数和连续函数。

MATLAB对这两种函数数据的可视化都提供了相应的指令。

仔细阅读教材【例5.1-1】的实现代码, 运行并保存结果；并改用stem函数, 画出【例5.1-1】的序列图。

仔细阅读教材【例5.1-2】的实现代码, 运行并保存结果；并分别使用描点和连折线方式, 画出连续函数y=xcosx的近似图形（采样点数自定, 要求画出的图尽量接近原连续函数的图）。

仔细阅读【例5.2-2】的实现代码, 理解plot指令画多条曲线的运用方法, 运行并保存结果；并使用plot函数和legend函数, 在同一个图形窗口上画出y=sint和y=sin(2t)在[0,2pi]区间上的图形, 并标出图例。

仔细阅读【例5.2-4】的实现代码, 理解图形标识选项的运用方法, 运行并保存结果；并修改代码, 把“sin(t)”字体改为正体, 大小改为20, “极大值”改为宋体。

阅读【例5.2-6】, 理解使用hold on指令画多幅图的方法, 运行并保存结果。

阅读【例5.2-8】, 理解使用subplot函数画多个子图的方法, 运行并保存结果。

（1）综合实验: 阅读以下关于通过绘制二阶系统阶跃响应综合演示图形标识的示例, 理解示例中所有图形标识指令的作用, 掌握各个图形标识指令的运用方法, 并在原指令上改动以实现以下功能:（2）把横坐标范围改为0至5pi, 纵坐标范围改为0至2；（3）把图中的横轴的刻度改为从0开始到4pi, 中间各点间隔为pi/2；纵轴刻度改为从0开始到1.5, 中间各点间隔为0.3；（4）把图中的α改为σ。

国家开放大学《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 图形基础班级：11电信一班姓名：何得中学号：实验目的：1.掌握MATLAB图形绘制功能；2.学会使用MATLAB提供的绘图函数；3.了解MATLAB对图形线性、色彩、光线、视角等的制定和处理。

实验仪器：MA TLAB软件，电脑实验数据：绘制fracreturnl变量同一图轴下绘制两条曲线增加子图分割子图删除图轴后仍保留两个子图效果添加图轴信息设定曲线样式在交互式绘图工具中增加图例>> rand('state',27)>> startprice=50;>>fracreturns1=0.0015*randn(200,1)+0.0003; >> x1=[startprice;1+fracreturns1];>> prices1=cumprod(x1);>> t=(1:length(prices1))';>> randn('state',7)>>fracreturns2=0.0015*randn(200,1)+0.0003; >> x2=[startprice;1+fracreturns2];>> prices2=cumprod(x2);>> whosName Size Bytes ClassL 51x51 20808 double arrayfracreturns1 200x1 1600 double arrayfracreturns2 200x1 1600 double arrayl1 1x1 8 double arrayl2 1x1 8 double arrayprices1 201x1 1608 double arrayprices2 201x1 1608 double arrays 1x1 8 double arraystartprice 1x1 8 double arrayt 201x1 1608 double arrayx1 201x1 1608 double arrayx2 201x1 1608 double arrayGrand total is 4010 elements using 32080 bytes>> plot (prices2, 'DisplayName', 'prices2', 'YDataSource', 'prices2'); figure(gcf) >> x=0:pi/1000:2*pi;>> y=exp(-0.5*x).*sin(2*pi*x+pi/4);>> plot(x,y)plot(x,y,x,y+1,x,y-1)>> B=pascal(6)B = 1 1 1 1 1 11 2 3 4 5 61 3 6 10 15 211 4 10 20 35 561 5 15 35 70 1261 6 21 56 126 252 >> plot(B)>> x=(0:pi/100:2*pi)';>> y1=2*exp(-0.5*x)*[1,-1];>> y2=2*exp(-0.5*x).*sin(2*pi*x);>> x1=(0:12)/2;>> y3=2*exp(-0.5*x1).*sin(2*pi*x1);>> plot(x,y1,'g:',x,y2,'b--',x1,y3,'rp');>> x1=0:pi/100:2*pi;>> x2=0:pi/100:3*pi;>> y1=exp(-0.5*x1).*sin(2*pi*x1);>> y2=1.5*exp(-0.1*x2).*sin(x2);>> plotyy(x1,y1,x2,y2);>> x=(0:pi/100:2*pi)';y1=2*exp(-0.5*x)*[1,-1];y2=2*exp(-0.5*x).*sin(2*pi*x);x1=(0:12)/2;y3=2*exp(-0.5*x1).*sin(2*pi*x1);plot(x,y1,'g:',x,y2,'b--',x1,y3,'rp');>> title('曲线及其包络线'); %加图形标题>> xlabel('independent variable X'); %加X轴说明>> ylabel('independent variable Y'); %加y轴说明>> text(2.8,0.5,'包络线'); %在指定位置添加图形说明>> text(0.5,0.5,'曲线y');>> text(1.4,0.1,'离散数据点');>> legend('包络线','包络线','曲线y','离散数据点') %加图例>> x=(0:pi/100:2*pi)';>> y1=2*exp(-0.5*x)*[1,-1];>>y2=2*exp(-0.5*x).*sin(2*pi*x);>> plot(x,y1,'b:');>> axis([0,2*pi,-2,2]);>> hold on;>> plot(x,y2,'k');>> grid on;>> box off;>> hold off; >> x=(0:pi/100:2*pi)';>> y1=2*exp(-0.5*x)*[1,-1];>> y2=2*exp(-0.5*x).*sin(2*pi*x);>> plot(x,y1,'b:');>> hold on;>> plot(x,y2,'k');>> grid on;>> box off;>> hold off>>alpha=-0.5;>>beta=3;>>A=50;>>t=0:0.01:10;>>y=A*exp(alpha*t).*sin(beta*t);>>plot(t,y);>>title('\fontname{隶书}\fontsize{16}{隶书}\fontname{Impact}{Impact}')>>xlabel('^{上标}and_{下标}')>>ylabel('Some\bf粗体\rm and some\it{斜体}')>>txt={'y={\itAe}^{\alphax}sin(\beta\itt)',... >> ['\itA\rm','='.num2str(A)],...>> ['\alpha='num2str(alpha)],...>> ['\beta=',num2str(beta)];>> text(2,22,txt)x=linspace(0,2*pi,60);y=sin(x);z=cos(x);t=sin(x)./(cos(x)+eps);ct=cos(x)./(sin(x)+eps); subplot(2,2,1);plot(x,y);title('sin(x)');axis([0,2*pi,-1,1]); subplot(2,2,2);plot(x,z);title('cos(x)');axis([0,2*pi,-1,1]);subplot(2,2,3);plot(x,t);title('tangent(x)');axis([0,2*pi,-40,40] );subplot(2,2,4);plot(x,ct);title('cotangent(x)');axis([0,2*pi,-40,40]);x=linspace(0,2*pi,60);y=sin(x);z=cos(x);t=sin(x)./(cos(x)+eps);ct=cos(x)./(sin(x)+eps); subplot(2,2,1);stairs(x,y);title('sin(x)-1');axis([0,2*pi,-1,1]); subplot(2,1,2);stem(x,y);title('sin(x)-2');axis([0,2*pi,-1,1]); subplot(4,4,3);plot(x,y);title('sin(x)');axis([0,2*pi,-1,1]); subplot(4,4,4);plot(x,z);title('cos(x)');axis([0,2*pi,-1,1]); subplot(4,4,7);plot(x,t);title('tangent(x)');axis([0,2*pi,-40,40] );subplot(4,4,8);plot(x,ct);title('cotangent(x)');axis([0,2*pi,-40, 40]);x=0:0.35:7;y=2*exp(-0.4*x);subplot(2,2,1);bar(x,y,'g');title('bar(x,y,''g'')');axis([0,7,0,2]);subplot(2,2,2);fill(x,y,'r');title('fill(x,y,''r'')');axis([0,7,0,2]);subplot(2,3,4);stairs(x,y,'b');title('stairs(x,y,''b'')');axis([0,7,0,2]); subplot(2,3,5);stem(x,y,'k');title('stem(x,y,k'')');axis([0,7,0,2]);subplot(2,3,6);area(x,y);title('area(x,y)');axis([0,7,0,2]);theta=0:0.01:2*pi;rho=sin(3*theta).*cos(3*theta);polar(theta,rho,'r');x=0:0.1:10;y=12*x.*x;subplot(2,2,1);plot(x,y);title('plot(x,y)');grid on;subplot(2,2,2);semilogx(x,y);title('semilogx(x, y)');grid on;subplot(2,2,3);semilogx(x,y);title('semipogx(x ,y)');grid on;subplot(2,2,4);loglog(x,y);title('loglog(x,y)');grid on;function y=myf(x);y=cos(tan(pi*x));fplot('myf',[-0.2,1.2],1e-4);subplot(1,2,1);pie([7,17,23,19,5]);title('饼图');legend('优秀','良好','中等','及格','不及格');subplot(1,2,2);compass([3+2i,4.5-i,-1.5+5i]);title('相量图');t=0:pi/50:2*pi;x=8*cos(t);y=4*sqrt(2)*sin(t);z=-4*sqrt(2)*si n(t);plot3(x,y,z,'p');title('Line in 3-D space');text(0,0,0,'origin'); xlabel('X'),ylabel('Y'),zlabel('Z');grid;x=5:29;y=14:35;[x,y]=meshgrid(x,y);z=2*x+5*y;k=find(z==126);x(k)',y(k)'>> bans =8 13 18 23 28ans =22 20 18 16 14x=0:0.1:2*pi;[x,y]=meshgrid(x);z=sin(y).*cos (x);mesh(x,y,z);xlabel('x-axis'),ylabel('y-axis'),zla bel('z-axis');title('mesh');x=0:0.1:2*pi;[x,y]=meshgrid(x);z=sin(y).*cos (x);surf(x,y,z);xlabel('x-asis'),ylabel('y-axis'),zlabe l('z-axis');title('surf');x=0:0.1:2*pi;[x,y]=meshgrid(x);z=sin(y).*cos (x);plot3(x,y,z);xlabel('x-axis'),ylabel('y-axis'),zla bel('z-axis');title('plot3');grid;[x,y]=meshgrid(-10:0.5:10);z=sin(sqrt(x.^3+y.^2))./sqrt(x.^2+y.^2+eps); subplot(2,2,1);meshc(x,y,z);title('meshc(x,y,x)');subplot(2,2,2);meshz(x,y,z);title('meshz(x,y,x)');subplot(2,2,3);surfc(x,y,z);title('surfc(x,y,z)');subplot(2,2,4);surfl(x,y,z);title('surfl(x,y,z)');t=0:pi/20:2*pi;[x,y,z]=cylinder(2+sin(t),30);subplot(1,3,1);surf(x,y,z);subplot(1,3,2);[x,y,z]=sphere;surf(x,y,z);subplot(1,3,3);[x,y,z]=peaks(300);meshz(x,y,z);subplot(2,2,1);bar3(magic(5));subplot(2,2,2);y=3*sin(0:pi/10:2*pi);stem3(y);subplot(2,2,3);pie3([2347,1827,2043,3025]);subplot(2,2,4);fill3(rand(3,6),rand(3,6),rand(3,6),'r');subplot(1,2,1);[x,y,z]=peaks;waterfall(x,y,z);xlabel('x-axis'),ylabel('y-axis'),zlabel('z-axis'); title('waterfall of peaks');subplot(1,2,2);contour3(x,y,z,12,'k');xlabel('x-axis'),ylabel('y-axis'),zlabel('z-axis');title('contour3 of peaks');subplot(2,2,1);mesh(peaks);view(-37.5,30);title('azimuth=-37.5,elevation=30')subplot(2,2,2);mesh(peaks);view(0,90);title('azimuth=0,elevation=90')subplot(2,2,3);mesh(peaks);view(90,0);title('azimuth=90,elevation=0')subplot(2,2,4);mesh(peaks);view(-7,-10);title('azimuth=-7,elevation=-10')z=peaks(20);colormap(copper);subplot(1,3,1);surf(z);subplot(1,3,2);surf(z);shading flat;subplot(1,3,3);surf(z);shading interp;x=0:0.1:2*pi;[x,y]=meshgrid(x);z=sin(y).*cos(x);[I,J]=find(z>0.2);for ii=1:length(I)z(I(ii),J(ii))=NaN;endsurf(x,y,z);。