繼續(xù)講解！之前文章里面講過(guò)一種等效方式對(duì)于是否采用滾降成型都實(shí)用！這么神奇？之前在QPSK調(diào)制解調(diào)器里面提及了，本文就來(lái)展示程序，在程序中讓大家感受道這個(gè)計(jì)算的細(xì)節(jié)。本文由于含有程序，長(zhǎng)達(dá)三千多字，請(qǐng)大家耐心看！

不管采用什么加噪方式，也不管你設(shè)置了多少dB的信噪比，按照這種方式計(jì)算出來(lái)的就是真實(shí)的信噪比值！！！怎么在程序中體現(xiàn)呢？那就趕快看下面的仿真程序吧！

%%% 仿真環(huán)境

% 軟件版本：R2021a

clear all;

%%****************** 程序主體 *******************%%

%%%%%%%%%% 發(fā)射機(jī)部分 %%%%%%%%%

cnr_number = 1; % 仿真循環(huán)次數(shù)

%－－－－－－－－－ 幀頭 20 位 －－－－－－－－

Iframe_head = [ 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1];

%－－－－－－－－ 信息內(nèi)容 150000 位 －－－－－－

indisplay1 = randi([0 1],1,150000);

% indisplay1 = ones(1,150000);

information1 = 2*indisplay1-1;

%－－－－－－－－起始內(nèi)容是BTR位 －－－－－－－－－－－

Ibit_flow = [Iframe_head information1];

% information Iframe_head information Iframe_head information];

%－－－－－－－－－－ 調(diào)制 －－－－－－－－－－－－－

len = length(Ibit_flow);

inf_phase = zeros(1,len);

inf_phase(1) = -pi/4;

%－－－－－－－－ BPSK －－－－－－－－－－－

Tra_baseI = Ibit_flow;

%－－－－－－－ 加入采樣頻偏 －－－－－－－

sr = 2.4e3;

offset = 0; % 不加采樣頻偏

ins_value = 8;

ideal_samplefre = ins_value*sr;

samplefre= ideal_samplefre+offset;

time1= ceil(sr/ideal_samplefre.*(1:length(Tra_baseI)*8));

%－－－－－ 生成以8*2.4e3采樣頻率采樣的數(shù)據(jù) －－－－－－－

Isignalsample = Tra_baseI(time1);

%－－－－－－－ 生成有載波調(diào)制的數(shù)據(jù)－－－－－－－－－－－

fcarrier = 3*sr;

foffset = 30; % Hz 頻偏值

ophase = pi/3;

time = (0:length(Isignalsample)-1)/samplefre;

tra_IFsignal = Isignalsample.*cos(2*pi*(fcarrier+foffset).*time+ophase);

power_tra_IFsignal = sum(tra_IFsignal.^2)/len*2;

% 計(jì)算信號(hào)功率 為什么要乘以2呢？

%%%%%%%%%%%%%%%% %%%%%%%%%%

%%%%%%%%% 信道部分 %%%%%%%%%

%%%%%%%%%%%%%%%%%%%% %%%%%%

for snr = 36

%－－－－經(jīng)過(guò)加性高斯白噪聲信道－－－

% spow=sum(tra_IFsignal.^2)/len; % sum: built in function

% attn=0.5*spow*1/2*10.^(-ebn0/10);

% attn=sqrt(attn); % sqrt: built in function

% inoise = randn(1,length(tra_IFsignal)).*attn;

IFsignal = awgn(tra_IFsignal, snr, 'measured'); %設(shè)置載噪比

power_noise = sum((IFsignal-tra_IFsignal).^2)/len;

審核編輯 ：李倩