Ray_qdata=qdata1.*ray_qch; Ray_s=Ray_idata+Ray_qdata;
% 高斯信道
s1=awgn(s,SNR); s11=abs(fft(s1)); s111=s1-s;
Awgn_s=awgn(Ray_s,SNR); % 通过高斯信道再通过瑞利衰落信道
%----------------------------- QPSK 解调部分
% 解调部分（高斯信道） idata2=s1.*a; qdata2=s1.*b;
idata3=zeros(1,nb/2); qdata3=zeros(1,nb/2);
% 抽样判决的过程，与0作比较，data>=0,则置1，否则置0 for n=1:nb/2
if sum(idata2((n-1)/delta_T+1:n/delta_T))>=0 idata3(n)=1; else idata3(n)=0; end end end
% 为了显示星座图,将信号进行处理 idata4=zeros(1,nb/2); qdata4=zeros(1,nb/2); for n=1:nb/2
Awgn_ichsum(n)=sum(idata2((n-1)/delta_T+1:n/delta_T))*delta_T; if Awgn_ichsum(n)>=0 idata4(n)=1;
31
else idata4(n)=0; end
Awgn_qchsum(n)=sum(qdata2((n-1)/delta_T+1:n/delta_T))*delta_T; if Awgn_qchsum(n)>=0 qdata4(n)=1; else qdata4(n)=0; end end
% 将判决之后的数据存放进数组 demodata=zeros(1,nb); demodata(1:ml:(nb-1))=idata3; demodata(2:ml:nb)=qdata3;
%为了显示，将它变成波形信号（即传输一个1代表单位宽度的高电平） demodata1=zeros(1,nb/delta_T); for q=1:nb
demodata1((q-1)/delta_T+1:q/delta_T)=demodata(q); end % 累计误码数
% abs(demodata-data)求接收端和发射端 % 数据差的绝对值，累计之后就是误码个数 Awgn_num_BER=sum(abs(demodata-data))
% 解调部分（瑞利+高斯）
Ray_idata2=Ray_s.*a; Ray_qdata2=Ray_s.*b;
% 为了显示星座图,将信号进行处理 Ray_idata4=zeros(1,nb/2); Ray_qdata4=zeros(1,nb/2); for n=1:nb/2
Ray_ichsum(n)=sum(idata2((n-1)/delta_T+1:n/delta_T))*delta_T; if Ray_ichsum(n)>=0 Ray_idata4(n)=1; else Ray_idata4(n)=0; end
Ray_qchsum(n)=sum(qdata2((n-1)/delta_T+1:n/delta_T))*delta_T;
32
if Ray_qchsum(n)>=0 Ray_qdata4(n)=1; else Ray_qdata4(n)=0; end end
% 将判决之后的数据存放进数组 Ray_demodata=zeros(1,nb);
Ray_demodata(1:ml:(nb-1))=Ray_idata4; Ray_demodata(2:ml:nb)=Ray_qdata4;
%为了显示，将它变成波形信号（即传输一个1代表单位宽度的高电平） Ray_demodata1=zeros(1,nb/delta_T); for q=1:nb
Ray_demodata1((q-1)/delta_T+1:q/delta_T)=Ray_demodata(q); end % 累计误码数
% abs(demodata-data)求接收端和发射端 % 数据差的绝对值，累计之后就是误码个数 Ray_num_BER=sum(abs(Ray_demodata-data))
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% 误码率计算
%% 调用了cm_sm32();和cm_sm33(）函数 %%声明：函数声明在另外俩个M文件中 %%作用：cm_sm32()用于瑞利信道误码率的计算 %% cm_sm33()用于高斯信道误码率的计算
%% ecoh on/off 作用在于决定是否显示指令内容
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
SNRindB1=0:1:6; SNRindB2=0:0.1:6; % 瑞利衰落信道
for i=1:length(SNRindB1),
[pb,ps]=cm_sm32(SNRindB1(i)); % 比特误码率 smld_bit_ray_err_prb(i)=pb;
33
smld_symbol_ray_err_prb(i)=ps; disp([ps,pb]); echo off; end; % 高斯信道 echo on;
for i=1:length(SNRindB1),
[pb1,ps1]=cm_sm33(SNRindB1(i)); smld_bit_awgn_err_prb(i)=pb1; smld_symbol_awgn_err_prb(i)=ps1; disp([ps1,pb1]); echo off; end; % 理论曲线 echo on;
for i=1:length(SNRindB2),
SNR=exp(SNRindB2(i)*log(10)/10); theo_err_awgn_prb(i)=0.5*erfc(sqrt(SNR)); theo_err_ray_prb(i)=0.5*(1-1/sqrt(1+1/SNR)); echo off; end;
h = spectrum.welch;
%---------------------------- 输出显示部分
% 第一部分（理想） figure(1) subplot(3,2,1);
plot(data0),title('基带信号'); axis([0 20000 -2 2]); subplot(3,2,2);
psd(h,data1,'fs',fs),title('基带信号功率谱密度'); subplot(3,2,3); plot(s),title('调制信号'); axis([0 500 -3 3]);
34
subplot(3,2,4);
psd(h,s,'fs',fs),title('调制信号功率谱密度'); subplot(3,2,5);
plot(demodata1),title('解调输出'); axis([0 20000 -2 2]); subplot(3,2,6);
psd(h,demodata1,'fs',fs),title('解调输出功率谱密度');
% 通过高斯信道 figure(2) subplot(2,2,1);
plot(s1),title('调制信号(Awgn)'); axis([0 500 -5 5]); subplot(2,2,2);
psd(h,s1,'fs',fs),title('调制信号功率谱密度(Awgn)'); subplot(2,2,3);
plot(s111),title('高斯噪声曲线'); axis([0 2000 -5 5]); figure(3) for i=1:nb/2
plot(idata(i),qdata(i),'r+'),title('QPSK信号星座图（Awgn）');hold on; axis([-2 2 -2 2]);
plot(Awgn_ichsum(i),Awgn_qchsum(i),'*');hold on; legend('理论值（发射端）','实际值（接收端）'); end
%通过高斯信道再通过瑞利衰落信道 figure(4) subplot(2,2,1)
plot(Ray_s),title('调制信号(Ray+Awgn)'); axis([0 500 -5 5]); subplot(2,2,2);
psd(h,Ray_s,'fs',fs),title('调制信号功率谱密度(Ray)'); figure(5) for i=1:nb/2
plot(idata(i),qdata(i),'r+'),title('QPSK信号星座图（Awgn+Ray）');hold on; axis([-2 2 -2 2]);
35