Question

MPSK传输方式的仿真分析

各位大侠大家好。
我现做的题目是《MPSK传输方式的仿真分析》，要求用MATLAB做出仿真结果。
我想请问一下，这个题目的论文应该从几个方面写，我现在有点无从下手。
众所周知，MPSK的M可取2^n，一边取2,4,8,16做分析。我写论文的时候是否是取其中的两项做下分析？
问题2：传输包括调制和解调2部分，写论文的时候是否应该分2部分写？
如果大家能帮我把正文提纲和仿真结果给我就更好了。
如果答案很不错的话可以追加分数（10~300）
请大家踊跃发言！

Answer 1

8psk的源程序：

function[pb,ps]=cm_sm32(snr_in_dB)

N=10000;
E=1;
snr=10^(snr_in_dB/10);
sgma=sqrt(E/(3*2*snr));

s000=[1 0];
s001=[cos(pi/4) sin(pi/4)];
s011=[0 1];
s010=[cos(3*pi/4) sin(3*pi/4)];
s110=[-1 0];
s111=[cos(5*pi/4) sin(5*pi/4)];
s101=[0 -1];
s100=[cos(7*pi/4) sin(7*pi/4)];

for i=1:N,
temp=rand;

if(temp<0.125),
dsource1(i)=0;
dsource2(i)=0;
dsource3(i)=0;
elseif(temp<0.25),
dsource1(i)=0;
dsource2(i)=0;
dsource3(i)=1;
elseif(temp<0.375),
dsource1(i)=0;
dsource2(i)=1;
dsource3(i)=0;
elseif(temp<0.5),
dsource1(i)=0;
dsource2(i)=1;
dsource3(i)=1;
elseif(temp<0.625),
dsource1(i)=1;
dsource2(i)=0;
dsource3(i)=0;
elseif(temp<0.75),
dsource1(i)=1;
dsource2(i)=0;
dsource3(i)=1;
elseif(temp<0.875),
dsource1(i)=1;
dsource2(i)=1;
dsource3(i)=0;
else
dsource1(i)=1;
dsource2(i)=1;
dsource3(i)=1;
end;
end;

numofsymbolerror=0;
numofbiterror=0;

for i=1:N,
n(1)=gngauss(sgma);
n(2)=gngauss(sgma);

if((dsource1(i)==0)&(dsource2(i)==0)&(dsource3(i)==0)),
r=s000+n;
elseif((dsource1(i)==0)&(dsource2(i)==0)&(dsource3(i)==1)),
r=s001+n;
elseif((dsource1(i)==0)&(dsource2(i)==1)&(dsource3(i)==0)),
r=s010+n;
elseif((dsource1(i)==0)&(dsource2(i)==1)&(dsource3(i)==1)),
r=s011+n;
elseif((dsource1(i)==1)&(dsource2(i)==0)&(dsource3(i)==0)),
r=s100+n;
elseif((dsource1(i)==1)&(dsource2(i)==0)&(dsource3(i)==1)),
r=s101+n;
elseif((dsource1(i)==1)&(dsource2(i)==1)&(dsource3(i)==0)),
r=s110+n;
else
r=s111+n;
end;

c000=dot(r,s000);
c001=dot(r,s001);
c010=dot(r,s010);
c011=dot(r,s011);
c100=dot(r,s100);
c101=dot(r,s101);
c110=dot(r,s110);
c111=dot(r,s111);

c_max=max([c000 c001 c010 c011 c100 c101 c110 c111]);

if(c000==c_max),
decis1=0;decis2=0;decis3=0;
elseif(c001==c_max),
decis1=0;decis2=0;decis3=1;
elseif(c010==c_max),
decis1=0;decis2=1;decis3=0;
elseif(c011==c_max),
decis1=0;decis2=1;decis3=1;
elseif(c100==c_max),
decis1=1;decis2=0;decis3=0;
elseif(c101==c_max),
decis1=1;decis2=0;decis3=1;
elseif(c110==c_max),
decis1=1;decis2=1;decis3=0;
else
decis1=1;decis2=1;decis3=1;
end;

symbolerror=0;
if(decis1~=dsource1(i)),
numofbiterror=numofbiterror+1;
symbolerror=1;
end;

if(decis2~=dsource2(i)),
numofbiterror=numofbiterror+1;
symbolerror=1;
end;

if(decis3~=dsource3(i)),
numofbiterror=numofbiterror+1;
symbolerror=1;
end;

if(symbolerror==1),
numofsymbolerror=numofsymbolerror+1;
end;
end;

ps=numofsymbolerror/N;
pb=numofbiterror/(3*N);

4PSK的源程序：

% 对M=4的PSK通信系统进行蒙特卡罗仿真 %
echo on
SNRindB1=0:2:10;
SNRindB2=0:0.1:10;
for i=1:length (SNRindB1),
[pb , ps]=cm_sm32 (SNRindB1 (i)); %仿真比特和符号误码率
smld_bit_err_prb (i)=pb;
smld_symbol_err_prb (i)=ps;
end;
for i=1:length (SNRindB2),
SNR=exp(SNRindB2(i)*log(10)/10); %信噪比
theo_err_prb(i)=Qfunct(sqrt(2*SNR)); %理论比特误码率
end;
%随后绘图曲线
semilogy (SNRindB1,smld_bit_err_prb,'*');
hold
semilogy (SNRindB1,smld_symbol_err_prb,'O');
semilogy (SNRindB2,theo_err_prb);
grid
xlabel('Eb/No in dB')
ylabel('error probability')
title('4PSK通信系统的蒙特卡洛仿真')
gtext('（注：“—”理论误码率；“*”误比特率；“o”误符号率）'