电子科大移动通信原理课程设计报告
移动通信课程设计报告2篇
移动通信课程设计报告移动通信课程设计报告2篇移动通信课程设计报告2篇篇一：北京移动LTEBookRRU应用情况汇报北京移动LTE BookRRU 3235e应用情况汇报一、产品规格：BookRRU规格：体积/重量： 4L / 4.5kg 【体积小、重量轻，便于安装】频段： 2575~2635MHz【Band38/41，可三载波扩容】载波： 3*20 MHz发射功率： 2*10 W 【覆盖50~100米路段优选】天线增益： 10 dBi供电： AC 【交流供电、安装方便】天馈模式：内置 / 外置二、BooKRRU应用情况：北京移动城一分公司51、53网格应用，从测试情况来看，效果良好：RSRP、SINR等均明显提升。
【站点1】站名：朝阳新族大厦HL；覆盖路段：朝阳秀水南街BookRRU建设前该路段的覆盖情况：2016-9-2华为保密信息,未经授权禁止扩散第1页, 共5页RSRP： SINR:速率：BookRRU建设后该路段RSRP（覆盖图、做主服小区图）、SINR、速率情况:安装效果：RSRP: SINR:速率：2016-9-2华为保密信息,未经授权禁止扩散第2页, 共5页【站点2】站名：朝阳劲松七区713号楼HL、覆盖路段：劲松南路靠近2环路方向 BookRRU建设前该路段RSRP、SINR、速率情况RSRP: SINR:速率：BookRRU建设后该路段RSRP（覆盖图、做主服小区图）、SINR、速率情况安装效果：2016-9-2华为保密信息,未经授权禁止扩散第3页, 共5页RSRP: SINR:速率：【站点3】站名：朝阳垂杨柳南里；覆盖路段：广和街与朝阳垂杨柳中街交叉口处 BookRRU建设前该路段RSRP（覆盖图、做主服小区图）、SINR、速率情况 RSRP：SINR：2016-9-2华为保密信息,未经授权禁止扩散第4页, 共5页速率：BookRRU建设后该路段RSRP、SINR、速率情况： RSRP: SINR:速率：三、结论：1，BookRRU具有安装方便、隐蔽，不易激起民扰的特点；2，BookRRU对50~100米路段弱覆盖具有明显改善效果；3，BookRRU跟普通宏站之间切换顺畅，精准控制覆盖。
电子科大移动通信原理课程设计报告
移动通信原理课程设计报告一、题目描述仿真一：M=1,选定BPSK调制，AWGN和瑞利信道下的误符号率性能曲线（横坐标为符号信噪比Es/N0），并与相应的理论曲线比较。
仿真二：对2发1收的STBC-MIMO系统（Alamouti空时码），分析2发射天线分别受到独立瑞利信道下的误码率性能曲线，并与相同条件下单天线曲线进行对比分析。
二、系统设置三、仿真代码3.1算法说明1、信号产生：利用Matlab中的随机整数随机数产生函数randi.2、调制方法的实现：不同的调制方式对应唯一的一个星座图；通过输入序列找出星座图上的对应位置，即可输出调制结果。
3、信道模拟实现方法：AWGN信道用MATLAB自带函数randn实现，对应平均噪声功率为零；瑞利信道用randn+j*randn，对应平均噪声功率为零。
4、误码率性能曲线：发射信号序列长度设定130比特，仿真4000次,使信噪比在[0,30]每隔2取值，求平均误比特率。
5、收发系统的实现方法：对于单发单收的模型，只需将发送信号加噪声信号即为接收信号；对于二发一收的模型，因为发射天线是相互独立的，所以每根发射天线的接收信号与单发单收模型的接收信号计算方法相同，最后采用最大比合并得到接收信号。
6、调制方式：BPSK7、编码和译码方法：二发一收空时编码，最大似然译码。
8、误码率的计算：错误比特数/传输的总比特数。
3.2仿真代码代码一：调制函数function[mod_symbols,sym_table,M]=modulator(bitseq,b)N_bits=length(bitseq);if b==1 %BPSK调制sym_table=exp(1i*[0,-pi]);sym_table=sym_table([1 0]+1);inp=bitseq;mod_symbols=sym_table(inp+1);M=2;elseif b==2 %QPSK调制sym_table=exp(1i*pi/4*[-3 3 1 -1]);sym_table=sym_table([0 1 3 2]+1);inp=reshape(bitseq,b,N_bits/b);mod_symbols=sym_table([2 1]*inp+1);M=4;elseif b==3 %8PSK调制sym_table=exp(1i*pi/4*[0:7]);sym_table=sym_table([0 1 3 2 6 7 5 4]+1);inp=reshape(bitseq,b,N_bits/b);mod_symbols=sym_table([4 2 1]*inp+1);M=8;elseif b==4 %16QAM调制m=0;sq10=sqrt(10);for k=-3:2:3for l=-3:2:2m=m+1;sym_table(m)=(k+1i*l)/sq10;endendsym_table=sym_table(...[0 1 3 2 4 5 7 6 12 13 15 14 8 9 11 10]+1);inp=reshape(bitseq,b,N_bits/b);mod_symbols=sym_table([8 4 2 1]*inp+1);M=16;elseerror('unimplemented modulation');end代码二：单发单收系统在高斯信道和瑞利信道下的仿真clear allL_frame=130; N_packet=4000;b=1; % Set to 1/2/3/4 for BPSK/QPSK/8PSK/16QAMSNRdBs=[0:2:30]; sq2=sqrt(2);NT=1;NR=1;% SISOfor i_SNR=1:length(SNRdBs)SNRdB=SNRdBs(i_SNR);sigma=sqrt(0.5/(10^(SNRdB/10)));for i_packet=1:N_packetsymbol_data=randi([0 1],L_frame*b,NT);[temp,sym_tab,P]=modulator(symbol_data.',b);X=temp.‘; % frlg=length(X), X为调制后的信号序列Hr = (randn(L_frame,1)+1i*randn(L_frame,1))/sq2 ;%Rayleigh Channel，见原理说明6.1Ha=randn(L_frame,1); %AWGN channelZ1=0;R1=Hr.*X+ sigma*(randn(L_frame,1)+1i*randn(L_frame,1));Z1=Z1+R1.*conj(Hr);R2=X+sigma*Ha;for m=1:Pd1(:,m)=abs(Z1-sym_tab(m)).^2; %最大似然译码d2(:,m)=abs(R2-sym_tab(m)).^2;end[y1,i1] = min(d1,[],2);Xd=sym_tab(i1).';[y2,i2]=min(d2,[],2);Xa=sym_tab(i2).';temp1 = X>0; temp2 = Xd>0;temp3=Xa>0;noeb_p1(i_packet)=sum(sum(temp1~=temp2));noeb_p2(i_packet)=sum(sum(temp1~=temp3));endBER1(i_SNR)=sum(noeb_p1)/(N_packet*L_frame*b);BER2(i_SNR)=sum(noeb_p2)/(N_packet*L_frame*b);SNRw=10^(SNRdB/10);BER3(i_SNR)=1/2*erfc(sqrt(SNRw));%theoretical BER in AWGN channelBER4(i_SNR)=1/2*(1-sqrt(SNRw/(1+SNRw)));endsemilogy(SNRdBs,BER1,'-rx'), hold on, axis([SNRdBs([1 end]) 1e-6 1e0])semilogy(SNRdBs,BER2,'-ro'), hold on, axis([SNRdBs([1 end]) 1e-6 1e0])semilogy(SNRdBs,BER3,'-^'), hold on, axis([SNRdBs([1 end]) 1e-6 1e0])semilogy(SNRdBs,BER4,'-g*'), hold on, axis([SNRdBs([1 end]) 1e-6 1e0])title(‘BER perfoemancde of AWGN and Rayleigh channel'), xlabel('SNR[dB]'), ylabel('BER')grid on, set(gca,'fontsize',9)legend('Rayleigh practice','AWGN practice','AWGN theoretical','Rayleigh theoretical')代码三：二发一收系统空时编码仿真算法流程图：%Alamounti_scheme.mclear;N_frame=130;N_packets=4000;NT=2;NR=1; %two transmitter and one reciever diversityb=2;SNRdbs=[0:2:30];sq_NT=sqrt(NT);sq2=sqrt(2);for i_SNR=1:length(SNRdbs)SNRdb=SNRdbs(i_SNR);sigma=sqrt(0.5/(10^(SNRdb/10)));for i_packet=1:N_packetsmsg_symbol=randi([0 1],N_frame*b,NT);tx_bits=msg_symbol.';tmp=[];tmp1=[];for i=1:NT[tmp1,sym_tab,P]=modulator(tx_bits(i,:),b);tmp=[tmp;tmp1];endX=tmp.';%pay attention to the matrix dimension%space-time codingX1=X;X2=[-conj(X(:,2)) conj(X(:,1))];%channel known by receiverfor n=1:NTHr(n,:,:)=(randn(N_frame,NT)+1i*randn(N_frame,NT))/sq2;endH=reshape(Hr(n,:,:),N_frame,NT);%receive signalR1=sum(H.*X1,2)/sq_NT+sigma*(randn(N_frame,1)+1i*randn(N_frame,1));R2=sum(H.*X2,2)/sq_NT+sigma*(randn(N_frame,1)+1i*randn(N_frame,1));%MLD decoder，见原理6.2说明Z1=R1.*conj(H(:,1))+conj(R2).*H(:,2);Z2=R1.*conj(H(:,2))-conj(R2).*H(:,1);for m=1:Pd1(:,m)=abs(sum(Z1,2)-sym_tab(m)).^2;d2(:,m)=abs(sum(Z2,2)-sym_tab(m)).^2;end[y1,i1]=min(d1,[],2);S1d=sym_tab(i1).';clear d1[y2,i2]=min(d2,[],2);S2d=sym_tab(i2).';clear d2Xd=[S1d S2d];tmp1=X>0;tmp2=Xd>0;noeb_p(i_packet)=sum(sum(tmp1~=tmp2));endBER(i_SNR)=sum(noeb_p)/(N_packets*N_frame*b);endsemilogy(SNRdbs,BER,'-^');axis([SNRdbs([1 end]) 1e-6 1e0])grid on, hold onxlabel('SNR[db]'),ylabel('BER');四、仿真图图（1）图（2）五、仿真分析5.1题目一：由图（1）AWGN 和瑞利信道下的误符号率性能曲线与相应的理论曲线基本重合，且AWGN 性能强于瑞利信道。
移动通信工程课程设计报告
目录第一章设计目的 (1)第二章设计要求和设计指标 (2)2.1设计要求 (2)2.2 设计指标 (2)第三章设计内容 (3)3.1 语音测试—TEMS测试 (3)3.1.1 TEMS语音测试前准备 (4)3.1.2 进行测试 (4)3。
1.3测试数据的回放和处理 (5)3.1。
4 测试文件导出 (6)3。
1。
5测试数据统计 (7)3.2 优化调整方案 (8)第四章本课程设计改进和建议 (9)第五章总结 (10)参考文献 (11)附录 (12)第一章设计目的本课程设计练习移动通信的一般原理与组网技术,是一门实用性很强的课程。
设置本课程的目的是通过本课程设计之后,对移动通信的基本概念、基本原理和组网技术有较全面的了解和领会，应能应用移动通信的原理与技术分析阐释常见移动通信方式中信息传输的发送与接收原理，应能分析设计一些简单移动通信系统,为移动通信系统的管理维护、研究和开发打下必要的理论基础和技能.移动通信技术迅猛发展，移动通信网络也从第1代逐渐演进到了第4代（4G）。
国内各移动通信网络基础运营企业目前都已提供第3代(3G）的移动通信网络，并且移动通信应用规模在世界上位于前列，通过本课程设计要了解GSM无线网络的原理和功能，无线网络优化的原理，无线网络结构调整及频率优化，网络的测试方法及其指标，要通过数据证明本课题所得出的理论分析、网络结构调整和GSM无线参数的提取的正确性，并进一步说明了网络优化工作对各大运营商的重要性。
第二章设计要求和设计指标2.1设计要求在对数据进行详细采集、分析和研究后，常常会涉及到天馈系统的调整、基站的调测、频率规划的调整、系统参数的调整、话务均衡以及增加一些微蜂窝等优化方案实施活动.(1）对移动通信网络熟悉。
(2）要求优化方案正确合理。
(3）方案能解决一些问题。
（4)设计报告内容的正确性、全面性、逻辑性等。
2.2 设计指标对正式投入运行的GSM网络进行参数采集、数据分析、找出影响网络运行质量的原因，并且通过参数调整或采取某些技术手段使网络达到最佳运行状态,使现有网络资源获取最佳效益，同时也对GSM网络今后的维护及规划建设提出合理化建议。
移动通信课程设计报告
移动通信社会调查报告GSM网络优化班级电1005-1班姓名赵聪蕾学号 20102571一、调查目的学习移动通信这门课程之后，对移动通信的基本概念、基本原理和组网技术有较全面的了解和领会，应能应用移动通信的原理与技术分析阐释常见移动通信方式中信息传输的发送与接收原理，应能分析设计一些简单移动通信系统，为移动通信系统的管理维护、研究和开发打下必要的理论基础和实际技能。
我选择了GSM网络优化这个方面进行了调查，对正式投入运行的GSM网络进行参数采集、数据分析、找出影响网络运行质量的原因，并且通过参数调整或采取某些技术手段使网络达到最佳运行状态，使现有网络资源获取最佳效益，同时也对GSM网络今后的维护及规划建设提出合理化建议。
二、GSM网络基本原理1.GSM系统结构GSM（Global System for Mobile Communications；全球移动通信系统）主要分交换部分和无线部分。
其中交换部分和PSTN网很类似，而无线部分是GSM网络特有的由于无线特有的移动行，复杂性，以及传播条件恶劣所带来的衰落等原因，直接影响了无线通信的质量，所以无线部分是优化的重点对象。
一套完整的GSM蜂窝系统主要由：MS（移动台），BSS（基站子系统），NSS（交换网络子系统），OSS（操作支持子系统），这四大部分组成，GSM系统结构如图1-1所示。
图1-1 GSM 系统结构2.GSM网络组成分为交换系统（SS）和基站系统（BSS）。
另外，所有对网络的维护操作管理（OMC）是通过网管设备来完成的。
2.1 交换系统基本组成：MSC：Mobile services Switching Center，移动业务交换中心。
负责呼叫建立（也包括鉴权程序，呼叫控制，监视和计费。
短信发送。
GMSC：Gateway MSC，关口MSC。
主要用为移动网络和其他网络的接口局。
VLR：Visitor Location Register，拜访位置寄存器。
《移动通信原理与技术》实验报告
《移动通信原理与技术》实验报告实验一：TD-LTE硬件配置（1）实验名称TD-LTE硬件配置（2）实验目的1、熟练掌握移动通信系统的工作过程和工作原理，在移动通信实验教学中认识和了解通信网络和设备。
2、使用模拟现网的TD-LTE硬件平台和维护操作网络管理平台，使学生了解和掌握无线网络设备之中各个网元设备的工作配置原理，熟练掌握无线网络信令流程，理解无线网络对接数据的含义特征，提高学生对现网设备的安装、维护能力，提高学生对无线网络的开局能力。
（3）实验器材客户端、服务端、CCS2000U用户端程序、ZXSDR B8200 TL200设备物理接口、协议接口。
（4）实验原理实验原理图：ZXSDR B8200 TL200是一款支持多频段、多制式的基带单元，可同时支持GSM、UMTS及LTE等多种制式。
仅需进行软件配置和少量的硬件改动，即可将ZXSDR B8200 TL200配置为GERAN基站、UTRAN基站、LTE基站或者GUL多模基站。
ZXSDR B8200 TL200的软件结构分为SDR平台软件层、LTE适应软件层和LTE应用层。
SDR平台软件层：主要实现BSP、OSS和BRS的功能。
LTE 适应软件层：主要实现OAM和DBS的功能。
LTE应用层：实现LTE协议功能，包括控制面子系统、用户面子系统、调度器子系统、基带处理子系统等功能模块。
通过数据配置完成对两个E-UTRAN TDD小区的建立互通。
（5）实验方法1、进入WIN 2008操作系统。
数据配置前，首先打开网管服务器；2、创建子网，填写相关信息；3、创建网元，填写相关信息；4、运营商配置，填写相关运营商信息；5、填写PLMN信息，添加BBU侧设备（说明：各单板放置的位置要和实验室机柜中所用的BBU一致）；6、配置RRU，在机架图上点击图标添加RRU机架和单板，右键设备，点击添加RRU，会弹出RRU类型选择框，选中类型即可。
由于有2个RRU故需要增加2次（说明：RRU的类型必须与实际的硬件设备保持一致）；7.、时钟配置（默认配置即可）；8、光纤配置，是配置光接口板和RRU的拓扑关系（说明：2个RRU需要增加2条光纤）；9、物理层端口配置（说明：以太网方式配置参数直接手动改成1000）；10、以太网链路层配置；11、IP层配置；12、带宽配置；13、SCTP配置；14、业务与DSCP映射配置；15、静态路由配置；16、OMCB通道配置；17、创建无线网络；18、配置基带资源(说明:此处要配置2条基带资源，两次配置基带资源，若参考功率超出范围，要降低，不能都设为19.9）；19、S1AP配置；20、E-UTRAN TDD小区配置（说明：本网元有2个射频单元（2个RRU），需要再创建一个小区）；21、数据配置完成；22、测试网管与BBU是否建立连接，数据同步；23、验证数据配置是否正确，小区是否起来。
大专移动通信原理课程设计
大专移动通信原理课程设计一、课程目标知识目标：1. 理解并掌握移动通信的基本原理，包括信号传输、调制解调、多址技术等；2. 学习并了解移动通信系统的组成、结构及其工作原理；3. 掌握移动通信网络规划与优化的基本方法。
技能目标：1. 能够分析移动通信系统中信号传输的整个过程，并进行简单的计算；2. 学会使用相关软件工具进行移动通信网络的模拟与优化；3. 能够运用所学知识解决实际移动通信中的问题，具备一定的故障排查和优化能力。
情感态度价值观目标：1. 培养学生对移动通信技术的兴趣，激发学习热情，形成主动探索的精神；2. 增强学生的团队协作意识，培养良好的沟通能力；3. 使学生认识到移动通信技术在我国社会发展中的重要作用，培养学生的社会责任感和使命感。
课程性质：本课程为大专院校移动通信专业的基础课程，旨在让学生掌握移动通信的基本原理和实际应用。
学生特点：学生具备一定的电子技术和通信原理基础知识，但对移动通信系统的整体认识有限。
教学要求：结合学生特点和课程性质，注重理论与实践相结合，强调实际操作能力的培养，提高学生解决实际问题的能力。
在教学过程中，将课程目标分解为具体的学习成果，以便于后续的教学设计和评估。
二、教学内容1. 移动通信概述：介绍移动通信的发展历程、系统分类及未来发展趋势；- 教材章节：第一章 移动通信概述2. 信号传输原理：学习信号传输的基础知识，包括信号、信道、噪声等；- 教材章节：第二章 信号传输原理3. 调制解调技术：掌握各种调制解调技术及其在移动通信中的应用；- 教材章节：第三章 调制解调技术4. 多址技术：学习FDMA、TDMA、CDMA等多址技术的原理及其在移动通信系统中的应用；- 教材章节：第四章 多址技术5. 移动通信系统：介绍GSM、CDMA、TD-SCDMA、LTE等移动通信系统的结构、原理及关键技术；- 教材章节：第五章 移动通信系统6. 移动通信网络规划与优化：学习网络规划、基站选址、覆盖优化等方面的知识；- 教材章节：第六章 移动通信网络规划与优化7. 实践环节：结合理论教学，安排相应的实验和实践操作，提高学生的实际操作能力；- 教材章节：第七章 实践环节教学内容安排和进度：按照教材章节顺序，依次进行教学，共计14周。
移动通信课程设计
目录摘要 (1)一、绪论 (1)二、课程设计名称 (1)三、课程设计时间 (1)四、课程设计环境 (1)五、课程设计任务和要求 (2)5.1课程设计任务 (2)六、课程设计原理 (2)6.1 交织编码 (2)6.2 去交织硬件模型 (3)6.3去交织基本原理 (4)6.4单片机CPU2电路 (5)6.5去交织软件实现 (5)七、课程设计过程及调试、结果 (9)7.1 实验箱设置 (9)7.2 实时仿真方式开发程序 (9)7.3在系统编程（ISP）方式开发程序 (10)7.4 实验结果图 (11)八、课程设计体会 (15)参考文献 (16)附录（源程序） (17)摘要由于数字通信系统传输的是一个接一个按节拍传送的数字信号单元，即码元，因而在接收端必须按与发送端相同的节拍进行接收，否则，会因收发节拍不一致而导致接收性能变差。
此外，为了表述消息的内容，基带信号都是按消息内容进行编组的，因此，编组的规律在收发之间也必须一致。
在数字通信中，称节拍一致为“位同步”，称编组一致为“帧同步”。
在时分复用通信体统中，为了正确地传输信息，必须在信息码流中插入一定数量的帧同步码，它可以是一组特定的码组，也可以是特定宽度的脉冲，可以集中插入，也可以分散插入。
CDMA移动通信系统收端帧同步提取，我采用交织编码的方法。
交织编码是在实际移动通信环境下改善移动通信信号衰落的一种通信技术。
将造成数字信号传输的突发性差错，利用交织编码技术可离散并纠正这种突发性差错，改善移动通信的传输特性。
关键字：帧同步提取，帧同步，交织编码1、绪论数字移动通信系统使用数字信号传送信息，为第二代移动通信系统。
交织编码的目的是把一个较长的突发差错离散成随机差错，再用纠正随机差错的编码（FEC）技术消除随机差错。
交织深度越大，则离散度越大，抗突发差错能力也就越强。
但交织深度越大，交织编码处理时间越长，从而造成数据传输时延增大，也就是说，交织编码是以时间为代价的。
移动通信原理课程设计_实验报告_
电子科技大学通信抗干扰技术国家级重点实验室实验报告课程名称移动通信原理实验内容无线信道特性分析；BPSK/QPSK通信链路搭建与误码性能分析；SIMO系统性能仿真分析课程教师胡苏1，必做题目1.1无线信道特性分析1.1.1实验目的1)了解无线信道各种衰落特性；2)掌握各种描述无线信道特性参数的物理意义；3)利用MATLAB中的仿真工具模拟无线信道的衰落特性。
1.1.2实验内容1)基于simulink搭建一个QPSK发送链路，QPSK调制信号经过了瑞利衰落信道，观察信号经过衰落前后的星座图，观察信道特性。
仿真参数：信源比特速率为500kbps，多径相对时延为[0 4e-06 8e-06 1.2e-05]秒，相对平均功率为[0 -3 -6 -9]dB，最大多普勒频移为200Hz。
例如信道设置如下图所示：1.1.3实验仿真(1)实验框图（2）图表及说明图一：Before Rayleigh Fading1 #上图为QPSK相位图，由图可以看出2比特码元有四种。
图二：After Rayleigh Fading#从上图可以看出，信号通过瑞利信道后，满足瑞利分布，相位和幅度发生随机变化，所以图三中的相位不是集中在四点，而是在四个点附近随机分布。
图三：Impulse Response#从冲激响应的图可以看出相位在时间上发生了偏移。
图四：Impulse Response#从频率响应的图可以看出，信号的频率响应失真比较严重。
（3）实验结论根据题目中给出的参数，计算瑞利衰落信道的相干带宽和相干时间：相干带宽 410*2787.421==τπσc B Hz相干时间 005.01==mc f T s1.2 BPSK/QPSK 通信链路搭建与误码性能分析1.2.1 实验目的掌握基于simulink 的BPSK 、QPSK 典型通信系统的链路实现，仿真BPSK/QPSK 信号在AWGN 信道、单径瑞利衰落信道下的误码性能。
1.2.2 实验作业1.基于simulink 搭建BPSK/QPSK 通信链路，经过AWGN 信道，接收端相干解调，仿真并绘出BPSK 和QPSK 信号在0b E N 为0~10dB 时（间隔：1dB ）误码性能曲线。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
移动通信原理课程设计报告一、题目描述仿真一：M=1,选定BPSK调制，AWGN和瑞利信道下的误符号率性能曲线（横坐标为符号信噪比Es/N0），并与相应的理论曲线比较。
仿真二：对2发1收的STBC-MIMO系统（Alamouti空时码），分析2发射天线分别受到独立瑞利信道下的误码率性能曲线，并与相同条件下单天线曲线进行对比分析。
二、系统设置三、仿真代码3.1算法说明1、信号产生：利用Matlab中的随机整数随机数产生函数randi.2、调制方法的实现：不同的调制方式对应唯一的一个星座图；通过输入序列找出星座图上的对应位置，即可输出调制结果。
3、信道模拟实现方法：AWGN信道用MATLAB自带函数randn实现，对应平均噪声功率为零；瑞利信道用randn+j*randn，对应平均噪声功率为零。
4、误码率性能曲线：发射信号序列长度设定130比特，仿真4000次,使信噪比在[0,30]每隔2取值，求平均误比特率。
5、收发系统的实现方法：对于单发单收的模型，只需将发送信号加噪声信号即为接收信号；对于二发一收的模型，因为发射天线是相互独立的，所以每根发射天线的接收信号与单发单收模型的接收信号计算方法相同，最后采用最大比合并得到接收信号。
6、调制方式：BPSK7、编码和译码方法：二发一收空时编码，最大似然译码。
8、误码率的计算：错误比特数/传输的总比特数。
3.2仿真代码代码一：调制函数function[mod_symbols,sym_table,M]=modulator(bitseq,b)N_bits=length(bitseq);if b==1 %BPSK调制sym_table=exp(1i*[0,-pi]);sym_table=sym_table([1 0]+1);inp=bitseq;mod_symbols=sym_table(inp+1);M=2;elseif b==2 %QPSK调制sym_table=exp(1i*pi/4*[-3 3 1 -1]);sym_table=sym_table([0 1 3 2]+1);inp=reshape(bitseq,b,N_bits/b);mod_symbols=sym_table([2 1]*inp+1);M=4;elseif b==3 %8PSK调制sym_table=exp(1i*pi/4*[0:7]);sym_table=sym_table([0 1 3 2 6 7 5 4]+1);inp=reshape(bitseq,b,N_bits/b);mod_symbols=sym_table([4 2 1]*inp+1);M=8;elseif b==4 %16QAM调制m=0;sq10=sqrt(10);for k=-3:2:3for l=-3:2:2m=m+1;sym_table(m)=(k+1i*l)/sq10;endendsym_table=sym_table(...[0 1 3 2 4 5 7 6 12 13 15 14 8 9 11 10]+1);inp=reshape(bitseq,b,N_bits/b);mod_symbols=sym_table([8 4 2 1]*inp+1);M=16;elseerror('unimplemented modulation');end代码二：单发单收系统在高斯信道和瑞利信道下的仿真clear allL_frame=130; N_packet=4000;b=1; % Set to 1/2/3/4 for BPSK/QPSK/8PSK/16QAMSNRdBs=[0:2:30]; sq2=sqrt(2);NT=1;NR=1;% SISOfor i_SNR=1:length(SNRdBs)SNRdB=SNRdBs(i_SNR);sigma=sqrt(0.5/(10^(SNRdB/10)));for i_packet=1:N_packetsymbol_data=randi([0 1],L_frame*b,NT);[temp,sym_tab,P]=modulator(symbol_data.',b);X=temp.‘; % frlg=length(X), X为调制后的信号序列 Hr = (randn(L_frame,1)+1i*randn(L_frame,1))/sq2 ;%Rayleigh Channel，见原理说明6.1Ha=randn(L_frame,1); %AWGN channelZ1=0;R1=Hr.*X+ sigma*(randn(L_frame,1)+1i*randn(L_frame,1));Z1=Z1+R1.*conj(Hr);R2=X+sigma*Ha;for m=1:Pd1(:,m)=abs(Z1-sym_tab(m)).^2; %最大似然译码d2(:,m)=abs(R2-sym_tab(m)).^2;end[y1,i1] = min(d1,[],2);Xd=sym_tab(i1).';[y2,i2]=min(d2,[],2);Xa=sym_tab(i2).';temp1 = X>0; temp2 = Xd>0;temp3=Xa>0;noeb_p1(i_packet)=sum(sum(temp1~=temp2));noeb_p2(i_packet)=sum(sum(temp1~=temp3));endBER1(i_SNR)=sum(noeb_p1)/(N_packet*L_frame*b);BER2(i_SNR)=sum(noeb_p2)/(N_packet*L_frame*b);SNRw=10^(SNRdB/10);BER3(i_SNR)=1/2*erfc(sqrt(SNRw));%theoretical BER in AWGN channel BER4(i_SNR)=1/2*(1-sqrt(SNRw/(1+SNRw)));endsemilogy(SNRdBs,BER1,'-rx'), hold on, axis([SNRdBs([1 end]) 1e-6 1e0]) semilogy(SNRdBs,BER2,'-ro'), hold on, axis([SNRdBs([1 end]) 1e-6 1e0]) semilogy(SNRdBs,BER3,'-^'), hold on, axis([SNRdBs([1 end]) 1e-6 1e0]) semilogy(SNRdBs,BER4,'-g*'), hold on, axis([SNRdBs([1 end]) 1e-6 1e0]) title(‘BER perfoemancde of AWGN and Rayleigh channel'),xlabel('SNR[dB]'), ylabel('BER')grid on, set(gca,'fontsize',9)legend('Rayleigh practice','AWGN practice','AWGNtheoretical','Rayleigh theoretical')代码三：二发一收系统空时编码仿真算法流程图：%Alamounti_scheme.mclear;N_frame=130;N_packets=4000;NT=2;NR=1; %two transmitter and one reciever diversityb=2;SNRdbs=[0:2:30];sq_NT=sqrt(NT);sq2=sqrt(2);for i_SNR=1:length(SNRdbs)SNRdb=SNRdbs(i_SNR);sigma=sqrt(0.5/(10^(SNRdb/10)));for i_packet=1:N_packetsmsg_symbol=randi([0 1],N_frame*b,NT);tx_bits=msg_symbol.';tmp=[];tmp1=[];for i=1:NT[tmp1,sym_tab,P]=modulator(tx_bits(i,:),b);tmp=[tmp;tmp1];endX=tmp.';%pay attention to the matrix dimension%space-time codingX1=X;X2=[-conj(X(:,2)) conj(X(:,1))];%channel known by receiverfor n=1:NTHr(n,:,:)=(randn(N_frame,NT)+1i*randn(N_frame,NT))/sq2;endH=reshape(Hr(n,:,:),N_frame,NT);%receive signalR1=sum(H.*X1,2)/sq_NT+sigma*(randn(N_frame,1)+1i*randn(N_frame,1));R2=sum(H.*X2,2)/sq_NT+sigma*(randn(N_frame,1)+1i*randn(N_frame,1));%MLD decoder，见原理6.2说明Z1=R1.*conj(H(:,1))+conj(R2).*H(:,2);Z2=R1.*conj(H(:,2))-conj(R2).*H(:,1);for m=1:Pd1(:,m)=abs(sum(Z1,2)-sym_tab(m)).^2;d2(:,m)=abs(sum(Z2,2)-sym_tab(m)).^2;end[y1,i1]=min(d1,[],2);S1d=sym_tab(i1).';clear d1[y2,i2]=min(d2,[],2);S2d=sym_tab(i2).';clear d2Xd=[S1d S2d];tmp1=X>0;tmp2=Xd>0;noeb_p(i_packet)=sum(sum(tmp1~=tmp2));endBER(i_SNR)=sum(noeb_p)/(N_packets*N_frame*b); endsemilogy(SNRdbs,BER,'-^');axis([SNRdbs([1 end]) 1e-6 1e0]) grid on, hold onxlabel('SNR[db]'),ylabel('BER');四、仿真图图（1）图（2）五、仿真分析5.1题目一：由图（1）AWGN 和瑞利信道下的误符号率性能曲线与相应的理论曲线基本重合，且AWGN 性能强于瑞利信道。