MATLAB实现通信系统仿真实例

补充内容：模拟调制系统的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;

%%根据抽样时间间隔进行抽样

t=(0:ts:pi/2);

envelop=cos(2*pi*t);%%DSB 信号包络

y=cos(2*pi*t).*cos(6*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]);

注意观察：过零点“反相”现象。

2.信号的频谱和功率

)()(f F t f ???→←傅立叶变换

【例3】画出DSB 信号)10cos()2cos(t t y ππ=的频谱并计算已调信号的功率。

clf

%%计算抽样时间间隔

fh=1;%%调制信号带宽(Hz)

fs=20*fh;

ts=1/fs;

%%根据抽样时间间隔进行抽样

t=(0:ts:2*pi);%抽样时间间隔要足够小，要满足抽样定理。

y=cos(2*pi*t).*cos(10*pi*t);%已调信号

%计算已调信号的傅立叶变化

n=2^nextpow2(length(t));%傅立叶变换的序列长度

spectrum=fft(y,n);

df=fs/n;%频率分辨率

f=[0:df:df*(length(spectrum)-1)]-fs/2;%FFT 频率向量

%计算信号功率

p=(norm(y)^2)/length(y)

%画出已调信号频谱

plot(f,fftshift(abs(spectrum)),'b','LineWidth',3);

xlabel('frequency/Hz');title('Magnitude-spectrum of y');

p =

0.2520

3.AWGN 的产生及功率谱密度

已知AWGN 信道中噪声的比特信噪比为10dB ，信号功率0.25W ，求噪声功率并画出噪声的时域波形。 clf

%%计算抽样时间间隔

fh=1;%%调制信号带宽(Hz)

fs=40*fh;

ts=1/fs;

t=(0:ts:2*pi);%根据抽样间隔生成抽样时间

%生成AWGN 噪声序列

snr=10;%信噪比10dB

snr_lin=10^(snr/10);%分贝信噪比转换为线性信噪比

noise_power=p_DSB/snr_lin

noise_std=sqrt(noise_power);

noise=randn(1,length(y))*noise_std;

plot(t,noise,'LineWidth',2)

noise_power =

4.DSB 调制解调

【例4】画出基带信号)2cos()(t t m π=和DSB 信号)10cos(

)2cos(t t y ππ=的频谱并计算基带信号和已调信号的功功率。

clf

%%计算抽样时间间隔

fh=1;%%调制信号带宽(Hz)

fs=40*fh;

ts=1/fs;

%%根据抽样时间间隔进行抽样

t=(0:ts:pi);%抽样时间间隔要足够小，要满足抽样定理。

m=cos(2*pi*t);%基带信号

y=m.*cos(10*pi*t);%已调信号

%计算基带信号的傅立叶变换

n=2^nextpow2(length(t));%傅立叶变换的序列长度

spectrum_baseband=fft(m,n);

df=fs/n;%频率分辨率

f=[0:df:df*(length(spectrum_baseband)-1)]-fs/2;%FFT频率向量

%计算已调信号的频谱

n=2^nextpow2(length(t));%傅立叶变换的序列长度

spectrum_DSB=fft(y,n);

df=fs/n;%频率分辨率

f=[0:df:df*(length(spectrum_DSB)-1)]-fs/2;%FFT频率向量

%计算信号功率

p_baseband=(norm(m)^2)/length(m)

p_DSB=(norm(y)^2)/length(y)

%画出已调信号频谱

subplot(2,1,1);

plot(f,fftshift(abs(spectrum_baseband)),'b','LineWidth',3); xlabel('frequency/Hz');title('Magnitude-spectrum of baseband'); subplot(2,1,2);

plot(f,fftshift(abs(spectrum_DSB)),'b','LineWidth',3);

xlabel('frequency/Hz');title('Magnitude-spectrum of DSB');

p_baseband =

0.5145

p_DSB =

clf

snr=10;%信噪比10dB

snr_lin=10^(snr/10);%分贝信噪比转换为线性信噪比

noise_power=p_DSB/snr_lin;

noise_std=sqrt(noise_power);

noise=randn(1,length(y))*noise_std;

signal_in_channel=y+noise;

subplot(2,1,1);

plot(t,signal_in_channel,'b')

xlabel('t');ylabel('DSB signal plus noise');axis([0,pi,-1,1]); subplot(2,1,2);

plot(t,y,'r');

xlabel('t');ylabel('DSB signal');

5.AM 调制解调

【例4】画出基带信号)2cos()(t t m π=和AM 信号)10cos()1)2(cos(t t y ππ+=以及载波)10cos()(t t s π=的频谱并计算基带信号和已调信号的功率。并比较该AM 信号与例题3中的DSB 信号的时域波形。 clf

%%计算抽样时间间隔

fh=1;%%调制信号带宽(Hz)

fs=20*fh;

ts=1/fs;

%%根据抽样时间间隔进行抽样

t=(0:ts:pi);%抽样时间间隔要足够小，要满足抽样定理。

m=cos(2*pi*t);

c=cos(10*pi*t);

y=(1+m).*c;%已调信号

%计算基带信号的傅立叶变换

n=2^nextpow2(length(t));%傅立叶变换的序列长度

spectrum_baseband=fft(m,n);

df=fs/n;%频率分辨率

f=[0:df:df*(length(spectrum_baseband)-1)]-fs/2;%FFT 频率向量

%计算载波功率

n=2^nextpow2(length(t));%傅立叶变换的序列长度

spectrum_carrier=fft(c,n);

df=fs/n;%频率分辨率

f=[0:df:df*(length(spectrum_baseband)-1)]-fs/2;%FFT 频率向量

%计算已调信号的频谱

n=2^nextpow2(length(t));%傅立叶变换的序列长度

spectrum_AM=fft(y,n);

df=fs/n;%频率分辨率

f=[0:df:df*(length(spectrum_AM)-1)]-fs/2;%FFT 频率向量

%计算信号功率

p_baseband=(norm(m)^2)/length(m)

p_AM=(norm(y)^2)/length(y)

p_carrier=(norm(c)^2)/length(c)

%画出已调信号频谱

subplot(3,1,1);

plot(f,fftshift(abs(spectrum_baseband)),'b','LineWidth',3);

xlabel('frequency/Hz');title('Magnitude-spectrum of baseband');

subplot(3,1,2);

plot(f,fftshift(abs(spectrum_AM)),'b','LineWidth',3);

xlabel('frequency/Hz');title('Magnitude-spectrum of AM');

subplot(3,1,3);

plot(f,fftshift(abs(spectrum_carrier)),'b','LineWidth',3); xlabel('frequency/Hz');title('Magnitude-spectrum of carrier'); p_baseband =

0.5168

p_AM =

0.8297

p_carrier =

clf

DSB=m.*c;%DSB信号

envelop_DSB=m;%包络

envelop_AM=1+m;

subplot(2,1,1);

plot(t,DSB,'b')

hold on

plot(t,envelop_DSB,'r:');

plot(t,-envelop_DSB,'r:');

hold off

xlabel('t');

ylabel('DSB signal');

subplot(2,1,2);

plot(t,y,'b');

hold on;

plot(t,envelop_AM,'r:');

plot(t,-envelop_AM,'r:');

hold off

xlabel('t');

ylabel('AM signal');

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);

基于MATLAB的MIMO通信系统仿真(DOC)

目录 （一）基于MATLAB的MIMO通信系统仿真………………………… 一、基本原理……………………………………………………… 二、仿真…………………………………………………………… 三、仿真结果……………………………………………………… 四、仿真结果分析…………………………………………………（二）自选习题部分…………………………………………………（三）总结与体会……………………………………………………（四）参考文献…………………………………………………… 实训报告 （一）基于MATLAB的MIMO通信系统仿真 一、基本原理 二、仿真 三、仿真结果 四、仿真结果分析 OFDM技术通过将频率选择性多径衰落信道在频域内转换为平坦信道，减小了多径衰落的影响。OFDM技术如果要提高传输速率，则要增加带宽、发送功率、子载波数目，这对于频谱资源紧张的无线通信时不现实的。 MIMO能够在空间中产生独立并行信道同时传输多路数据流，即传输速率很高。这些增加的信道容量可以用来提高信息传输速率，也可以通过增加信息冗余来提高通信系统的传输可靠性。但是MIMO却不能够克服频率选择性深衰落。 所以OFDM和MIMO这一对互补的技术自然走到了一起，现在是3G，未来也是4G，以及新一代WLAN技术的核心。总之，是核心物理层技术之一。 1、MIMO系统理论：

核心思想：时间上空时信号处理同空间上分集结合。 时间上空时通过在发送端采用空时码实现： 空时分组、空时格码，分层空时码。 空间上分集通过增加空间上天线分布实现。此举可以把原来对用户来说是有害的无线电波多径传播转变为对用户有利。 2、MIMO 系统模型： 11h 12 h 21 h 22 h r n h 1r n h 21 R n h 2 R n h 1 n n R h 可以看到，MIMO 模型中有一个空时编码器，有多根天线，其系统模型和上述MIMO 系统理论一致。为什么说nt>nr ，因为一般来说，移动终端所支持的天线数目总是比基站端要少。 接收矢量为：y Hx n =+，即接收信号为信道衰落系数X 发射信号+接收端噪声 3、MIMO 系统容量分析： （附MIMO 系统容量分析程序） 香农公式的信道容量（即信息传送速率）为： 2log (1/)C B S N =+ 4、在MIMO 中计算信道容量分两种情况： 未知CSI 和已知CSI （CSI 即为信道状态信息），其公式推导较为复杂，推导结果为信道容量是信噪比与接收、发射天线的函数。 在推导已知CSI 中，常用的有waterfilling ，即著名的注水原理。但是，根据相关文献资料，通常情况下CSI 可以当做已知，因为发送，接收端会根据具体信道情况估算CSI 的相关参数。 在这里对注水原理做一个简单介绍：之所以成为注水原理是因为理想的注水原理是在噪声大的时候少分配功率，噪声小时多分配功率，最后噪声+功率=定值，这如果用图形来表示，则类似于给水池注水的时候，水池低的地方就多注水，也就是噪声小分配的功率就多，故称这种达到容量的功率分配方式叫做注水原理。通过给各个天线分配不同的发射功率，增加系统容量。核心思想就是上面所阐述的，信道条件好，则分配更多功率；信道条件差，则分配较少的功率。 在MIMO 的信道容量当中要注意几个问题：（下面说已知CSI 都是加入了估计CSI 的算法，并且采用了注水原理。） 1. 已知CSI 的情况下的信道容量要比发送端未知CSI 的情况下的信道容量高，这是 由于当发送端已知CSI 的时候，发送端可以优化发送信号的协方差矩阵。也就是

MATLAB通信建模实验仿真实验报告

实验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)

MATLAB通信系统仿真心得体会

MATLAB通信系统仿真心得体会 课程名称(中文) MATLAB通信系统仿真成绩姓名班级学号日期 学习MATLAB通信系统仿真心得体会 经过一学期的MATLAB通信系统仿真的学习，使我对通信原 理及仿真实践有了更深层次的理解。在学习过程当中，了解到了MATLAB的语言基础以及应用的界面环境，基本操作和语法，通信仿真工具箱的应用，simulink 仿真基础，信号系统分析等一系列的内容。我明白学好这门课程是非常的重要。 在学习当中，我首先明白了通信系统仿真的现实意义，系统模型是对实际系统的一种抽象，是对系统本质(或是系统的某种特性)的一种描述。模型可视为对真实世界中物体或过程的信息进行形式化的结果。模型具有与系统相似的特性，可以以各种形式给出我们所感兴趣的信息。知道了通信系统仿真的必要性，利用系统建模和软件仿真技术，我们几乎可以对所有的设计细节进行分层次的建模和评估。通过仿真技术和方法，我们可以有效地将数学分析模型和经验模型结合起来。利用系统仿真方法，可以迅速构建一个通信系统模型，提供一个便捷，高效和精确的评估平台。明白了MATLAB通信系统仿真课程重点就是系统仿真软件 Matlab / Simulink 在通信系统建模仿真和性能评估方面的应用原理，通信系统仿真的一般原理和方法。 MATLAB集成度高，使用方便，输入简捷，运算高效，内容丰富，并且很容易由用户自行扩展，与其它计算机语言相比， MATLAB有以下显著特点:1.MATLAB是一种解释性语言;2(变量的“多功能性”;3.运算符号的“多功能性”;4(人机界面适合科技人员;5(强大而简易的作图功能;6(智能化程度高;7(功能丰富，可扩展性强。在MATLAB的Communication Toolbox(通 信工具箱)中提供了许多仿真函数和模块，用于对通信系统进行仿真和分析。

MATLAB 2psk通信系统仿真报告

实验一 2PSK调制数字通信系统 一实验题目 设计一个采用2PSK调制的数字通信系统 设计系统整体框图及数学模型； 产生离散二进制信源，进行信道编码（汉明码），产生BPSK信号； 加入信道噪声（高斯白噪声）； BPSK信号相干解调，信道解码； 系统性能分析（信号波形、频谱，白噪声的波形、频谱，信道编解 二实验基本原理 数字信号的传输方式分为基带传输和带通传输，在实际应用中，大多数信道具有带通特性而不能直接传输基带信号。为了使数字信号在带通信道中传输，必须使用数字基带信号对载波进行调制，以使信号与信道的特性相匹配。这种用数字基带信号控制载波，把数字基带信号变换为数字带通信号的过程称为数字调制。 数字调制技术的两种方法：①利用模拟调制的方法去实现数字式调制，即把数字调制看成是模拟调制的一个特例，把数字基带信号当做模拟信号的特殊情况处理；②利用数字信号的离散取值特点通过开关键控载波，从而实现数字调制。这种方法通常称为键控法，比如对载波的相位进行键控，便可获得相移键控（PSK）基本的调制方式。 图1 相应的信号波形的示例 1 0 1 调制原理 数字调相：如果两个频率相同的载波同时开始振荡，这两个频率同时达到正最大值，同时达到零值，同时达到负最大值，它们应处于"同相"状态；如果其中一个开始得迟了一点，就可能不相同了。如果一个达到正最大值时，另一个达到负最大值，则称为"反相"。一般把信号振荡一次（一周）作为360度。如果一个波比另一个波相差半个周期，我们说两个波的

相位差180度，也就是反相。当传输数字信号时，"1"码控制发0度相位，"0"码控制发180度相位。载波的初始相位就有了移动，也就带上了信息。 相移键控是利用载波的相位变化来传递数字信息，而振幅和频率保持不变。在2PSK中，通常用初始相位0和π分别表示二进制“1”和“0”。因此，2PSK信号的时域表达式为(t)=Acos t+) 其中，表示第n个符号的绝对相位： = 因此，上式可以改写为 图2 2PSK信号波形 解调原理 2PSK信号的解调方法是相干解调法。由于PSK信号本身就是利用相位传递信息的，所以在接收端必须利用信号的相位信息来解调信号。下图2-3中给出了一种2PSK信号相干接收设备的原理框图。图中经过带通滤波的信号在相乘器中与本地载波相乘，然后用低通滤波器滤除高频分量，在进行抽样判决。判决器是按极性来判决的。即正抽样值判为1，负抽样值判为0. 2PSK信号相干解调各点时间波形如图 3 所示. 当恢复的相干载波产生180°倒相时,解调出的数字基带信号将与发送的数字基带信号正好是相反,解调器输出数字基带信号全部出错.

matlab与通信仿真实验指导书(上)

《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实现通信系统仿真实例

补充内容：模拟调制系统的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; %%根据抽样时间间隔进行抽样

MATLAB仿真实验报告

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/ 」- □歹

北邮dsp软件matlab仿真实验报告

题目: 数字信号处理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的通信系统仿真毕业论文

创新实践报告
报 告 题 目： 学 院 名 称： 姓 名：
基于 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的通信系统仿真

创新实践报告 报告题目: 基于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与通信仿真》实验指导书-(通信级)word版本

电子信息学院 《MATLAB与通信仿真》 实验指导书 刘紫燕编写 适用专业：通信工程 贵州大学

二O一四年二月

前言 本课程是通信工程专业的选修课程。课程内容包含MATLAB的基本操作、MATLAB程序设计、函数文件、MATLAB的图形和数据处理、SIMULNK的基本应用及其在通信工程中的应用等。 通过本课程的学习，掌握MATLAB软件使用和编程方法，验证和深化书本知识，从而加强基础知识，掌握基本技能，提高MATLAB软件的编程能力，并用MATLAB软件解决通信工程专业中的实际问题，同时，针对通信工程专业的特点，要求学生掌握使用MATLAB来研究和开发与本专业相关的系统的方法。本课程设置5个实验，均为设计性实验。建议实验学时为12学时。 实验一是MATLAB软件的基本操作；实验二是MATLAB程序设计；实验三是MATLAB的图形绘制；实验四是MATLAB的数据处理；实验五是MATLAB/Simulink在电路中的仿真设计。 每个实验2个小时左右，为了使学生更好的掌握实验内容，学生务必要做到以下几点： (1)实验前认真预习实验。明确实验目的，熟悉实验内容，理论分析实验结果， 编写相应的程序代码，并撰写出预习报告。 (2)实验过程中积极思考，深入分析命令、程序的执行过程和执行结果，对比理 论分析结果，分析评判实验结果，并把实验中出现的问题及解决方法记录下来。 (3)实验完成后，总结本次实验有哪些收获，还存在什么问题，撰写并提交最终 的实验报告。 本指导书实验项目和要求明确，学生容易着手实验并得出实验结果。 本实验指导书适用于通信工程专业。

目录 实验一MATLAB的基本操作 (1) 实验二MATLAB程序设计 (9) 实验三MATLAB的图形绘制 (11) 实验四MATLAB的数据处理 (13) 实验五MATLAB/Simulink在电路中的仿真设计 (15) 实验报告的基本内容及要求 (20) 贵州大学实验报告 (21)

基于matlab的通信系统仿真要点

创新实践报告 报告题目：基于matlab的通信系统仿真学院名称：信息工程学院 姓名： 班级学号： 指导老师： 二O一四年十月十五日

一、引言 现代社会发展要求通信系统功能越来越强，性能越来越高，构成越来越复杂；另一方面，要求通信系统技术研究和产品开发缩短周期，降低成本，提高水平。这样尖锐对立的两个方面的要求，只有通过使用强大的计算机辅助分析设计技术和工具才能实现。在这种迫切的需求之下，MATLAB应运而生。它使得通信系统仿真的设计和分析过程变得相对直观和便捷，由此也使得通信系统仿真技术得到了更快的发展。通信系统仿真贯穿着通信系统工程设计的全过程，对通信系统的发展起着举足轻重的作用。通信系统仿真具有广泛的适应性和极好的灵活性，有助于我们更好地研究通信系统性能。通信系统仿真的基本步骤如下图所示：

二、仿真分析与测试 （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 D0D2D1D3 + + M V1 V2 OUT 02 1V D D =⊕0123 2V D D D D =⊕⊕⊕

matlab通信仿真课程设计样本

《matlab通信仿真设计》课程设计指导书 11月

课程设计题目1: 调幅广播系统的仿真设计 模拟幅度调制是无线电最早期的远距离传输技术。在幅度调制中, 以声音信号控制高频率正弦信号的幅度, 并将幅度变化的高频率正弦信号放大后经过天线发射出去, 成为电磁波辐射。 波动的电信号要能够有效地从天线发送出去, 或者有效地从天线将信号接收回来, 需要天线的等效长度至少达到波长的1/4。声音转换为电信号后其波长约在15~1500km之间, 实际中不可能制造出这样长度和范围的天线进行有效信号收发。因此需要将声音这样的低频信号从低频率段搬移到较高频率段上去, 以便经过较短的天线发射出去。 人耳可闻的声音信号经过话筒转化为波动的电信号, 其频率范围为20~20KHz。大量实验发现, 人耳对语音的频率敏感区域约为300~3400Hz, 为了节约频率带宽资源, 国际标准中将电话通信的传输频带规定为300~3400Hz。调幅广播除了传输声音以外, 还要播送音乐节目, 这就需要更宽的频带。一般而言, 调幅广播的传输频率范围约为100~6000Hz。 任务一: 调幅广播系统的仿真。 采用接收滤波器Analog Filter Design模块, 在同一示波器上观察调幅信号在未加入噪声和加入噪声后经过滤波器后的波形。采用另外两个相同的接收滤波器模块, 分别对纯信号和纯噪声滤波, 利用统计模块计算输出信号功率和噪声功率, 继而计算输出信噪比, 用Disply显示结果。 实例1: 对中波调幅广播传输系统进行仿真, 模型参数指标如下。

1.基带信号: 音频, 最大幅度为1。基带测试信号频率在100~6000Hz 内可调。 2.载波: 给定幅度的正弦波, 为简单起见, 初相位设为0, 频率为550~1605Hz 内可调。 3.接收机选频放大滤波器带宽为12KHz, 中心频率为1000kHz 。 4.在信道中加入噪声。当调制度为0.3时, 设计接收机选频滤波器输出信噪比为20dB, 要求计算信道中应该加入噪声的方差, 并能够测量接收机选频滤波器实际输出信噪比。 仿真参数设计: 系统工作最高频率为调幅载波频率1605KHz, 设计仿真采样率为最高工作频率的10倍, 因此取仿真步长为 8max 1 6.2310(1-1)10step t s f -==? 相应的仿真带宽为仿真采样率的一半, 即 18025.7(1-2)2step W KHz t == 设基带测试正弦信号为m(t)=Acos2πFt, 载波为c(t)=cos2πf c t, 则调制度为m a 的调制输出信号s(t)为 ()(1cos 2)cos 2(1-3)a c s t m Ft f t ππ=+ 容易求出, s(t)的平均功率为 21(1-4)24a m P =+ 设信道无衰减, 其中加入的白噪声功率谱密度为N 0/2, 那么仿真带宽(-W, W)内噪声样值的方差为 2002(1-5)2N W N W σ=?=

MATLAB通信系统仿真实验报告

实验一、MATLAB的基本使用与数学运算 目的：学习MATLAB的基本操作，实现简单的数学运算程序。 内容： 1-1 要求在闭区间[0，2π]上产生具有10个等间距采样点的一维数组。试用两种不同的指令实现。 运行代码：x=[0:2*pi/9:2*pi] 运行结果： 1-2 用M文件建立大矩阵x x=[ 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3.1 3.2 3.3 3.4 3.5 3.6 3.7 3.8 3.9] 代码：x=[ 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3.1 3.2 3.3 3.4 3.5 3.6 3.7 3.8 3.9] m_mat 运行结果： 1-3已知A=[5，6;7，8]，B=[9，10;11,12]，试用MATLAB分别计算 A+B,A*B,A.*B,A^3，A.^3，A/B,A\B. 代码：A=[5 6;7 8] B=[9 10;11 12] x1=A+B X2=A-B X3=A*B X4=A.*B X5=A^3 X6=A.^3 X7=A/B X8=A\B 运行结果： 1-4任意建立矩阵A，然后找出在[10,20]区间的元素位置。 程序代码及运行结果： 代码：A=[12 52 22 14 17;11 10 24 03 0;55 23 15 86 5 ] c=A>=10&A<=20

Matlab与通信仿真课程设计报告

《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与系统仿真实验手册

MATLAB与控制系统仿真实验指导书 河北大学电子信息工程学院 2006年9月

目录 实验一MATLAB基本操作与基本运算 (1) 实验二M文件及数值积分仿真方法设计 (3) 实验三MATLAB 的图形绘制 (4) 实验四函数文件设计和控制系统模型的描述 (6) 实验五控制系统的分析与设计 (7) 实验六连续系统离散化仿真方法设计 (8) 实验七SIMULINK 仿真 (9) 实验八SIMULINK 应用进阶 (10) 附录MATLAB常用函数 (12)

实验一MATLAB基本操作与基本运算 一、实验目的及要求： 1．熟悉MATLAB6.5的开发环境； 2．掌握MATLAB6.5的一些常用命令； 3．掌握矩阵、变量、表达式的输入方法及各种基本运算。 二、实验内容： 1.熟悉MATLAB6的开发环境： ①MATLAB的各种窗口： 命令窗口、命令历史窗口、启动平台窗口、工作空间窗口、当前路径窗口。 图1 MA TLAB界面窗口 ②路径的设置： 建立自己的文件夹，加入到MA TLAB路径中，并保存。 设置当前路径，以方便文件管理。 ③改变命令窗口数据的显示格式 >> format short >> format long

然后键入特殊变量：pi （圆周率），比较显示结果。 2.掌握MATLAB 常用命令 >> who %列出工作空间中变量 >> whos %列出工作空间中变量，同时包括变量详细信息 >>save test %将工作空间中变量存储到test.mat 文件中 >>load test %从test.mat 文件中读取变量到工作空间中 >>clear %清除工作空间中变量 >>help 函数名 %对所选函数的功能、调用格式及相关函数给出说明 >>lookfor %查找具有某种功能的函数但却不知道该函数的准确名称 如： lookfor Lyapunov 可列出与Lyapunov 有关的所有函数。 3. 在MATLAB 的命令窗口计算: 1) )2sin(π 2) 5.4)4.05589(÷?+ 4. 设计M 文件计算: 已知 求出： 1) a+b a-b a+b*5 a-b+I （单位阵） 2) a*b a.*b a/b 3) a^2 a.^2 注意：点运算 . 的功能，比较结果。 5. 设计M 文件计算: x=0:0.1:10 当sum>1000时停止运算，并显示求和结果及计算次数。 三、实验报告要求： 1）体会1、2的用法； 2）对3、4、5写出程序及上机的结果。 ?? ??? ?????=987654321a ?? ?? ? ?????=300120101b ) 2(100 2 i i i x x sum -=∑=

通信原理MATLAB仿真

小学期报告 实习题目通信原理Matlab仿真专业通信与信息工程 班级 学号 学生姓名 实习成绩 指导教师 2010年

通信原理Matlab仿真 目录 一、实验目的------------------------------------------------------------------------------------------------2 二、实验题目------------------------------------------------------------------------------------------------2 三、正弦信号波形及频谱仿真------------------------------------------------------------------------2 （一）通信原理知识--------------------------------------------------------------------------------------2 （二）仿真原理及思路--------------------------------------------------------------------------------------2 （三）程序流程图------------------------------------------------------------------------------------------- 3 （四）仿真程序及运行结果------------------------------------------------------------------------------3 （五）实验结果分析---------------------------------------------------------------------------------------5 四、单极性归零波形及其功率谱密度仿真--------------------------------------------------------5 （一）通信原理知识--------------------------------------------------------------------------------------6 （二）仿真原理及思路------------------------------------------------------------------------------ -------6 （三）程序流程图-------------------------------------------------------------------------------------------6 （四）仿真程序及运行结果--------------------------------------------------------------------------------6 （五）实验结果分析-------------------------------------------------------------------------------- -------6 五、升余弦滚降波形的眼图及功率谱密度仿真-------------------------------------------------8 （一）通信原理知识--------------------------------------------------------------------------------------8 （二）仿真原理及思路------------------------------------------------------------------------------ -------9 （三）程序流程图------------------------------------------------------------------------------- -----------9 （四）仿真程序及运行结果------------------------------------------------------------------------------10 （五）实验结果分析---------------------------------------------------------------------------------------11 六、PCM编码及解码仿真-----------------------------------------------------------------------------12 （一）通信原理知识---------------------------------------------------------------------------------- ---12 （二）仿真原理及思路------------------------------------------------------------------------------ ------ 13 （三）程序流程图------------------------------------------------------------------------------- -----------14 （四）仿真程序及运行结果------------------------------------------------------------------------------15 （五）实验结果分析---------------------------------------------------------------------------------------18 七、实验心得---------------------------------------------------------------------------- -------------------18