4PSK.rar_4PSK MATLAB_4PSK解调_4psk调制simulink_MATLAB 4PSK_PSK

function [f,Ep,p]=cm_sm32(snr_in_dB) E=1; SNR=10^(snr_in_dB/10); %信噪比 sgma=sqrt(E/SNR)/2; %噪声的标准偏差 M=20000; %码元数 L=32; %每个码元采样点数 N=M*L*2; %采样点数 Rb=2; %码速率 Ts=1/(2*Rb); %码元间隔 dt=Ts/L; %时域采样间隔 df=1/(N*dt); %频域采样间隔 d=16; %半个码元采样点数 T=N*dt; %截短时宽 Bs=N*df/2; %频带宽度 t=[-T/2+dt/2:dt:T/2]; %时域坐标 t1=[-T/2+dt/2:32*dt:T/2]; %时域坐标 f=[-Bs+df/2:df:Bs]; %频域坐标 alpha=0.5; %滚降系数 g1=sin(pi*t/Ts)./(pi*t/Ts); g2=cos(alpha*pi*t/Ts)./(1-(2*alpha*t/Ts).^2); g=g1.*g2; G=t2f(g,dt); %发送和接收滤波器频谱 Cf=1-0.5*exp(-2*i*pi*f*d); %信道传输函数 Gf=G.*Cf; %系统总的频谱 %求系统后的码元样本序列 aa=zeros(1,N); for ii=1:M temp=rand; if(temp<0.25) a1_temp(ii)=0;a2_temp(ii)=0;a1(ii)=1;a2(ii)=0; elseif(temp<0.5) a1_temp(ii)=0;a2_temp(ii)=1;a1(ii)=0;a2(ii)=1; elseif(temp<0.75) a1_temp(ii)=1;a2_temp(ii)=0;a1(ii)=0;a2(ii)=-1; else a1_temp(ii)=1;a2_temp(ii)=1;a1(ii)=-1;a2(ii)=0; end for jj=1:L %构造4PSK信号 aa(2*L*(ii-1)+jj)=a1(ii); aa(2*L*(ii-1)+jj+L)=a2(ii); end end AA=t2f(aa,dt); Ep=(AA.*Gf).*conj(AA.*Gf)/T; %信号的功率谱密度 S=AA.*Gf; s=f2t(S,dt); s1=zeros(1,N); %此段程序将经傅立叶变换得到的数字序列 s1(N/2+1:N)=s(1:N/2);%在中轴反转,如111-1-11,经变换变成-1-11111 s1(1:N/2)=s(N/2+1:N);%因此应把它反过来,这主要是由 s1=real(s1); %f2t和t2f函数引起的. for ii=1:M %对接收到的信号进行采样，采样偏差为Ts/4 ra1=s1(8:64:N); ra2=s1(40:64:N); end %随后为判决和误码概率的计算 numoferr=0; for ii=1:M, %随后为判决 n(1)=gngauss(sgma); r=[ra1(ii),ra2(ii)]+n; c00=dot(r,[1,0]); %采用最佳判决，判断接收到的信号 c01=dot(r,[0,1]); c10=dot(r,[-1,0]); c11=dot(r,[0,-1]); c_max=max([c00,c01,c10,c11]); if(c00==c_max) decis1=0;decis2=0; elseif(c01==c_max) decis1=0;decis2=1; elseif(c11==c_max) decis1=1;decis2=0; else decis1=1;decis2=1; end if(decis1~=a1_temp(ii)) numoferr=numoferr+1; %如果出错，计数器加1 end; if(decis2~=a2_temp(ii)) numoferr=numoferr+1; end end; p=numoferr/(2*M); %误差概率估算