香农误码率与信噪比 matlab,画QPSK的信噪比与误码率之间的线性关系图

最近编写了一个程序，现在想看看信噪比OSNR和误码率之间的线性关系，并画出图，可是一直都不对，拜托大神们能给予帮助，帮修改一下程序。现附上所写的代码。

clear

Nfft=256;%fft步长

Nc=127;%有效载波

Ncp=Nfft/8;%循环前缀

Ns=Nfft+Ncp;%一个完整的OFDM符号长度

Nd=300;%ofdm符号的个数

Nts=200;%训练序列个数

Ntotal=Nts+Nd;%总符号个数

sr=40*10^6;%传输比特率

vr=sr/(Ns*2);%OFDM符号传输速率。QPSK是2bit

M=2;%QPSK为2bit

OSNR=5;%信噪比

%%%%%训练符号产生(和数据产生形式相同)%%%%%

train=floor(rand(1,M*Nc*Nts)*2);%训练序列产生

%%QPSK调制%%

train1=1-2*train;

qpsktrain=reshape(train1,2,length(train1)/2);

for i=1:length(train1)/2

qpsktrain1(i)=qpsktrain(1,i)+1i*qpsktrain(2,i);

end

trainqpsk=qpsktrain1;

trainqpsk1=reshape(trainqpsk,Nc,Nts);

train2=zeros(Nfft,Nts);

train2(2:Nfft/2,:)=trainqpsk