练习题基于Matlab的模拟通信系统的仿真设计.doc
matlab模拟设计一个通信系统
1.2 主要步骤
本通信传输系统的 MATLAB 仿真包括以下内容：单频正弦波模拟信号经过抽 样（抽样频率为 Fs=20F） 、增量调制、二进制自然编码、AMI 码型变换基带传输， 循环编码 ASK 调制后发送到衰落信道，然后经过 ASK 解调，循环码译码、AMI 码 译码、数模转换后恢复出单频正弦波模拟信号。 信源：采用单频率正弦波； 数字化方式：增量调制（DM） ； 基带码：AMI 码； 信道码：循环码； 调制方式：ASK； 信道类型：衰落信道。
通信系统课群综合训练课程设计说明书
目录
摘要.............................................................................................................................................1 Abstract...................................................................................................................................... 2 1 设计任务..................................................................................................................................3 1.1 系统简述......................................................................................................................3 1.2 主要步骤......................................................................................................................3 1.3 系统框图.....................................................................................................................3 2 系统原理及实现......................................................................................................................5 2.1 增量调制......................................................................................................................5 2.2 基带传输 AMI 码.........................................................................................................6 2.3 信道传输码循环码......................................................................................................7 2.4 ASK 调制......................................................................................................................9 2.5 衰落信道....................................................................................................................10 3 系统仿真................................................................................................................................11 3.1 增量调制仿真............................................................................................................11 3.2 基带传输 AMI 码.......................................................................................................12 3.3 信道传输码循环码....................................................................................................12 3.4 ASK 调制....................................................................................................................14 3.5 衰落信道....................................................................................................................14 3.6 解码过程....................................................................................................................15 3.7 数模转化及失真分析................................................................................................16 4 小结与体会............................................................................................................................18 5 参考文献................................................................................................................................19 6 附录程序................................................................................................................................20
基于matlab的通信系统仿真
创新实践报告报告题目：基于matlab的通信系统仿真学院名称：信息工程学院姓名：班级学号：指导老师：二O一四年十月十五日一、引言现代社会发展要求通信系统功能越来越强，性能越来越高，构成越来越复杂；另一方面，要求通信系统技术研究和产品开发缩短周期，降低成本，提高水平。
这样尖锐对立的两个方面的要求，只有通过使用强大的计算机辅助分析设计技术和工具才能实现。
在这种迫切的需求之下，MATLAB应运而生。
它使得通信系统仿真的设计和分析过程变得相对直观和便捷，由此也使得通信系统仿真技术得到了更快的发展。
通信系统仿真贯穿着通信系统工程设计的全过程，对通信系统的发展起着举足轻重的作用。
通信系统仿真具有广泛的适应性和极好的灵活性，有助于我们更好地研究通信系统性能。
通信系统仿真的基本步骤如下图所示：二、仿真分析与测试（1）随机信号的生成利用Matlab 中自带的函数randsrc 来产生0、1等概分布的随机信号。
源代码如下所示：global NN=300;global pp=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、译码原理卷积码译码方法主要有两类：代数译码和概率译码。
通信原理基于matlab的计算机仿真
例1-1%周期信号（方波）的展开,fb_jinshi.mclose all;clear all;N=100; %取展开式的项数为2N＋1项T=1;fs=1/T;N_sample=128; %为了画出波形，设置每个周期的采样点数dt = T/N_sample;t=0:dt:10*T-dt;n=-N:N;Fn = sinc(n/2).*exp(-j*n*pi/2);Fn(N+1)=0;ft = zeros(1,length(t));for m=-N:Nft = ft + Fn(m+N+1)*exp(j*2*pi*m*fs*t);endplot(t,ft)例 1-2利用FFT计算信号的频谱并与信号的真实频谱的抽样比较。
脚本文件T2F.m定义了函数T2F，计算信号的傅立叶变换。
function [f,sf]= T2F(t,st)%This is a function using the FFT function to calculate a signal's Fourier %Translation%Input is the time and the signal vectors,the length of time must greater %than 2%Output is the frequency and the signal spectrumdt = t(2)-t(1);T=t(end);df = 1/T;N = length(st);f=-N/2*df:df:N/2*df-df;sf = fft(st);sf = T/N*fftshift(sf);脚本文件F2T.m定义了函数F2T，计算信号的反傅立叶变换。
function [t st]=F2T(f,sf)%This function calculate the time signal using ifft function for the input %signal's spectrumdf = f(2)-f(1);Fmx = ( f(end)-f(1) +df);dt = 1/Fmx;N = length(sf);T = dt*N;%t=-T/2:dt:T/2-dt;t = 0:dt:T-dt;sff = fftshift(sf);st = Fmx*ifft(sff);另写脚本文件fb_spec.m如下：%方波的傅氏变换, fb_spec.mclear all;close all;T=1;N_sample = 128;dt=T/N_sample;t=0:dt:T-dt;st=[ones(1,N_sample/2), -ones(1,N_sample/2)]; %方波一个周期subplot(211);plot(t,st);axis([0 1 -2 2]);xlabel('t'); ylabel('s(t)');subplot(212);[f sf]=T2F(t,st); %方波频谱plot(f,abs(sf)); hold on;axis([-10 10 0 1]);xlabel('f');ylabel('|S(f)|');%根据傅氏变换计算得到的信号频谱相应位置的抽样值sff= T^2*j*pi*f*0.5.*exp(-j*2*pi*f*T).*sinc(f*T*0.5).*sinc(f*T*0.5);plot(f,abs(sff),'r-')例1-3%信号的能量计算或功率计算,sig_pow.mclear all;close all;dt = 0.01;t = 0:dt:5;s1 = exp(-5*t).*cos(20*pi*t);s2 = cos(20*pi*t);E1 = sum(s1.*s1)*dt; %s1(t)的信号能量P2 = sum(s2.*s2)*dt/(length(t)*dt); %s2(t)的信号功率s[f1 s1f]= T2F(t,s1);[f2 s2f]= T2F(t,s2);df = f1(2)-f1(1);E1_f = sum(abs(s1f).^2)*df; %s1(t)的能量,用频域方式计算df = f2(2)-f2(1);T = t(end);P2_f = sum(abs(s2f).^2)*df/T; %s2(t)的功率，用频域方式计算figure(1)subplot(211)plot(t,s1);xlabel('t'); ylabel('s1(t)');subplot(212)plot(t,s2)xlabel('t'); ylabel('s2(t)');例1-4%方波的傅氏变换，sig_band.mclear all;close all;T=1;N_sample = 128;dt=1/N_sample;t=0:dt:T-dt;st=[ones(1,N_sample/2) -ones(1,N_sample/2)];df=0.1/T;Fx = 1/dt;f=-Fx:df:Fx-df;%根据傅氏变换计算得到的信号频谱sff= T^2*j*pi*f*0.5.*exp(-j*2*pi*f*T).*sinc(f*T*0.5).*sinc(f*T*0.5);plot(f,abs(sff),'r-')axis([-10 10 0 1]);hold on;sf_max = max(abs(sff));line([f(1) f(end)],[sf_max sf_max]);line([f(1) f(end)],[sf_max/sqrt(2) sf_max/sqrt(2)]); %交点处为信号功率下降3dB处Bw_eq = sum(abs(sff).^2)*df/T/sf_max.^2; %信号的等效带宽例 1-5%带通信号经过带通系统的等效基带表示，sig_bandpass.mclear all;close all;dt = 0.01;t = 0:dt:5;s1 = exp(-t).*cos(20*pi*t); %输入信号[f1 s1f]= T2F(t,s1); %输入信号的频谱s1_lowpass = hilbert(s1).*exp(-j*2*pi*10*t); %输入信号的等效基带信号[f2 s2f]=T2F(t,s1_lowpass); %输入等效基带信号的频谱h2f = zeros(1,length(s2f));[a b]=find( abs(s1f)==max(abs(s1f)) ); %找到带通信号的中心频率h2f( 201-25:201+25 )= 1;h2f( 301-25:301+25) = 1;h2f = h2f.*exp(-j*2*pi*f2); %加入线性相位，[t1 h1] = F2T(f2,h2f); %带通系统的冲激响应h1_lowpass = hilbert(h1).*exp(-j*2*pi*10*t1); %等效基带系统的冲激响应figure(1)subplot(521);plot(t,s1);xlabel('t'); ylabel('s1(t)'); title('带通信号');subplot(523);plot(f1,abs(s1f));xlabel('f'); ylabel('|S1(f)|'); title('带通信号幅度谱');subplot(522)plot(t,real(s1_lowpass));xlabel('t');ylabel('Re[s_l(t)]');title('等效基带信号的实部');subplot(524)plot(f2,abs(s2f));xlabel('f');ylabel('|S_l(f)|');title('等效基带信号的幅度谱');%画带通系统及其等效基带的图subplot(525)plot(f2,abs(h2f));xlabel('f');ylabel('|H(f)|');title('带通系统的传输响应幅度谱'); subplot(527)plot(t1,h1);xlabel('t');ylabel('h(t)');title('带通系统的冲激响应');subplot(526)[f3 hlf]=T2F(t1,h1_lowpass);plot(f3,abs(hlf));xlabel('f');ylabel('|H_l(f)|');title('带通系统的等效基带幅度谱');subplot(528)plot(t1,h1_lowpass);xlabel('t');ylabel('h_l(t)');title('带通系统的等效基带冲激响应');%画出带通信号经过带通系统的响应及等效基带信号经过等效基带系统的响应tt = 0:dt:t1(end)+t(end);yt = conv(s1,h1);subplot(529)plot(tt,yt);xlabel('t');ylabel('y(t)');title('带通信号与带通系统响应的卷积')ytl = conv(s1_lowpass,h1_lowpass).*exp(j*2*pi*10*tt);subplot(5,2,10)plot(tt,real(yt));xlabel('t');ylabel('y_l(t)cos(20*pi*t');title('等效基带与等效基带系统响应的卷积×中心频率载波')例 1-6%例：窄带高斯过程，文件 zdpw.mclear all; close all;N0=1; %双边功率谱密度fc=10; %中心频率B=1; %带宽dt=0.01;T=100;t=0:dt:T-dt;%产生功率为N0*B的高斯白噪声P = N0*B;st = sqrt(P)*randn(1,length(t));%将上述白噪声经过窄带带通系统，[f,sf] = T2F(t,st); %高斯信号频谱figure(1)plot(f,abs(sf)); %高斯信号的幅频特性[tt gt]=bpf(f,sf,fc-B/2,fc+B/2); %高斯信号经过带通系统glt = hilbert(real(gt)); %窄带信号的解析信号，调用hilbert函数得到解析信号glt = glt.*exp(-j*2*pi*fc*tt);[ff,glf]=T2F( tt, glt );figure(2)plot(ff,abs(glf));xlabel('频率(Hz)'); ylabel('窄带高斯过程样本的幅频特性')figure(3)subplot(411);plot(tt,real(gt));title('窄带高斯过程样本')subplot(412)plot(tt,real(glt).*cos(2*pi*fc*tt)-imag(glt).*sin(2*pi*fc*tt))title('由等效基带重构的窄带高斯过程样本')subplot(413)plot(tt,real(glt));title('窄带高斯过程样本的同相分量')subplot(414)plot(tt,imag(glt));xlabel('时间t(秒)'); title('窄带高斯过程样本的正交分量')%求窄带高斯信号功率;注：由于样本的功率近似等于随机过程的功率，因此可能出现一些偏差P_gt=sum(real(gt).^2)/T;P_glt_real = sum(real(glt).^2)/T;P_glt_imag = sum(imag(glt).^2)/T;%验证窄带高斯过程的同相分量、正交分量的正交性a = real(glt)*(imag(glt))'/T;用到的子函数function [t,st]=bpf(f,sf,B1,B2)%This function filter an input at frequency domain by an ideal bandpass filter %Inputs:% f: frequency samples% sf: input data spectrum samples% B1: bandpass's lower frequency% B2: bandpass's higher frequency%Outputs:% t: frequency samples% st: output data's time samplesdf = f(2)-f(1);T = 1/df;hf = zeros(1,length(f));bf = [floor( B1/df ): floor( B2/df )] ;bf1 = floor( length(f)/2 ) + bf;bf2 = floor( length(f)/2 ) - bf;hf(bf1)=1/sqrt(2*(B2-B1));hf(bf2)=1/sqrt(2*(B2-B1));yf=hf.*sf.*exp(-j*2*pi*f*0.1*T);[t,st]=F2T(f,yf);例 1-7%显示模拟调制的波形及解调方法DSB，文件mdsb.m%信源close all;clear all;dt = 0.001; %时间采样间隔fm=1; %信源最高频率fc=10; %载波中心频率T=5; %信号时长t = 0:dt:T;mt = sqrt(2)*cos(2*pi*fm*t); %信源%N0 = 0.01; %白噪单边功率谱密度%DSB modulations_dsb = mt.*cos(2*pi*fc*t);B=2*fm;%noise = noise_nb(fc,B,N0,t);%s_dsb=s_dsb+noise;figure(1)subplot(311)plot(t,s_dsb);hold on; %画出DSB信号波形plot(t,mt,'r--'); %标示mt的波形title('DSB调制信号');xlabel('t');%DSB demodulationrt = s_dsb.*cos(2*pi*fc*t);rt = rt-mean(rt);[f,rf] = T2F(t,rt);[t,rt] = lpf(f,rf,2*fm);subplot(312)plot(t,rt); hold on;plot(t,mt/2,'r--');title('相干解调后的信号波形与输入信号的比较');xlabel('t')subplot(313)[f,sf]=T2F(t,s_dsb);psf = (abs(sf).^2)/T;plot(f,psf);axis([-2*fc 2*fc 0 max(psf)]);title('DSB信号功率谱');xlabel('f');function [t st]=lpf(f,sf,B)%This function filter an input data using a lowpass filter %Inputs: f: frequency samples% sf: input data spectrum samples% B: lowpass's bandwidth with a rectangle lowpass%Outputs: t: time samples% st: output data's time samplesdf = f(2)-f(1);T = 1/df;hf = zeros(1,length(f));bf = [-floor( B/df ): floor( B/df )] + floor( length(f)/2 ); hf(bf)=1;yf=hf.*sf;[t,st]=F2T(f,yf);st = real(st);例1-8%显示模拟调制的波形及解调方法AM,文件mam.m%信源close all;clear all;dt = 0.001; %时间采样间隔fm=1; %信源最高频率fc=10; %载波中心频率T=5; %信号时长t = 0:dt:T;mt = sqrt(2)*cos(2*pi*fm*t); %信源%N0 = 0.01; %白噪单边功率谱密度%AM modulationA=2;s_am = (A+mt).*cos(2*pi*fc*t);B = 2*fm; %带通滤波器带宽%noise = noise_nb(fc,B,N0,t); %窄带高斯噪声产生%s_am = s_am + noise;figure(1)subplot(311)plot(t,s_am);hold on; %画出AM信号波形plot(t,A+mt,'r--'); %标示AM的包络title('AM调制信号及其包络');xlabel('t');%AM demodulationrt = s_am.*cos(2*pi*fc*t); %相干解调rt = rt-mean(rt);[f,rf] = T2F(t,rt);[t,rt] = lpf(f,rf,2*fm); %低通滤波subplot(312)plot(t,rt); hold on;plot(t,mt/2,'r--');title('相干解调后的信号波形与输入信号的比较'); xlabel('t')subplot(313)[f,sf]=T2F(t,s_am);psf = (abs(sf).^2)/T;plot(f,psf);axis([-2*fc 2*fc 0 max(psf)]);title('AM信号功率谱');xlabel('f');例 1-9%显示模拟调制的波形及解调方法SSB，文件mssb.m%信源close all;clear all;dt = 0.001; %时间采样间隔fm=1; %信源最高频率fc=10; %载波中心频率T=5; %信号时长t = 0:dt:T;mt = sqrt(2)*cos(2*pi*fm*t); %信源%N0 = 0.01; %白噪单边功率谱密度%SSB modulations_ssb = real( hilbert(mt).*exp(j*2*pi*fc*t) );B=fm;%noise = noise_nb(fc,B,N0,t);%s_ssb=s_ssb+noise;figure(1)subplot(311)plot(t,s_ssb);hold on; %画出SSB信号波形plot(t,mt,'r--'); %标示mt的波形title('SSB调制信号');xlabel('t');%SSB demodulationrt = s_ssb.*cos(2*pi*fc*t);rt = rt-mean(rt);[f,rf] = T2F(t,rt);[t,rt] = lpf(f,rf,2*fm);subplot(312)plot(t,rt); hold on;plot(t,mt/2,'r--');title('相干解调后的信号波形与输入信号的比较');xlabel('t')subplot(313)[f,sf]=T2F(t,s_ssb);psf = (abs(sf).^2)/T;plot(f,psf);axis([-2*fc 2*fc 0 max(psf)]);title('SSB信号功率谱');xlabel('f');例 2-0%显示模拟调制的波形及解调方法VSB，文件mvsb.m%信源close all;clear all;dt = 0.001; %时间采样间隔fm=5; %信源最高频率fc=20; %载波中心频率T=5; %信号时长t = 0:dt:T;mt = sqrt(2)*( cos(2*pi*fm*t)+sin(2*pi*0.5*fm*t) ); %信源%VSB modulations_vsb = mt.*cos(2*pi*fc*t);B=1.2*fm;[f,sf] = T2F(t,s_vsb);[t,s_vsb] = vsbpf(f,sf,0.2*fm,1.2*fm,fc);figure(1)subplot(311)plot(t,s_vsb);hold on; %画出VSB信号波形plot(t,mt,'r--'); %标示mt的波形title('VSB调制信号');xlabel('t');%VSB demodulation[f,rf] = T2F(t,rt);[t,rt] = lpf(f,rf,2*fm);subplot(312)plot(t,rt); hold on;plot(t,mt/2,'r--');title('相干解调后的信号波形与输入信号的比较');xlabel('t')subplot(313)[f,sf]=T2F(t,s_vsb);psf = (abs(sf).^2)/T;plot(f,psf);axis([-2*fc 2*fc 0 max(psf)]);title('VSB信号功率谱');xlabel('f');function [t,st]=vsbpf(f,sf,B1,B2,fc)%This function filter an input by an residual bandpass filter %Inputs: f: frequency samples% sf: input data spectrum samples% B1: residual bandwidth% B2: highest freq of the basedband signal%Outputs: t: frequency samples% st: output data's time samplesdf = f(2)-f(1);T = 1/df;hf = zeros(1,length(f));bf1 = [floor( (fc-B1)/df ): floor( (fc+B1)/df )] ;bf2 = [floor( (fc+B1)/df )+1: floor( (fc+B2)/df )];f1 = bf1 + floor( length(f)/2 ) ;f2 = bf2 + floor( length(f)/2 ) ;stepf = 1/length(f1);hf(f1)=0:stepf:1-stepf;hf(f2)=1;f3 = -bf1 + floor( length(f)/2 ) ;f4 = -bf2 + floor( length(f)/2) ;hf(f3)=0:stepf:(1-stepf);hf(f4)=1;yf=hf.*sf;[t,st]=F2T(f,yf);st = real(st);例 2-1%显示模拟调制的波形及解调方法AM、DSB、SSB,%信源close all;clear all;dt = 0.001;fm=1;fc=10;t = 0:dt:5;N0 = 0.1;%AM modulationA=2;s_am = (A+mt).*cos(2*pi*fc*t);B = 2*fm;noise = noise_nb(fc,B,N0,t); s_am = s_am + noise;figure(1)subplot(321)plot(t,s_am);hold on;plot(t,A+mt,'r--');%AM demodulationrt = s_am.*cos(2*pi*fc*t);rt = rt-mean(rt);[f,rf] = T2F(t,rt);。
基于MATLAB的通信系统的仿真报告参考例文
摘要Simulink是Mathworks公司推出的基于Matlab平台的著名仿真环境Simulin作为一种专业和功能强大且操作简单的仿真工具，目前已被越来越多的工程技术人员所青睐，它搭建积木式的建模仿真方式既简单又直观，而且已经在各个领域得到了广泛的应用。
本文主要是以simulink为基础平台，对2ASK、2FSK、2PSK信号的仿真。
文章第一章内容是对simulink的简单介绍和通信技术的目前发展和未来展望；第二章是对2ASK、2FSK和2PSK信号调制及解调原理的详细说明；第三章是本文的主体也是这个课题所要表现的主要内容，第三章是2ASK、2FSK和2PSK信号的仿真部分，调制和解调都是simulink建模的的方法，在解调部分各信号都是采用相干解调的方法，而且在解调的过程中都对整个系统的误码率在display模块中有所显示本文的主要目的是对simulink的熟悉和对数字通信理论的更加深化和理解。
关键词：2ASK、2FSK、2PSK，simulink,调制，相干解调目录第一章绪论 (31)1.1 MATLAB/Smulink的简介 (31)1.2 通信发展简史........................................ 错误!未定义书签。
1 1.3 通信技术的现状和发展趋势............................ 错误!未定义书签。
4 第二章 2ASK、2FSK、2PSK和2DPSK的基本原理和实现....... 错误!未定义书签。
7 2.1 2ASK的基本原理和调制解调实现...................... 错误!未定义书签。
8 2.2 2FSK的基本原理和调制解调实现..................... 错误!未定义书签。
11 2.3 2PSK的基本原理和调制解调实现.................... 错误!未定义书签。
14 2. 2DPSK的基本原理和调制解调实现.................... 错误!未定义书签。
基于Matlab 下的通信系统仿真
IX
.-.. .. ..- /--.. ./.... .
2.程序或仿真模型 2.1 设计思想（流程图）
上图中，假如连续时间信号是一个带限信号，其频率是 - m ~ m ，抽样脉冲为理想单位 冲激串，其数学表达式为、
由图可见，模拟信号 X(t)经抽样后，得到已知抽样信号 Xs(t);
X
.-.. .. ..- /--.. ./.... .
中包含的有效内容，也即信息(Information) 。消息是具体的、表面的，而信息是抽象的、 本质的，且消息中包含的信息的多少可以用信息量来度量。 通信系统就是传递信息所需
要的一切技术设备和传输媒质的总和，包括信息源、发送设备、信道、接收设备和信宿(受 信者) ，它的一般模型如图 3-1。
1 模拟通信系统模型和数字通信系统模型 1.1 模拟通信系统模型 在模拟通信系统中，信源（信息源，也称发终端）的作用是把
I
.-.. .. ..- /--.. ./.... .
上述的实际抽样过程中，很容易用简单的数学公式来描述，设连续时间信号用 X(t)表示，抽 样周期为 Ts ，抽样频率为 Ws ，则以抽样信号的数学表达式为
I
.-.. .. ..- /--.. ./.... .
Matlab 原理
设计项目
基于 Matlab 下的通信系统仿真（信号的抽样） 基于 Simulink 下的通信系统仿真（信号的抽样）
姓名：许美茹 学号：1467119128 学院：信工院通信一班
II
.-.. .. ..- /--.. ./.... .
号称为模拟信号。在用户线上传输模拟信号的通信方式称为“模拟通信”。 数字信号与模拟信号不同，它是一种离散的、脉冲有无的组合形式，是负载数字信息的信 号。电报信号就属于数字信号。现在最常见的数字信号是幅度取值只有两种(用 0 和 1 代表) 的波形，称为“二进制信号”。“数字通信”是指用数字信号作为载体来传输信息，或者用数字 信号对载波进行数字调制后再传输的通信方式。 数字通信与模拟通信相比具有明显的优 点：首先是抗干扰能力强。数字通信中的信息是包含在脉冲的有无之中的，只要噪声绝对 值不超过某一门限值，接收端便可判别脉冲的有无，以保证通信的可靠性。其次是远距离 传输仍能保证质量。因为数字通信是采用再生中继方式，能够消除噪音，再生的数字信号 和原来的数字信号一样，可继续传输下去，这样通信质量便不受距离的影响，可高质量地 进行远距离通信。此外，它还便于采用大规模集成电路，便于实现加密处理，便于实现通 信网的计算机管理等优点。 实现数字通信，必须使发送端发出的模拟信号变为数字信号， 这个过程称为“模数变换”。模拟信号数字化最基本的方法有三个过程，第一步是“抽样”，就 是对连续的模拟信号进行离散化处理，通常是以相等的时间间隔来抽取模拟信号的样值。 第二步是“量化”，将模拟信号样值变换到最接近的数字值。因抽样后的样值在时间上虽是离 散的，但在幅度上仍是连续的，量化过程就是把幅度上连续的抽样也变为离散的。第三步 是“编码”，就是把量化后的样值信号用一组二进制数字代码来表示，最终完成模拟信号的数 字化。数字信号送入数字网进行传输。接收端则是一个还原过程，把收到的数字信号变为 模拟信号，即“数据摸变换”，从而再现声音或图像。 如果发送端发出的信号本来就是数字 信号，则用不着进行模数变换过程，数字信号可直接进入数字网进行传输。 区别在于调
实验一 模拟通信的MATLAB仿真..
实验一 模拟通信的MATLAB 仿真姓名：左立刚 学号：031040522简要说明：实验报告注意包括AM ，DSB ，SSB ，VSB ，FM 五种调制与解调方式的实验原理，程序流程图，程序运行波形图，simulink 仿真模型及波形，心得体会，最后在附录中给出了m 语言的源程序代码。
一．实验原理1．幅度调制（AM ）幅度调制（AM ）是指用调制信号去控制高频载波的幅度，使其随调制信号呈线性变化的过程。
AM 信号的数学模型如图3-1所示。
图2-1 AM 信号的数学模型为了分析问题的方便，令δ=0，1.1 AM 信号的时域和频域表达式()t S AM=[A 0+m ()t ]cos t cω （2-1）()t S AM =A 0π[()()ωωωωδC C ++-]+()()[]ωωωωc c M M ++-21（2-2）AM 信号的带宽2=BAMfH（2-3）式中，fH为调制信号的最高频率。
2.1.3 AM 信号的功率P AM 与调制效率ηAMP AM=()2222t m A +=PP mc + （2-4）式中，P C=2A为不携带信息的载波功率；()22t m P m=为携带信息的边带功率。
()()t t m A m PP AMCAM222+==η（2-5） AM 调制的优点是可用包络检波法解调，不需要本地同步载波信号，设备简单。
AM 调制的最大缺点是调制效率低。
2.2、双边带调制（DSB ）如果将在AM 信号中载波抑制，只需在图3-1中将直流 A 0去掉，即可输出抑制载波双边带信号。
2.2.1 DSB 信号的时域和频域表达式()()t t m t cDSB S ωcos= （2-6）()()()[]ωωωωωC C DSBM M S ++-=21 （2-7） DSB 信号的带宽fB BHAM DSB2== （2-8）DSB 信号的功率及调制效率由于不再包含载波成分，因此，DSB 信号的功率就等于边带功率，是调制信号功率的一半，即()()t t m PS P CDSB DSB 2221=== （2-9） 显然，DSB 信号的调制效率为100%。
基于matlab的FM通信系统仿真设计与实现_课程设计报告
题目：基于matlab地FM通信系统仿真设计与实现学生姓名：杨丽君学生学号： 1008030317系别：电气信息工程学院专业：电子信息工程届别： 14届指导教师：马立宪电气信息工程学院制2013年5月基于matlab地FM通信系统仿真设计与实现学生：杨丽君指导教师：马立宪电气信息工程学院电子信息工程1课程设计地任务与要求1.1 课程设计地任务(1)熟悉MATLAB文件中M文件地使用方法，包括函数、原理和方法地应用.(2)加深对FM信号调制原理地理解.(3)画出基于MATLAB地FM通信系统仿真设计与实现设计地原理图.1.2 课程设计地要求(1) 学会MATLAB软件地安装.(2)在做完FM调制仿真之后，在今后遇到类似地问题，学会对所面对地问题进行系统地分析，并能从多个层面进行比较.(3) 熟练并且掌握对MA TLAB软件地使用，学会输入程序并且加以运行.1.3 课程设计地研究基础通信地目地是传输信号.通信系统地作用是将信息从信息源发送到一个或者多个目地地.模拟信号是时间和幅值上都连续地信号.调制是用原始信号即调制信号去控制高频载波信号地某一参数，是指随着原始信号幅度地变化而变化.而FM频率调制是高频载波信号地频率随着原始信号幅度变化而变化.解调是将已调制地信号恢复成原始信号即基带调制信号.以下是通信系统地一般模型：图1通信系统地一般模型（1)信息源信息源（简称信源）地作用是把各种消息转换成电信号.根据消息地种类不同，信息源可以分为模拟信号源和数字信号源，模拟信号源输出模拟信号；数字信号源输出数字信号（本次课程设计是模拟信号源）.（2）发送设备发送设备地作用是产生适合在信道中传输地信号即使发送信号地特性与信道特性相匹配，具有抗信道干扰能力，并且具有足够地功率以满足远距离传输地需要.因此，发送设备涵盖地内容很多，包含变换、放大、滤波、编码、调制地过程.（3）信道信道是一种物理媒介，用于将来自发送设备地信号发送到接收端.信道分为无线信道和有线信道.在有线信道中可以是明线、电缆、光纤.在无线信道中，信道可以是自由空间.信道地固有特性及引入地干扰与噪声直接关系到通信地质量.（4）噪声源噪声源是信道中及分布在系统中地其他各处噪声集中表示，噪声是随机地、形式是多样地，它地出现直接干扰信号地传输.（5）接收设备接收设备地功能是将信号放大及反变换（如译码、解调等）.目地是从受到减损地信号中正确恢复出原始信号，减少在传输过程中噪声与干扰所带来地影响.2 FM通信系统方案制定2.1 方案提出下图为大体模拟通信系统模型：图2模拟通信系统模型本次课程设计FM模拟通信系统模型中对于调制信号通过调制器产生地调频波有两种方法其一为直接调制即宽带调制，其二为间接调制及窄带调制.解调器对应地解调方法也有两种其一为相干解调另外一种为非相干解调.所以据此我提出了四种方案：方案一：模型中调制器中地调制方法为直接调制即宽带调制，解调器对应地解调方法为相干解调.方案二：模型中调制器中地调制方法为直接调制即宽带调制，解调器对应地解调方法为非相干解调.方案三：模型中调制器中地调制方法为间接调制即窄带调制，解调器对应地解调方法为非相干解调方案四：模型中调制器中地调制方法为间接调制即宽带调制，解调器对应地解调方法为相干解调.2.2 方案论证窄带调频地应用更广泛与宽带调频，我们对此也更为熟悉，技术也更为成熟.此外，它地最大频率偏移较小，占据带宽较窄、抗干扰性能更好等.所以本次课程设计选择窄带调频.一般情况下，相干解调法较适用于窄带调频.所以在以上选择地前提下，本次课程设计地解调方法选择相干解调针对以上地分析选择地最佳方案为方案四.3 FM通信系统方案设计3.1 FM通信系统模型设计图3 模拟通信系统模型设计3.2 FM 通信系统各部分地功能调制器: 使信号与信道相匹配, 有利于信号在信道中传输.发滤波器: 滤除调制器输出地无用信号.收滤波器: 滤除信号频带以外地噪声.一般设N(t)为高斯白噪声,则Ni(t)为窄带白噪声.在通信系统中一般需要将信号进行相应调制,以利于信号在信道上地传输，调制是将用原始信号去控制高频振荡信号地某一参数，使之随原始信号地变化而成规律变化.调制可分为线性调制和非线性调制.线性调制有AM 、DSB 等，非线性调制有FM 、PM 等，这里主要讨论FM 调制通信系统 3.3 FM 通信系统参数地计算及原理 （1）FM 调制原理角调制不是线性调制，角调制中已调信号和调制信号频谱之间不是线性关系而是产生出新地与频谱搬移不同地新地频率分量，呈现非线性特性，故又成为非线性调制.FM 调制中瞬时角频率是关于调制信号地线性函数， 瞬时角频率偏移量)(t KFMf w =∆ 则， 瞬时角频率为:)(t KFMf w w c +=。
MATLAB通信原理仿真五章作业
通信原理仿真作业第五章 模拟调制1. AM 、DSB 调制及解调用matlab 产生一个频率为1Hz ，功率为1的余弦信源()m t ，设载波频率10c Hz ω=，02m =,试画出：● AM 及DSB 调制信号的时域波形；● 采用相干解调后的AM 及DSB 信号波形；● AM 及DSB 已调信号的功率谱；● 调整载波频率及m0，观察分的AM 的过调与DSB 反相点现象。
● 在接收端带通后加上窄带高斯噪声，单边功率谱密度00.1n =，重新解调。
2. SSB 调制及解调用matlab 产生一个频率为1Hz ，功率为1的余弦信源，设载波频率10c Hz ω=，,试画出：● SSB 调制信号的时域波形；● 采用相干解调后的SSB 信号波形；● SSB 已调信号的功率谱；● 在接收端带通后加上窄带高斯噪声，单边功率谱密度00.1n =，重新解调。
3. FM 调制及解调设输入信号为()cos 2m t t π=，载波中心频率为10c z f H =，VCO 的压控振荡系数为5/z H V ，载波平均功率为1W 。
试画出：● 已调信号的时域波形；● 已调信号的振幅谱；● 用鉴频器解调该信号，并与输入信号比较。
加入相同功率的高斯白噪声信号，DSB信号解调较AM信号解调误差较大。
二代码clc;close all;clear all;fm=1;fc=10;T=5;t=linspace(-5,5,1024);%采样频率a=sqrt(2);mt=a*cos(2*pi*fm*t);%带通滤波器A=2;s_db=mt.*cos(2*pi*fc*t);s_am=(A+mt).*cos(2*pi*fc*t); figure(1);subplot(2,1,1);0.511.522.533.544.55-11t/smt)AM相干解调0.511.522.533.544.55-11t/sm(t)DSB相干解调plot(t,s_am); hold on;plot(t,A+mt,'r--'); title('AM调制信号');xlabel('t/s'); ylabel('幅度');subplot(2,1,2);plot(t,s_db);title('DSB调制信号');xlabel('t/s');ylabel('幅度');hold on;plot(t,mt,'r--');sp=s_am.*cos(2*pi*fc*t);sp2=s_db.*cos(2*pi*fc*t);fs1=0.3;fp1=0.5;fpu=1.8;fsu=2.0;Fs=50;Rp=1;Rs=42;ws1=2*pi*fs1/Fs;wp1=2*pi*fp1/Fs;wpu=2*pi*fpu/Fs;wsu=2 *pi*fsu/Fs;wd1=0.5*(wp1+ws1);wd2=0.5*(wpu+wsu);bt=min(abs(ws1-wp1),abs(wsu-wpu));n=ceil(6.6*pi/bt);wn=[wd1/pi,wd2/pi];hn=fir1(n-1,wn,'bandpass',hamming(n));y=conv(sp,hn);y2=conv(sp2,hn);y=2*y;y2=2*y2;figure(2);subplot(2,1,1);plot(t,y(1:1024));hold on;plot(t,mt,'r--');axis([0 5 -1.5 1.5]);xlabel('t/s'),ylabel('m(t)'),title('AM相干解调');subplot(2,1,2);plot(t,y2(1:1024));hold on;plot(t,mt,'r--');axis([0 5 -1.5 1.5]);xlabel('t/s'),ylabel('m(t)'),title('DSB相干解调、DSB调制及解调用matlab产生一个频率为1Hz，功率为1的余弦信源()mt，设载波频率10cHzω=，2m=,试画出：●AM及DSB调制信号的时域波形；●采用相干解调后的AM及信号波形；●AM及DSB已调信号的功率谱；●调整载波频率及m0，观察分的AM的过调与DSB反相点现象。
通信系统MATLAB仿真实验2010FALL
通信系统MATLAB 仿真实验实验一 模拟调制系统实验目的：● 熟悉matlab 仿真环境。
● 掌握基于matlab 的通信系统建模方法。
●验证AM 、DSB 、SSB 、FM 系统的调制解调方法。
实验内容：1. 建立AM 信号的发送端模型,分别用包络检波法和相干解调法解调AM 信号。
2. 搭建DSB 系统调制解调模型。
3. 搭建SSB 系统调制解调模型。
4. 搭建FM 系统调制解调模型。
1 AM 信号调制的调制解调：1.1 AM 信号调制的调制AM 调制器模型如图1：图1-1AM 调制器原理框图s AM (t)=［A 0+m(t)］cos2πf c t (1-1)设调制信号为单一频率的正弦波，频率f =10Hz ，载波频率为100Hz 。
matlab 代码如下：A0=2; %直流偏置电压c S Afm=10; % 信号10Hz fc=100; %载波频率100Hzfs=800; %取样频率800Hz,采样频率选取原则是满足系统中可能出现的信号最高频率的2倍以上。
t=0:1/fs:0.1; %仿真时间长度0.1s mt=sin(2*pi*fm*t);%基带信号,可选多音调制比如mt=sin(2*pi*fm*ts)+ cos(2*pi*3*fm*ts) s_am=(mt+A0).*sin(2*pi*fc*t); 观察调制信号波形plot(t,s_am)观察调制信号频谱periodogram(s_am,[],'twosided',[],1/ts)；%有关此函数的用法参考matlab 随机帮助文档 思考题：看到的谱线如何？如何观察到更清晰的谱线？1.2 信号在信道中传输仅考虑加性高斯白噪声信道。
信号在高斯白噪声信道中传输时会叠加高斯白噪声，这里介绍两种叠加高斯白噪声的方法。
利用randn 或wgn 函数产生与信号长度相同的白噪声序列函数randn (1,N)产生1行N 列互不相关的、均值为零、标准差为1的高斯随机数，N 等于输入信号的长度，也意味着采样频率与信号的采样频率相同，如果信号的采样频率为f s ,则噪声单边功率谱密度snnnf B N 2202σσ==,其中B n =f s /2为噪声带宽,这里需要区别所加噪声功率与解调器输入噪声功率.解调器输入噪声功率是经过带通滤波以后的噪声功率N i =N 0B,B 为信号带宽,B<Bn.噪声功率就是其方差，改变噪声功率可以通过给所有的采样值乘以一个系数（标准差来完成），例如：产生方差为1的高斯噪声：x= randn (1,N)； 产生方差为4的高斯噪声：x=2* randn (1,N)； %可以利用matlab 函数var(x)来验证x 的方差. 利用Matlab 提供的函数awgn awgn 函数的用法为: y = awgn(x,snr) ;x.为输入信号, snr 为信噪比,单位是 dB. 如果 x 是复数, awgn 加复噪声. 这种用法假定 x 的功率是 0 dBW.(1w)信号功率x （单位w ）与dBW 的关系是：dBW=10log 10(x/1w). y = awgn(x,snr,sigpower) ；x,snr 同上, sigpower 指定信号功率,单位是dBW. y = awgn(x,snr,'measured') ； 让Matlab 自动测量信号功率。
matlab通信仿真课程设计
《matlab通信仿真设计》课程设计指导书2009年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 ，要求计算信道中应该加入噪声的方差，并能够测量接收机选频滤波器实际输出信噪比。
基于MATLAB的模拟通信系统的仿真与实现
• 139•针对通信原理课程的教学特点和传统实验教学存在的问题，讨论了将Matlab软件引入到通信原理课程教学的必要性。
以模拟调制系统为例，利用Matlab的工具箱和Simulink界面对通信系统进行可视化教学，并给出了仿真结果。
实践证明，不仅在课堂教学中以更加直观的方式进行讲解，而且补充和完善传统实验的不足，提高学生学习积极性，教学效果得到较大提升。
随着5G通信的到来，通信技术在人们日常生活中是无处不在，现代通信技术取得了显著进展。
通信原理作为高校通信工程和电子信息等本科专业课程体系中重要的专业基础课，系统阐述了模拟和数字通信系统的基本概念、基本原理和基本分析方法，为学生学习后续课程储备专业素养（王海华，Matlab/Simulink仿真在“通信原理”教学中的应用研究：湖北理工学院学报，2015）。
然而这门课程理论内容丰富，系统模型抽象，数学公式多，推理过程繁琐，学生普遍感到枯燥难懂，抓不住重点，学习吃力，不能顺利学好本课程（基于Matlab_Simulink的通信原理虚拟仿真实验教学方法研究：现代电子技术，2015；邵玉斌，Matlab/Simulink通信系统建模与仿真实例分析：清华大学出版社，2008）。
为此，在教学过程中引入Matlab仿真技术，理论联系实践开展教学工作，通过simulink界面搭建系统模型，调整参数，观察通信系统性能，激发学生的学习积极性，提升教学质量，实现良好的教学模式。
1 Matlab软件介绍Matlab在工程数值运算和系统仿真方面具有强大的功能，主要包括数值分析、仿真建模、系统控制和优化等功能（牛磊，赵正平，郭博，Matlab仿真在通信原理教学中的应用：阜阳师范学院学报，2014）。
在Matlab的Communication Toolbox（通信工具箱）中提供了许多仿真函数和模块，用于对通信系统进行仿真和分析。
Simulink平台是Matlab中一种可视化仿真工具，提供了建立模型方框图的图形用户界面（GUI），可以将图形化的系统模块连接起来，从而建立直观、功能丰富的动态系统模型（黄琳，曹杉杉，熊旭辉.基于Matlab的通信原理实验课程设计：湖北师范大学学报，2017）。
基于matlab的通信系统仿真实验
程序clear all;close all;fs=4000000; %抽样频率k=20000; %设定数字基带信号的频率fc=200000; %正弦载波频率t=0:1/fs:9525/fs; %仿真时间范围be=[];xin=0; %误信率for mm=0:9;msg=randint(1,32,2);generator=[1 1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 1]; c=[1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];multip=conv(c,msg);[divid,remainder]=deconv(multip,generator);remainder=mod(remainder(end-15:end),2);s=[msg,remainder]; %产生生成码m=s(ceil(k*t+0.01)); %将基带生成时域信号figure(1)subplot(311)plot(t,m);axis([0 25e-4 -0.2 1.2]);grid on;title('数字基带信号');q=randint(1,48,2);%将生成的基带转换为差分码for i=1:48if (i==1)if (s(i)==0)q(i)=0;elseq(i)=1;endelseif (s(i)==q(i-1))q(i)=0;elseq(i)=1;endendn=q(ceil(k*t+0.01)); %将差分码生成时域信号subplot(312)plot(t,n);axis([0 25e-4 -0.2 1.2]);grid on;title('差分码')x=(n-0.5).*2;car=cos(2*pi*fc*t); %定义载波dpsk=x.*car; %信号的载波调制subplot(313);plot(t,car);axis([0 25e-5 -1.2 1.2]);title('余弦载波');figure(2)subplot(211);plot(t,dpsk);axis([0 10e-4 -1.2 1.2]);title('2DPSK信号');grid on;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%mss=[];ps=22;dpskn=awgn(dpsk,ps);%高斯白噪声subplot(212);plot(t,dpskn);axis([0 10e-4 -1.2 1.2]);title('加噪后信号');grid on;fBW=40e3; %带通滤波器f=(0:3e3:4e5);w=2*pi*f/fs;z=exp(w*1i);%j为虚数BW=2*pi*fBW/fs;a=.8547;p=(1i^2*a^2);gain=.135;Hz=gain*(z+1).*(z-1)./(z.^2-(p));Hz(Hz==0)=10^(8);a=[1 0 0.7305];b=[0.135 0 -0.135];dait=filter(b,a,dpskn);dait=dait.*10;figure(3)subplot(311);plot(t,dait);axis([0 25e-4 -1.5 1.5]);title('通过带通滤波后输出');grid on;cm=dait.*car; %2dpsk相干解调subplot(312);plot(t,cm);axis([0 25e-4 -1.5 1.5]);grid on;title('通过相乘器后输出');p=0.72; %低通滤波器gain1=0.14;Hz1=gain1*(z+1)./(z-(p));a1=[1 -0.72];b1=[0.14 0.14];dit=filter(b1,a1,cm);dit=dit-mean(dit);subplot(313);plot(t,dit);axis([0 25e-4 -0.5 0.5]);title('通过低通滤波器后输出');grid on;H=1; %抽样判决器L=0;Z=0;len=length(dit);Vs=(dit);for ii=1:lenif dit(ii)>=ZVs(ii)=H;elseVs(ii)=L;endendfigure(4)subplot(311)plot(t,Vs)title('解调后差分信号')axis([0 2.5e-4 -0.2 1.2])grid on;c=randint(1,48,2); %产生解调后的差分码元for f=0:46c(f+1)=fix(Vs(f*200+50)+0.2);%fix:向0靠拢取整函数将连续信号转为21位矩阵endd=randint(1,48,2); %定义差分译码后的码元for l=1:48 %得到差分译码后的码元if (l==1)if (c(1)==0)d(1)=0;elsed(1)=1;endelseif (c(l)==q(l-1))d(l)=0;elsed(l)=1;endend[divid, remainder]=deconv(d,generator); %crc译码remainder=mod(remainder(end-15:end),2);if isequal(remainder,[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]); fprintf('校验正确' );elsefprintf('校验错误');xin=xin+1;endgrid on;y=d(ceil(k*t+0.01));subplot(313);plot(t,y);axis([0 25e-4 -0.2 1.2]);title('输出信号'); %基带信号与解调后的信号对比subplot(312)plot(t,m);axis([0 25e-4 -0.2 1.2]);title('原始基带信号');%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%计算误码率num2=0;for n=1:(length(d))if(s(n)~=d(n))num2=num2+1;disp(('出错的码元是：'));disp(n);endendberror=num2/(length(d));disp(('误码率为：'));disp(berror);endwuxin=xin/10;disp(('误信率是：'));disp(wuxin);。
基于MATLABSimulink的模拟通信系统的仿真课程设计
陕西理工学院通信原理课程设计题目基于SIMULINK的通信系统仿真毕业设计（论文）原创性声明和使用授权说明原创性声明本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及取得的成果。
尽我所知，除文中特别加以标注和致谢的地方外，不包含其他人或组织已经发表或公布过的研究成果，也不包含我为获得及其它教育机构的学位或学历而使用过的材料。
对本研究提供过帮助和做出过贡献的个人或集体，均已在文中作了明确的说明并表示了谢意。
作者签名：日期：指导教师签名：日期：使用授权说明本人完全了解大学关于收集、保存、使用毕业设计（论文）的规定，即：按照学校要求提交毕业设计（论文）的印刷本和电子版本；学校有权保存毕业设计（论文）的印刷本和电子版，并提供目录检索与阅览服务；学校可以采用影印、缩印、数字化或其它复制手段保存论文；在不以赢利为目的前提下，学校可以公布论文的部分或全部内容。
作者签名：日期：学位论文原创性声明本人郑重声明：所呈交的论文是本人在导师的指导下独立进行研究所取得的研究成果。
除了文中特别加以标注引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写的成果作品。
对本文的研究做出重要贡献的个人和集体，均已在文中以明确方式标明。
本人完全意识到本声明的法律后果由本人承担。
作者签名：日期：年月日学位论文版权使用授权书本学位论文作者完全了解学校有关保留、使用学位论文的规定，同意学校保留并向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅和借阅。
本人授权大学可以将本学位论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。
涉密论文按学校规定处理。
作者签名：日期：年月日导师签名：日期：年月日II题目基于SIMULINK的模拟通信系统的仿真（线性调制）摘要在模拟通信系统中，由模拟信源产生的携带信息的消息经过传感器转换成电信号，模拟基带信号在经过调制将低通频谱搬移到载波频率上适应信道，最终解调还原成电信号；在数字传输系统中，数字信号对高频载波进行调制，变为频带信号，通过信道传输，在接收端解调后恢复成数字信号。
基于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加运算特点可以得知奇数个1模2运算后结果仍是1，偶数个1模2运算后结果是0。
021V D D =⊕01232V D D D D=⊕⊕⊕2、译码原理卷积码译码方法主要有两类：代数译码和概率译码。
MATLAB仿真实例(通信原理)
值。仿真时我们用这个样值集合来表示信号 s( t ) 。显然 t 反映了仿真系统对信号波形的分 辨率， t 越小则仿真的精确度越高。据通信原理所学，信号被取样以后的频谱是频率的周 1 1 期函数，其重复周期是 。如果信号的最高频率为 f H ，那么必须有 f H 才能保证不 t 2 t 1 发生混叠失真。我们称 Bs 为仿真系统的系统带宽。如果我们的仿真程序中设定的采 2 t 样间隔是 t ，那么我们不能用此仿真程序来研究带宽大于 Bs 的信号或系统。 此外，信号 s( t ) 的频谱 S f 通常来说也是定义在频率区间 , 上的连续函数，所以仿 真频域特性时，我们也必须把 S f 截短并取样。考虑到系统带宽是 Bs ，我们把频域的截短
三、 实验内容
3.1 傅里叶变换与傅里叶反变换 对于确定信号 f (t ) ，其傅里叶变换为：
F( f )
F ( f ) 傅里叶反变换为：
f (t )e 2 j ft dt
f (t ) F ( f )e2 j ft df
在通信原理仿真中，傅里叶变换与傅里叶反变换会经常用到，我们可以利用 MATLAB 的快速傅里叶变换函数 fft 与快速傅里叶反变换函数 ifft 编写傅里叶变换子程序与傅里叶反 变换子程序。其 程 序 代 码 如 下 ：
傅里叶变换子程序：
%傅里叶变换子程序 function X=t2f(x) global dt df N t f T %X=t2f(x) %x 为时域的取样值矢量 %X 为 x 的傅氏变换 %X 与 x 长度相同,并为 2 的整幂。
基于Matlab的模拟通信系统的仿真设计.pptx
学海无 涯
应用。 3）.掌握参数设置方法和性能分析方法。 4）.通过实验中波形的变换，学会分析实验现象。 5）.学习掌握如何生成 GUI 界面以及如何使用它来显示相应的波形。
第二章 通信系统的调制与解调
1、通信系统的概念
通信是为了传输信息，通信系统就是将信息从信源发送到一个或多个目的 地，对于电通信来说，首先要把消息转变成电信号，然后经过发送设备，将信号 送入信道，在接收端利用接受设备对接收信号作相应的处理后，送给信宿再转换 为原来的消息，这一过程可利用图 1 所示的通信系统一般模型来概括。模拟信息源调制器噪信声
道
源
3
学海无 涯
受信者
解调器
图 2 模拟通信系统模型型
2.调制和解调的概念
调制：把信号转换成适合在信道中传输的形式的一种过程，广义的调制分为基带 调制和带通调制（也称为载波调制）。在无线通信中和其他大多数场合，调制均 指载波调制。
载波调制，就是用调制信号去控制载波的参数的过程，使载波的某一个或 某几个参数按照调制信号的规律而变化。调制信号是指来自信源的消息信号（ 基 带信号），这些信号可以是模拟的，也可以是数字的。未受调制的周期性振荡 信 号称为载波，它可以是正弦波，也可以是非正弦波（如周期性脉冲序列）。载 波 调制后称为已调信号，它含有调制信号的全部特征。
本次课程设计首先介绍了通信系统的概念，进而引出调制和解调，然 后介绍了我们常用的几种调制和解调的方法。由于 MATLAB 具有的强大功能 所以详细介绍了 MATLAB 通信系统工具箱，并给出了基于 MATLAB 的通信系统 的调制与解调的实现，运用 MATLAB 仿真软件进行仿真。
第一章 课程设计内容及要求
学海无涯
目录
摘要------------------------------------------------------4 第一章 课程设计内容及要求--------------------------------4
matlab通信系统仿真[1]
![matlab通信系统仿真[1]](https://img.taocdn.com/s3/m/c3856f661ed9ad51f01df2b1.png)
大作业三通信系统仿真电子系无71李乾2007011252FDMA的Simulink仿真：1，连接图如下，参数设置同教材。
2，其中六个滤波器的作用，前边三个是将调制后的信号滤出来，便于将其加在一起进行传输。
后边三个是将加在一起的信号，根据其频率不同，进行选频，并且滤除噪声干扰。
将选出的信号交给后边的部件进行解调。
参数设置为：都是butterworth型带通滤波器。
前三个低频截止频率和高频截止频率分别是：Filter1:30*2*pi~50*2*piFilter2:50*2*pi~70*2*piFilter3:70*2*pi~90*2*pi后三个对应相同。
三个scope的波形依次为：第一路正弦信号：第二路方波信号：第三路锯齿波信号：三路相加的信号频谱：没深入考虑。
)TDMA的Simulink仿真：1，设计如下：pulse Genrtator单元用来产生三分之一方波，参数设置如下：Unit Delay用来产生三分之一周期的延迟，参数设置如下：产生的方波将信号分别取值并放到各自的时间域。
然后送入传输线路。
2，接受单元和时分多址复用单元完全相同。
其门控脉冲要同时发出，且周期占空比都相同。
即达到完全同步。
3，仿真结果如下：CDMA的Simulink仿真1，initial seed设置成不同值，就可以产生三路不同的信号。
第1到6号Relay单元设置如下：第7到9号Relay单元设置如下：两个Integer Delay分别设置如下：PN序列产生器：滤波器设置如下：最后得到的三路信号和输出的信号：（这部分的主要参数设置还是改sample time，不过有些单元的设置还是有些不明就里，大量参照了别人的设置。
）2，语音信号由workspace导入，并且转化成二进制码流，输入到CDMA的三个输入端，输出的二进制信号再转化成语音信号导入到workspace中。
语音信号可能会好些。
不过遗憾的是，没有听到远近效应。
三种多址方式的比较：1，关于频带资源占用，FDMA占用最大，TDMA和CDMA较小。
基于Matlab的模拟通信系统的仿真设计
目录摘要------------------------------------------------------4 第一章课程设计内容及要求--------------------------------41、课程设计的内容-----------------------------------42、课程设计的要求-----------------------------------4 第二章通信系统的调制与解调------------------------------51、通信系统的概念----------------------------------52、调制和解调的概念--------------------------------6 第三章 MATLAB软件及功能介绍------------------------------71、MATLAB软件简介-----------------------------------72、GUI功能简介--------------------------------------73、基于MATLAB相关函数介绍---------------------------8 第四章四种模拟信号的调制解调---------------------------101、AM的调制与解调---------------------------------102、DSB的调制与解调--------------------------------133、SSB的调制与解调--------------------------------164、FM的调制与解调---------------------------------195、GUI界面的设计----------------------------------23 第五章总结与结束语-------------------------------------251、各调制解调方式性能分析总结----------------------252、结束语------------------------------------------26参考文献-------------------------------------------------26摘要：通信系统是用以完成信息传输过程的技术系统的总称，作用是将信息从信源发送到一个或多个目的地。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
目录摘要------------------------------------------------------4 第一章课程设计内容及要求--------------------------------41、课程设计的内容-----------------------------------42、课程设计的要求-----------------------------------4 第二章通信系统的调制与解调------------------------------51、通信系统的概念----------------------------------52、调制和解调的概念--------------------------------6 第三章 MATLAB软件及功能介绍------------------------------71、MATLAB软件简介-----------------------------------72、GUI功能简介--------------------------------------73、基于MATLAB相关函数介绍---------------------------8 第四章四种模拟信号的调制解调---------------------------101、AM的调制与解调---------------------------------102、DSB的调制与解调--------------------------------133、SSB的调制与解调--------------------------------164、FM的调制与解调---------------------------------195、GUI界面的设计----------------------------------23 第五章总结与结束语-------------------------------------251、各调制解调方式性能分析总结----------------------252、结束语------------------------------------------26参考文献-------------------------------------------------26摘要：通信系统是用以完成信息传输过程的技术系统的总称，作用是将信息从信源发送到一个或多个目的地。
调制与解调在信息的传输过程中占据着重要的地位，是不可或缺的，因此研究系统的调制和解调过程就极为重要。
MATLAB是集数值计算、图形绘制、图像处理及系统仿真等强大功能于一体的科学计算语言，它强大的矩阵运算和图形可视化的功能以及丰富的工具箱，为通信系统的调制和解调过程的分析提供了极大的方便。
本次课程设计首先介绍了通信系统的概念，进而引出调制和解调，然后介绍了我们常用的几种调制和解调的方法。
由于MATLAB具有的强大功能所以详细介绍了MATLAB通信系统工具箱，并给出了基于MATLAB的通信系统的调制与解调的实现，运用MATLAB仿真软件进行仿真。
第一章课程设计内容及要求1、课程设计的内容1）.编写 MATLAB 程序实现 AM的调制与解调；2）.编写 MATLAB 程序实现 DSB的调制与解调；3）.编写 MATLAB 程序实现 SSB 的调制与解调；4）.编写 MATLAB 程序实现 FM 的调制与解调；5）.调用GUI实现以上相应的波形。
2、课程设计的要求1）.掌握线性幅度(AM、DSB、SSB）的调制和解调原理，以及非线性角度（FM）的调制与解调原理。
2）.学会Matlab仿真软件在幅度调制与解调和角度调制与解调中的应用。
3）.掌握参数设置方法和性能分析方法。
4）.通过实验中波形的变换，学会分析实验现象。
5）.学习掌握如何生成GUI界面以及如何使用它来显示相应的波形。
第二章通信系统的调制与解调1、通信系统的概念通信是为了传输信息，通信系统就是将信息从信源发送到一个或多个目的地，对于电通信来说，首先要把消息转变成电信号，然后经过发送设备，将信号送入信道，在接收端利用接受设备对接收信号作相应的处理后，送给信宿再转换为原来的消息，这一过程可利用图1所示的通信系统一般模型来概括。
图1 通信系统的一般模型通信(Communication)传输的消息是多种多样的，可以分成两大类：一类称为连续消息；另一类称为离散消息。
消息的传递是通过它的物理载体电信号来实现的，按信号参量的取值不同，可以把信号分为两类：模拟信号和数字信号。
通常按照信道中传输的是模拟信号还是数字信号，相应的可以把通信系统分为模拟通信系统和数字通信系统。
1）.模拟通信系统模拟通信系统是利用模拟信号来传递信息的通信系统，其模型如图2所示，其中包含两种重要变换。
第一种变换是，在发送端把连续消息变换成原始电信号，在接收端进行相反的变换，这种变换由信源和信宿来完成，通常称为原始电信号为基带信号，基带的意思是指信号的频谱从零频附近开始。
有些信道可以直接传输基带信号，而以自由空间作为信道的无线电传输却无法直接传输这些信号。
因此，模拟通信系统中常常需要进行第二种变换：把基带信号变换成适合在信道中传输的信号，并在接收端进行反变换。
完成这种变换和反变换的通常是调制器和解调器[1]。
图2 模拟通信系统模型型2.调制和解调的概念调制：把信号转换成适合在信道中传输的形式的一种过程，广义的调制分为基带调制和带通调制（也称为载波调制）。
在无线通信中和其他大多数场合，调制均指载波调制。
载波调制，就是用调制信号去控制载波的参数的过程，使载波的某一个或某几个参数按照调制信号的规律而变化。
调制信号是指来自信源的消息信号（基带信号），这些信号可以是模拟的，也可以是数字的。
未受调制的周期性振荡信号称为载波，它可以是正弦波，也可以是非正弦波（如周期性脉冲序列）。
载波调制后称为已调信号，它含有调制信号的全部特征。
解调：将已调信号中的调制信号恢复出来，是调制的逆过程。
调制方式不同，解调方法也不一样。
解调可以分为正弦波解调（有时也称为连续波解调）。
正弦波解调还可再分为幅度解调、频率解调和相位解调，此外还有一些变种如单边带信号解调、残留边带信号解调等。
同样，脉冲波解调也可分为脉冲幅度解调、脉冲相位解调、脉冲宽度解调和脉冲编码解调等。
调制方式有很多，根据调制信号是模拟信号还是数字信号，载波是连续波（通常是正弦波）还是脉冲序列，相应的调制方式有模拟连续波调制（简称模拟调制）、数字连续波调制（简称数字调制）、模拟脉冲调制和数字脉冲调制等。
最重要和最常用的模拟调制方式是用正弦波作为载波的幅度调制和角度调制。
常见的调幅、双边带、单边带和残留边带等调制就是幅度调制的几个典型实例。
解调的方法可以分为两类：相干解调和非相干解调（包络检波）。
解调过程与采用何种解调方式有关，对于常规幅度调制，一般用包络检波进行解调，由于在这种解调方式中，接收机对载波频率和相应精度的了解是无关紧要的，所以解调过程相对简单。
对于DSB调制和SSB调制，用相干解调的方法，它要求在接收机中有一个与载波同频同相的信号，接收机中产生所需要的正弦波振荡器，为本地振荡器。
数字通信系统中，采用键控调制方式，这里不一一介绍。
第三章 MATLAB软件及功能介绍1、MATLAB软件简介MATLAB是matrix&laboratory两个词的组合，意为矩阵工厂（矩阵实验室）。
是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。
它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。
MATLAB和Mathematica、Maple并称为三大数学软件。
它在数学类科技应用软件中在数值计算方面首屈一指。
MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。
编程创造的功能：在命令行或者在手稿里确定单个行的函数的匿名的函数，当任何表示是真实的时，有条件的休息点，让你停止块意见出于记分员可发表意见的一个代码的整个部分数学整数算术，让你处理更大的整数数据集合单精度运算，线性代数，FFT和过滤，使你能够处理更大的单精度数据集合更坚固的计算几何学程序使用Qhull.1，给更大的对算法选择的控制linsolve 功能，使你能够迅速通过指定系数的基体的结构解决线性方程序的系统和multipoint边值问题编程环境：MATLAB由一系列工具组成。
这些工具方便用户使用MATLAB的函数和文件，其中许多工具采用的是图形用户界面。
包括MATLAB桌面和命令窗口、历史命令窗口、编辑器和调试器、路径搜索和用于用户浏览帮助、工作空间、文件的浏览器。
随着MATLAB的商业化以及软件本身的不断升级，MATLAB的用户界面也越来越精致，更加接近Windows的标准界面，人机交互性更强，操作更简单。
而且新版本的MATLAB提供了完整的联机查询、帮助系统，极大的方便了用户的使用。
简单的编程环境提供了比较完备的调试系统，程序不必经过编译就可以直接运行，而且能够及时地报告出现的错误及进行出错原因分析。
2、GUI功能简介GUI 是 Graphical User Interface 的简称，即图形用户界面，通常人机交互图形化用户界面设计经常读做“goo-ee”准确来说 GUI 就是屏幕产品的视觉体验和互动操作部分。
GUI 是一种结合计算机科学、美学、心理学、行为学，及各商业领域需求分析的人机系统工程，强调人—机—环境三者作为一个系统进行总体设计。
这种面向客户的系统工程设计其目的是优化产品的性能，使操作更人性化，减轻使用者的认知负担，使其更适合用户的操作需求，直接提升产品的市场竞争力。
主要组成部分：多文件界面在一个窗口之内进行多个数据管理的方式。
这种情况下，窗口的管理简单化，但是操作变为双重管理。
标签多文件界面的数据管理方式中使用的一种界面，将数据的标题在窗口中并排，通过选择标签标题显示必要的数据，这样使得接入数据方式变得更为便捷。
上述中，多文件界面主要是微软视窗系统采用。
而在其他环境中，通常多是单文件界面，所以无所谓单一/多文件界面的称呼问题。
菜单将系统可以执行的命令以阶层的方式显示出来的一个界面。