%譯碼算法
clear?all
Nt=4;%發(fā)射天線數(shù)
Nr=1;%接收天線數(shù)
num_S=4;%每次發(fā)射的符號數(shù)量
num_slot=4;%時隙數(shù)
nit=1000;%采樣點(diǎn)數(shù)
M_qam=4;%調(diào)制方式
num_bit_per_sym=log2（M_qam）;%每個符號含有的比特數(shù)
snr_min=0;%最小信噪比
snr_max=20;%最大信噪比

graph_inf_bit=zeros（snr_max-snr_min+12）;
for?SNR=snr_min:1:snr_max
????clc
????disp（‘Wait?Until?SNR=‘）;%disp文本形式顯示矩陣
????disp（snr_max）;
????SNR??%顯示仿真進(jìn)度
????graph_inf_bit（SNR-snr_min+11）=SNR;%橫坐標(biāo)取值；當(dāng)前SNR存儲到graph_inf_bit的第一列中的SNR-snr_min+1位置
????num_error_sym=0;%初始化誤符號數(shù)
????num_error_bit=0;%初始化誤比特數(shù)
????for?num_sym=1:nit
????????s0=randint（num_Snum_bit_per_sym）;%產(chǎn)生隨機(jī)數(shù)列；rand隨機(jī)產(chǎn)生0到1之間的數(shù)，randint隨機(jī)產(chǎn)生0到1之間的整數(shù)，即0或1
????????s=zeros（num_S1）;
????????for?i=1:num_S
????????????if?s0（i:）==[0?0]
????????????????s（i1）=-1+j;
????????????elseif?s0（i:）==[0?1]?
???????????????s（i1）=-1-j;
????????????elseif?s0（i:）==[1?0]
???????????????s（i1）=1+j;
????????????else
???????????????s（i1）=1-j;???
????????????end
????????end??%調(diào)制過程；以上是4QAM映射，s0（i:）表示s0矩陣的第i行，所有列?
????????
????????h=（randn（NtNt）+j*randn（NtNt））/sqrt（2）;%路徑增益?滿足正態(tài)分布的偽隨機(jī)數(shù)4*1
????????de_snr=10^（SNR/10）;%將信噪比從db數(shù)換算成十進(jìn)制
????????sig=sqrt（Nt/de_snr）;
????????noise=sig*（randn（num_slotNr）+j*randn（num_slotNr））/sqrt（2）;%歸一化處理
????????r=h*s+noise;???%接收到的信號4*1
???????
%----------------------part?1-----MMSE檢測算法--------------

W_mmse=inv（h‘*h+eye（4）./de_snr）*h‘;%inv表示取逆
ss=W_mmse*r;
%Signal?demodulation
re_sym=zeros（num_Snum_bit_per_sym）;%存放譯碼結(jié)果，num_bit_per_sym每個符號含有的比特數(shù)

%----------------------part?2-----成對譯碼-------------
?qam_space_2=[-1+j
????????-1-j
????????1+j
????????1-j];
????for?i=1:num_S
????????for?ii=1:M_qam
????????????f_1_3=（abs（ss（i1）-qam_space_2））.^2;
????????end
????????[M1?N1]=min（f_1_3）;%M1為（f_1_3）的最小值，N1為對應(yīng)的s3
????????re_sym（i1）=qam_space_2（N1）
????end


????????for?i=1:length（re_sym）
?????????????if?re_sym（i）==-1+j
????????????????re_sym1（i:）=[0?0]?;???
????????????elseif?re_sym（i）==-1-j
????????????????re_sym1（i:）=[0?1]?;
????????????elseif?re_sym（i）==1-j
????????????????re_sym1（i:）=[1?1]?;
?????????????else?re_sym1（i:）=[1?0]?;
????????????end
????????end??????????????????????????????????
????????de_sym=bi2de（re_sym1）;%二進(jìn)制轉(zhuǎn)化為十進(jìn)制
????????de_s0=bi2de（s0）;
????????for?rp=1:num_S
????????????if?de_sym（rp1）~=de_s0（rp1）%發(fā)射矩陣與譯碼接收的矩陣比較
????????????????num_error_sym=num_error_sym+1;
????????????????for?mm=1:num_bit_per_sym
????????????????????if?re_sym1（rpmm）~=s0（rpmm）
????????????????????????num_error_bit=?num_error_bit+1;
????????????????????end
????????????????end
????????????end
????????end
????end
????per_bit=?num_error_bit/（num_S*num_bit_per_sym*nit）;%誤比特率的計算，誤比特率=錯誤的比特數(shù)/（符號長度*錯誤符號個數(shù)*符號數(shù)）
????graph_inf_bit（SNR-snr_min+12）=per_bit;
?end?
x=graph_inf_bit（:1）;
y=graph_inf_bit（:2）;
semilogy（xy‘k-o‘）;
grid?on
hold?on

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????3095??2017-12-02?14:00??MMSE.m

-----------?---------??----------?-----??----

?????????????????3095????????????????????1