clear
clc
%輸入?yún)?shù)變量的值
p?=?0.1;????%序列{+1-1}中+1出現(xiàn)的概率為p，－1出現(xiàn)的概率為1-p
signal_power?=?1;???%信號(hào)功率為1
noise_power?=?0.0;???%噪聲功率為0
noise_power_times?=?50;?%信噪比的次數(shù)
test_times?=?100;??%重復(fù)實(shí)驗(yàn)的次數(shù)
length?=?1000;??%序列的長(zhǎng)度
SNR_dB?=?zeros（noise_power_times1）;??%保留每次實(shí)驗(yàn)的SNR
BER_mean_MAP?=?zeros（noise_power_times1）;??%保留每次MAP實(shí)驗(yàn)BER的均值
BER_mean_ML?=?zeros（noise_power_times1）;??%保留每次ML實(shí)驗(yàn)BER的均值
BER_all_MAP?=?zeros（test_times1）;??%保留每次MAP實(shí)驗(yàn)的BER
BER_all_ML?=?zeros（test_times1）;??%保留每次ML實(shí)驗(yàn)的BER
%運(yùn)算過(guò)程
for?i?=?1:?noise_power_times????%信噪比的個(gè)數(shù)
????noise_power?=?noise_power+i/10;?%噪聲功率N
????noise_mean?=?0;?%高斯白噪聲的均值
????noise_var?=?noise_power*0.5;????%高斯白噪聲的方差????
????for?j?=?1:test_times????%進(jìn)行重復(fù)實(shí)驗(yàn)
????????%信號(hào)序列
????????signal?=?randsrc（1length[-1?1;1-p?p]）;???%randsrc（）產(chǎn)生-1和1的貝努利隨機(jī)序列，概率均為?0.5
????????%噪聲序列
????????noise?=?normrnd（noise_meannoise_var[1length]）;???%normrnd（）產(chǎn)生N（noise_meannoise_var）的高斯白噪聲
????????%加噪聲的信號(hào)序列
????????mixed_signal?=?signal?+?noise;
????????judged_signal_MAP?=?zeros（size（signal））;???%zeros（）初始化序列的值為0，size（）獲得矩陣的行數(shù)和列數(shù)
????????judged_signal_ML?=?zeros（size（signal））;???%zeros（）初始化序列的值為0，size（）獲得矩陣的行數(shù)和列數(shù)
????????for?k?=?1:length
????????????%最佳MAP判決的信號(hào)序列
????????????i