电子科大移动通信原理课程设计报告

《电子科大移动通信原理课程设计报告》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、移动通信原理课程设计报告一、题目描述仿真一：M=1,选定BPSK调制，AWGN和瑞利信道下的误符号率性能曲线（横坐标为符号信噪比Es/N0），并与相应的理论曲线比较。仿真二：对2发1收的STBC-MIMO系统（Alamouti空时码），分析2发射天线分别受到独立瑞利信道下的误码率性能曲线，并与相同条件下单天线曲线进行对比分析。二、系统设置系统设置如下表所示：调制方式收发天线数信道条件编码方式仿真一系统一BPSK一发一收高斯信道无系统二一发一收瑞利信道无仿真二系统三二发一收瑞利信道Alamouti空时码

2、三、仿真代码3.1算法说明1、信号产生：利用Matlab中的随机整数随机数产生函数randi.2、调制方法的实现：不同的调制方式对应唯一的一个星座图；通过输入序列找出星座图上的对应位置，即可输出调制结果。3、信道模拟实现方法：AWGN信道用MATLAB自带函数randn实现，对应平均噪声功率为零；瑞利信道用randn+j*randn，对应平均噪声功率为零。4、误码率性能曲线：发射信号序列长度设定130比特，仿真4000次,使信噪比在[0,30]每隔2取值，求平均误比特率。5、收发系统的实现方法：对于单

3、发单收的模型，只需将发送信号加噪声信号即为接收信号；对于二发一收的模型，因为发射天线是相互独立的，所以每根发射天线的接收信号与单发单收模型的接收信号计算方法相同，最后采用最大比合并得到接收信号。6、调制方式：BPSK7、编码和译码方法：二发一收空时编码，最大似然译码。8、误码率的计算：错误比特数/传输的总比特数。3.2仿真代码代码一：调制函数function[mod_symbols,sym_table,M]=modulator(bitseq,b)N_bits=length(bitseq);ifb==1

4、%BPSK调制sym_table=exp(1i*[0,-pi]);sym_table=sym_table([10]+1);inp=bitseq;mod_symbols=sym_table(inp+1);M=2;elseifb==2%QPSK调制sym_table=exp(1i*pi/4*[-331-1]);sym_table=sym_table([0132]+1);inp=reshape(bitseq,b,N_bits/b);mod_symbols=sym_table([21]*inp+1);M=4;

5、elseifb==3%8PSK调制sym_table=exp(1i*pi/4*[0:7]);sym_table=sym_table([01326754]+1);inp=reshape(bitseq,b,N_bits/b);mod_symbols=sym_table([421]*inp+1);M=8;elseifb==4%16QAM调制m=0;sq10=sqrt(10);fork=-3:2:3forl=-3:2:2m=m+1;sym_table(m)=(k+1i*l)/sq10;endendsym_ta

6、ble=sym_table(...[0132457612131514891110]+1);inp=reshape(bitseq,b,N_bits/b);mod_symbols=sym_table([8421]*inp+1);M=16;elseerror('unimplementedmodulation');end代码二：单发单收系统在高斯信道和瑞利信道下的仿真clearallL_frame=130;N_packet=4000;b=1;%Setto1/2/3/4forBPSK/QPSK/8PSK/16Q

7、AMSNRdBs=[0:2:30];sq2=sqrt(2);NT=1;NR=1;%SISOfori_SNR=1:length(SNRdBs)SNRdB=SNRdBs(i_SNR);sigma=sqrt(0.5/(10^(SNRdB/10)));fori_packet=1:N_packetsymbol_data=randi([01],L_frame*b,NT);[temp,sym_tab,P]=modulator(symbol_data.',b);X=temp.‘;%frlg=length(X),X为调

8、制后的信号序列Hr=(randn(L_frame,1)+1i*randn(L_frame,1))/sq2;%RayleighChannel，见原理说明6.1Ha=randn(L_frame,1);%AWGNchannelZ1=0;R1=Hr.*X+sigma*(randn(L_frame,1)+1i*randn(L_frame,1));Z1=Z1+R1.*conj(Hr);R2=X+sigma*Ha;form=1:Pd1(:,m)=abs(Z1-sym_