MAAB实验报告Word文档下载推荐.docx

1、减.数组左除*矩阵乘/矩阵右除.*数组乘./数组右除矩阵乘方矩阵转置.数组乘方.数组转置表3MATLAB关系运算符=等于=不等于大于小于=大于等于小于等于表4MATLAB逻辑运算符逻辑运算符逻辑运算说明&And逻辑与|Or逻辑或Not逻辑非Xor逻辑异或表5MATLAB特殊运算符号功能说明示例：1:4;2:11.；分隔行.，分隔列（）%注释构成向量、矩阵！调用操作系统命令构成单元数组用于赋值4.MATLAB的一维、二维数组的寻访表6子数组访问与赋值常用的相关指令格式指令格式指令功能A（r,c）数组A中r指定行、c指定列之元素组成的子数组A（r,:）数组A中r指定行对应的所有列之元素组成的子数组

2、A（:,c）数组A中c指定列对应的所有行之元素组成的子数组数组A中各列元素首尾相连组成的“一维长列”子数组A（i）一维长列子数组中的第i个元素A（r,c）=Sa数组A中r指定行、c指定列之元素组成的子数组的赋值）=D（:数组全元素赋值，保持A的行宽、列长不变，A、D两组元素总和应相同5.MATLAB的基本运算表7两种运算指令形式和实质内涵的异同表数组运算矩阵运算指令含义A非共轭转置A共轭转置A=s把标量s赋给A的每个元素s+B标量s分别与B元素之和s-B,B-s标量s分别与B元素之差s.*A标量s分别与A元素之积s*A标量s分别与A每个元素之积s./B,B.s标量s分别被B的元素除S*inv（

3、B）B阵的逆乘sA.nA的每个元素自乘n次AnA阵为方阵，自乘n次A.p对A各元素分别求非整数幂Ap方阵A的非整数乘方A+B对应元素相加矩阵相加A-B对应元素相减矩阵相减A.*B对应元素相乘A*B内维相同矩阵相乘A./BA的元素别B的对应元素除A/BA右除BB.A与上相同BAA左除Bexp（A）以自然数e为底，分别以A的元素为指数，求幂expm（A）A的矩阵指数函数log（A）对A的各元素求对数logm（A）A的矩阵对数函数sqrt（A）对A的各元素求平方根sqrtm（A）A的矩阵平方根函数6.MATLAB的常用函数表8标准数组生成函数diag产生对角形数组（对高维不适用）rand产生均匀分布

4、随机数组eye产生单位数组（对高维不适用）randn产生正态分布随机数组magic产生魔方数组（对高维不适用）zeros产生全0数组ones产生全1数组size（）返回指定矩阵的行数和列数表9数组操作函数提取对角线元素，或生成对角阵flipud以数组“水平中线”为对称轴，交换上下对称位置上的数组元素fliplr以数组“垂直中线”为对称轴，交换左右对称位置上的数组元素reshape在总元素数不变的前提下，改变数组的“行数、列数”rot90矩阵逆时针旋转90度det方阵的行列式值rank矩阵的秩三、实验内容1、学习安装MATLAB软件。2、学习使用help命令，例如在命令窗口输入helpeye，然

5、后根据帮助说明，学习使用指令eye（其它不会用的指令，依照此方法类推）3、学习使用clc、clear，观察commandwindow、commandhistory和workspace等窗口的变化结果。4、初步程序的编写练习，新建M-file，保存（自己设定文件名，例如exerc1、exerc2、exerc3），学习使用MATLAB的基本运算符、数组寻访指令、标准数组生成函数和数组操作函数。注意：每一次M-file的修改后，都要存盘。四、实验结果练习A：（1） helprand，然后随机生成一个26的数组，观察commandwindow、commandhistory和workspace等窗口的变

6、化结果。（2）学习使用clc、clear，了解其功能和作用。（3）输入C=1:20，则C（i）表示什么？其中i=1,2,3,10。（4）输入A=715；256；315，B=111;222;333，在命令窗口中执行下列表达式，掌握其含义：A（2,3）A（:,2）A（3,:）A（:,1:3）A（:,3）.*B（:,2）A（:,3）*B（2,:）A*BA.*BA2A.2B/AB./A（5）二维数组的创建和寻访，创建一个二维数组（48）A，查询数组A第2行、第3列的元素，查询数组A第2行的所有元素，查询数组A第6列的所有元素。（6）两种运算指令形式和实质内涵的比较。设有3个二维数组A24，B24，C2

7、2，写出所有由2个数组参与的合法的数组运算和矩阵指令。（7）学习使用表8列的常用函数（通过help方法）（8）学习使用表9数组操作函数。练习B（9）用reshape指令生成下列矩阵（不超过3次操作），并取出方框内的数组元素。五、心得体会这次实验是一次上机进行，所以只有一步一步按照实验的步骤做，通过这次实验更加明白了一些知识和一些用法，光从书本上还是无法真的知道各种用法的区别的。实验二MATLAB数值及符号运算一、 实验目的：1、掌握矩阵的基本运算2、掌握矩阵的数组运算3、掌握多项式的基本运算4、会求解代数方程5、掌握创建符号表达式和矩阵的方法6、掌握符号表达式的微分和积分运算二、 实验基本知识

8、：1、多项式运算poly产生特征多项式系数向量roots求多项式的根p=poly2str（c,x）（以习惯方式显示多项式）conv,convs多项式乘运算deconv多项式除运算多项式微分polyder（p）:求p的微分polyder（a,b）:求多项式a,b乘积的微分2、代数方程组求解1）.恰定方程组的解方程ax=b（a为非奇异）两种求解方法：x=inv（a）?b采用求逆运算解方程x=ab采用左除运算解方程2）超定方程组的解方程ax=b,m0.5）-1产生，该输入信号的起始作用时刻为0。试用直杆图（提示：用stem指令）画出分别显示该系统输入、输出信号的两张子图。四、心得体会这次实验进一步提

9、升自己的能力，虽然刚开始实验时没有把所有的实验做完，回来后又一个人把没做完的补完，学习的过程是很美好的。希望一次比一次有进步。 实验三MATLAB语言的程序设计1、熟悉MATLAB程序编辑与设计环境2、掌握各种编程语句语法规则及程序设计方法3、函数文件的编写和设计4、了解和熟悉跨空间变量传递和赋值1、程序流程控制语句for循环结构语法：fori初值：增量：终值语句1语句nend说明：1i初值：终值，则增量为1。2初值、增量、终值可正可负，可以是整数，也可以是小数，只须符合数学逻辑。while循环结构while逻辑表达式循环体语句1、whiIe结构依据逻辑表达式的值判断是否执行循环体语勾。若表达

10、式的值为真，执行循环体语句一次、在反复执行时，每次都要进行判断。若表达式的值为假，则程序执行end之后的语句。2、为了避免因逻辑上的失误，而陷入死循环，建议在循环体语句的适当位置加break语句、以便程序能正常执行。（执行循环体的次数不确定；每一次执行循环体后，一定会改变while后面所跟关系式的值。3、while循环也可以嵌套、其结构如下：while逻辑表达式1循环体语句1while逻辑表达式2循环体语句2循环体语句3if-else-end分支结构if表达式1elseif表达式2（可选）语句2else（可选）语句31.if结构是一个条件分支语句，若满足表达式的条件，则往下执行；若不满足，则跳

11、出if结构。2elseif表达式2与else为可选项，这两条语句可依据具体情况取舍。3注意：每一个if都对应一个end，即有几个if，记就应有几个end。switch-case结构switch表达式case常量表达式1语句组1case常量表达式2语句组2otherwise语句组n1switch后面的表达式可以是任何类型，如数字、字符串等。2当表达式的值与case后面常量表达式的值相等时，就执行这个case后面的语句组如果所有的常量表达式的值都与这个表达式的值不相等时，则执行otherwise后的执行语句。程序流程控制指令break、return、pause2、函数文件的结构和编写方法函数只能在

12、函数体内对变量进行操作，也就是只能访问函数本身工作空间中的变量。M函数文件的结构（1）函数定义行（function）（2）H1行（函数帮助文本的第一行）（3）函数帮助文本（4）函数体（5）注释函数文件编写后，保存时，其文件名必须与函数名相同。函数名不要与MATLAB自身的函数命令相同。练习A2、用for循环语句实现求1100的和。?3、用for循环语句实现编写一个求n阶乘的函数文件。4、一行100列的Fibonacci数组a，元素满足Fibonacci规则：a（1）=a（2）=1，a（i）=a（i-1）+a（i-2）（i=3,4,），求该数组中第一个大于10000的元素。练习B思考题4、已知一

13、维数组A=2,4,5,8,10,B=4,6,9,3,4，用for循环语句实现求和函数可用sum（）。5、利用for循环求1！+2！+3！+?+5！的值。6、编写M函数表示曲线y2=e-t/3sin（3t）以及它的包络线y1=e-t/3，并从命令窗口输入命令语句绘制曲线，t的取值范围是0,4。7、设，编写一个M函数文件，使得调用f（x）时，x可用矩阵代入，得出的f（x）为同阶矩阵。这次实验，是我第一次练习用m文件写程序，然后再运行。用m文件要比直接在窗口上写要好一些，因为m文件是可以保存下来，这样自己当时没有截取结果，以后也是可以在自己的电脑上打开。实验四采用SIMULINK的系统仿真1、熟悉S

14、IMULINK工作环境及特点2、掌握线性系统仿真常用基本模块的用法3、掌握SIMULINK的建模与仿真方法1了解SIMULINK模块库中各子模块基本功能2.SIMULINK的建模与仿真方法（1）打开模块库，找出相应的模块。鼠标左键点击相应模块，拖拽到模型窗口中即可。（2）创建子系统：当模型大而复杂时，可创建子系统。（3）设置仿真控制参数。三、实验内容及结果练习SIMULINK仿真实际应用.1. 利用SIMULINK求解微分方程，方程的初始条件为x（0）=1，在增益模块Gain取值分别为2和100的情况下（即运行，给出运行结果。2. 已经系统框图，求传递函数，并求阶跃响应。此次实验，我们是需要用

15、SIMULINK来完成。刚开始用SIMULINK时，有些慌慌张张，因为不清楚一些模块的位置。后来通过老师的提点知道位置并成功把图画出来了，同时自己也掌握SIMULINK的建模与仿真方法。实验五采用SIMULINK的系统仿真（二）4、子系统的创建和封装设计三、实验内容及实验结果1、用两种滤波器设计方法实现下列问题：构建一个低通滤波系统的SIMULINK模型。输入信号是一个受正态噪声干扰的采样信号在此，而（秒）采用10阶Butterworth低通滤波器，以便从输入信号中过滤获得10Hz的输出信号。2、利用Simulink求解在区间0,1的积分，并求出积分值I（1）.3、时变系统的仿真：对时变受控对

16、象模型考虑一个控制系统模型如图所示，其中均为常数，试分析闭环系统的阶跃响应曲线。这次的实验相比于上一次的实验难度要大一些，因为用到的模块多而且比较复杂也不好寻找。但是却让我们学习了很多，解决一个问题的方法有很多，题目的解答也不是只有一种。通过多种方法的解答，更能学习到很多。实验六MATLAB的图形绘制1、学习MATLAB图形绘制的基本方法；2、熟悉和了解MATLAB图形绘制程序编辑的基本指令；3、熟悉掌握利用MATLAB图形编辑窗口编辑和修改图形界面，并添加图形的各种标注；4、掌握plot、subplot的指令格式和语法。 1、plot（x,y）-绘制由x,y所确定的曲线；2、多组变量绘图：p

17、lot（x1,y1,选项1,x2,y2,选项2,）；3、图形窗口的分割；4、图形编辑窗口的使用。二维曲线绘制的基本指令plot 基本调用格式plot（x,y,s特别注意：s是离散点形、连续线型和点线色彩三种设置数据构成的字符串，三种设置数据可以组合使用,顺序任意。 plot的衍生调用格式：1. 单色或多色绘制多条曲线：plot（X,Y,s）：用s指定的点形线型色彩绘制多条曲线。plot（X,Y）：采用默认的色彩次序用细实线绘制多条曲线。2. 多三元组绘制多条曲线plot（X1,Y1,s1,X2,Y2,s2,.,Xn,Yn,sn3. 单输入量绘线plot（Y） 当X、Y均为（mn）数组时，将绘制

18、出n条曲线。每条曲线的几何位置由X、Y对应的列确定。 当X、Y两个输入量中有一个是一维数组，且该数组的长度与另一个输入量的“行数”（或“列数”）相等时，将绘制出“列数”（或“行数”）条曲线。 plot（X,Y,s）：只能用s指定的同一色彩绘制多条曲线。【1】二维曲线绘图基本指令演示。本例运作后，再试验plot（t）,plot（Y）,plot（Y,t），以观察产生图形的不同。t=（0:pi/50:2*pi）;k=0.4:0.1:1;Y=cos（t）*k;plot（t,Y）【2】用图形表示连续调制波形Y=sin（t）sin（9t）及其包络线。t=0:pi/1000:pi;y1=（t）sin（t）.

19、*sin（9*t）;y2=sin（t）;plot（t,y1（t）,t,y2*1,-1,r-holdont0=linspace（0,pi,10）;fori=1:length（t0）t00=fzero（y1,t0（i）;plot（t00,0,oplot（pi,0,【3】在一个图形窗口绘制正弦和余弦曲线，要求给图形加标题“正弦和余弦曲线”，X轴Y轴分别标注为“时间t”和“正弦、余弦”，在图形的某个位置标注“sin（t）”“cos（t）”，并加图例，显示网格，坐标为正方形坐标系。pi/100:2*pi;y1=sin（t）;y2=cos（t）;plot（t,y1,r.,t,y2,g.）;legend（s

20、in（t）cos（t）axissquare;title（sintcosttext（1,-0.4,cos（t）rightarrowtext（2.8,0.4,leftarrowsin（t）xlabel（时间ylabel（正弦、余弦这次是通过写程序得到绘制的图形，突然觉得MATLAB这个软件非常的神奇，先自己写，再运行，然后再对比结果进行修改。这个过程看到了波形格式各样的变换很精彩。实验七MATLAB在信号与系统中的应用1掌握利用MATLAB求解LTI系统的方法2掌握用MATLAB中滤波分析的函数命令。3.熟悉MATLAB信号处理工具箱。1、绘制三维网格曲面图指令mesh（）。2.对于线性定常离散系

21、统的响应函数为filter（）。函数调用格式：filter（B,A,x）:对输入为向量x进行滤波处理。1、已知二阶系统传递函数为，当时，试计算阻尼比从0.11时的二阶系统的阶跃响应，并绘制阶跃响应三维网格曲面图。2、某LTI系统的单位冲激响应，输入，初始条件为0，求系统响应y（t）.3、已知滤波器fc=30,fs=1000,为10阶的butterworth滤波器，对信号输入x（n）=0.8n（n=0:31）的响应y（n）.画出输入x（n）、输出y（n）波形。此次实验学习了很多，掌握利用MATLAB求解LTI系统的方法，同时也掌握用MATLAB中滤波分析的函数命令。实验八：频率采样法设计FIR数

22、字滤波器数字滤波器从滤波器功能上可以分为低通、高通、带通、带阻以及全通滤波器；根据系统的单位冲击响应的特性，又可以分为有限长单位冲激响应滤波器（FIR）和无限长单位冲激响应滤波器（IIR）。 与IIR滤波器相比，FIR滤波器在保证幅度特性满足技术要求的同时，很容易做到有严格的线性相位特性。设FIR滤波器单位脉冲响应h（n）的长度为N,其系统函数为：当滤波器的系数N满足一定的对称条件时，就可以获得线性相位。线性相位FIR滤波器共分为四种类型，分别为：1） 类型，系数对称，即,N为奇数。2） 类型，N为偶数。3） 类型，系数反对称，即，N为奇数。4） 类型1）频率采样法设计滤波器时，通常给出的是幅频特性的技术指标要求，可直接在频域进行处理，按照理想的频率特性，在到2之间等间隔采样N点，得到：然后用H（k）的傅里叶逆变换作为滤波器的系数：构成一个系统传递函数为的实际的FI