实验一单边带调幅系统的建模仿真
班级：姓名：学号：
一、实验目的
1.了解单边带调幅系统的工作原理
2.掌握单边带调幅系统的Matlab和Simulink建模过程
二、实验内容
1、Matlab设计一个单边带发信机、带通信道和相应的接收机，参数要
求如下。
（1）输入话音信号为一个话音信号，采样率8000Hz。话音输入后首先
进行预滤波，预滤波器是一个频率范围在[300,3400]Hz的带通滤波器。
其目的是将话音频谱限制在3400Hz以下。单边带调制的载波频率设计
为10KHz，调制输出上边带。要求观测单边带调制前后的信号功率谱。
（2）信道是一个带限高斯噪声信道，其通带频率范围是[10000,
13500]Hz。要求能够根据信噪比SNR要求加入高斯噪声。
（3）接收机采用相干解调方式。为了模拟载波频率误差对解调话音音
质的影响，设本地载波频率为9.8KHz，与发信机载波频率相差200Hz。
解调滤波器设计为300Hz到3400Hz的带通滤波器。
程序框图：
设计思想：
程序分为三部分：
一：SSB调制模块
首先从计算机中读入音频信号，作为原信号，读入完成后，对源信号进行参数采集和与滤波处理。进行与滤波之后，对信号进行希尔伯特变换，将原来的信号和载波相乘，将希尔伯特变换后的信号和载波进行希尔伯特变换后的信号相乘之后两者想减，得到SSB调制后的信号。
二：信道加噪声模块
通过信道，通过设置信道的信噪比来加入相应的噪声
三：解调模块：
将SSB调制后的信号通过信道加入噪声以后得到新的信号，并将信号和本地载波相乘进行想干解调，得到输出信号，并通过语句输出到相应的目录下。
2、用Simulink方式设计一个单边带传输系统并通过声卡输出接收机解调的结果声音。系统参数参照实例5.9，系统仿真参数设置为50KH
系统设计：
总体仿真框图
信号输入：
clc;
clear;
[wavs,fs]=wavread('GDGvoice8000.wav');
t_end=1/fs*length(wavs);
t=(1/fs:1/fs:t_end)';
source=[t wavs];
通过workspace编写程序对模块进行信号输入，程序如上
系统参数设置：
首先将系统的仿真步进时长设置为1/50000，即仿真采样频率为50Khz，对信号进行预滤波，这里采用butter带通滤波器频率设置为【300 3400】hz。
由于希尔伯特变换需要输入信号离散，所以将信号滤波之后进入一个零阶保持器，将信号离散化，保持器采样时间与系统采样时间相同。同样的，由于本地载波也需要希尔伯特变换，所以对本地载波也要进行离散化变化参数设置同上。本地载波频率设置为（10KHz）初相位设置为pi/2（为了得到cos 信号），经过希尔伯特变化后进入一个实部和虚部的分离器，然后将二者实部虚部分别相乘之后再相加就可以得到SSB调制信号。因为输出时需要的是离散化的信号，所以要将乘法器的采样时间设置为1/50000.。得到调制信号之后就可以进行信道的加入噪声，将信号经过高斯白噪声信道（SNR可以自由设置，越大则对调制信号的影响越小）这里设置为50。然后对信号进行想干解调，为了能对比结果将想干解调的本地载波设置与源信号的本地载波有所差距，这里设置为9800Hz。最后通过一个与前一个相同的带通滤波器，最后将解调后的信号输出。
输出：
wavwrite(demod_out,50000,'linkSSbDemod_OUT.wav');
这里输出也同样利用编程将信号输出到指定目录下，同时保存文件名为linkSSbDemod_OUT.wav 音频文件
显示结果
（1）能观察音频信号和SSB加载后的信号
输入音频信号SSB加载后的信号（2）能观察音频信号频谱和SSB加载后的信号频谱
（3）解调结果放到.wav音频文件，改变信道信噪比听解调的结果
当高斯白信道噪声为50时，解调后声音会有少量的噪声，但是基本清楚，当信噪比增加时信号会越来越清楚，当信噪比减小时，噪声会盖
过源信号的声音，到达很小时就只能听到噪声。
实验总结：
经过本次实验，我了解了单边带调幅系统的工作原理和单边带系统的matlab和Simulink建模的过程，虽然现在还对这些不太懂，但还是学到了一些。
实验二数字通信系统的建模仿真
一、实验目的
1.了解数字通信系统的建模过程
2.了解数字通信系统的仿真过程，并掌握对建模的好坏进行分析
二、实验内容
设计框图：
框图解析：
实验中要求用仿真模块来产生二进制数据源，得到数据源后和本地伪随机码相乘来实现对源信号的扩频，完成之后对信号进行极性变换，然后通过BPsk调制经过信道加入噪声之后，再和本地的载波信号相乘实现对信号的解扩之后进行解调得到一个新的信号，并且可以和原始信号进行比较，计算调制解调过程中产生的误码。
实验结果：
1）建立一个直接序列扩频体制(默认M序列)
观察收发数据波形；并用频谱仪观察原信号、直接序列扩频后的信号频谱、加噪声的信号频谱、解扩后的信号频谱；测试误码率
系统设计
参数设置：
本地二进制产生码源频率设置为100hz，本地M序列产生器设置为2000hz，极性转换设置为2，经过相乘器扩频后到达调制信道，Bpsk每一个码源采样个数设置为400，因为此处的频率是2000.经过高斯白噪声信道加入噪声，SNR可以自由设置（此处设置为20）。和本地一个频率为300hz，采样时间为2000的单频正弦波相加后进入解调信道之后得到解调后的信号与原信号相比，计算误码率。实验中各个测量点的图谱如下所示：
图一扩频前图二扩频后
图三加入噪声
图四解扩后图五解调后
图六误码率
从上面刘附图可以看出，原始信号经过扩频后频率范围增加，再加入信道噪声和干扰后，在图像波形中出现脉冲干扰，因为这里的SNR设置的比较大（20），所以这里的误码率接近为零，当减小SNR的时候，信道噪声增加，误码率就会增加，除此之外，信道时延的大小对误码率的影响也很大。
2）对比以Gold序列、m序列（已做）以及随机整数发生器Random Integer Generator作为直接序列扩频码源的传输性能，观察波形、频谱(扩频后、加噪声后、解扩后、解调后）、误码率，比较分析传输性能
这里系统设计大致与1）中的系统相同，只是将本地载波序列的产生模块作调整，m序列已经完成，首先介绍GOLD序列
2.1）GOLD序列
Gold Sequence Generator用于产生gold的序列，其他参数与M序列时相同，只要将Gold模块的采样时间改为1/2000就好了。得到结果如下：
扩频前扩频后
加入噪声
解扩后
解调后误码率
2.2）随机数序列
Random Integer Generator模块产生随机数，其他参数也相同，随机数产生序列改为二进制码，采样时间设置为1/2000.运行后结果如下：
扩频前扩频后
加入噪声
解扩后解调
后
误码率
比较上面三个不同的伪随机模块产生方法很明显的看出，Gold序列扩频后的信号毛刺很多，而且表面不是很光滑，解扩后的效果也受到影响，幅度变化比其他两者都要明显，输出后信号的是真更高。随机数虽然表面有毛刺，但是轮廓上比较光滑，相对Gold序列效果还是相对比较好的。三者当中GOLD序列相对最好，扩
频解扩后得到的信号与原信号相差不大，效果最好。
3、建立并测试一个直接序列扩频体制的码分多址传输系统
观察两路CDMA码源的收发数据波形；观察误码率；观察频谱（各自扩频后、加噪声后、各自解扩后、各自解调后）。码分多址和单序列的区别在于，码分多址是在一个信道上面传递两个以上的信号，因此这里的程序设计框图的设计和之前的单序列扩频的相同，基本设计思想也相同，只是在之前的基础之上加了一路信号。系统的设计如下所示：
框图中的两个CDMA子系统就是单扩频序列中扩频和调制的部分，后面的模块和单扩频中的模块是一样的。五个频谱图的观测点分别是扩频后的两路信号，加了噪声的一路信号和解调后的两路信号，这里运行后对所有的结果进行对比。
系统参数设置
因为是基于单序列扩频的基础之上，所以各个模块的参数设置单序列扩频是一样的，有所区别的是，这里是在一个信道上同时传递两路信号。此外，为了区别对比两路信号，将两路信号的产生种子数设置为61和67，这样在显示接收到的信号时就会有不同的波形出现。
3.1）M序列作为伪随机序列
扩频后对比
解调后对比
加噪声的扩频信号
接收到的信号
4、对比以Gold序列、m序列以及随机整数发生器Random Integer Generator作为直接序列扩频码源的码分多址传输系统
M序列已经完成，GOLD序列和随机整数序列，总体的设计相同，只是将扩频序列稍作改动
1.gold序列
扩频后对比
解调后对比
加入噪声接收到的两路信号
2.伪随机整数序列
扩频后对比
解调后对比
加入噪声之后的信号接收到的两路信号
结果分析：
首先看扩频之后的两路信号，可以看出，gold序列作为随即发生器时，产生的频谱波形整体上的光滑度比随机整数和M序列都要差得多。其次，看解调之后，解调之后三个频谱图看着都很相近，没有很大的区别。加入信道噪声之后，可以看出Gold和随机整数的表面都有很明显的毛刺出现，说明噪声对二者的影响还是很大的。最后看结果图片，每个接收到的波形都有不同程度上的误码。
5）调频扩频传输系统仿真模型
系统模型框图
参数设置：
数据率为100bps，数据调制采用2FSK调制方式，频率间隔为100hz，调频频点为32个，调频频率间隔为50hz。所以讲二进制数据源采样时间设置为0.01s，然后经过Mfsk模块，参数设置为：调制码元为2，频率间隔100hz，每个符号采样点为40.pN序列模块采样时间设置为1/250s，并且设置按帧输出，每帧5个码片。调频的MFSK则设置为调制元数32，频率间隔为50hz，每个符号采样点数为80.扩频经过高斯白噪声信道并且加入150hz的单品正弦波作为干扰源。输出结果如
下图所示：
从左到右从上到下分别是1调频扩频前的信号2加入信道噪声的信号3调频解跳后的信号4跳频扩频后的信号
总结和心得体会：
经过本次实验，了解了通信系统的建模过程和仿真系统并学会了对其的分析，总的来说学的不深，但学会了一些基础。
课题一： 连续时间信号和系统时域分析及MATLAB实现 课题要求： 深入研究连续时间信号和系统时域分析的理论知识。利用MATLAB强大的图形处理功能、符号运算功能以及数值计算功能，实现连续时间信号和系统时域分析的仿真波形。 课题内容： 一、用MATLAB实现常用连续时间信号的时域波形（通过改变参数，分析其时域特性）。 1、单位阶跃信号， 2、单位冲激信号， 3、正弦信号， 4、实指数信号， 5、虚指数信号， 6、复指数信号。 二、用MATLAB实现信号的时域运算 1、相加， 2、相乘， 3、数乘， 4、微分， 5、积分 三、用MATLAB实现信号的时域变换（参数变化，分析波形变化） 1、反转， 2、使移（超时，延时）， 3、展缩， 4、倒相， 5、综合变化 四、用MATLAB实现信号简单的时域分解 1、信号的交直流分解， 2、信号的奇偶分解 五、用MATLAB实现连续时间系统的卷积积分的仿真波形 给出几个典型例子，对每个例子，要求画出对应波形。 六、用MATLAB实现连续时间系统的冲激响应、阶跃响应的仿真波形。 给出几个典型例子，四种调用格式。 七、利用MATLAB实现连续时间系统对正弦信号、实指数信号的零状态响应的仿真波形。 给出几个典型例子，要求可以改变激励的参数，分析波形的变化。 课题二： 离散时间信号和系统时域分析及MATLAB实现。 课题要求： 深入研究离散时间信号和系统时域分析的理论知识。利用MATLAB强大的图
形处理功能、符号运算功能以及数值计算功能，实现离散时间信号和系统时域分析的仿真波形。 课题内容： 一、用MATLAB绘制常用信号的时域波形（通过改变参数分析其时域特性） 1、单位序列， 2、单位阶跃序列， 3、正弦序列， 4、离散时间实指数序列， 5、离散时间虚指数序列， 6、离散时间复指数序列。 二、用MATLAB实现信号的时域运算 1、相加， 2、相乘， 3、数乘。 三、用MATLAB实现信号的时域变换（参数变化，分析波形的变化） 1、反转， 2、时移（超时，延时）， 3、展缩， 4、倒相。 四、用MATLAB实现离散时间系统卷积和仿真波形 给出几个典型例子，对每个例子要求画出e(k)，h(k)，e(i)，h(i)，h(-i)，Rzs(k)波形。 五、用MATLAB实现离散时间系统的单位响应，阶跃响应的仿真波形 给出几个典型例子，四中调用格式。 六、用MATLAB实现离散时间系统对实指数序列信号的零状态响应的仿真波形 给出几个典型例子，要求可以改变激励的参数，分析波形的变化。 课题三： 连续时间信号傅里叶级数分析及MATLAB实现。 课题要求： 深入研究连续时间信号傅里叶级数分析的理论知识，利用MATLAB强大的图形处理功能，符号运算功能以及数值计算功能，实现连续时间周期信号频域分析的仿真波形。 课题内容： 一、用MATLAB实现周期信号的傅里叶级数分解与综合 以周期矩形波信号为例，绘出包含不同谐波次数的合成波形，观察合成波形与原矩形 波形之间的关系及吉布斯现象。
Matlab通信原理仿真 学号： 2142402 姓名：圣斌
实验一Matlab 基本语法与信号系统分析 一、实验目的： 1、掌握MATLAB的基本绘图方法； 2、实现绘制复指数信号的时域波形。 二、实验设备与软件环境： 1、实验设备：计算机 2、软件环境：MATLAB R2009a 三、实验内容： 1、MATLAB为用户提供了结果可视化功能，只要在命令行窗口输入相应的命令，结果就会用图形直接表示出来。 MATLAB程序如下： x = -pi::pi; y1 = sin(x); y2 = cos(x); %准备绘图数据 figure(1); %打开图形窗口 subplot(2,1,1); %确定第一幅图绘图窗口 plot(x,y1); %以x，y1绘图 title('plot(x,y1)'); %为第一幅图取名为’plot(x,y1)’ grid on; %为第一幅图绘制网格线 subplot(2,1,2) %确定第二幅图绘图窗口 plot(x,y2); %以x，y2绘图 xlabel('time'),ylabel('y') %第二幅图横坐标为’time’,纵坐标为’y’运行结果如下图： 2、上例中的图形使用的是默认的颜色和线型，MATLAB中提供了多种颜色和线型，并且可以绘制出脉冲图、误差条形图等多种形式图： MATLAB程序如下： x=-pi:.1:pi; y1=sin (x); y2=cos (x); figure (1); %subplot (2,1,1); plot (x,y1); title ('plot (x,y1)'); grid on %subplot (2,1,2); plot (x,y2);
课程设计任务书 学生姓名：董航专业班级：电信1006班 指导教师：阙大顺，李景松工作单位：信息工程学院 课程设计名称：Matlab应用课程设计 课程设计题目：Matlab运算与应用设计5 初始条件： 1.Matlab6.5以上版本软件； 2.课程设计辅导资料：“Matlab语言基础及使用入门”、“Matlab及在电子信息课程中的应 用”、线性代数及相关书籍等； 3.先修课程：高等数学、线性代数、电路、Matlab应用实践及信号处理类相关课程等。 要求完成的主要任务：（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求） 1.课程设计内容：根据指导老师给定的7套题目，按规定选择其中1套完成； 2.本课程设计统一技术要求：研读辅导资料对应章节，对选定的设计题目进行理论分析， 针对具体设计部分的原理分析、建模、必要的推导和可行性分析，画出程序设计框图，编写程序代码（含注释），上机调试运行程序，记录实验结果（含计算结果和图表），并对实验结果进行分析和总结。具体设计要求包括： ①初步了解Matlab、熟悉Matlab界面、进行简单操作； ②MATLAB的数值计算：创建矩阵矩阵运算、多项式运算、线性方程组、数值统计； ③基本绘图函数：plot, plot3, mesh, surf等，要求掌握以上绘图函数的用法、简单图形 标注、简单颜色设定等； ④使用文本编辑器编辑m文件，函数调用； ⑤能进行简单的信号处理Matlab编程； ⑥按要求参加课程设计实验演示和答辩等。 3.课程设计说明书按学校“课程设计工作规范”中的“统一书写格式”撰写，具体包括： ①目录； ②与设计题目相关的理论分析、归纳和总结； ③与设计内容相关的原理分析、建模、推导、可行性分析； ④程序设计框图、程序代码（含注释）、程序运行结果和图表、实验结果分析和总结； ⑤课程设计的心得体会（至少500字）； ⑥参考文献（不少于5篇）； ⑦其它必要内容等。 时间安排：1.5周（分散进行） 参考文献： [1](美)穆尔,高会生,刘童娜,李聪聪．MA TLAB实用教程(第二版) . 电子工业出版社,2010. [2]王正林,刘明．精通MATLAB(升级版) .电子工业出版社,2011. [3]陈杰. MA TLAB宝典(第3版) . 电子工业出版社,2011. [4]刘保柱,苏彦华,张宏林. MATLAB 7.0从入门到精通(修订版) . 人民邮电出版社,2010. 指导教师签名：年月日 系主任（或责任教师）签名：年月日
实验1：上采样与内插 一、实验目的 1、了解上采样与内插的基本原理和方法。 2、掌握上采样与内插的matlab程序的设计方法。 二、实验原理 上采样提高采样频率。上采样使得周期降低M倍，即新采样周期Tu和原有采样周期Ts的关系是T u=T s/M,根据对应的连续信号x(t),上采样过程从原有采样值x(kT s)生成新采样值x(kT u)=x(kT s/M)。操作的结果是在每两个采样值之间放入M-1个零值样点。 更实用的内插器是线性内插器，线性内插器的脉冲响应定义如下： 上采样值x(kT u)=x(kT s/M)通过与线性内插器的脉冲响应的卷积来完成内插。 三、实验内容 仿真正弦波采样和内插，通过基本采样x（k）,用M=6产生上采样x u(k)，由M=6线性内插得到样点序列x i（k）。 四、实验程序 % File: c3_upsampex.m M = 6; % upsample factor h = c3_lininterp(M); % imp response of linear interpolator t = 0:10; % time vector tu = 0:60; % upsampled time vector x = sin(2*pi*t/10); % original samples xu = c3_upsamp(x,M); % upsampled sequence subplot(3,1,1) stem(t,x,'k.') ylabel('x') subplot(3,1,2) stem(tu,xu,'k.') ylabel('xu') xi = conv(h,xu); subplot(3,1,3) stem(xi,'k.') ylabel('xi') % End of script file. % File: c3_upsample.m function out=c3_upsamp(in,M)
信源函数 randerr 产生比特误差样本 randint 产生均匀分布的随机整数矩阵 randsrc 根据给定的数字表产生随机矩阵 wgn 产生高斯白噪声 信号分析函数 biterr 计算比特误差数和比特误差率 eyediagram 绘制眼图 scatterplot 绘制分布图 symerr 计算符号误差数和符号误差率 信源编码 compand mu律/A律压缩/扩张 dpcmdeco DPCM（差分脉冲编码调制）解码dpcmenco DPCM编码 dpcmopt 优化DPCM参数 lloyds Lloyd法则优化量化器参数 quantiz 给出量化后的级和输出值 误差控制编码 bchpoly 给出二进制BCH码的性能参数和产生多项式convenc 产生卷积码 cyclgen 产生循环码的奇偶校验阵和生成矩阵cyclpoly 产生循环码的生成多项式 decode 分组码解码器 encode 分组码编码器 gen2par 将奇偶校验阵和生成矩阵互相转换gfweight 计算线性分组码的最小距离 hammgen 产生汉明码的奇偶校验阵和生成矩阵rsdecof 对Reed-Solomon编码的ASCII文件解码rsencof 用Reed-Solomon码对ASCII文件编码rspoly 给出Reed-Solomon码的生成多项式
syndtable 产生伴随解码表 vitdec 用Viterbi法则解卷积码 （误差控制编码的低级函数） bchdeco BCH解码器 bchenco BCH编码器 rsdeco Reed-Solomon解码器 rsdecode 用指数形式进行Reed-Solomon解码 rsenco Reed-Solomon编码器 rsencode 用指数形式进行Reed-Solomon编码 调制与解调 ademod 模拟通带解调器 ademodce 模拟基带解调器 amod 模拟通带调制器 amodce 模拟基带调制器 apkconst 绘制圆形的复合ASK-PSK星座图 ddemod 数字通带解调器 ddemodce 数字基带解调器 demodmap 解调后的模拟信号星座图反映射到数字信号dmod 数字通带调制器 dmodce 数字基带调制器 modmap 把数字信号映射到模拟信号星座图（以供调制）qaskdeco 从方形的QASK星座图反映射到数字信号qaskenco 把数字信号映射到方形的QASK星座图 专用滤波器 hank2sys 把一个Hankel矩阵转换成一个线性系统模型hilbiir 设计一个希尔伯特变换IIR滤波器 rcosflt 升余弦滤波器 rcosine 设计一个升余弦滤波器 （专用滤波器的低级函数） rcosfir 设计一个升余弦FIR滤波器 rcosiir 设计一个升余弦IIR滤波器 信道函数
自控系统仿真软件课程设计报告 MATLAB 设计题目：牛顿摆球 姓名： 学号： 院系： 班级：1203 指导教师： 2014年12月20日
一．课程设计目的 1、熟悉课程设计的基本流程； 2、掌握MATLAB语法结构及调试方法； 3、熟悉MATLAB函数调用，熟练二维画图; 4、掌握MATLAB语言在控制方面的运用； 5、学会用MATLAB进行基本仿真； 6、掌握MATLAB编程技巧，提高编程水平。 二．系统分析 1.题目的描述： （1）牛顿摆球原理描述 五个质量相同的球体由吊绳固定，彼此紧密排列。当摆动最右侧的球并在回摆时碰撞紧密排列的另外四个球，最左边的球将被弹出，并仅有最左边的球被弹出。当然此过程也是可逆的，当摆动最左侧的球撞击其它球时，最右侧的球会被弹出。当最右侧的两个球同时摆动并撞击其他球时，最左侧的两个球会被弹出。同理相反方向同样可行，并适用于更多的球。 为了更接近现实，在这里我将考虑重力及空气阻力的影响，摆球将不会永无止境的运动下去，由于外界因素的影响，摆球运动一段时间后将回归静止状态。（2）通过MATLAB动画程序制作软件，实现下述过程 当运行程序时，把最右边的小球拉到一定的高度放下，让其碰撞其余四个小球，仅让最左边的小球被弹出，当最左边小球回摆碰撞其它球时，最右边小球又被弹出，如此循环。由于是非理想条件下，摆球的摆动幅度会随摆动次数的增加越来越小，直到静止。 时间停顿两秒，把右边两小球一起拉到一定高度放下，让其碰撞其余三个球，同样仅让左边两球被弹出，当球回摆再次碰撞时，最右边两球又被同时弹出，如此循环，因为外界因素的影响，最终五个球都会静止下来。 （3）整个实验看似简单，但要在MATLAB上完成这样一个动画过程，还是需要下点功夫，克服困难的。经过自己的努力，终于实现了整个过程，这也是一种不小的收获。 2.设计要求： （1）能够实现有阻尼摆动，即摆幅随摆动次数增加越来越小，直到静止。（2）能够让摆球弧线摆动。 三．系统设计 1.系统设计过程 （1）通过函数axis建立坐标系 （2）在坐标系范围内通过函数line画各个支架 （3）通过函数title添加标题“动量守恒实验”、函数text添加标注“牛顿摆球” （4）通过函数line画出五个球，并设定其初始位置，颜色,大小,线条的擦拭方式
《Matlab与通信仿真》实验指导书(上) 刘毓杨辉徐健和煦黄庆东吉利萍编著 通信与信息工程学院 2011-1
目录 第一章 MALTAB基础知识 (1) 1.1MATLAB基础知识 (1) 1.2MATLAB基本运算 (2) 1.3MATLAB程序设计 (7) 第二章 MATLAB计算结果可视化和确知信号分析 (13) 2.1计算结果可视化 (13) 2.2确知信号分析 (17) 第三章随机信号与数字基带仿真 (23) 3.1基本原理 (23) 3.2蒙特卡罗算法 (30) 第四章模拟调制MATLAB实现 (34) 4.1模拟调制 (34) 4.2信道加性高斯白噪声 (35) 4.3AM调制解调的MATLAB实现 (36) 第五章模拟信号的数字传输 (45) 5.1脉冲编码调制 (45) 5.2低通抽样定理 (45) 5.3均匀量化原理 (46) 5.4非均匀量化 (48) 第六章数字频带传输系统 (52) 6.1数字频带传输原理 (52) 6.2数字频带传输系统的MATLAB实现 (53) 第七章通信系统仿真综合实验 (67) 7.1基本原理 (67) 7.2实验内容 (67)
第一章 MALTAB基础知识 本章目标 ●了解MATLAB 程序设计语言的基本特点，熟悉MATLAB软件运行环境 ●掌握创建、保存、打开m文件及函数的方法 ●掌握变量等有关概念，具备初步的将一般数学问题转化为对应的计算机模型并进行处理的能力 1.1 MATLAB基础知识 1.1.1 MATLAB程序设计语言简介 MATLAB，Matrix Laboratory的缩写，是由MathWorks公司开发的一套用于科学工程计算的可视化高性能语言，具有强大的矩阵运算能力。与大家常用的Fortran和C等高级语言相比，MATLAB的语法规则更简单，更贴近人的思维方方式，被称为“草稿纸式的语言”。MATLAB软件主要由主包、仿真系统（simulink）和工具箱（toolbox）三大部分组成。 1.1.2 MATLAB界面及帮助 MATLAB基本界面如图1-1所示，命令窗口包含标题栏、菜单栏、工具栏、命令行区、状态栏、垂直和水平波动条等区域。 图1-1 MATLAB基本界面 （1）菜单栏
补充内容：模拟调制系统的MATLAB 仿真 1.抽样定理 为了用实验的手段对连续信号分析，需要先对信号进行抽样（时间上的离散化），把连续数据转变为离散数据分析。抽样（时间离散化）是模拟信号数字化的第一步。 Nyquist 抽样定律：要无失真地恢复出抽样前的信号，要求抽样频率要大于等于两倍基带信号带宽。 抽样定理建立了模拟信号和离散信号之间的关系，在Matlab 中对模拟信号的实验仿真都是通过先抽样，转变成离散信号，然后用该离散信号近似替代原来的模拟信号进行分析的。 【例1】用图形表示DSB 调制波形)4cos()2cos(t t y ππ= 及其包络线。 clf %%计算抽样时间间隔 fh=1;%%调制信号带宽(Hz) fs=100*fh;%%一般选取的抽样频率要远大于基带信号频率，即抽样时间间隔要尽可能短。 ts=1/fs; %%根据抽样时间间隔进行抽样，并计算出信号和包络 t=(0:ts:pi/2)';%抽样时间间隔要足够小，要满足抽样定理。 envelop=cos(2*pi*t);%%DSB 信号包络 y=cos(2*pi*t).*cos(4*pi*t);%已调信号 %画出已调信号包络线 plot(t,envelop,'r:','LineWidth',3); hold on plot(t,-envelop,'r:','LineWidth',3); %画出已调信号波形 plot(t,y,'b','LineWidth',3); axis([0,pi/2,-1,1])% hold off% xlabel('t'); %写出图例 【例2】用图形表示DSB 调制波形)6cos()2cos(t t y ππ= 及其包络线。 clf %%计算抽样时间间隔 fh=1;%%调制信号带宽(Hz) fs=100*fh;%抽样时间间隔要足够小，要满足抽样定理。 ts=1/fs; %%根据抽样时间间隔进行抽样
第一类：单位转换 1.长度单位换算的设计与实现 2.面积单位换算的设计与实现 3.体积单位换算的设计与实现 4.容积单位换算的设计与实现 5.质量单位换算的设计与实现 6.时间单位换算的设计与实现 7.温度单位换算的设计与实现 7.压强单位换算的设计与实现 8.角度单位换算的设计与实现 8.功率单位换算的设计与实现 第二类：曲线绘制 1.直线的自动绘制和相关计算 2.椭圆的自动绘制和相关计算 3.双曲线的自动绘制和相关计算 4.抛物线的自动绘制和相关计算 5.心脏线的自动绘制和相关计算 6.渐开线的自动绘制和相关计算 7.滚圆线的自动绘制和相关计算 8.三叶玫瑰线的自动绘制和相关计算9.四叶玫瑰线的自动绘制和相关计 10.阿基米德螺线的自动绘制和相关计算第三类：曲面绘制 1.球面的自动绘制和相关计算 2.椭球面的自动绘制和相关计算 3.单叶双曲面的自动绘制和相关计算 4.双叶双曲面的自动绘制和相关计算 5.抛物面的自动绘制和相关计算 6.双曲抛物面的自动绘制和相关计算 7.双曲柱面的自动绘制和相关计算 8.椭圆柱面的自动绘制和相关计算 9.抛物柱面的自动绘制和相关计算 10.圆锥面的自动绘制和相关计算 第四类：线性回归 1.男士身高体重相关计算经验公式 2.女士身高体重相关计算经验公式 3.男士胖瘦等级的确定 4.女士胖瘦等级的确定 5.男士身高脚长相关计算经验公式 6.女士身高脚长相关计算经验公式 7.父子身高相关性研究 8.母子身高相关性研究 9.父女身高相关性研究 10.母女身高相关性研究 第五类：学习成绩 1.期末总评自动计算的设计与实现 2.成绩等级自动评定的设计与实现 3.成绩分段自动统计的设计与实现 4.成绩分布折线自动绘制的设计与实现 5.成绩自动统计分析的设计与实现 6.试卷分布自动分析的设计与实现 7.试卷难度自动分析的设计与实现 8.考试成绩名次自动生成的设计与实现
Matalab课后作业 学院：电气信息工程及其自动化 班级： 学号： 姓名： 完成日期： 2012年12月23日
1、 matlab 软件主要功能是什么？电气工程及其自动化专业本科生主要用到哪 些工具箱，各有什么功能？ 答：（1）主要功能：工业研究与开发； 数学教学，特别是线性代数；数值分析和科学计算方面的教学与研究；电子学、控制理论和物理学等工程和科学学科方面的教学与研究； 经济学、化学和生物学等计算问题的所有其他领域中的教学与研究；符号计算功能；优化工具；数据分析和可视化功能；“活”笔记本功能；工具箱；非线性动态系统建模和仿真功能。 （2）常用工具箱： （a ） MATLAB 主工具箱：扩充matlab 的数值计算、符号运算功能、图形建模仿真功能、文字处理功能以及与硬件实时交互功能。 （b ）符号数学工具箱：符号表达式、符号矩阵的创建；符号可变精度求解；因式分解、展开和简化；符号代数方程求解；符号微积分；符号微分方程。 （c ） SIMULINK 仿真工具箱： Simulink 是用于动态系统和嵌入式系统的多领域仿真和基于模型的设计工具。对各种时变系统，包括通讯、控制、信号处理、视频处理和图像处理系统，Simulink 提供了交互式图形化环境和可定制模块库来对其进行设计、仿真、执行和测试。 （d ）信号处理工具箱：数字和模拟滤波器设计、应用及仿真；谱分析和估计；FFT 、DCT 等 变换；参数化模型。 （e ）控制系统工具箱：连续系统设计和离散系统设计；状态空间和传递函数以及模型转换；时域响应（脉冲响应、阶跃响应、斜坡响应）；频域响应（Bode 图、Nyquist 图）；根轨迹、极点配置。 2、设y=23e t 4-sin(43t+3 ),要求以0.01秒为间隔，求出y 的151个点，并求出其导数的值和曲线。 程序如下： clc clear x=0:0.01:1.5; y=sqrt(3)/2*exp(-4*x).*sin(4*sqrt(3)*x+pi/3); y1=diff(y); subplot(2,1,1) plot(x,y) subplot(2,1,2) plot(x(1:150),y1) 曲线如下图所示：
《MATLAB与通信仿真》课程设计指导老师: 水英、汪泓 班级:07通信(1)班 学号:E07680104 :林哲妮
目录 目的和要求 (1) 实验环境 (1) 具体容及要求 (1) 实验容 题目一 (4) 题目容 流程图 程序代码 仿真框图 各个参数设置 结果运行 结果分析 题目二 (8) 题目容 流程图 程序代码 仿真框图 各个参数设置 结果运行 结果分析 题目三 (17) 题目容 流程图 程序代码 仿真框图 各个参数设置 结果运行 结果分析 题目四 (33) 题目容 流程图 程序代码 仿真框图 各个参数设置 结果运行 结果分析 心得与体会 (52)
目的和要求 通过课程设计，巩固本学期相关课程MATLAB与通信仿真所学知识的理解，增强动手能力和通信系统仿真的技能。在强调基本原理的同时，更突出设计过程的锻炼。强化学生的实践创新能力和独立进行科研工作的能力。 要求学生在熟练掌握MATLAB和simulink仿真使用的基础上，学会通信仿真系统的基本设计与调试。并结合通信原理的知识，对通信仿真系统进行性能分析。 实验环境 PC机、Matlab/Simulink 具体容及要求 基于MATLAB编程语言和SIMULINK通信模块库，研究如下问题： (1)研究BFSK在加性高斯白噪声信道下（无突发干扰）的误码率性能与信噪比之间的关系； (2)研究BFSK在加性高斯白噪声信道下（有突发干扰）的误码率性能与信噪比之间的关系； 分析突发干扰的持续时间对误码率性能的影响。 (3)研究BFSK+信道编码（取BCH码和汉明码）在加性高斯白噪声信道下（无突发干扰） 的误码率性能与信噪比之间的关系；分析不同码率对误码率性能的影响。比较不同信道编码方式的编码增益性能。 (4)研究BFSK+信道编码（取BCH码和汉明码）在加性高斯白噪声信道下（有突发干扰） 的误码率性能与信噪比之间的关系；分析突发干扰的持续时间对误码率性能的影响。分析不同码率对误码率性能的影响。比较不同信道编码方式的编码增益性能。
创新实践报告
报 告 题 目： 学 院 名 称： 姓 名：
基于 matlab 的通信系统仿真 信息工程学院 余盛泽
班 级 学 号： 指 导 老 师： 温 靖
二 O 一四年十月十五日
目录
一、引言........................................................................................................................ 3 二、仿真分析与测试 ................................................................................................... 4
2.1 随机信号的生成 ............................................................................................................... 4 2.2 信道编译码 ........................................................................................................................ 4 2.2.1 卷积码的原理 ........................................................................................................ 4 2.2.2 译码原理 ................................................................................................................ 5 2.3 调制与解调 ....................................................................................................................... 5 2.3.1 BPSK 的调制原理 .................................................................................................. 5 2.3.2 BPSK 解调原理 ...................................................................................................... 6 2.3.3 QPSK 调制与解调 ................................................................................................. 7 2.4 信道 .................................................................................................................................... 8
《MATLAB与通信仿真》课程设计指导老师: 张水英、汪泓 班级:07通信(1)班 学号:E07680104 姓名:林哲妮
目录 目的和要求 (1) 实验环境 (1) 具体内容及要求 (1) 实验内容 题目一 (4) 题目内容 流程图 程序代码 仿真框图 各个参数设置 结果运行 结果分析 题目二 (8) 题目内容 流程图 程序代码 仿真框图 各个参数设置 结果运行 结果分析 题目三 (17) 题目内容 流程图 程序代码 仿真框图 各个参数设置 结果运行 结果分析 题目四 (33) 题目内容 流程图 程序代码 仿真框图 各个参数设置 结果运行 结果分析 心得与体会 (52)
目的和要求 通过课程设计，巩固本学期相关课程MATLAB与通信仿真所学知识的理解，增强动手能力和通信系统仿真的技能。在强调基本原理的同时，更突出设计过程的锻炼。强化学生的实践创新能力和独立进行科研工作的能力。 要求学生在熟练掌握MATLAB和simulink仿真使用的基础上，学会通信仿真系统的基本设计与调试。并结合通信原理的知识，对通信仿真系统进行性能分析。 实验环境 PC机、Matlab/Simulink 具体内容及要求 基于MATLAB编程语言和SIMULINK通信模块库，研究如下问题： (1)研究BFSK在加性高斯白噪声信道下（无突发干扰）的误码率性能与信噪比之间的关系； (2)研究BFSK在加性高斯白噪声信道下（有突发干扰）的误码率性能与信噪比之间的关系； 分析突发干扰的持续时间对误码率性能的影响。 (3)研究BFSK+信道编码（取BCH码和汉明码）在加性高斯白噪声信道下（无突发干扰） 的误码率性能与信噪比之间的关系；分析不同码率对误码率性能的影响。比较不同信道编码方式的编码增益性能。 (4)研究BFSK+信道编码（取BCH码和汉明码）在加性高斯白噪声信道下（有突发干扰） 的误码率性能与信噪比之间的关系；分析突发干扰的持续时间对误码率性能的影响。分析不同码率对误码率性能的影响。比较不同信道编码方式的编码增益性能。
MATLA仿真实验报告 学院：计算机与信息学院 课程：—随机信号分析 姓名： 学号： 班级： 指导老师： 实验一
题目：编写一个产生均值为1,方差为4的高斯随机分布函数程序, 求最大值，最小值，均值和方差，并于理论值比较。 解：具体的文件如下，相应的绘图结果如下图所示 G仁random( 'Normal' ,0,4,1,1024); y=max(G1) x=mi n(G1) m=mea n(G1) d=var(G1) plot(G1);
实验二 题目：编写一个产生协方差函数为CC)=4e":的平稳高斯过程的程序，产生样本函数。估计所产生样本的时间自相关函数和功率谱密度，并求统计自相关函数和功率谱密度，最后将结果与理论值比较。 解：具体的文件如下，相应的绘图结果如下图所示。 N=10000; Ts=0.001; sigma=2; beta=2; a=exp(-beta*Ts); b=sigma*sqrt(1-a*a); w=normrnd(0,1,[1,N]); x=zeros(1,N); x(1)=sigma*w(1); for i=2:N x(i)=a*x(i-1)+b*w(i); end %polt(x); Rxx=xcorr(x0)/N; m=[-N+1:N-1]; Rxx0=(sigma A2)*exp(-beta*abs(m*Ts)); y=filter(b,a,x) plot(m*Ts,RxxO, 'b.' ,m*Ts,Rxx, 'r');
periodogram(y,[],N,1/Ts); 文件旧硯化）插入（1〕 ZMCD 克闻〔D ］窗口曲） Frequency (Hz) 50 100 150 200 250 300 350 400 450 500 NH---.HP)&UO 二 balj/ 」- □歹
% ch3example1A.m clear; f_p=2400; f_s=5000; R_p=3; R_s=25; % 设计要求指标 [n, fn]=buttord(f_p,f_s,R_p,R_s, 's'); % 计算阶数和截止频率 Wn=2*pi*fn; % 转换为角频率 [b,a]=butter(n, Wn, 's'); % 计算H(s) f=0:100:10000; % 计算频率点和频率范围 s=j*2*pi*f; % s=jw=j*2*pi*f H_s=polyval(b,s)./polyval(a,s); % 计算相应频率点处H(s)的值 figure(1); subplot(2,1,1); plot(f, 20*log10(abs(H_s))); % 幅频特性 axis([0 10000 -40 1]); xlabel('频率Hz');ylabel('幅度dB'); subplot(2,1,2); plot(f, angle(H_s)); % 相频特性 xlabel('频率Hz');ylabel('相角rad'); figure(2); freqs(b,a); % 也可用指令freqs直接画出H(s)的频率响应曲线。 % ch3example1B.m clear; f_p=2400; f_s=5000; R_p=3; R_s=25; % 设计要求指标 [n, fn]=ellipord(f_p,f_s,R_p,R_s,'s'); % 计算阶数和截止频率 Wn=2*pi*fn; % 转换为角频率 [b,a]=ellip(n,R_p,R_s,Wn,'s'); % 计算H(s) f=0:100:10000; % 计算频率点和频率范围 s=j*2*pi*f; % s=jw=j*2*pi*f H_s=polyval(b,s)./polyval(a,s); % 计算相应频率点处H(s)的值 figure(1); subplot(2,1,1); plot(f, 20*log10(abs(H_s))); % 幅频特性 axis([0 10000 -40 1]); xlabel('频率Hz');ylabel('幅度dB'); subplot(2,1,2); plot(f, angle(H_s)); % 相频特性 xlabel('频率Hz');ylabel('相角rad'); figure(2); freqs(b,a); % 也可用指令freqs直接画出H(s)的频率响应曲线。 % ch3example2A.m f_N=8000; % 采样率 f_p=2100; f_s=2500; R_p=3; R_s=25; % 设计要求指标 Ws=f_s/(f_N/2); Wp=f_p/(f_N/2); % 计算归一化频率 [n, Wn]=buttord(Wp,Ws,R_p,R_s); % 计算阶数和截止频率 [b,a]=butter(n, Wn); % 计算H(z) figure(1); freqz(b,a, 1000, 8000) % 作出H(z)的幅频相频图, freqz(b,a, 计算点数, 采样率)
实验一单边带调幅系统的建模仿真 班级：姓名：学号： 一、实验目的 1.了解单边带调幅系统的工作原理 2.掌握单边带调幅系统的Matlab和Simulink建模过程 二、实验内容 1、Matlab设计一个单边带发信机、带通信道和相应的接收机，参数要 求如下。 （1）输入话音信号为一个话音信号，采样率8000Hz。话音输入后首先 进行预滤波，预滤波器是一个频率范围在[300,3400]Hz的带通滤波器。 其目的是将话音频谱限制在3400Hz以下。单边带调制的载波频率设计 为10KHz，调制输出上边带。要求观测单边带调制前后的信号功率谱。 （2）信道是一个带限高斯噪声信道，其通带频率范围是[10000, 13500]Hz。要求能够根据信噪比SNR要求加入高斯噪声。 （3）接收机采用相干解调方式。为了模拟载波频率误差对解调话音音 质的影响，设本地载波频率为9.8KHz，与发信机载波频率相差200Hz。 解调滤波器设计为300Hz到3400Hz的带通滤波器。 程序框图：
设计思想： 程序分为三部分： 一：SSB调制模块 首先从计算机中读入音频信号，作为原信号，读入完成后，对源信号进行参数采集和与滤波处理。进行与滤波之后，对信号进行希尔伯特变换，将原来的信号和载波相乘，将希尔伯特变换后的信号和载波进行希尔伯特变换后的信号相乘之后两者想减，得到SSB调制后的信号。 二：信道加噪声模块 通过信道，通过设置信道的信噪比来加入相应的噪声 三：解调模块： 将SSB调制后的信号通过信道加入噪声以后得到新的信号，并将信号和本地载波相乘进行想干解调，得到输出信号，并通过语句输出到相应的目录下。 2、用Simulink方式设计一个单边带传输系统并通过声卡输出接收机解调的结果声音。系统参数参照实例5.9，系统仿真参数设置为50KH 系统设计：
题目: 数字信号处理MATLAB仿真实验 姓名 学院 专业 班级 学号 班内序号
实验一：数字信号的 FFT 分析 1、实验内容及要求 (1) 离散信号的频谱分析： 设信号 此信号的0.3pi 和 0.302pi 两根谱线相距很近，谱线 0.45pi 的幅度很小，请选择合适的序列长度 N 和窗函数，用 DFT 分析其频谱，要求得到清楚的三根谱线。 (2) DTMF 信号频谱分析 用计算机声卡采用一段通信系统中电话双音多频（DTMF ）拨号数字 0～9的数据，采用快速傅立叶变换（FFT ）分析这10个号码DTMF 拨号时的频谱。 2、实验目的 通过本次实验，应该掌握： (a) 用傅立叶变换进行信号分析时基本参数的选择。 (b) 经过离散时间傅立叶变换（DTFT ）和有限长度离散傅立叶变换（DFT ） 后信号频谱上的区别，前者 DTFT 时间域是离散信号，频率域还是连续的，而 DFT 在两个域中都是离散的。 (c) 离散傅立叶变换的基本原理、特性，以及经典的快速算法（基2时间抽选法），体会快速算法的效率。 (d) 获得一个高密度频谱和高分辨率频谱的概念和方法，建立频率分辨率和时间分辨率的概念，为将来进一步进行时频分析（例如小波）的学习和研究打下基础。 (e) 建立 DFT 从整体上可看成是由窄带相邻滤波器组成的滤波器组的概念，此概念的一个典型应用是数字音频压缩中的分析滤波器，例如 DVD AC3 和MPEG Audio 。 3、程序代码 (1) N=5000; n=1:1:N; x=0.001*cos(0.45*pi*n)+sin(0.3*pi*n)-cos(0.302*pi*n-pi/4); y=fft(x,N); magy=abs(y(1:1:N/2+1)); k=0:1:N/2; w=2*pi/N*k; stem(w/pi,magy) axis([0.25,0.5,0,50]) （2） column=[1209,1336,1477,1633]; line=[697,770,852,941]; fs=10000; N=1024; 00010450303024().*cos(.)sin(.)cos(.)x n n n n ππππ=+--
创新实践报告 报告题目: 基于matlab的通信系统仿真学院名称: 信息工程学院 姓名: 班级学号: 指导老师: 二O一四年十月十五日
一、引言 现代社会发展要求通信系统功能越来越强,性能越来越高,构成越来越复杂;另一方面,要求通信系统技术研究与产品开发缩短周期,降低成本,提高水平。这样尖锐对立的两个方面的要求,只有通过使用强大的计算机辅助分析设计技术与工具才能实现。在这种迫切的需求之下,MA TLAB应运而生。它使得通信系统仿真的设计与分析过程变得相对直观与便捷,由此也使得通信系统仿真技术得到了更快的发展。通信系统仿真贯穿着通信系统工程设计的全过程,对通信系统的发展起着举足轻重的作用。通信系统仿真具有广泛的适应性与极好的灵活性,有助于我们更好地研究通信系统性能。通信系统仿真的基本步骤如下图所示: 二、仿真分析与测试 (1)随机信号的生成 利用Matlab中自带的函数randsrc来产生0、1等概分布的随机信号。源代码如下所示: global N N=300; global p
p=0、5; source=randsrc(1,N,[1,0;p,1-p]); (2)信道编译码 1、卷积码的原理 卷积码(convolutional code)就是由伊利亚斯(p 、Elias)发明的一种非分组码。在前向纠错系统中,卷积码在实际应用中的性能优于分组码,并且运算较简单。 卷积码在编码时将k 比特的信息段编成n 个比特的码组,监督码元不仅与当前的k 比特信息段有关,而且还同前面m=(N-1)个信息段有关。 通常将N 称为编码约束长度,将nN 称为编码约束长度。一般来说,卷积码中k 与n 的值就是比较小的整数。将卷积码记作(n,k,N)。卷积码的编码流程如下所示。 可以瞧出:输出的数据位V1,V2与寄存器D0,D1,D2,D3之间的关系。根据模2加运算特点可以得知奇数个1模2运算后结果仍就是1,偶数个1模2运算后结果就是0。 2、译码原理 卷积码译码方法主要有两类:代数译码与概率译码。代数译码主要根据码本身的代数特性进行译码,而信道的统计特性并没有考虑在内。目前,代数译码的主要代表就是大数逻辑解码。该译码方法对于约束长度较短的卷积码有较好的效果,并且设备较简单。概率译码,又称最大似然译码,就是基于信道的统计特性与卷积 码的特点进行计算。在现代通信系统中,维特比译码就是目前使用最广泛的概率 译码方法。 02 1V D D =⊕01232V D D D D =⊕⊕⊕
课程设计任务书 学生姓名：专业班级：电信 指导教师：工作单位：信息工程学院 题目:MATLAB运算与应用设计2 初始条件： 1.MATLAB6.5以上版本软件； 2.课程设计辅导资料：“MATLAB语言基础及使用入门”、“MATLAB及在电子信息课程中的 应用”等； 3.先修课程：信号与系统、数字信号处理、MATLAB应用实践及信号处理类课程等。 要求完成的主要任务: （包括课程设计工作量及其技术要求，以及说 明书撰写等具体要求） （1）选择一本《MATLAB教程》，学习该教程的全部内容，包括使用方法、数组运算、矩阵运算、数学运算、程序设计、符号计算、图形绘制、GUI设计等内容； （2）对该套综合题的10道题，进行理论分析，针对具体设计部分的原理分析、建模、必要的推导和可行性分析，画出程序设计框图，编写程序代码（含注释），上机调试运行程序，记录实验结果（含计算结果和图表）。 （3）对实验结果进行分析和总结； （4）要求阅读相关参考文献不少于5篇； （5）根据课程设计有关规范，按时、独立完成课程设计说明书。 时间安排: (1) 布置课程设计任务,查阅资料,学习《MATLAB教程》十周； (2) 进行编程设计一周； (3) 完成课程设计报告书一周； 指导教师签名：年月日 系主任（或责任教师）签名：年月日
目录 1 MATLAB概述 (3) 1.1MATLAB简介 (4) 1.2MATLAB的功能 (4) 1.3MATLAB 的典型应用 (6) 2设计题目：MATLAB运算与应用设计套题二 (6) 3设计内容 (8) 3.1 题一 (8) 3.2 题二 (8) 3.3 题三 (9) 3.4 题四 (10) 3.5 题五 (15) 3.6 题六 (15) 3.7 题七 (15) 3.8 题八 (16) 3.9 题九 (17) 3.10题十 (18) 4 课程设计心得 (20) 5参考文献 (21) 6 本科生课程设计成绩评定表 (22)