《Simulink系统仿真课程设计.pdf》由会员分享，可在线阅读，更多相关《Simulink系统仿真课程设计.pdf（23页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、信息系统仿真课程设计信息系统仿真课程设计课程设计报告课程设计报告题目：信息系统课程设计仿真院（系）：信息科学及技术工程学院专业班级：学生姓名：学号：指导教师：信息系统仿真课程设计任务书信息系统仿真课程设计任务书一、设计（调查报告一、设计（调查报告/论文）题目论文）题目MATLAB 仿真设计：（1）自编程序实现动态卷积；（2）用双线性变换法设计 IIR 数字滤波器。Simulink 仿真设计：（1）低通抽样定理的实现；（2）抽样量化编码器的设计。二、设计（调查报告二、设计（调查报告/论文）主要内容论文）主要内容MATLAB 仿真设计：（1）自编程序实现动态卷积：动态演示序列f1=u(n)(0n1
2、1)和f2=e0.1n(0n16)的卷积和的过程。（2）用双线性变换法设计IIR 数字滤波器：采用双线性变换法设计一个巴特沃斯数字低通滤波器，要求：wp0.2，Rp1 dB；ws0.3，As15 dB，滤波器采样频率 Fs100 Hz。自编程序实现上述滤波器的设计，满足各个参数要求。Simulink 仿真设计：（1）低通抽样定理的实现：建立 Simulink 动态仿真模型，设置模块参数及系统仿真参数，进行系统仿真。对仿真结果进行调试、分析，得出结论。（2）抽样量化编码器的设计：设计一个PCM 编码器，当输入一个电平，通过该编码器得到相应的编码输出。三、原始资料三、原始资料电子课件，仿真参考例子
3、。四、要求的设计（调查四、要求的设计（调查/论文）成果论文）成果对于 MATLAB 仿真设计部分：（1）自编卷积程序要出现动态卷积的过程，并且结果正确；（2）双线性变换法设计的 IIR 滤波器应满足参数要求，设计结果合理、正确。对于 Simulink 仿真设计部分：（1）建立 Simulink 仿真模型完成低通抽样定理的实现；（2）建立 Simulink 仿真模型仿真完成 PCM 编码器的设计，调试。当输入一个电平，应得出正确的编码输出。五、进程安排五、进程安排讲述 MATLAB 语言及 Simulink 仿真技术。（2.5 天）布置课设任务。（0.5 天）Simulink 仿真设计。（2.5
4、 天）MATLAB 仿真设计。（2.5 天）验收课程设计结果。（1 天）撰写课程设计报告。（1 天）六、主要参考资料六、主要参考资料1李贺冰等.Simulink 通信仿真教程.北京：国防工业出版社，2006.2刘舒帆等.数字信号处理实验（MATLAB 版）.西安：西安电子科技大学出版社，2008.3党宏社.信号及系统实验（MATLAB 版）.西安：西安电子科技大学出版社，2007.4郭文彬等.通信原理基于 Matlab 的计算机仿真.北京：北京邮电大学出版社，2006.指导教师（签名）指导教师（签名）：20 20年年月月日日目录目录摘要摘要 .5 .5一、一、SimulinkSimulink
5、仿真设计仿真设计 .6 .6 1.1 1.1 低通抽样定理低通抽样定理 .6 .6 1.2 1.2抽样量化编码抽样量化编码 .9 .9二、二、MATLABMATLAB 仿真设计仿真设计 .12 .12 2.1 2.1、自编程序实现动态卷积、自编程序实现动态卷积 .12 .12 2.1.1 2.1.1 编程分析编程分析 .12 .12 2.1.2 2.1.2 自编自编 matlabmatlab 程序：程序：.13 .13 2.1.3 2.1.3 仿真图形仿真图形 .13 .132.1.42.1.4 仿真结果分析仿真结果分析 .15 .15 2.2 2.2 用双线性变换法设计用双线性变换法设计 I
6、IRIIR 数字滤波器数字滤波器 .15 .15 2.2.1 2.2.1 双线性变换法的基本知识双线性变换法的基本知识 .15 .15 2.2.2 2.2.2 采用双线性变换法设计一个巴特沃斯数字低通滤波器采用双线性变换法设计一个巴特沃斯数字低通滤波器.16.16 2.2.3 2.2.3 自编自编 matlabmatlab 程序程序 .16 .16 2.2.4 2.2.4 仿真波形仿真波形 .17 .17 2.2.5 2.2.5 仿真结果分析仿真结果分析 .17 .17三、总结三、总结 .19 .19四、参考文献四、参考文献 .19 .19五、课程设计成绩五、课程设计成绩 .20 .20摘要摘
7、要 Matlab 是一种广泛应用于工程设计及数值分析领域的高级仿真平台。它功能强大、简单易学、编程效率高，目前已发展成为由 MATLAB 语言、MATLAB 工作环境、MATLAB 图形处理系统、MATLAB 数学函数库和 MATLAB 应用程序接口五大部分组成的集数值计算、图形处理、程序开发为一体的功能强大的系统。本次课程设计主要包括 MATLAB 和 SIMULINKL 两个部分。首先利用 SIMULINKL 实现了连续信号的采样及重构，通过改变抽样频率来实现过采样、等采样、欠采样三种情况来验证低通抽样定理，绘出原始信号、采样信号、重构信号的时域波形图。然后利用 SIMULINKL 实现抽
8、样量化编码，首先用一连续信号通过一个抽样量化编码器按照 A 律 13 折线进量化行，观察其产生的量化误差，其次利用折线近似的 PCM编码器对一连续信号进行编码。最后利用 MATLAB 进行仿真设计，通过编程，在编程环境中对程序进行调试，实现动态卷积以及双线性变换法设计 IIR 数字滤波器。本次课程设计加深理解和巩固通信原理、数字信号处理课上所学的有关基本概念、基本理论和基本方法，并锻炼分析问题和解决问题的能力。一、Simulink 仿真设计1.1 低通抽样定理输入信号为频率为 10Hz 的正弦波，观察对于同一输入信号在不同的抽样频率时，恢复信号的不同波形形态。(1)当抽样频率大于信号频率的两倍
9、，（eg:30Hz）建立模型；图 1-1 抽样仿真框图图 1-3参数设置图 1-2 F=30HZ 抽样及恢复波形(2)改 变 抽 样 频 率 为 信 号 频 率 的 两（eg:20Hz）图 1-4 F=20HZ 的抽样及恢复波形(3)改变抽样频率小于信号频率的两倍（eg:5Hz）图 1-5 F=5HZ 的抽样及恢复波形分析：分析：1 1、由仿真图可知当抽样频率大于或等于被抽样信号频率的、由仿真图可知当抽样频率大于或等于被抽样信号频率的 2 2 倍时，抽样倍时，抽样输出通过模拟低通滤波器能够恢复出被抽样信号；输出通过模拟低通滤波器能够恢复出被抽样信号；当抽样频率小于被抽样信号频当抽样频率小于被抽
10、样信号频率的率的 2 2 倍，模拟低通滤波器的输出波形的形状已失真，即不能恢复出原始信号。倍，模拟低通滤波器的输出波形的形状已失真，即不能恢复出原始信号。2 2、验证抽样定理，即一个频带限制在、验证抽样定理，即一个频带限制在 0f0f 内的连续信号内的连续信号 m(t)m(t)，如果取样速，如果取样速率率FsFs大于或等于大于或等于 2f2f，则可以由抽样值序列则可以由抽样值序列m(nTs)m(nTs)无失真地重建原始信号无失真地重建原始信号 m(t)m(t)。即一个频带限制在（即一个频带限制在（0 0，fH）内的时间连续信号）内的时间连续信号m(t)，如果以，如果以T T1秒的间隔对秒的间隔
11、对2 fH它进行等间隔抽样，则它进行等间隔抽样，则m(t)将被所得到的抽样值完全确定。将被所得到的抽样值完全确定。3 3、采样信号通过模拟低通滤波器后，将其截止频率设置为、采样信号通过模拟低通滤波器后，将其截止频率设置为 80HZ80HZ，为了接近，为了接近理想低通滤波器，将其滤波阶数设置的高一点，信号通过低通滤波器后，除了被理想低通滤波器，将其滤波阶数设置的高一点，信号通过低通滤波器后，除了被滤除阻带上的高频分量外，幅度还会有衰减，因此还需通过一个增益放大器。滤除阻带上的高频分量外，幅度还会有衰减，因此还需通过一个增益放大器。4 4、低通巴特沃兹滤波器的阶次要选择适宜，低通巴特沃兹滤波器的阶
12、次要选择适宜，阶次过高或过，阶次过高或过，输出的信号都有输出的信号都有失真另外，低截止频率不能过大，否则会使多余的谐波通过，也会导致复原的信失真另外，低截止频率不能过大，否则会使多余的谐波通过，也会导致复原的信号失真。号失真。5 5、当输入的连续信号仿真时长也会影响信号的波形，当输入的连续信号仿真时长也会影响信号的波形，仿真是长越短越好，仿真是长越短越好，输输出的波形越好。出的波形越好。由于计算机处理的信号本质上不是模拟信号，由于计算机处理的信号本质上不是模拟信号，为了更接近模拟信为了更接近模拟信号，就将仿真时长尽可能的减小，但不可能减少到零，减少到零后的波形也有一号，就将仿真时长尽可能的减小
13、，但不可能减少到零，减少到零后的波形也有一定程度上的失真。定程度上的失真。1.2抽样量化编码 1）用一个正弦信号通过一个抽样量化编码器后按照 A 律 13 折线产生量化输出信号，从示波器（Scope）上可以观察到产生的量化误差。正弦信号为幅度为 1，频率为 1Hz 的连续时间信号。建立模型；图 1-6抽样量化编码仿真模块框图 1-7参数设置图图 1-81-8 参数设置参数设置图 1-9量化波形图 1-10量化误差分析：分析：1 1、正弦信号通过一个抽样量化编码器后按照、正弦信号通过一个抽样量化编码器后按照 A A 律律 1313 折线产生量化输出信折线产生量化输出信号，从示波器（号，从示波器（
14、scopescope）上可以观察到产生的量化误差，为了比较量化之前和之）上可以观察到产生的量化误差，为了比较量化之前和之后的正弦信号，后的正弦信号，正弦信号产生器和抽样量化编码器的第正弦信号产生器和抽样量化编码器的第 2 2 个输出端口的输出信号个输出端口的输出信号通过一个复用器连接到示波器通过一个复用器连接到示波器 1 1（scope1scope1）2 2、抽样量化编码器用于产生、抽样量化编码器用于产生 1A1A 律律 3 3 折线，折线，它把正弦信号产生器产生的正弦信号它把正弦信号产生器产生的正弦信号转换成量化信号，转换成量化信号，并且计算这个过程中产生的量化噪声。并且计算这个过程中产生的
15、量化噪声。抽样量化编码的参数设抽样量化编码的参数设置如下：置如下：quantization partitionquantization partition 设置为设置为-1/2-1/4-1/8-1/16-1/32-1/64-1/128-1/2-1/4-1/8-1/16-1/32-1/64-1/1280 1/128 1/64 1/32 1/16 1/8 1/4 1/2 10 1/128 1/64 1/32 1/16 1/8 1/4 1/2 13 3、示波器运行结果，图中平滑曲线是表示抽样之前的正弦信号，折线表示通过、示波器运行结果，图中平滑曲线是表示抽样之前的正弦信号，折线表示通过抽样量化编码之
16、后的信号。抽样量化编码之后的信号。可以看到，可以看到，抽样量化之后的信号及原来的连续信号之抽样量化之后的信号及原来的连续信号之间存在一定的量化误差，间存在一定的量化误差，示波器示波器 1 1 运行结果，运行结果，从中可以看出抽样量化产生的量化从中可以看出抽样量化产生的量化误差。误差。1.3 设计一个 13 折线近似的 PCM 编码器，设计一个 13 折线近似的 PCM 编码器使它能够对取值在-1，1内的归一化信号样值进行编码。（eg：当输入为 843，输出编码应为 11101010）图 1-11 13 折线 PCM 编码器模型模型图1-12 常数发生器模块参数图1-13限幅器模块参数图 1-1
17、4 查表模块参数图 1-15 增益参数图 1-16 延时器参数图 1-17 绝对值参数图 1-18 量化器参数图 1-19 十进制转二进制参数分析：分析：1 1 由常数发生器模块产生一个对由常数发生器模块产生一个对 20482048 个量化台阶的归一化的电压，连续个量化台阶的归一化的电压，连续信号通过限幅器将输入信号的幅度限制在信号通过限幅器将输入信号的幅度限制在 PCMPCM 编码定义的范围编码定义的范围-1-1，11内，内，将延将延迟的门限设置为迟的门限设置为 0 0，用来判断输入电平的极性，其输出作为，用来判断输入电平的极性，其输出作为 PCMPCM 编码的极性码。编码的极性码。样值取绝
18、对值后，通过查表模块进行样值取绝对值后，通过查表模块进行 1313 折线压缩，由于以取绝对值，表格中的折线压缩，由于以取绝对值，表格中的值只取值只取 1313 折线第一象限中的九个分段点的坐标，查表方法为内插法，以便能精折线第一象限中的九个分段点的坐标，查表方法为内插法，以便能精确近似，确近似，然后用增益模块将样值放大然后用增益模块将样值放大 128128 倍，倍，用间距为用间距为 1 1 的的 QuantierQuantier 进行四舍进行四舍五入量化，最后将整数编码为五入量化，最后将整数编码为 7 7 位二进制序列，作为位二进制序列，作为 PCMPCM 编码的低编码的低 7 7 位。位。2
19、 2、在设置增益时，增益为、在设置增益时，增益为 127127 或或 128128 时输出的量化电平正确，时输出的量化电平正确，二、MATLAB 仿真设计2.1、自编程序实现动态卷积2.1.1 编程分析1）掌握离散卷积过程：序列倒置移位相乘求和；2）自编程序实现两个序列的动态卷积。题目：动态演示序列 f1=u(n)(0n11)和 f2=e0.1n(0n16)的卷积和的过程。2.1.2 自编 matlab 程序：clear;clc;%清除数据区，指令区nf1=1:10;lf1=length(nf1);%f1 序列长度f1=ones(1,lf1);%f1 序列nf2=1:15;lf2=length
20、(nf2);%f2 长度f2=exp(-0.1.*nf2);%f2 序列f3=fliplr(f2);%f2 翻转h=zeros(1,lf2+1),f2,zeros(1,lf1-1);%f2 补零u=zeros(1,lf2+1),f1,zeros(1,lf2-1);%f1 补零nf3=-lf2:lf1+lf2-1;%求卷积for k=0:lf1+lf2%卷积的动态显示过程p=zeros(1,k),f3,zeros(1,lf1+lf2-k);%序列补零并且进行循环移位y1=p.*u;%求乘积yk=sum(y1);%各点之和，即一个点卷积结果y(k+1)=yk;%将卷积结果存入 y 序列z=zero
21、s(1,lf2-1),y,zeros(1,lf2+lf1-k);%y 序列进行补零后的卷积序列%绘图subplot(5,1,1);stem(nf3,u,filled);title(u(n);ylabel(u(n);subplot(5,1,2);stem(nf3,h,filled);title(h(n);ylabel(h(n);subplot(5,1,3);stem(nf3,p,filled);title();ylabel(h(n-k);subplot(5,1,4);stem(nf3,y1,filled);title();ylabel(u(n)*h(n);subplot(5,1,5);stem(
22、nf3,z,filled);title(卷积结果);ylabel(y(n);if k=0 disp(开始卷积),pause;End乘积翻转移位If k=1f1+1 disp/(暂停),pause;Else pause(0.5)Endend2.1.3 仿真图形图 2-2 卷积动态2.1.4仿真结果分析:1 1、卷积后的结果卷积后的结果0 0 0.90484 1.7236 2.4644 3.1347 3.7412 4.290 0 0.90484 1.7236 2.4644 3.1347 3.7412 4.294.7866 5.236 5.6425 6.0104 5.4384 4.9209 3.64
23、55 3.0967 2.60014.7866 5.236 5.6425 6.0104 5.4384 4.9209 3.6455 3.0967 2.6001 2.1508 1.7442 1.3763 1.0435 0/74226 0.46973 0.22313 2.1508 1.7442 1.3763 1.0435 0/74226 0.46973 0.223132 2、u u（n n）的长度为）的长度为 1010，h h（n n）的长度为）的长度为 1515，卷积后的长度，卷积后的长度 Ny=Nu+Nh-1=24.Ny=Nu+Nh-1=24.有图有图分析知，分析知，y y（n n）的长度为）的长
24、度为 2424 点，满足卷积定理。点，满足卷积定理。2.2 用双线性变换法设计 IIR 数字滤波器 2.2.1 双线性变换法的基本知识1）双线性变换法是将整个 s 平面映射到整个 z 平面，其映射关系为2 1 z1s T 1 z1双线性变换法克服了脉冲响应不变法从 s 平面到 z 平面的多值映射的缺点，消除了频谱混叠现象。但其在变换过程中产生了非线性的畸变，在设计 IIR 数字滤波器的过程中需要进行一定的预修正。2）设计中有关双线性变换法的子函数。Bilinear功能：双线性变换将 s 域(模拟域)映射到 z 域(数字域)的标准方法，将模拟滤波器变换成离散等效滤波器。调用格式：numd，den
25、dbilinear(num，den，Fs)；将模拟域传递函数变换为数字域传递函数，Fs 为取样频率。numd，dendbilinear(num，den，Fs，Fp)；将模拟域传递函数变换为数字域传递函数，Fs 为取样频率，Fp 为通带截止频率。zd，pd，kdbilinear(z，p，k，Fs)；将模拟域零极点增益系数变换到数字域，Fs 为取样频率。zd，pd，kdbilinear(z，p，k，Fs，Fp)；将模拟域零极点增益系数变换到数字域，Fs 为取样频率，Fp 为通带截止频率。Ad，Bd，Cd，Ddbilinear(A，B，C，D，Fs)；将模拟域状态变量系数变换到数字域，Fs 为取样频
26、率。3）用 MATLAB 双线性变换法进行 IIR 数字滤波器设计的步骤：(1)输入给定的数字滤波器设计指标；2(2)根据公式 tan()进行预修正，将数字滤波器指标转换成模拟滤波器T2设计指标；(3)确定模拟滤波器的最小阶数和截止频率；(4)计算模拟低通原型滤波器的系统传递函数；(5)利用模拟域频率变换法，求解实际模拟滤波器的系统传递函数；(6)用双线性变换法将模拟滤波器转换为数字滤波器。2.2.2 采用双线性变换法设计一个巴特沃斯数字低通滤波器，要求：wp0.2，Rp1 dB；ws0.3，As15 dB，滤波器采样频率 Fs100 Hz。2.2.3 自编matlab程序wp=0.2*pi;
27、%滤波器的通带截止频率ws=0.3*pi;%滤波器的阻带截止频率Rp=1;As=15;%输入滤波器的通阻带衰减指标ripple=10(-Rp/20);Attn=10(-As/20);Fs=100;T=1/Fs;Omgp=(2/T)*tan(wp/2);%原型通带频率预修正Omgs=(2/T)*tan(ws/2);%原型阻带频率预修正n,Omgc=buttord(Omgp,Omgs,Rp,As,s);%计算阶数n和3dB截止频率z0,p0,k0=buttap(n);%归一化原型设计ba=k0*real(poly(z0);%求原型滤波器系数baa=real(poly(p0);%求原型滤波器系数ab
28、a1,aa1=lp2lp(ba,aa,Omgc);%变换为模拟低通滤波器系数b,abd,ad=bilinear(ba1,aa1,Fs)%用双线性变换法求数字滤波器系数b,asos,g=tf2sos(bd,ad);%由直接型变换为级联型H,w=freqz(bd,ad);dbH=20*log10(abs(H)+eps)/max(abs(H);%化为分贝值subplot(2,2,1),plot(w/pi,abs(H);ylabel(|H|);title(幅度响应);axis(0,1,0,1.1);set(gca,XTickMode,manual,XTick,0,0.25,0.4,1);set(gca
29、,YTickMode,manual,YTick,0,Attn,ripple,1);gridsubplot(2,2,2),plot(w/pi,angle(H)/pi);ylabel(phi);title(相位响应);axis(0,1,-1,1);set(gca,XTickMode,manual,XTick,0,0.25,0.4,1);set(gca,YTickMode,manual,YTick,-1,0,1);gridsubplot(2,2,3),plot(w/pi,dbH);title(幅度响应(dB);ylabel(dB);xlabel(频率(pi);axis(0,1,-40,5);set(
30、gca,XTickMode,manual,XTick,0,0.25,0.4,1);set(gca,YTickMode,manual,YTick,-50,-15,-1,0);gridsubplot(2,2,4),zplane(bd,ad);2.2.4 仿真波形2.2.52.2.5仿真结果分析仿真结果分析1 1、滤波器阶数、滤波器阶数n=6,3dBn=6,3dB截止频率截止频率Omgc=76.623.Omgc=76.623.分子系数分子系数bdbd=0.0007 0.0044 0.0111 0.0148 0.0111 0.0044=0.0007 0.0044 0.0111 0.0148 0.011
31、1 0.0044 0.0007 0.0007分母系数分母系数 ad 1.0000 -3.1836 4.6222 -3.7795 1.8136 -0.4800ad 1.0000 -3.1836 4.6222 -3.7795 1.8136 -0.4800 0.0544 0.0544二阶环节级联结构的分子分母系数矩阵二阶环节级联结构的分子分母系数矩阵sossosSos=1 2.0004 1 1 -0.90437 0.21552Sos=1 2.0004 1 1 -0.90437 0.21552 1 2.018 1.0184 1 -1.0106 0.35827 1 2.018 1.0184 1 -1.0
32、106 0.35827 1 1.9816 0.98194 1 -1.2686 0.70513 1 1.9816 0.98194 1 -1.2686 0.70513增益增益g=0.00073782g=0.000737822 2、由频率特性曲线可知，设计结果在通带截止频率处满足由频率特性曲线可知，设计结果在通带截止频率处满足 R Rp p1dB1dB，AsAs15dB15dB 的的设计指标要求，系统的极点全部在单位圆内，是一个稳定的系统，由设计指标要求，系统的极点全部在单位圆内，是一个稳定的系统，由 N=6N=6 可知，设可知，设计的巴特沃斯数字低通滤波器是一个计的巴特沃斯数字低通滤波器是一个6
33、6 阶的系统，原型阶的系统，原型H H（s s）在）在s=-s=-处有六个零处有六个零点，映射到点，映射到z=-1z=-1 处。处。3 3、由幅度相应图可看出，通带有富余量，这是由于用最小阻带频率指标计算的结由幅度相应图可看出，通带有富余量，这是由于用最小阻带频率指标计算的结果。果。4 4、计算过程：、计算过程：W Wp p=0.2=0.2 R Rp p=1dB Ws=0.3=1dB Ws=0.3 As=15dB As=15dBOmgp=(2/T)*tan(wp/2)=200tan0.1Omgp=(2/T)*tan(wp/2)=200tan0.1=64.9839=64.9839Omgs=(2/
34、T)*tan(ws/2)=200tan0.15Omgs=(2/T)*tan(ws/2)=200tan0.15=101.9051=101.9051spsp=Omgs/Omgp=1.568=Omgs/Omgp=1.568K Kspsp=(10as/10-1)/(10ap/10-1)=10.875=10.875N=lgkN=lgkspsp/lg/lgspsp=5.3056=5.3056N N 取取 6 6Omgc=OmgsOmgc=Omgs（10100.1as0.1as-1)-1)-1/2N-1/2N=76.6229=76.6229G Gp p=1/(p=1/(p2 2+0.5176p+1)(p+0
35、.5176p+1)(p2 2+1.4142p+1)(p+1.4142p+1)(p2 2+1.9319p+1)+1.9319p+1)P=S/Omgc s=2/T(1-ZP=S/Omgc s=2/T(1-Z)(1+Z)(1+Z)H(z)=H(z)=(0.0007+0.0007+0.0044z0.0044z-1-1+0.0111z0.0111z-2-2+0.0148z0.0148z-3-3+0.0111z0.0111z-4-4+0.0044z0.0044z-5-5+0.0007z+0.0007z-6-6)/)/(1.0000-3.1836z(1.0000-3.1836z-1-1+4.6222z+4.
36、6222z-2-2-3.7795z-3.7795z-3-3+1.8136z+1.8136z-4-4-0.4800z-0.4800z-5-5+0.0544z+0.0544z-6-6)本题所设计的滤波器是一个六阶滤波器，由于本题所设计的滤波器是一个六阶滤波器，由于 MATLABMATLAB 精度有限，表现的精度有限，表现的 sossos 有有误差，把它们调整到精确值后，滤波器的级连型系统函数应该为误差，把它们调整到精确值后，滤波器的级连型系统函数应该为H H（z z）=0.00073782=0.00073782（1+z1+z-1-1)-6-6/（1-0.904371-0.90437 z z-1-1
37、0.21552z0.21552z-2-2）(1-1.0106(1-1.0106 z z-1-1+0.35827+0.35827z z-2-2)(1-1.2686 z)(1-1.2686 z-1-1+0.70513 z+0.70513 z-2-2)-1-11 1四、总结通过本次课程设计，我学会了用 Simulink 建模的方法和调试过程，验证了低通抽样定理，实现 PCM 编码及 13 折线量化，通过 MATLAB 实现卷积过程的动态，实现用双线性变换法设计 IIR 数字滤波器。学会了用 simulink 建立模型来完成仿真过程，全新的接触了另外一种结合书本学习的方法，通过建立模型后对其参数进行设
38、置，得到我们所需的结果，在接受 MATLAB 后，另外一种编程完成仿真过程，通过对程序的不断调试和更改，最后完成了本次课设的任务。了解了仿真技术及通信仿仿真的基本概念，知道了仿真的基本步骤和应该注意的问题，对 simulink 中的模块的基本了解及其作用，然后对 MATLAB 工作环境的认知，本次课程设计需要掌握通信原理和数字信号处理的理论知识，通过老师的讲解及到图书馆查阅各种有关的资料，掌握了利用 simulink 来建模实现抽样定理及 13 折线编码，建模锻炼了我的动手和独立思考的能力。学会了利用 MATLAB进行双线性变换法设计 FIR 数字低通滤波器，在进行编程和调试的过程会遇到各种各
39、样的问题，通过自己的思考，向同学老师请教逐次解决。通过本次课程设计，加深理解和巩固通信原理、数字信号处理课上所学的有关基本概念、基本理论和基本方法，锻炼了我的分析问题和解决问题的能力，同时培养良好的独立工作习惯和科学素质，初步掌握了 MATALB 的使用为今后的学习打下良好的基础。五、参考文献1李贺冰等.Simulink 通信仿真教程.北京：国防工业出版社，2006.2刘舒帆等.数字信号处理实验（MATLAB 版）.西安：西安电子科技大学出版社，2008.3党宏社.信号及系统实验（MATLAB 版）.西安：西安电子科技大学出版社，2007.4郭文彬等.通信原理基于 Matlab 的计算机仿真.北京：北京邮电大学出版社，2006.课程设计成绩评定表课程设计成绩评定表项项目目比例比例得得分分成成绩绩平时成绩（百分制记分）平时成绩（百分制记分）30%业务考核成绩（百分制记分）业务考核成绩（百分制记分）评评定定总评成绩（百分制记分）总评成绩（百分制记分）70%100%评定等级评定等级指导教师（签名）指导教师（签名）：优优良良中中及格及格不及格不及格2020年年月月日日