《电子科大移动通信原理课程设计报告》由会员分享，可在线阅读，更多相关《电子科大移动通信原理课程设计报告（6页珍藏版）》请在金锄头文库上搜索。
1、-移动通信原理课程设计报告一、题目描述仿真一：M=1,选定BPSK调制，AWGN和瑞利信道下的误符号率性能曲线（横坐标为符号信噪比Es/N0），并与相应的理论曲线比较。仿真二：对2发1收的STBC-MIMO系统（Alamouti空时码），分析2发射天线分别受到独立瑞利信道下的误码率性能曲线，并与相同条件下单天线曲线进行对比分析。二、系统设置系统设置如下表所示：调制方式收发天线数信道条件编码方式仿真一系统一BPSK一发一收高斯信道无系统二一发一收瑞利信道无仿真二系统三二发一收瑞利信道Alamouti空时码三、仿真代码3.1算法说明1、信号产生：利用Matlab中的随机整数随机数产生函数randi
2、.2、调制方法的实现：不同的调制方式对应唯一的一个星座图；通过输入序列找出星座图上的对应位置，即可输出调制结果。3、信道模拟实现方法：AWGN信道用MATLAB自带函数randn实现，对应平均噪声功率为零；瑞利信道用randn+j*randn，对应平均噪声功率为零。4、误码率性能曲线：发射信号序列长度设定130比特，仿真4000次,使信噪比在0,30每隔2取值，求平均误比特率。5、收发系统的实现方法：对于单发单收的模型，只需将发送信号加噪声信号即为接收信号；对于二发一收的模型，因为发射天线是相互独立的，所以每根发射天线的接收信号与单发单收模型的接收信号计算方法相同，最后采用最大比合并得到接收信
3、号。6、调制方式：BPSK7、编码和译码方法：二发一收空时编码，最大似然译码。8、误码率的计算：错误比特数/传输的总比特数。3.2仿真代码代码一：调制函数functionmod_symbols,sym_table,M=modulator(bitseq,b)N_bits=length(bitseq);if b=1 %BPSK调制 sym_table=e*p(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=e*p(1i*
4、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=e*p(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=
5、0; sq10=sqrt(10); for k=-3:2:3 for l=-3:2:2 m=m+1; sym_table(m)=(k+1i*l)/sq10; end endsym_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;else error(unimplemented modulation);end代码二：单发单收系统在高斯信道和瑞利信道下的仿真clear allL_f
6、rame=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_packet symbol_data=randi(0 1,L_frame*b,NT); temp,sym_tab,P=modulator(symbol_data.,b); *=te
7、mp.; % frlg=length(*), *为调制后的信号序列 Hr = (randn(L_frame,1)+1i*randn(L_frame,1)/sq2 ;%Rayleigh Channel，见原理说明6.1 Ha=randn(L_frame,1); %AWGN channelZ1=0; R1=Hr.*+ sigma*(randn(L_frame,1)+1i*randn(L_frame,1); Z1=Z1+R1.*conj(Hr); R2=*+sigma*Ha; for m=1:P d1(:,m)=abs(Z1-sym_tab(m).2; %最大似然译码 d2(:,m)=abs(R2-
8、sym_tab(m).2; end y1,i1 = min(d1,2); *d=sym_tab(i1).; y2,i2=min(d2,2); *a=sym_tab(i2).; temp1 = *0; temp2 = *d0; temp3=*a0; 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);
9、 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,-r*), hold on, a*is(SNRdBs(1 end) 1e-6 1e0)semilogy(SNRdBs,BER2,-ro), hold on, a*is(SNRdBs(1 end) 1e-6 1e0)semilogy(SNRdBs,BER3,-), hold on, a*is(SNRdBs(
10、1 end) 1e-6 1e0)semilogy(SNRdBs,BER4,-g*), hold on, a*is(SNRdBs(1 end) 1e-6 1e0)title(BER perfoemancde of AWGN and Rayleigh channel), *label(SNRdB), ylabel(BER) grid on, set(gca,fontsize,9)legend(Rayleigh practice,AWGN practice,AWGN theoretical,Rayleigh theoretical)代码三：二发一收系统空时编码仿真算法流程图：%Alamounti_s
11、cheme.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_packets msg_symbol=randi(0 1,N_frame*b,NT); t*_bits=msg_symbol.; t
12、mp=; tmp1=; for i=1:NT tmp1,sym_tab,P=modulator(t*_bits(i,:),b); tmp=tmp;tmp1; end *=tmp.;%pay attention to the matri* dimension %space-time coding *1=*; *2=-conj(*(:,2) conj(*(:,1); %channel known by receiver for n=1:NT Hr(n,:,:)=(randn(N_frame,NT)+1i*randn(N_frame,NT)/sq2; end H=reshape(Hr(n,:,:),
13、N_frame,NT); %receive signal R1=sum(H.*1,2)/sq_NT+sigma*(randn(N_frame,1)+1i*randn(N_frame,1); R2=sum(H.*2,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:P d1(:,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)
