通信原理大作业 (2).pdf

上传人:修**** 文档编号:75980992 上传时间:2023-03-06 格式:PDF 页数:15 大小:1.33MB
返回 下载 相关 举报
通信原理大作业 (2).pdf_第1页
第1页 / 共15页
通信原理大作业 (2).pdf_第2页
第2页 / 共15页
点击查看更多>>
资源描述

《通信原理大作业 (2).pdf》由会员分享,可在线阅读,更多相关《通信原理大作业 (2).pdf(15页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、通信原理大作业通信原理大作业年年级级20102010 级级班班级级0210*0210*姓名学号姓名学号张张*专专业业电子信息工程电子信息工程学学院院电子工程学院电子工程学院西安电子科技大学2013 年 6 月选题一:选题一:2FSK2FSK 信号传输仿真信号传输仿真一、一、题目要求:题目要求:按照 2FSK 产生模型和解调模型分别产生 2FSK 信号和高斯白噪声,经过信道传输后进行解调。对调制解调过程中的波形进行时域和频域观察,并且对解调结果进行误码率测量。2FSK 信号的解调可以选用包络解调或者相干解调法。二、二、2FSK2FSK 信号传输介绍信号传输介绍数字信号的传输方式分为基带传输和带通

2、传输。然而,实际中的大多数信道因具有带通特性而不能直接传送基带信号。为了使数字信号在带通系统中传输,必须用数字基带信号对载波进行调制,以使信号与信道的特性相匹配。这种用数字基带信号控制载波,把数字基带信号变换为数字带通信号的过程称为数字调制。在接收端,通过解调器把带通信号还原为数字基带信号的过程称为数字解调。移频键控(FSK)是数据通信中最常用的一种调制方式。FSK 方法简单,易于实现,并且解调不需要恢复本地载波,可以异步传输,抗噪声和抗衰落性能较强。缺点是占用频带较宽,频带利用不够经济。FSK 主要应用于低中速数据传输,以及衰落信道和频带较宽的信道中。三、三、2FSK2FSK 信号传输实现方

3、法信号传输实现方法1.调制过程FSK 即 Frequency Shift Key 叫做移频键控或频移键控。调制的方法一般有两种,一种叫直接调频法,另一种叫键控法。所谓直接调频法,就是将输入的基带脉冲去控制一个振荡器的某种参数,而达到改变振荡频率的目的。键控法就是利用矩形脉冲序列控制的开关电路,对两个不同的独立频率源进行选通。一般来说,键控法采用两个独立的振荡器,得到的是相位不连续的FSK 信号;而且直接调频法f1,f2 由同一个谐振电路产生,则得到相位连续的FSK 信号。2FSK 信号便是 0 符号对应于载频 f1,1 符号对应于载频 f2(与 f1 不同的另一个载频)的一调制波形,而 f1

4、与 f2 的改变是瞬间完成的。2.解调过程2FSK 的常用解调方法是采用相干解调(如下图)。解调原理是将 2FSK 信号分解为上下两路 2ASK 信号分别进行解调,然后进行判决。判决规则与调制规则相呼应,调制时若规定“0”符号对应载波频率 f1,则接收时上支路的样值较大,应判为“0”,反之则为“1”。1带 通 滤波器相乘器低 通 滤波器抽样脉冲低 通 滤波器抽 样 判决器输入输出Cos1t带 通 滤波器2相乘器Cos2t四、四、MatlabMatlab 仿真思路仿真思路1.首先设定仿真的码元个数、载波频率 f1,f2、采样频率和波特率。2.利用 Matlab 函数 rand 和 round 设

5、定初始码元,并计算得出基带信号及其反码。3.产生载波信号,利用模拟相乘得到 FSK 信号,再利用 awgn 函数在信号中加入高斯白噪声得到已调信号表达式 y_2FSK_n。4.在 2FSK 的解调过程中,根据原理图,信号首先通过带通滤波器,设置带通滤波器参数(firl 函数),后利用一维数字滤波器 filter 对已调信号进行滤波处理。输出经过带通滤波器后的信号波形。由于已调信号中有两个不同的载波(1,2),则经过两个不同频率的带通滤波器后输出两个不同的信号波形 z1,z2。5.经过带通滤波器后的 2FSK 信号再经过相乘器(cos1,cos2),输出得到相乘后的两个不同的 2FSK 波 st

6、1_1 和 st2_1。6.经过相乘器输出的波形再通过低通滤波器,设置低通滤波器的参数(函数 firl),用一维数字滤波韩式 filter 对信号的数据进行新的一轮的滤波处理。输出经过低通滤波器后的两个波形(st1_2,st2_2)。7.将信号(st1_2,st2_2)同时经过抽样判决器,输出code_s_1=1。code_s_1。对抽样判决器经定义一个时间变量长度 i,当 st1_2(i)=st2_2(i)时,则 code_s_1=0,否则8.计算误码率,将 code_s_1 进行还原,还原办法是将一个码元周期内的 0,1 个数比较,取较多的那一个。将结果与码元 code 比较,将传输有误的

7、码元个数除以总码元个数即可。五、五、MatlabMatlab 仿真结果仿真结果在仿真过程中,设定码元个数 10、载波频率 fc1 为 30khz,fc2 为 60khz、采样频率 fc2*5、信号波特率800。以下图片是在数据码流为1100011100 时仿真的结果。1.基带信号,2FSK 信号以及 2FSK 信号混高斯白噪声的时域波形2.2FSk 信号频谱分析3.2FSK 信号经过带通滤波器后频域波形和时域波形4.两路信号经低通滤波器后的频谱和时域波形5.解调后误码率计算经过计算,误码率为 0六、六、仿真过程中问题仿真过程中问题1.2FSK 信号传输的误码率跟载波频率和采样频率有关,载波频率

8、和采样频率均越低,出现误码的可能性越高。2.2FSK 信号传输的误码率跟波特率有关,波特率越高,误码率越高。3.经过分析,上述两个因素综合起来就是码元周期,经过仿真与分析,码元周期越长越好4.同时也更加深入的了解了为何国际电信联盟会建议数据率低于1200b/s 时采用2FSK 体制选题二:信道噪声特性仿真选题二:信道噪声特性仿真一、一、题目要求:题目要求:产生信道高斯白噪声,设计信道带通滤波器对高斯白噪声进行滤波,得到窄带高斯噪声。对信道带通滤波器的输入输出的噪声的时域、频域特性进行统计和分析,画出其时域和频域的图形。二、二、信道噪声高斯白噪声介绍信道噪声高斯白噪声介绍在分析通信系统的抗噪声性

9、能时,常用高斯白噪声作为通信信道中的噪声模型。这是因为,通信系统中常见的热噪声近似为白噪声,且热噪声的取值恰好服从高斯分布。另外,实际信道或滤波器的带宽存在一定限制,白噪声通过后,其结果是带限噪声,若其谱密度在通带范围内仍具有白色特性,则称其为带限白噪声,他又可以分为低通白噪声和带通白噪声。如果白噪声取值的概率分布服从高斯分布,则称之为高斯白噪声,我们常用它作为通信信道中的噪声模型。高斯白噪声在任意两个不同时刻上的随机变量之间,不仅是互不相关的,而且还是统计独立的。三、三、滤波器的实现滤波器的实现滤波器是一种对信号有处理作用的器件或电路。分为有源滤波器和无源滤波器。主要作用是让有用信号尽可能无

10、衰减的通过,对无用信号尽可能大的反射。滤波器一般有两个端口,一个输入信号、一个输出信号,利用这个特性可以选通通过滤波器的一个方波群或复合噪波,而得到一个特定频率的正弦波。滤波器的响应种类有巴特沃斯响应,贝塞尔响应,切比雪夫响应。这里我采用的是巴特沃斯滤波器。巴特沃斯响应能够最大化滤波器的通带平坦度。该响应非常平坦,接近 DC 信号,然后慢慢衰减至截止频率点为-3dB,最终逼近-20ndB/decade的衰减率,其中n 为滤波器的阶数。巴特沃斯滤波器特别适用于低频应用,其对于维护增益的平坦性来说非常重要。四、四、MatlabMatlab 仿真思路仿真思路1.用 wgn 函数产生高斯白噪声信号 y

11、1,频谱信号为 ya。2.用 butter 函数产生巴特沃斯滤波器,滤波器的系统函数为 filter_h。3.使信号 ya 通过滤波器产生信号 yb,通过反傅里叶变换得到时域信号 y2。4.对时域信号进行均值,最大值,最小值,标准差,方差,自相关函数的分析。5.对频域信号进行频谱,功率谱密度(自相关函数的频谱)分析五、五、MatlabMatlab 仿真结果仿真结果使用采样 1000 个点且采样频率为 1khz 的高斯白噪声模拟信道中的噪声。1.高斯白噪声的时域图和频谱图2.滤波器的幅频特性3.对高斯白噪声和窄带高斯白噪声时域的分析波均值形滤0.0134波器前滤-9.1e-19+i-1.3e-1

12、波9器后-0.001-i0.0004-5.9888+i0.31440.51730.2675-3.09962.8918极小值极大值标准方差差1.03051.06204.通过滤波器前后的噪声频谱图对比5.通过滤波器前后的功率谱密度图对比六、六、仿真中遇到的问题仿真中遇到的问题功率谱密度理论上来说应该在频带上均匀分布,而仿真结果却不是。附录附录一、一、2FSK2FSK 信号传输代码信号传输代码%2FSK 信号传输仿真,省略画图代码clc;close all;clear allcodn=10;%仿真的码元个数fc1=30000;fc2=60000;%设定载波频率 Fc1=30khz,Fc2=60Khz

13、fs=fc2*5;%数据采样率bode=800;%信号波特率code=round(rand(1,codn);code_len=round(1/bode/(1/fs);%得到一个码元周期的数据长度%产生基带信号st1=0;for i=1:codnif code(i)=1for m=(i-1)*code_len+1:i*code_lenst1(m)=0;endelsefor m=(i-1)*code_len+1:i*code_lenst1(m)=1;endendendn=length(st1);%产生基带信号的反码st2=0;for i=1:nif st1(i)=1st2(i)=0;elsest2

14、(i)=1;endendt=(1:n)/fs;s1=cos(2*pi*fc1*t);%产生频率为 Fc1 和 Fc2 的载波 s1,s2s2=cos(2*pi*fc2*t);y_2fsk=st1.*s1+st2.*s2;%采用模拟相乘法产生 FSK 信号y_2fsk_n=awgn(y_2fsk,5);%生成高斯白噪声ay=abs(fft(y_2fsk_n);%进行幅频特性分析f=0:fs/length(y_2fsk_n):fs/2;ay=ay(1:length(f);%低通滤波b1=fir1(101,0.18 0.22);b2=fir1(101,0.38 0.42);%带通滤波器系数z1=fi

15、lter(b1,1,y_2fsk_n);%经过带通滤波器的信号z2=filter(b2,1,y_2fsk_n);ay_1_1=abs(fft(z1);%进行幅频特性分析ay_1_1=ay_1_1(1:length(f);ay_1_2=abs(fft(z2);ay_1_2=ay_1_2(1:length(f);st1_1=z1.*s1;%相乘器st2_1=z2.*s2;bn=fir1(101,0.04 0.16);%低通滤波器st1_2=filter(bn,1,st1_1);st2_2=filter(bn,1,st2_1);ay_2_1=abs(fft(st1_2);%进行幅频特性分析ay_2_

16、1=ay_2_1(1:length(f);ay_2_2=abs(fft(st2_2);ay_2_2=ay_2_2(1:length(f);%判决器for i=1:nif st1_2(i)=st2_2(i)code_s_1(i)=0;elsecode_s_1(i)=1;endend%计算误码率for i=1:codntag=code_s_1(i-1)*code_len+1:i*code_len);j=length(find(tag=0);if j=(code_len/2)code_s(i)=1;elsecode_s(i)=0;endendP1=sum(abs(code_s-code)/codn;

17、%误码率二、二、信道噪声特性仿真代码信道噪声特性仿真代码%产生高斯白噪声省略画图代码t=0:0.001:0.999;%确定取点数和采样频率y1=wgn(1,1000,0.5);%高斯白噪声ya=abs(fft(y1);f1=(0:length(ya)-1)*1/length(ya);yb=ya.*filter_h;y2=ifft(yb);%滤波器ws=0.35,0.65;%ws 为归一化阻带边界频率wp=0.36,0.64;%wp 为归一化通带边界频率%Ap 单位为 dB 的通带波形Ap=3;As=15;%As 单位为 dB 的最小阻带衰减n1,wn1=buttord(wp,ws,Ap,As)

18、;%n1为数字滤波器的阶数%wn1 为归一化截止频率num,den=butter(n1,wn1);%设计数字滤波器,num 系统函数的分子%计算单位为 dB 的增益函数den 系统函数的分母999:pi;%单位园上将函数转化为999 等分w=0:pi/%计算相位filter_h=freqz(num,den,w);%频率响应计算0,pig=20*log10(abs(filter_h);%计算幅度pha=angle(filter_h);%相角计算%时域分析a,b=xcorr(y1);%高斯白噪声的自相关性c,d=xcorr(y2);%窄带高斯白噪声的自相关性m1=mean(y1);disp(m1)

19、;%计算平均值mi1=min(y1);disp(mi1);%极小值mx1=max(y1);disp(mx1);%极大值st1=std(y1);disp(st1);%标准差fc1=st1.2;%方差m2=mean(y2);disp(m2);%计算平均值mi2=min(y2);disp(mi2);%极小值mx2=max(y2);disp(mx2);%极大值st2=std(y2);disp(st2);%标准差fc2=st2.2;%方差%频域分析fa=(0:length(ya)-1)*1/length(ya);%换算出对应时域的频域横坐标值fb=(0:length(yb)-1)*1/length(yb);%功率谱密度nfft=1500;CXk1=fft(a,nfft);Pxx1=abs(CXk1);index=0:round(nfft-1);k=index*1000/nfft;plot_Pxx1=10*log10(Pxx1(index+1);CXk2=fft(c,nfft);Pxx2=abs(CXk2);plot_Pxx2=10*log10(Pxx2(index+1);

展开阅读全文
相关资源
相关搜索

当前位置:首页 > 教育专区 > 高考资料

本站为文档C TO C交易模式,本站只提供存储空间、用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。本站仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知淘文阁网,我们立即给予删除!客服QQ:136780468 微信:18945177775 电话:18904686070

工信部备案号:黑ICP备15003705号© 2020-2023 www.taowenge.com 淘文阁