1、1Matlab 语言及其在电子信息科学中的应用实验指导书实验一：Matlab 操作环境熟悉一、实验目的1初步了解 Matlab 操作环境。2学习使用图形函数计算器命令 funtool 及其环境。二、实验内容熟悉 Matlab 操作环境，认识命令窗口、内存工作区窗口、历史命令窗口；学会使用 format 命令调整命令窗口的数据显示格式；学会使用变量和矩阵的输入，并进行简单的计算；学会使用 who 和 whos 命令查看内存变量信息；学会使用图形函数计算器funtool，并进行下列计算：1单函数运算操作。 求下列函数的符号导数(1) y=sin(x);(2) y=(1+x)3*(2-x); 求下列
2、函数的符号积分(1) y=cos(x);(2) y=1/(1+x2);(3) y=1/sqrt(1-x2);(4) y=(x-1)/(x+1)/(x+2); 求反函数(1) y=(x-1)/(2*x+3);(2) y=exp(x);(3) y=log(x+sqrt(1+x2); 代数式的化简(1) (x+1)*(x-1)*(x-2)/(x-3)/(x-4);(2) sin(x)2+cos(x)2;(3) x+sin(x)+2*x-3*cos(x)+4*x*sin(x);2函数与参数的运算操作。 从 y=x2 通过参数的选择去观察下列函数的图形变化(1) y1=(x+1)2(2) y2=(x+2
3、)2(3) y3=2*x2(4) y4=x2+2(5) y5=x4(6) y6=x2/23两个函数之间的操作 求和(1) sin(x)+cos(x)(2) 1+x+x2+x3+x4+x5 乘积(1) exp(-x)*sin(x)2Matlab 语言及其在电子信息科学中的应用实验指导书(2) sin(x)*x 商(1) sin(x)/cos(x);(2) x/(1+x2);(3) 1/(x-1)/(x-2); 求复合函数(1) y=exp(u) u=sin(x)(2) y=sqrt(u) u=1+exp(x2)(3) y=sin(u) u=asin(x)(4) y=sinh(u) u=-x三、设
4、计提示1初次接触 Matlab 应该注意函数表达式的文本式描述。2在使用图形函数计算器 funtool 时，注意观察 1 号和 2 号窗口中函数的图形。四、实验报告要求1针对图形函数计算器 funtool，对每一类型计算记录其中一个图形的曲线。2书写实验报告时要结构合理，层次分明，在分析描述的时候，需要注意语言的流畅。3Matlab 语言及其在电子信息科学中的应用实验指导书实验二：M 文件和 Mat 文件操作一、实验目的1定制自己的工作环境。2编写简单的 M 文件。3保存内存工作区中的变量到.mat 文件。4学会只用 Matlab 帮助。二、实验内容1使用 format 命令和 File|Pe
5、ferences 菜单定制自己的工作环境。2编写如下 M 文件，试调整参数 a 的大小，观察并记录 y1、y2 的波形特征。%example1.mt=0:pi/100:4*pi;a=3;y2=exp(-t/a);y1=y2.*sin(a*t);plot(t,y1,-r,t,y2,:b,t,-y2,:b);3保存内存工作区变量 a、t、y1、y2 到 example1.mat 文件；关闭 Matlab，再重新启动；观察内存工作区；重新根据.mat 文件恢复原来的工作区变量。4在命令窗口中查看 exp 函数的帮助；运行 helpwin 查看超文本格式的帮助文件，试翻译并记录下信号处理工具箱（Sig
6、nal Processing Toolbox）中的函数分类（Functions - Categorical List） 。三、设计提示1可以用命令语句、菜单或按钮等多种方式执行命令。2用于编辑 M 文件的文本编辑器还可以执行和调试程序。3不同的工具箱可能包含同名的函数，查看帮助时应注意在左侧栏选择相应的工具箱类别。四、实验报告要求1对实验内容 2，说明参数 a 的大小对 y1、y2 波形特征的影响。2翻译命令窗口中的 exp 函数的帮助信息。3运行 helpwin，试翻译并记录下信号处理工具箱（Signal Processing Toolbox）中的函数分类（Functions - Categ
7、orical List） 。4书写实验报告时要结构合理，层次分明，在分析描述的时候，需要注意语言的流畅。4Matlab 语言及其在电子信息科学中的应用实验指导书实验三：矩阵运算与元素群运算一、实验目的1掌握数组与矩阵的创建。2掌握矩阵运算与数组运算。3掌握基本元素群运算。4掌握向量与矩阵的特殊处理。二、实验内容1 “：”号的用法。用“：”号生成行向量 a=1 2 3 4 5 6 7 8 9 10、b=5 3 1 -1 -3 -5；2用线性等分命令 linspace 重新生成上述的 a 和 b 向量。3在 100 和 10000 之间用对数等分命令 logspace 生成 10 维的向量 c。4
8、生成范围在0，10 、均值为 5 的 35 维的均匀分布随机数矩阵 D。5利用 magic 函数生成 55 维的魔方矩阵，取其对角向量 e，并根据向量 e 生成一个对角矩阵 E。 （所谓魔方矩阵就是各行、各列、各对角线元素之和相等。 ）6另 AA 是 33 维魔方矩阵，BB 是由 A 旋转 180得到。CC 是一个复数矩阵，其实部为 AA，虚部为 BB。DD 是 CC 的转置，EE 是 CC 的共轭。分别计算CC 和 EE 的模和幅角。7是一个首项为 20，公比为 0.5 的 10 维等比数列；g 是一个首项为 1，公差为3 的 10 维等差数列。试计算向量 f 和 g 的内积 s。8生成一个
9、 99 维的魔方矩阵，提取其中心的 33 维子矩阵 M，利用 sum 函数检验其各行和各列的和是否相等。9已知 ，利用函数生成左上三角矩阵 。1234567T 12340T三、设计提示1等比数列可利用首项和公比的元素群的幂乘积生成。2提取子矩阵，可灵活应用“：”号或空阵 。3尽量用 Matlab 函数生成上述矩阵或向量，不要用手工逐个输入。四、实验报告要求1编写实现第二节实验内容中所使用的函数命令，并记录相应的生成结果。2思考题：是否存在 22 维的魔方矩阵？。3书写实验报告时要结构合理，层次分明，在分析描述的时候，需要注意语言的流畅。五，参考答案1.代码：a=1:10 5Matlab 语言及
10、其在电子信息科学中的应用实验指导书b=5:-2:-5 a = 1 2 3 4 5 6 7 8 9 10b =5 3 1 -1 -3 -51. 代码:linspace(1,10,10)linspace(5,-5,6)ans =1 2 3 4 5 6 7 8 9 10ans =5 3 1 -1 -3 -52. 代码：logspace(2,4,10)ans =1.0e+004 *0.0100 0.0167 0.0278 0.0464 0.0774 0.1292 0.2154 0.3594 0.5995 1.00003. 代码:D=10*rand(3,5)结果:D =0.1527 9.3181 8.4
11、622 6.7214 6.81287.4679 4.6599 5.2515 8.3812 3.79484.4510 4.1865 2.0265 0.1964 8.31805.代码:A=magic(5);e=diag(A);6Matlab 语言及其在电子信息科学中的应用实验指导书E=diag(e)结果；E =17 0 0 0 00 5 0 0 00 0 13 0 00 0 0 21 00 0 0 0 96.代码:AA=magic(3);BB=rot90(AA,2);CC=AA+j*BB;DD=CC.;EE=conj(CC);A=abs(CC)B=angle(CC)C=abs(EE)D=angle
12、(EE)结果:A =8.2462 9.0554 7.21117.6158 7.0711 7.61587.2111 9.0554 8.2462B =0.2450 1.4601 0.58801.1659 0.7854 0.40490.9828 0.1107 1.3258C =8.2462 9.0554 7.21117.6158 7.0711 7.61587.2111 9.0554 8.24627Matlab 语言及其在电子信息科学中的应用实验指导书D =-0.2450 -1.4601 -0.5880-1.1659 -0.7854 -0.4049-0.9828 -0.1107 -1.32587.代码
13、:f=20*0.5.(0:9);g=1:3:28;s=dot(f,g)结果:s =158.67198.代码：A=magic(9);M=A(4:6,4:6);sum(A,1)sum(A,2)结果：ans =369 369 369 369 369 369 369 369 369ans =3693693693693693693693693699代码：T=1,2,3,4;2,3,4,5;3,4,5,6;4,5,6,7A=rot90(T);B=tril(A);T1=rot90(B,-1)8Matlab 语言及其在电子信息科学中的应用实验指导书结果：T =1 2 3 42 3 4 53 4 5 64 5
14、6 7T1 =1 2 3 42 3 4 03 4 0 04 0 0 09Matlab 语言及其在电子信息科学中的应用实验指导书实验四：线性方程组的求解一、实验目的1掌握恰定方程组的解法。2了解欠定方程组、超定方程组的解法。3掌握左除法求解线性方程组。4学会测试程序的运行时间。二、实验内容1用两种方法求下列方程组的解，并比较两种方法执行的时间。 1234512345123457 49103562 73710 5xxxxxx2判定下列方程是恰定方程组、欠定方程组还是超定方程组，并求其解。 12345123456916847294216930xxxxxx3用网孔电流法求如下电路的各支路电流。 VIa
15、bcd10Matlab 语言及其在电子信息科学中的应用实验指导书4用结点电压法求如下电路的结点电压 un1、u n2。2k130A+-g=. (S)三、设计提示1在计算程序的执行时间之前，应注意用 clear 命令将内存变量清空。2求得线性方程组的解之后，代入原方程验证是否正确。四、实验报告要求1编写实现第二节实验内容中所使用的函数命令，并记录相应的生成结果。2对于电路的求解，应列出相应的网孔方程和结点方程，并注意方向。3书写实验报告时要结构合理，层次分明，在分析描述的时候，需要注意语言的流畅。五、参考答案第一题：A=7,14,-9,-2,5;3,-15,-13,-6,-4;-11,-9,-2
16、,5,7;5,7,14,16,-2;-2,5,12,-11,-4;B=100;200;300;400;500;X=AB结果:X =370.9455-224.7276238.3333-138.1891503.637811Matlab 语言及其在电子信息科学中的应用实验指导书第二题：A=6,9,14,-11,5;1,14,-7,-15,-6;-2,1,-7,12,-1;6,11,11,-9,-13;B=68;294;-441;103;X=AB结果：X =0-26.3759-5.4126-38.3960-8.2392第三题：A=1,1,-1,1;0,20,40,0;60,0,40,0;0,0,40,
17、40;B=0;10;50;-40;X=ABX =0.78570.35710.0714-1.071412Matlab 语言及其在电子信息科学中的应用实验指导书实验五：函数编写与程序设计一、实验目的1掌握函数的编写规则。2掌握函数的调用。3会用 Matlab 程序设计实现一些工程算法问题。二、实验内容1编写一个y,y1,y2=mwave(f1,m1,f2,m2)函数，实现以下功能，并绘出y1、y2、y 在 t0，2区间 500 个样点的图形。 （其中调用参数 2 f1、f2 20 HZ；0.5 m1、m2 2）1122sin();ymft2程序设计：相传古代印度国王要褒奖他的聪明能干的宰相达依尔（
18、国际象棋发明者） ，问他要什么？达依尔回答：“陛下只要在国际象棋棋盘的第一个格子上放一粒麦子，第二个格子上放二粒麦子，以后每个格子的麦子数都按前一格的两倍计算。如果陛下按此法给我 64 格的麦子，就感激不尽，其他什么也不要了。 ”国王想：“这还不容易！”让人扛了一袋麦子，但很快用光了，再扛出一袋还不够，请你为国王算一下共要给达依尔多少小麦？（1 袋小麦约1.4108 粒） 。3程序设计：公元前五世纪我国古代数学家张丘建在算经一书中提出了“百鸡问题”：鸡翁一值钱五，鸡母一值钱三，鸡雏三值钱一。百钱买百鸡，问鸡翁、母、雏各几何？三、设计提示1函数名和函数文件名应相同；对调用参数的取值范围要检验是否
19、符合要求，如不符合要求，应给出出错信息（用 if 和 error 函数实现） 。2程序设计“百鸡问题”答案不唯一。提示：设 x：鸡翁数，则 x 的范围：019；y：鸡母数，则 y 的范围：033；z：鸡雏数，则 z 的范围：0100。四、实验报告要求1编写实现第二节实验内容（1）所使用的函数文件。2程序设计用 M 文件编写，并记录执行结果。 “百鸡问题”答案不唯一，要给出所有答案。3书写实验报告时要结构合理，层次分明，在分析描述的时候，需要注意语言的流畅。13Matlab 语言及其在电子信息科学中的应用实验指导书五、参考答案第 1 题：function y,y1,y2=mwave(f1,m1,
20、f2,m2)% 根据给定的频率和幅值计算标准正弦函数 y1、y2 及其叠加 y 的波形。if (f120) error(f1 超出范围！ ), return, endif (f220) error(f2 超出范围！ ), return, endif (m12) error(m1 超出范围！), return, endif (m22) error(m2 超出范围！), return, end% -t=0:2*pi/(500-1):2*pi;y1=m1*sin(2*pi*f1*t);y2=m2*sin(2*pi*f2*t);y=y1+y2;figuresubplot(311); plot(t,y1
21、); title(y1 波形 );subplot(312); plot(t,y2); title(y2 波形 );subplot(313); plot(t,y); title(y=y1+y2 波形);% =第 2 题：a=1;s=0for i=1:64s=s+a;a=2*a;endn=s/1.4/108运行后得：n=1.3176e+011第 3 题：解 设 x:鸡翁数，则 x 的范围：019y:鸡母数，则 y 的范围：033z:鸡雏数，则 z 的范围：0100则： x+y+z=1005x+3y+z/3=100这是一个欠定方程。for x=0:19for y=0:33for z=0:100if
22、(x+y+z=100)figure(1)subplot(221);plot(t,r);title(连线图)subplot(222);stem(t,y);title(脉冲图);subplot(223);stairs(t,b);title(阶梯图);subplot(224);bar(t,g);title(条形图);17Matlab 语言及其在电子信息科学中的应用实验指导书0 5 1000.20.40.60.81 冲冲冲0 5 1000.20.40.60.81 冲冲冲0 5 1000.51 冲冲冲1 2 3 4 5 6 7 8 91000.51 冲冲冲第二题：x=1:1:50;y1=2.x;y2=(
23、1/2).x;figure(1)plot(x,y1,r,x,y2,g);gtext(y1=2.x);gtext(y2=(1/2).x);18Matlab 语言及其在电子信息科学中的应用实验指导书0 5 10 15 20 25 30 35 40 45 50024681012x 1014y1=2.xy2=(1/2).x第 3 题：function y,t=mcircle(r)% 根据给定的半径 r，以原点为圆心画一个红色空心圆。t=0:2*pi/64:2*pi;y=r*ones(size(t);subplot(121), polar(t,y,*r)X,Y=pol2cart(t,y);% 或者采用如
24、下方法转换% X=r*cos(t);% Y=r*sin(t);subplot(122), plot(X,Y,*r)axis equal;axis square;第 4 题：t=0:pi/180:2*pi*5;r1=10/2;x1=r1*cos(t);y1=r1*sin(t);z=t/(2*pi);subplot(121), plot3(x1,y1,z)grid on19Matlab 语言及其在电子信息科学中的应用实验指导书%-r2=linspace(5,0,length(t);x2=r2.*cos(t);y2=r2.*sin(t);subplot(122), plot3(x2,y2,z)gri
25、d on20Matlab 语言及其在电子信息科学中的应用实验指导书实验七：Matlab 多项式和符号运算一、实验目的1掌握 Matlab 多项式的运算。2了解符号运算。二、实验内容（边做实验，边将生成结果和图形拷贝到 Word 文档中）1将多项式 化为 x 的降幂排列。()2)(37)(1Pxx2求一元高次方程的根。 9876543253016820576-80x x3求一元高次方程的根，并画出左边多项式函数在 区间内的曲线。,42x4求多项式 和 的乘积 ；并321()57fxx322()864fx()fx求 的商和余式。2()f5求 的符号导数。52tan(4)3yx6用符号运算求实验内容
26、 4 中的 的表达式。()fx三、设计提示1关于多项式运算的函数有 poly、roots 等。2多项式做加减运算时要注意等长度。3符号表达式的输入可以用字符串方式，也可以用 sym 函数。四、实验报告要求1编写实现第二节实验内容中所使用的函数文件，并记录相应的生成结果和图形。2对于多项式的结果应以多项式向量和多项式表达式两种方式记录。3书写实验报告时要结构合理，层次分明，在分析描述的时候，需要注意语言的流畅。五、参考答案第 1 题：P=poly(2,-3,7,-1)结果：21Matlab 语言及其在电子信息科学中的应用实验指导书P =1 -5 -19 29 42即： 432()5194xx第
27、2 题：P1=1 -5 -30 150 273 -1365 -820 4100 576 -2880;x1=roots(P1)结果：x1 =5.0000-4.00004.0000-3.00003.0000-2.0000-1.00002.00001.0000第 3 题：P2=1 0 -2 0 1;x2=roots(P2)n=1;for x=-2:0.01:2y(n)=sum(P2.*(x.(length(P2)-1):-1:0);% 或者 y(n)=x4-2*x2+1;n=n+1;endx=-2:0.01:2;plot(x,y)结果：(有重根！ )x2 =1.0000 + 0.0000i1.000
28、0 - 0.0000i-1.0000 + 0.0000i-1.0000 - 0.0000i第 4 题：f1=1 3 5 7;f2=8 -6 4 -2;f=conv(f1,f2)f11=zeros(1,length(f)-length(f1),f1 % 补 0，与 f 同维q,r=deconv(f-f11,f2)结果：f =8 18 26 36 -28 18 -14即： 65432()126814xxxxf11 =0 0 0 1 3 5 722Matlab 语言及其在电子信息科学中的应用实验指导书q =1.0000 3.0000 5.0000 6.8750r =0 0 0 0 -3.7500 -
29、4.5000 -7.2500第 5 题：y=x5+tan(4*x2)+3;diff(y)结果：ans =5*x4+8*(1+tan(4*x2)2)*x第 6 题：f1=sym(x3+3*x2+5*x+7);f2=sym(8*x3-6*x2+4*x-2);f=f1*f2collect(f)(f-f1)/f2collect(ans)结果：f =(x3+3*x2+5*x+7)*(8*x3-6*x2+4*x-2)f =8*x6+18*x5+26*x4+36*x3-28*x2+18*x-14h =(-21+8*x6+18*x5+26*x4+35*x3-31*x2+13*x)/(8*x3-6*x2+4*x
30、-2)h =(-21+8*x6+18*x5+26*x4+35*x3-31*x2+13*x)/(8*x3-6*x2+4*x-2)23Matlab 语言及其在电子信息科学中的应用实验指导书实验八：线性时不变系统的时域响应一、实验目的1掌握线性时不变系统的三种描述形式传递函数描述法、零极点增益描述法、状态空间描述法。2掌握三种描述形式之间的转换。3掌握连续和离散系统频率响应的求解。二、实验内容（边做实验，边将生成结果和图形拷贝到 Word 文档中）1生成 20 个点的单位脉冲信号、单位阶跃信号，并记录下函数命令和波形。2生成占空比为 30的矩形波。3将连续系统 转化为传递函数的形式，并显示其表(1)
31、3()0.524)sHs达式。4将离散系统 转化为零极点增益的描述121343().6.0.9.5zz z形式，并显示其表达式。5分别求实验内容 3 和 4 的频率响应（对离散系统取 256 样点，采样频率取8000Hz） 。6分别求实验内容 3 和 4 的单位冲激响应（对离散系统，作 60 样点图） 。三、设计提示1显示传递函数模型用 tf(b,a)；显示零极点增益模型用 zpk(z,p,k)。注意：z、p 为列向量。2连续系统频率响应用 freqs 函数；离散系统用 freqz 函数。3连续系统冲激响应用 impulse 函数；离散系统用 impz 函数。四、实验报告要求1编写实现第二节实
32、验内容中所使用的函数文件，并记录相应的生成结果。2书写实验报告时要结构合理，层次分明，在分析描述的时候，需要注意语言的流畅。五、参考答案第一题：k=-4:15x =zeros(1,4),1,zeros(1,15)24Matlab 语言及其在电子信息科学中的应用实验指导书subplot (1,2,1)stem(k,x)y=zeros(1,4),ones(1,16)subplot (1,2,2)stem(k,y)结果：-5 0 5 10 1500.10.20.30.40.50.60.70.80.91-5 0 5 10 1500.10.20.30.40.50.60.70.80.91第二题f=2;t=
33、0:0.001:2w=2*pi*f*ty=square(w,30)axis(0 2 0 2)25Matlab 语言及其在电子信息科学中的应用实验指导书hold onplot(t,y)结果：0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 200.20.40.60.811.21.41.61.82第 3 题：k=0.5;z=1,-3;p=-1,-2,-4;sys_zpk=zpk(z,p,k)b,a=zp2tf(z,p,k);sys_tf=tf(b,a)结果：Zero/pole/gain:0.5 (s-1) (s+3)-(s+1) (s+2) (s+4)Transfer fun
34、ction:0.5 s2 + s - 1.526Matlab 语言及其在电子信息科学中的应用实验指导书-s3 + 7 s2 + 14 s + 8第 4 题：b=3 5 2;a=1 -1.6 1.3 -.9 .5;systf=tf(b,a,variable,z-1)%显示传递函数z,p,k=tf2zp(b,a);syszpk=zpk(z,p,k,variable,z-1)%显示零点增益结果：Transfer function:3 + 5 z-1 + 2 z-2-1 - 1.6 z-1 + 1.3 z-2 - 0.9 z-3 + 0.5 z-4Sampling time: unspecifiedZ
35、ero/pole/gain:3 z-2 (1+z-1) (1+0.6667z-1)-(1 - 1.685z-1 + 0.8654z-2) (1 + 0.08497z-1 + 0.5778z-2)Sampling time: unspecified第 5 题：k=0.5;z=1,-3;p=-1,-2,-4;b,a=zp2tf(z,p,k);freqs(b,a)结果：27Matlab 语言及其在电子信息科学中的应用实验指导书b=3 5 2;a=1 -1.6 1.3 -.9 .5;freqz(b,a,256,8000) % Fs=8000结果：第 6 题：k=0.5;z=1,-3;p=-1,-2,-
36、4;sys3=zpk(z,p,k);impulse(sys3)%频率响应结果：28Matlab 语言及其在电子信息科学中的应用实验指导书b=3 5 2;a=1 -1.6 1.3 -.9 .5;impz(b,a,60)29Matlab 语言及其在电子信息科学中的应用实验指导书实验九：数字滤波器设计与语音信号处理一、实验目的1掌握 Butterworth 滤波器、 Chebyshev 型、型滤波器的设计，并能用于对语音信号子频带的提取。2了解 FIR 滤波器的设计。二、实验内容（边做实验，边将生成结果和图形拷贝到 Word 文档中）1将 Test1.wav 和 Test2.wav 两个语音文件复制
37、到 C:Matlab6p5work目录下。2设计一个 4 阶的 Chebyshev 型模拟高通滤波器，要求截止频率 ，阻1kHzcf带纹波不超过-30dB。取采样频率 。要求结果给出该滤波器的传递10kHzsF函数表达式，并画出该滤波器的频率响应曲线图。3设计一个 6 阶的 Butterworth 数字带通滤波器，要求截止频率 ，10cf，取采样频率 。要求结果给出该滤波器的传递函数表20Hzcfs达式，并画出该滤波器的频率响应曲线图。4设计一个 6 阶的 Butterworth 数字低通滤波器，从一段含噪声语音信号Test2.wav（04kHz，Fs=8kHz）中提取出 500Hz 以下的子
38、带信号。在一个图形窗口中分上下两个子窗口分别显示原始语音信号波形和滤波后的 500Hz 以下子带信号波形。5设计一个 65 阶的 FIR 数字滤波器，通带（归一化）：00.45、0.650.85。画出该滤波器的频率响应曲线图。三、设计提示1设计滤波器时，注意频率的归一化。2Chebyshev 型模拟高通滤波器设计，用函数 cheby2(n,Rs,Wn,high,s)。其中Rs=30（即表示：阻带纹波不超过 -30dB） 。3设计 IIR 带通或带阻滤波器时，若给滤波器函数的参数为 n，则设计出来的是2n 阶的滤波器。4设计高通滤波器，用high关键字；带通用bandpass关键字。5查看模拟滤
39、波器的频率响应曲线用 freqs 函数，查看数字滤波器的频率响应曲线用 freqz 函数。6对于 FIR 滤波器的设计，第一通带如果从 0 开始，则应加关键字DC-1。四、实验报告要求1编写实现第二节实验内容所使用的 M 文件，并在 Word 文档中记录相应的生成结果和图形。2书写实验报告时要结构合理，层次分明，在分析描述的时候，需要注意语言的流畅。五、参考答案30Matlab 语言及其在电子信息科学中的应用实验指导书第 2 题：fc=1000;Fs=10000;Wn=fc/(Fs/2); %归一化n=4;Rs=30;b,a=cheby2(n,Rs,Wn,high,s);sys=tf(b,a)freqs(b,a);结果：s4 - 1.42e-017 s3 + 0.04 s2 - 2.938e-018 s + 0.0002-s4 + 0.6442 s3 + 0.2475 s2 + 0.05581 s + 0.006325第 3 题：fc1=1000;fc2=3000;Fs=10000;Wn=fc1,fc2/(Fs/2); %归一化n=6;b,a=butter(n/2,Wn,bandpass);sys=tf(b,a)freqz(b,a,512);结果：Transfer function:0.09853 - 0.2956 z-2 + 0.2956 z-4 - 0.09853 z-6
