matlab实验报告.doc

通信系统仿真实验报告姓名： 杨仕浩 专业： 通信工程 学号： 2014111347 实验一 MATLAB环境的熟悉与基本运算 一、实验目的： 1熟悉MATLAB开发环境2掌握矩阵、变量、表达式的各种基本运算二、实验基本知识： 1.熟悉MATLAB环境: MATLAB桌面和命令窗口、命令历史窗口、帮助信息浏览器、工作空间浏览器文件和搜索路径浏览器。2.掌握MATLAB常用命令clc 清除命令窗口中内容clear 清除工作空间中变量help 对所选函数的功能、调用格式及相关函数给出说明lookfor 查找具有某种功能的函数但却不知道该函数的准确名称who查询工作空间中的变量信息3.MATLAB变量与运算符变量命名规则如下： （1） 变量名可以由英语字母、数字和下划线组成（2） 变量名应以英文字母开头（3） 长度不大于31个（4） 区分大小写MATLAB中设置了一些特殊的变量与常量，列于下表。表1 MATLAB的特殊变量与常量变量名 功能说明 变量名 功能说明 ANS 默认变量名，以应答最近一次操作运算结果realmin 最小的正实数i或j 虚数单位INF(inf) 无穷大pi 圆周率NAN(nan) 不定值（0/0） eps 浮点数的相对误差nargin 函数实际输入参数个数realmax 最大的正实数nargout 函数实际输出参数个数 MATLAB运算符，通过下面几个表来说明MATLAB的各种常用运算符表2 MATLAB算术运算符操作符 功能说明 操作符 功能说明 + 加 矩阵左除- 减. 数组左除* 矩阵乘/ 矩阵右除.* 数组乘./ 数组右除 矩阵乘方 矩阵转置. 数组乘方. 数组转置表3 MATLAB关系运算符操作符 功能说明 = 等于= 不等于 大于= 大于等于= 小于等于表4 MATLAB逻辑运算符逻辑运算符 逻辑运算 说明 & And 逻辑与| Or 逻辑或 Not 逻辑非Xor 逻辑异或表5 MATLAB特殊运算符号 功能说明示例 符号 功能说明示例 ： 1:1:4;1:2:11 . ； 分隔行 . ， 分隔列 （） % 注释 构成向量、矩阵 ！ 调用操作系统命令 构成单元数组= 用于赋值 4.MATLAB的一维、二维数组的寻访表6 子数组访问与赋值常用的相关指令格式指令格式 指令功能 A(r,c) 数组A中r指定行、c指定列之元素组成的子数组A(r,:) 数组A中r指定行对应的所有列之元素组成的子数组A(:,c) 数组A中c指定列对应的所有行之元素组成的子数组A(:) 数组A中各列元素首尾相连组成的“一维长列”子数组A(i) 一维长列子数组中的第i个元素A(r,c)=Sa 数组A中r指定行、c指定列之元素组成的子数组的赋值A(:)=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(B)B阵的逆乘s A.n A的每个元素自乘n次An A阵为方阵，自乘n次A.p 对A各元素分别求非整数幂Ap 方阵A的非整数乘方A+B 对应元素相加A+B 矩阵相加A-B 对应元素相减A-B 矩阵相减A.*B 对应元素相乘A*B 内维相同矩阵相乘A./B A的元素别B的对应元素除A/B A右除B B.A 与上相同BA A左除B exp(A) 以自然数e为底，分别以A的元素为指数，求幂expm(A)A的矩阵指数函数log(A) 对A的各元素求对数logm(A) A的矩阵对数函数sqrt(A) 对A的各元素求平方根sqrtm(A)A的矩阵平方根函数 6.MATLAB的常用函数表8 标准数组生成函数指令 含义 指令 含义 diag 产生对角形数组（对高维不适用）rand 产生均匀分布随机数组eye 产生单位数组（对高维不适用） randn产生正态分布随机数组magic 产生魔方数组（对高维不适用） zeros产生全0数组ones 产生全1数组size()返回指定矩阵的行数和列数表9 数组操作函数指令 含义 diag 提取对角线元素，或生成对角阵flipud以数组“水平中线”为对称轴，交换上下对称位置上的数组元素fliplr 以数组“垂直中线”为对称轴，交换左右对称位置上的数组元素reshape在总元素数不变的前提下，改变数组的“行数、列数” rot90 矩阵逆时针旋转90度det 方阵的行列式值rank 矩阵的秩三、实验内容 1、学习安装MATLAB软件。2、学习使用help命令，例如在命令窗口输入help eye，然后根据帮助说明，学习使用指令eye（其它不会用的指令，依照此方法类推） 3、学习使用clc、clear，观察command window、command history和workspace等窗口的变化结果。4、初步程序的编写练习，新建M-file，保存（自己设定文件名，例如exerc1、exerc2、 exerc3），学习使用MATLAB的基本运算符、数组寻访指令、标准数组生成函数和数组操作函数。注意：每一次M-file的修改后，都要存盘。练习A： （1） help rand，然后随机生成一个26的数组，观察command window、command history和workspace等窗口的变化结果。（2）学习使用clc、clear，了解其功能和作用。（3）输入C=1:2:20，则C（i）表示什么？其中i=1,2,3,10。（4）输入A=7 1 5；2 5 6；3 1 5，B=1 1 1; 2 2 2; 3 3 3，在命令窗口中执行下列表达式，掌握其含义： A(2, 3) A(:,2) A(3,:) A(:,1:2:3) A(:,3).*B(:,2) A(:,3)*B(2,:) A*B A.*B A2 A.2 B/A B./A（5）二维数组的创建和寻访，创建一个二维数组（48）A，查询数组A第2行、第3列的元素，查询数组A第2行的所有元素，查询数组A第6列的所有元素。（6） 两种运算指令形式和实质内涵的比较。设有3个二维数组A24，B24，C22，写出所有由2个数组参与的合法的数组运算和矩阵指令。（7） 学习使用表8列的常用函数（通过help方法） （8） 学习使用表9数组操作函数。 练习B（9）用reshape指令生成下列矩阵（不超过3次操作），并取出方框内的数组元素。四、试验结果（1）（2） clc:清除命令窗口所有内容，数值不变；clear:初始化变量的值。（3）C(i)表示C中的第i个的数值。（4） （5）（6） （9）实验二 MATLAB数值及符号运算一、 实验目的：1、掌握矩阵的基本运算2、掌握矩阵的数组运算3、掌握多项式的基本运算4、会求解代数方程5、掌握创建符号表达式和矩阵的方法6、掌握符号表达式的微分和积分运算二、 实验基本知识：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指令）画出分别显示该系统输入、输出信号的两张子图四、实验结果（1） （2） （3） （4）（5）（7） （6）实验三 MATLAB语言的程序设计一、实验目的： 1、熟悉MATLAB程序编辑与设计环境2、掌握各种编程语句语法规则及程序设计方法 3、函数文件的编写和设计 4、了解和熟悉跨空间变量传递和赋值 二、实验基本知识： 1、程序流程控制语句 for循环结构语法：for i初值：增量：终值 语句1 语句n end 说明：1i初值：终值，则增量为1。 2初值、增量、终值可正可负，可以是整数，也可以是小数，只须符合数学逻辑。while 循环结构 语法：while 逻辑表达式 循环体语句 end 说明：1、whiIe结构依据逻辑表达式的值判断是否执行循环体语勾。若表达式的值为真，执行循环体语句一次、在反复执行时，每次都要进行判断。若表达式的值为假，则程序执行end之后的语句。2、为了避免因逻辑上的失误，而陷入死循环，建议在循环体语句的适当位置加break语句、以便程序能正常执行。（执行循环体的次数不确定；每一次执行循环体后，一定会改变while后面所跟关系式的值。） 3、while循环也可以嵌套、其结构如下： while逻辑表达式1 循环体语句1 while逻辑表达式2 循环体语句2 end 循环体语句3 end if-else-end分支结构if 表达式 1 语句1 else if 表达式 2(可选) 语句2 else(可选) 语句3 end end 说明：1.if结构是一个条件分支语句，若满足表达式的条件，则往下执行；若不满足，则跳出if结构。2else if表达式2与else为可选项，这两条语句可依据具体情况取舍。3注意：每一个if都对应一个end，即有几个if，记就应有几个end。switch-case结构语法：switch表达式case常量表达式1 语句组1 case常量表达式2 语句组2 otherwise 语句组n end 说明：1switch后面的表达式可以是任何类型，如数字、字符串等。2当表达式的值与case后面常量表达式的值相等时，就执行这个case后面的语句组如果所有的常量表达式的值都与这个表达式的值不相等时，则执行otherwise后的执行语句。程序流程控制指令break、return、pause 2、函数文件的结构和编写方法 函数只能在函数体内对变量进行操作，也就是只能访问函数本身工作空间中的变量。M函数文件的结构（1） 函数定义行（function） （2） H1行（函数帮助文本的第一行） （3） 函数帮助文本（4） 函数体（5） 注释函数文件编写后，保存时，其文件名必须与函数名相同。注意：函数名不要与MATLAB自身的函数命令相同。三、实验内容： 练习A 1、熟悉MATLAB程序编辑与设计环境2、用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、已知一维数组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) 为同阶矩阵。四、实验结果（1）（2）（3） （4） （5）（6）（7） 实验四 采用SIMULINK的系统仿真 一、实验目的： 1、熟悉SIMULINK 工作环境及特点2、掌握线性系统仿真常用基本模块的用法3、掌握SIMULINK 的建模与仿真方法二、实验基本知识： 1了解SIMULINK模块库中各子模块基本功能2. SIMULINK 的建模与仿真方法（1）打开模块库，找出相应的模块。鼠标左键点击相应模块，拖拽到模型窗口中即可。（2）创建子系统：当模型大而复杂时，可创建子系统。（3）设置仿真控制参数。三、实验内容 练习 SIMULINK仿真实际应用.1. 利用SIMULINK求解微分方程 ，方程的初始条件为x(0)=1，。在增益模块Gain取值分别为2和100的情况下（即运行，给出运行结果。2. 已经系统框图，求传递函数，并求阶跃响应。3、用两种滤波器设计方法实现下列问题：构建一个低通滤波系统的SIMULINK模型。输入信号是一个受正态噪声干扰的采样信号在此 ，而 （秒）采用10阶Butterworth低通滤波器，以便从输入信号中过滤获得10Hz的输出信号。四、实验结果（1） （2） （3） 幅频，相频响应曲线： ， 冲击响应：实验五 采用SIMULINK的系统仿真（二） 一、实验目的： 1、熟悉SIMULINK 工作环境及特点2、掌握线性系统仿真常用基本模块的用法3、掌握SIMULINK 的建模与仿真方法4、子系统的创建和封装设计二、实验基本知识： 1了解SIMULINK模块库中各子模块基本功能2. SIMULINK 的建模与仿真方法（1）打开模块库，找出相应的模块。鼠标左键点击相应模块，拖拽到模型窗口中即可。（2）创建子系统：当模型大而复杂时，可创建子系统。（3）设置仿真控制参数。三、实验内容 1、用两种滤波器设计方法实现下列问题：构建一个低通滤波系统的SIMULINK模型。输入信号是一个受正态噪声干扰的采样信号在此 ，而 （秒）采用10阶Butterworth低通滤波器，以便从输入信号中过滤获得10Hz的输出信号。2、利用Simulink求解在区间0,1的积分，并求出积分值I(1).)()()62sin()(52.022tutytedtdyedttydtt=+-3、时变系统的仿真：对时变受控对象模型 考虑一个控制系统模型如图所示，其中 均为常数，试分析闭环系统的阶跃响应曲线。 实验六 MATLAB的图形绘制 一、实验目的： 1、学习MATLAB图形绘制的基本方法； 2、熟悉和了解MATLAB图形绘制程序编辑的基本指令； 3、熟悉掌握利用MATLAB图形编辑窗口编辑和修改图形界面，并添加图形的各种标注； 4、掌握plot、subplot的指令格式和语法。二、实验基本知识： 1、plot(x,y)-绘制由x,y所确定的曲线； 2、多组变量绘图：plot(x1, y1, 选项1, x2, y2, 选项2, )； 3、图形窗口的分割； 4、图形编辑窗口的使用。 三、实验内容 练习A 二维曲线绘制的基本指令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,sn) 3. 单输入量绘线plot(Y) 当X、Y均为（mn）数组时，将绘制出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(Y,t) 【2】用图形表示连续调制波形Y=sin(t)sin(9t) 及其包络线。 t=(0:pi/100:pi);y1=sin(t)*-1,1;y2=sin(t).*sin(9*t);t3=pi*(0:9)/9;y3=sin(t3).*sin(9*t3); plot(t,y1,r:,t,y2,b,t3,y3,bo) axis(0,pi,-1,1)【3】在一个图形窗口绘制正弦和余弦曲线，要求给图形加标题“正弦和余弦曲线”，X轴Y轴分别标注为“时间t”和“正弦、余弦”，在图形的某个位置标注“sin(t)”“cos(t)”，并加图例，显示网格，坐标为正方形坐标系。t=0:pi/100:2*pi;y1=sin(t);y2=cos(t);plot(t,y1,r -,t,y2,k -); legend(sin(t),cos(t); axis square;grid on title(sint cost) text(1,-0.4,cos(t)rightarrow) text(2.8,0.4,leftarrowsin(t) xlabel(t);ylabel(正弦、余弦)实验七 MATLAB在信号与系统中的应用一、实验目的：1、掌握利用MATLAB求解LTI系统的方法； 2、掌握用MATLAB中滤波分析的函数命令； 3、熟悉MATLAB信号处理工具箱； 二、实验内容：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)波形。三、实验结果：（1）（2）（3）实验八：频率采样法设计FIR数字滤波器数字滤波器从滤波器功能上可以分为低通、高通、带通、带阻以及全通滤波器；根据系统的单位冲击响应的特性，又可以分为有限长单位冲激响应滤波器（FIR）和无限长单位冲激响应滤波器（IIR）。与IIR滤波器相比，FIR滤波器在保证幅度特性满足技术要求的同时，很容易做到有严格的线性相位特性。设FIR滤波器单位脉冲响应h(n)的长度为N,其系统函数为：Hz= n=0N-1h(n)z-n当滤波器的系数N满足一定的对称条件时，就可以获得线性相位。线性相位FIR滤波器共分为四种类型，分别为：1) 类型I，系数对称，即hn=h(N-1-n),N为奇数。2) 类型II，系数对称，即hn=h(N-1-n)，N为偶数。3) 类型III，系数反对称，即hn=-h(N-1-n)，N为奇数。4) 类型IV，系数反对称，即hn=-h(N-1-n)，N为偶数。1）频率采样法设计滤波器时，通常给出的是幅频特性的技术指标要求，可直接在频域进行处理，按照理想的频率特性H(ej)，在=0到2之间等间隔采样N点，得到：Hk=Hej | =2kN， k=0,1,2,N-1 然后用H(k)的傅里叶逆变换作为滤波器的系数：bn=hn=IDFTH(k) 构成一个系统传递函数为Hz=n=0N-1h(n)z-n的实际的FIR数字滤波器，这种设计方法称为频率采样法。 在使用频率采样法设计FIR数字滤波器时，应注意下列问题： （1）根据频域抽样定理，被采样的理想频率特性其采样点数N与滤波器的长度M应满足NM，否则将造成混叠。本实验取N=M。 （2）为保证滤波器的系数为实序列，作为复数序列的理想频率特性应具有共轭对称性，幅度特性应为偶函数，相位特性应为奇函数。注意：必须在02的全频段上才能观察到其对称图形。习惯上，我们一般利用其对称性，只做0频段上的图形。 （3）理想频率特性的相位特性应该与频率成线性关系，即满足线性相位的条件。2）频率采样法的优化设计为了提高阻带的衰减，减小通带的波动，可以采用频率采样的优化设计法，即在间断点区间内插入一个或几个过渡带采样点。例1. 用频率采样法设计一个FIR数字低通滤波器，3dB截止频率c=0.4，采样点数分别取N=21点和N=61点，分别显示其幅频特性和脉冲响应曲线，观察采样点数对滤波器特性的影响。参考程序如下：N = 21; %采样点数n = 0:N-1;wc = 0.4*pi; %截止频率N1 = fix(wc/(2*pi/N); %样点间隔为2*pi/N，N1为wc的样点数N2 = N - 2*N1 - 1; %N2为阻带样点数A = ones(1,N1+1),zeros(1,N2),ones(1,N1); %建立幅度特性样本序列theta = -pi*(N-1)/N * 0 : N-1; %建立相位特性样本序列wa = 0:N-1/N*2; %为作图建立对应的频率向量%画图subplot(211);plot(wa,A,.-);axis(0,2,-0.2,1.2);title(理想幅频响应及样点序列（N=21）);xlabel(频率（单位：pi）);ylabel(H(ejomega);subplot(212);plot(wa,theta,.-);axis(0,2,-80,0);title(理想相频响应及样点序列（N=21）);xlabel(频率（单位：pi）);ylabel(phi(omega);由图形应可见：该滤波器特性符合幅度特性为偶函数，相位特性为奇函数的特点。后面部分示例不再使用02的全频段，仅在0频段上显示特性曲线。将理想的幅频特性和实际的幅频特性在同一图形中显示，同时显示其脉冲响应曲线，观察采样点数对滤波器特性的影响。此时可通过对上面的程序进行修改实现目的：%N = input(N=); %由使用者输入NN = 41;n = 0:N-1;wc = 0.4*pi; %截止频率N1 = fix(wc/(2*pi/N); %样点间隔为2*pi/N,N1为wc的样点数N2 = N - 2*N1 -1; %N2为阻带样点数A = ones(1,N1+1),zeros(1,N2),ones(1,N1); %建立幅度特性样本序列theta = -pi*(N-1)/N *0: N-1; %建立相位特性样本序列Hk = A.*exp(j*theta); %建立频域特性样本序列h = real(ifft(Hk); %由反变换求脉冲序列h(n)，并取实部H,w = freqz(h,1,1000,whole); %计算滤波器的频率响应H1 = (H(1:1000); %对矩阵H取转置w1 = (w(1:1000);mag = abs(H1);db = 20*log(mag+eps)/max(mag); %将单位转换为dBpha = angle(H);wa = 0:N-1/N*2;%画图subplot(311);plot(wa,A,.-,w1/pi,mag);axis(0,2,-0.1,1.2);title(理想幅频、样点序列及实际滤波器幅频响应);xlabel(频率（单位：pi）);ylabel(H(ejomega);subplot(312);stem(n,h);title(滤波器脉冲响应);xlabel(n);ylabel(h(n);subplot(313);plot(w1/pi,pha),title(实际滤波器的相频响应);N = 21N = 61例2：在满足上面显示理想/实际滤波器幅频响应等的设计要求下，在过渡带中增加一个样点T1，取值0.42。要求显示其幅频特性曲线，观察增加了过渡带采样点后对滤波器特性的影响。提示：只需将建立幅频特性样本序列的一句程序修改结果：在过渡带增加了一点之后，通带波动迅速减小，阻带衰减也得到了很大的提高。增加过渡带采样点前后对比例3：与例2类似，从增加一个过渡带采样点变为增加两个采样点，T1=0.58，T2=0.12（这两个值可以自己尝试修改，并观察效果）。要求显示幅频特性曲线，观察增加两个过渡带采样点后对滤波器特性的影响。提示：只需将建立幅频特性样本序列的一句程序修改结果：在过渡带增加了两点之后，通带波动明显下降，阻带衰减明显提高，用较少的采样点数已经能够实现较高的滤波性能指标要求。例4：用频率采样法设计一个FIR数字高通滤波器（截止频率自己定，0）（可选择添加过渡带采样点）；设计一个FIR数字带通滤波器，要求：下阻带截止频率ws1=0.3,（As=35dB）;通带低端起始频率wp1=0.45，通带高端截止频率wp2=0.65，（Rp=1.5dB）;上阻带起始频率ws2=0.8.描绘实际滤波器的幅频响应、脉冲响应曲线。（可选择添加过渡带采样点）；【扩展：检验最小阻带衰减、检验通带波动例：As=-round(max(db(1:fix(ws/dw)+1)Rp=-(min(db(fix(wp/dw)+1:501)】例5：用fir2设计一个N=41（可修改）的FIR数字低通滤波器，要求通带起始频率ws=0.5(可修改)，描绘理想和实际滤波器的幅频响应曲线。参考程序如下：wc = 0.5;N = 41;n = 0:N-1;f = 0 wc wc 1; %建立理想幅频特性频率向量m = 1 1 0 0; %建立理想幅频特性幅度向量b = fir2(N-1,f,m) %计算滤波器系统函数的系数H,w = freqz(b,1); %求解频率特性dbH = 20*log10(abs(H)+eps)/max(abs(H); %单位化为分贝subplot(2,2,1);plot(f,m,w/pi,abs(H),title(理想和实际的高通滤波器幅频特性);xlabel(频率（单位：pi）);ylabel(H(ejomega);subplot(222);stem(n,b),title(滤波器脉冲响应);xlabel(n);ylabel(h(n);subplot(223);plot(w/pi,dbH),title(实际的高通滤波器幅频特性（dB）);xlabel(频率（单位：pi）);ylabel(G(dB);subplot(224);plot(w/pi,angle(H),title(实际的高通滤波器相频特性);xlabel(频率（单位：pi）);ylabel(phi(omega);扩展：1）将高通滤波器改为高通滤波器；2）将设计指标改为：通带截止频率ws=0.4，阻带起始频率为0.6，描绘理想和实际滤波器的幅频响应曲线。3）设计一个N=61的FIR数字带阻滤波器，要求：通带低端截止频率wc1=0.4，通带高端起始频率wc2=0.6.描绘理想和实际滤波器的幅频响应曲线。4）将上面的带阻滤波器设计指标改为：下通带截止频率wp1=0.3,阻带低端起始频率ws1=0.4，阻带高端截止频率ws2=0.6，上通带起始频率wp2=0.7。且在过渡带增加一点T1=0.38.附：matlab基础1. Matlab桌面系统简介：l 主界面的布局：Current folder, command window, workspace, command history, editor； l 常用工具栏工具;常用快捷键(ctrl+s, ctrl+z, clear; clc; clear all; ctrl+a; ctrl+I; ctrl+r; ctrl+t; ctrl+d; ctrl+w; close all等)；l 如何快速找到文件所在位置/打开指定位置的文件；当前打开的文件所在目录与目前指向路径current folder不一致时需要将指向路径变换到当前文件所在路径；l 新建脚本文件script/子函数function，脚本文件要和其需要调用的函数文件在同一目录下，子函数的函数名与文件名要一致。2. matlab帮助系统：l 帮助命令：help + 函数名（edit 函数名，可以查看该库函数的代码）；如果不知道具体的函数名，则直接在command window输入help，在弹出的列表中根据所需要的函数所属类型选择查看帮助内容（很多时候也可以直接百度）； (例：reshape，linspace，zeros，abs，find，length。)Lookfor +关键字 ，关键字搜索；l 导航浏览器帮助：doc + 函数名；l 演示系统帮助：demos （视频，脚本文件等例程）。3. 调试功能M程序调试相对简单，因为它没有指针的概念，这样就避免了一大类难以查找的错误。当运行脚本文件出错时，程序会中断执行，这时在command window窗口中会提示程序哪里有误，并有对错误的简单描述（都是英文提示，但是很多时候对调试很有帮助，能告知我们是什么类型的错误，最好养成看错误提示的习惯）。例如下图所示：上图中的test_Viterbi表示出错的函数，单击test_Viterbi会弹出对该函数的描述，后面的line 9表示出错的具体位置（行数），单击它会跳转到程序出错的那一行，从而便于我们找出错误的原因并对程序进行修改。l 断点设置断点可以使程序调试变得更为简单高效。 设置或清除断点：在你想添加断点的那一行程序前的短横线处单击即可 或 使用快捷键F12 （个人觉得单击添加断点更方便）（注意：只有程序行才可以添加断点，空白行或注释行无法添加断点） 执行点击Run按钮 或 用快捷键F5。当设置了断点后，程序若运行到断点处，Run按钮会变成continue 。点击continue，程序继续执行。若未变成continue，则说明断点前面有错误，可以根据command window中的提示去修改。4. M文件 脚本M文件：处理较为简单且重复性少的问题；共用matlab基本空间 函数M文件：处理重复性较多的问题，文件的第一有效行为函数声明行；开辟一个专用的临时工作空间上图为一个脚本M文件下图为一个函数