通信原理课程设计报告

遥光农林哥牧大学通信系统课程设计报告单边带电台仿真TOC\o"1-5"\h\z院 （系）:专业年级（班）:学 生:学 号:指导教师:完成时间:摘要这篇文章介绍了使用MATLAB软件编制程序，以实现单边带信号的调制和解调。首先，利用相移法从语音信号的双边带信号得到单边带信号，再编写MATLAB程序，使单边带信号得到调制和调解。分析调制前后的时域和频域波形图，以更加深入理解单边带信号的调制和解调的原理。关键词：单边带；调制和解调；MATLABSinglesidebandradioSimulationAbstract：ThisarticleintroducesthewaytomodulateanddemodulatethesinglesidebandwiththesoftprogramfromMATLAB.First,getthesinglesidebandsignalfromthespeechsignal'scouplesidesband,thenwrighttheMATLABprogramtomodulateanddemodulatethesignal.Analyzethepictures,andunderstandthetheoryfurtherly.Keywords：singlesideband,modulationanddemodulation,MATLAB目录TOC\o"1-5"\h\z设计任务与要求 -1-开发环境 -1-三•电台设计原理及方案 -1 -\o"CurrentDocument"3.1发送模块的设计方案 -1 -3.2信道模块的设计方案 -1 -\o"CurrentDocument"3.3接收模块的设计方案 -2 -四.电台详细设计 - 2 -\o"CurrentDocument"4.1对音频输入与调制的仿真 -2 -\o"CurrentDocument"4.2对指定信噪比信道的仿真 -2 -4.3话音信号的解调 -3 -\o"CurrentDocument"载频9.8KHz解调 -3 -\o"CurrentDocument"载频9.9KHz解调 -4 -载频lOKHz解调 -6 -4.4原信号和解调后的信号对比 -7 -\o"CurrentDocument"4.4.1原信号和解调后信号波形对比 -7 -4.4.2原信号和解调后信号频率对比 -8 -五•设计总结 -8 -参考文献： - 9 -\o"CurrentDocument"附录： -10-.设计任务与要求设计任务与要求这个仿真以真实的音频信号作为输入，设计一个单边带发信机。将基带信号调制为SSB信号后送入带通型高斯噪声信道，加入给定功率的噪声之后，再送入单边带接收机。单边带接收机将型号解调下来，通过计算机声卡将解调信号播放出来试听效果，从而对信道信噪比与解调音质之间的关系进行主观测试。现设计一个单边带发信机、带通信道和相应的接收机，参数定位如下值。1•输入信号为一个话音信号，采样率为8000Hz。话音输入后首先进行预滤波，滤波器是一个频率范围在[300,3400]Hz的带通滤波器，其目的是将话音频频谱限制在3400Hz以下。单边带调制的载波频率设计为10kHz，调制输出上边带。要求观测单边带调制前后的信号功率谱。2•信道是一个带限高斯噪声信道，其通带频率范围是[10000,13500]Hz。能够根据信噪比SNR的要求加入高斯噪声。3•接收机采用相干解调方式。为了模拟载波频率误差对解调话音音质的影响，相干解调，设本地载波频率分别为为9.8kHz，9.9kHz,10kHz。解调滤波器设计为300〜3400Hz的带通滤波器。开发环境操作系统：MicrosoftWindows7Ultimate（64位）交互工具：键盘/鼠标开发工具：MATLAB2010a电台设计原理及方案3.1发送模块的设计方案有限语音信号可以采用直接用程序读取，当读取一次之后，可以把音频文件保存起来，当用到没有麦克风的电脑时可以备用，发送的语音信号首先进行希尔比特变化，然后用函数可以进行单边带调制，保存单边带语音信号。3.2信道模块的设计方案读取发送模块保存的语音信号，然后加入高斯白噪声，然后改变信噪比，通常情况下加入的信噪比为20,然后滤波，模拟现实当中的信道传输，并保存得到的语音信号。3.3接收模块的设计方案首先读取信道传输过来的语音信号，加入不同频率的载波，然后解调，并观察分析不同载波解调出来的语音信号和原语音信号的区别，分析失真度。电台详细设计4.1对音频输入与调制的仿真根据设置的参数，系统中信号最高频率约为14kHz。为了较好地显示调制波形，系统仿真采样率设为50kHz，满足取样定理。由于话音信号的采样率为8000Hz，与系统仿真采样率不等，因此，在进行信号处理之前，必须将话音的采样率提高到50kHz，用插值函数来做这一任务。先编写程序将基带音频信号读入，进行[300，3400]Hz的带通滤波，并将信号采样率提高到50kHz，进行单边带调制之后，将调制输出结果保存为wav文件，文件名为SSB_OUT.wav。 'Figure1 |=|回FileEditViewInsertToolsDesktopWindowHelp浊基带信号波形0.50-0.5基带信号功率谱基带信号波形0.50-0.5基带信号功率谱100-10-2001002时间（样值数）日信号波形01002时间（样值数）日信号波形010002000频率FH巨SSB信号功率谱100-10-20图1基带信号和SSB的波形、功率谱分析调制前后的波形频率，可以得到SSB调制只是把原来波形搬移到1OKHz位置，并四.电台详细设计错误！未指定书签。没有改变原来波形的形状和频率范围，说明调制的效果很好。4.2对指定信噪比信道的仿真信道是一个带限高斯噪声信道，其带通频率范围是[10000,13500]Hz.,目前输入的信噪比SNR为20。仿真指定信噪比信道计算信噪比为20dB时的信道输出，将结果保存为Chanel_out.wav文件。方便以后的实用。运行结束：计算出的信噪比结果SNR_dB=19.91264.3话音信号的解调4.3.1载频9.8KHz解调接收机采用想干解调方式。为了模拟载波频率误差对解调话音音质的影响，社本地载波的频率为9.8KHz，与发信机载波频率相差200Hz。解调滤波器设计为300Hz至【」3400Hz。QFigure2 —I口-|區FileEditViewInsertToolsDesktopWindowHelp$nq凤饨脅⑥晏銘丁凰□匡| >>105）100200时问（样倩埶）解调信号液形0.30.20.10-0.1-0.2mphB&FI-20解调信号功率105）100200时问（样倩埶）解调信号液形0.30.20.10-0.1-0.2mphB&FI-20解调信号功率i晋10002000OOJ—■5图29.8KHz解调后波形和功率谱解调输出信号被保存为音频文件SSBDemo_0UT9.8KHz.wav,并由sound函数播放。聆听播放解调输出信号的声音可知，在20dB信道信噪比条件下，即使解调本地载波频率误差达到200Hz，声音仍然是清晰可懂的。mphgs槪悴O2OOJ—■mphgs槪悴O2OOJ—■图39.8KHz解调后波形和功率谱解调输出信号被保存为音频文件SSBDemo_0UT9.8KHz.wav,并由sound函数播放。聆听播放解调输出信号的声音可知，在10dB信道信噪比条件下，即使解调本地载波频率误差达到200Hz，声音仍然清晰地，但是也略有失真，效果明显比信噪比为20dB差。4.3.2载频9.9KHz解调接收机采用想干解调方式。为了模拟载波频率误差对解调话音音质的影响，社本地载波的频率为9.9KHz，与发信机载波频率相差100Hz。解调滤波器设计为300Hz到3400Hz。将得到的波形和波形的功率谱输出。四.电台详细设计错误！未指定书签。FileEditView0.3r0.2-0.1-0--0.1-0.2■解调信号液形10oJI■FileEditView0.3r0.2-0.1-0--0.1-0.2■解调信号液形10oJI■-55解调信号功率谱10002000-201002000时间InsertToolsDesktopWindowHelp勺□a 陰電凤鋼®晏摇，凰|□匡！|图49.9KHz解调后波形和功率谱解调输出信号被保存为音频文件SSBDemo_0UT9.9KHz.wav,并由sound函数播放。聆听播放解调输出信号的声音可知，在20dB信道信噪比条件下，即使解调本地载波频率误差达到100Hz，声音仍然是清晰可懂的。Figure3 k=I回I—FileEditViewInsertToolsDesktopWindowHelp□aqte電鱼巴⑧晏銘〒凰|口匡1w106100200鮮调信号汲形0.30.20.10-0.1-0.2mp)、B»£I106100200鮮调信号汲形0.30.20.10-0.1-0.2mp)、B»£InuJ—--20鮮调信号功率诸10002000图59.9KHz解调后波形和功率谱解调输出信号被保存为音频文件SSBDemo_0UT9.9KHz.wav,并由sound函数播放。聆听播放解调输出信号的声音可知，在10dB信道信噪比条件下，即使解调本地载波频率误差达到100Hz，声音仍然是清晰可懂的，但是效果比信噪比为20dB效果要差一点,不过基本分不出来。4.3.3载频10KHz解调接收机采用想干解调方式。为了与上述不同载波解调出的结果相对比，这次的在载波和原来调制的载波相同。解调滤波器设计为300Hz到3400Hz。将得到的波形和波形的功率谱输出Figure4 〔=I回1—^]FileEditViewInsertToolsDesktopWindowHelp理解调信号滋形】100200时问（存倩埶）0.30.20.10-0.1-0.210解调信号功率谱解调信号滋形】100200时问（存倩埶）0.30.20.10-0.1-0.210解调信号功率谱-150100020005OOJ—■-20图610KHz解调后波形和功率谱解调输出信号被保存为音频文件SSBDemo_OUT10KHz.wav，并由sound函数播放。聆听播放解调输出信号的声音可知，在解调波和载波频率相同，20dB信道信噪比条件下,声音是很清晰的。四.电台详细设计错误！未指定书签。图7lOKHz解调后波形和功率谱解调输出信号被保存为音频文件SSBDemo_0UT10KHz.wav,并由sound函数播放。聆听播放解调输出信号的声音可知，在解调波和载波频率相同，1OdB信道信噪比条件下,声音是很清晰的，并且和信噪比为20dB基本没有差别，耳朵几乎分辨不出。4.4原信号和解调后的信号对比4.4.1原信号和解调后信号波形对比图8原信号波形与解调后信号波形对比从波形上来看，三个频率解调出来的音频信号波形整体差别不大，要是要仔细看,在低频部分最后一个频率解调出来的波形效果最好。4.4.2原信号和解调后信号频率对比MFigure5FileEditViewInsertToolsDesktopWindowHelp勺

dHa心k\n®汇銘戸凰口 -基带信号功率谱频率/吃9£KHz解调信号功率谱9一眼七解调信号功率谱频率『Hz基带信号功率谱频率/吃9£KHz解调信号功率谱9一眼七解调信号功率谱频率『Hz频率/He■WKHe解调信号功率谱010002000频率FHz图9原信号波形与解调后信号频谱对比从频率上来看，三个频率解调出来的音频信号波形整体差别不大，在细微处还是有差别的，尤其是在高频和低频处，有部分的失真，但是对音频的整体效果影响不大，最后一个借条出来的频谱和原信号的频谱最接近，因此也是最好的解调频率。五•设计总结通过这门实验使我学习掌握了许多知识。首先是对matlab有了一个全新的认识，其次是对matlab的更多操作和命令的使用有了更高的掌握，最重要的事对matlab的处理能力有了一个更高的飞跃尤其是对相关函数的使用及相关问题的处理。就对matlab相关的命令操作而言，通过这次实验的亲身操作和实践，学习掌握了参考文献：错误！未指定书签。许多原本不知道的或者不太熟悉的命令。比如说相关m文件的建立，画图用到的标注，配色，坐标控制，同一张图里画几幅不同的图像，相关参数的设置以及相关函数的调用格式等等。就拿建立一个数学方程而言，通过设置不同的参数达到所需要的要求和结果，而且还可以在不同的窗口建立不同的函数而达到相同的效果。而自己对于矩阵及闭环传递函数的建立原本所掌握的知识几乎为零，而通过这次实验使我彻底的掌握了相关的命令操作和处理的方法，在这里我们不仅可以通过建立函数和参数来达到目标效果，而且还可以通过可视化的编程达到更快更方便，更简洁的效果。就拿可视化编程而言原本根本就只是听说而已罢了，从来就没有亲身去尝试过，然而现在自己却可以和容易的通过搭建不同功能木块来实现相关的函数及功能。这些在原本根本就不敢相信，然而通过学习和实验亲身操作这些原本看似不可能的操作在此就变的轻而易举的事了。通过对同一个模块分析其对应的不同的参数分析图的建立去分析和解释其对应的相关功能和技术指标和性能分析是非常重要的，我们不可能只需要建立相关的模块和功能就说自己掌握了所有的相关知识和技术，真正的技术和知识是怎么去分析和解释相关的技术指标和功能参数才是重中之重。就此而言，我坦诚的说自己所掌握的还是十分的有限的，但是老师给我们介绍的相关方法和技巧还是十分有效果的，如果自己真的想在这方面有什么建树对自己以后的要求还是需要更改的要求的，万不可以就此止步不前，自命不凡，我们还需掌握和了解还有许多许多，我们真正所掌握的只是皮毛，要想取得更大的成绩就得不断的去努力学习和汲取相关的知识和技巧。万不可自以为傲，裹足不前,matlab真的是个非常强大和有用的工具我们真正的能把它学懂学透的话还是需要下非常大的功夫和努力的。参考文献:张渭滨，陈方，苏武浔等.Chen-Mobius多路数字通信系统的 Matlab仿真[J].系统仿真学报，2007,19(8):1886-1889.D0I:10.3969/j.issn.l004-731X.2007.08.058.刘素心，王汝芳，张广森等•基于Matlab的通信系统的实验仿真设计[J].实验室科学，2008,(3):101-103,105.D0I:10.3969/j.issn.l672-4305.2008.03.043.王峰，丁金林.基于MATLAB/GUI的数字通信系统虚拟实验平台的设计[J].南通职业大学学报，2013,27(l):96-100.D0I:10.3969/j.issn.l008-5327.2013.01.026.陈蕾，姚远程，秦明伟等.自适应抗干扰通信系统中频谱感知技术研究[J].电视技术，2014,38(5):101-104.梅志红,邓文华.基于Matlab与VB集成技术的通信系统仿真软件的开发 [J].计算机仿真，2004,21(6):195-199.D0I:10.3969/j.issn.l006-9348.2004.06.056.范伟，翟传润，战兴群等•基于MATLAB的扩频通信系统仿真研究[J].微计算机信息，2006,22(19):242-244.D0I:10.3969/j.issn.l008-0570.2006.19.086.周晓兰，张杰.MATLAB在通信系统仿真中的应用[J].计算机技术与发展，2006,16(9):166-168.D0I:10.3969/j.issn.l673-629X.2006.09.057.解伟俊，胡修林，张蕴玉等.基于MATLAB的仿真通信语音库实现[J].华中理工大学学报，2000,28(2):77-79.D0I:10.3321/j.issn:1671-4512.2000.02.027.附录：%FileName:ForSSB.mclc;clearall;%功能，采样点数40000,采样率为8000%jilu=wavrecord(5*8000,8000,'double');%wavwrite(jilu,'GDGvoice8000.wav');[wav,fs]=wavreadCGDGvoice8000.wav');%计算声音时间长度t_end=1/fs*length(wav);%仿真系统采样时间点Fs=50000;t=1/Fs:1/Fs:t_end;%设计300~3400hz的带通滤波器[fenzi,fenmu]=butter(3,[3003400]/(fs/2));%对音频信号进行滤波wav=filter(fenzi,fenmu,wav);%输出滤波后的声音wavwrite(wav,'LVBO_OUT.wav');%利用插值函数将音频的采样率提升为50khzwav=interp1([1/fs:1/fs:t_end],wav,t,'spline');%音频信号的希尔伯特变换wav_hilbert=imag(hilbert(wav));%载波频率fc=10000;%单边带调制附录：错误！未指定书签。SSB_OUT=wav.*cos(2*pi*fc*t)-wav_hilbert.*sin(2*pi*fc*t);%输出波形figure(l);subplot(2,2,1);plot(wav(53550:53750));axis([0200-0.50.5]);title('基带信号波形');xlabel('时间(样值数)');subplot(2,2,2);psd(wav,10000,Fs);axis([02500-2010]);title('基带信号功率谱');xlabel('频率/Hz');ylabel('功率谱/(dB)');subplot(2,2,3);plot(SSB_OUT(53550:53750));axis([0200-0.50.5]);title('SSB信号波形');xlabel('时间(样值数)');subplot(2,2,4);psd(SSB_OUT,10000,Fs);axis([025000-2010]);title('SSB信号功率谱');xlabel('频率/Hz');ylabel('功率谱/(dB)');%将SSB调制输出存盘备用wavwrite(0.5*SSB_OUT,Fs,'SSB_OUT.wav');%FileName:ForSNB.m%功能;计算实际信噪比%clear;[in,Fs]=wavreadCSSB_OUT.wav');%设定信噪比为20时的新到输出SNR_db=20;out=ChanelSimulink(in,SNR_db);wavwrite(out,Fs,'Chanel_OUT.wav');%FileName:ChanelSimulink.mfunctionout=ChanelSimulink(in,SNR_db)%SNR_db设定信噪比%in输入信号序列%out信道输出序列%系统采样率Fs=50000;Power_of_in=var(in);Power_of_noise=Power_of_in/(10."(SNR_db/10));%信道带宽bandwidth=13500-10000;%噪声功率谱密度值W/HzNO=Power_of_noise/bandwidth;Gause_noise=sqrt(N0*Fs/2).*randn(size(in));%噪声通道10〜13.5kHz[num,den]=butter(4,[1000013500]/(Fs/2));signal_of_filter_out=filter(num,den,in);noise_of_filter_out=filter(num,den,Gause_noise);SNR_dB=10*logl0(var(signal_of_filter_out)/var(noise_of_filter_out));%测量得出信噪比%信道输出out=signal_of_filter_out+noise_of_filter_out;SNR_dB%FileName:Forjietiao9.8K.m%clear;Fs=50000;%读入信道输出信号数据[recvsignal,Fs]=wavread('Chanel_OUT.wav');t=(1/Fs:l/Fs:length(recvsignal)/Fs)';%本地载波频率fc_locall=10000-200;%fc_local=10000-100;%fc_local=10000;%本地载波local_carrier1=cos(2*pi*fc_local1.*t);%相干解调xianggan_out1=recvsignal.*local_carrier1;%设计300〜3400Hz的带通滤波器[fenzi,fenmu]=butter(3,[3003400]/(Fs/2));demod_out1=filter(fenzi,fenmu,xianggan_out1);%sound(demod_out1/max(demod_out1),Fs);wavwrite(demod_out1,Fs,'SSBDemod_0UT_9.8k.wav');figure(2);subplot(1,2,1);plot(5*demod_out1(53550:53750));axis([0200-0.30.3]);title('解调信号波形');xlabel('时间(样值数)');subplot(1,2,2);psd(5*demod_out1,10000,Fs);axis([02500-2010]);title('解调信号功率谱');xlabel('频率/Hz');ylabel('功率谱/(dB)');%FileName:Forjietiao9.9K.m%clear;Fs=50000;%读入信道输出信号数据[recvsignal,Fs]=wavread('Chanel_OUT.wav');t=(1/Fs:1/Fs:length(recvsignal)/Fs)';%本地载波频率%fc_local=10000-200;fc_local2=10000-100;附录：错误！未指定书签。%fc_local=10000;%本地载波local_carrier2=cos(2*pi*fc_local2.*t);%相干解调xianggan_out2=recvsignal.*local_carrier2;%设计300〜3400Hz的带通滤波器[fenzi,fenmu]=butter(3,[3003400]/(Fs/2));demod_out2=filter(fenzi,fenmu,xianggan_out2);%sound(demod_out2/max(demod_out2),Fs);wavwrite(demod_out2,Fs,'SSBDemod_0UT_9.9k.wav');figure(3);subplot(1,2,1);plot(5*demod_out2(53550:53750));axis([0200-0.30.3]);title('解调信号波形');xlabel('时间(样值数)');subplot(1,2,2);psd(5*demod_out2,10000,Fs);axis([02500-2010]);title('解调信号功率谱');xlabel('频率/Hz');ylabel('功率谱/(dB)');%FileName:Forjietiao10K.m%clear;Fs=50000;%读入信道输出信号数据[recvsignal,Fs]=wavread('Chanel_OUT.wav');t=(1/Fs:1/Fs:length(recvsignal)/Fs)';%本地载波频率%fc_local=10000-200;%fc_local=10000-100;fc_local3=10000;%本地载波local_carrier3=cos(2*pi*fc_local3.*t);%相干解调xianggan_out3=recvsignal.*local_carrier3;%设计300〜3400Hz的带通滤波器[fenzi,fenmu]=butter(3,[3003400]/(Fs/2));demod_out3=filter(fenzi,fenmu,xianggan_out3);%sound(demod_out/max(demod_out),Fs);wavwrite(demod_out3,Fs,'SSBDemod_OUT_10k.wav');figure(4);subplot(1,2,1);plot(5*demod_out3(53550:53750));axis([0200-0.30.3]);title('解调信号波形');xlabel('时间(样值数)');subplot(1,2,2);psd(5*demod_out3,10000,Fs);axis([02500-2010]);title('解调信号功率谱');xlabel('频率/Hz');ylabel('功率谱/(dB)');%插值函数interp1的用法hours=l:12;temps=[58915252