QPSK.rar_qpsk simulink

clear all; close all; %初始化 fc=0.01;%载波频率; num=20;%码元数 tnum=100;%码元长度 n=1:num; a=zeros(1,num); R=randi([0,1],1,num);%随机生成num个0,1 Base=[]; for i=1:1:num if(R(i)==0) base=-1*ones(1,tnum); else base=1*ones(1,tnum); end Base=[Base base]; end %基带信号 figure(1); subplot(411); t=1:1:num*tnum; plot(t,Base); axis([1 num*tnum -2 2]); set(gca,'Xtick',[0:200:2000],'Ytick',-1:1:1);%设置网格线坐标 xlabel('基带信号'); ylabel('幅值'); grid on; %串并转换 for i=1:1:num if(R(i)==0) R(i)=-1; else R(i)=1; end end R1=reshape(R,2,num/2); I=R1(1,:); Q=R1(2,:); %与载波相乘 I2=[]; for i=1:1:num/2 if(I(i)==-1) I1=-1*ones(1,tnum); else I1=1*ones(1,tnum); end I2=[I2,I1]; end Q2=[]; for i=1:1:num/2 if(Q(i)==-1) Q1=-1*ones(1,tnum); else Q1=1*ones(1,tnum); end Q2=[Q2,Q1]; end t=1:1:num*tnum/2; I3=I2.*sin(2*pi*fc*t); Q3=Q2.*cos(2*pi*fc*t); Qpsk=I3+Q3;%将两路信号相加得QPSK信号 %for i=1:1:num/2 %Qpsk(i*100)=NaN;%如何绘出不连续的点 %end figure(1); subplot(412); plot(t,Qpsk); axis([1 1000 -2 2]); set(gca,'Xtick',[0:100:1000],'Ytick',-1:1:1);%设置网格线坐标 xlabel('QPSK信号'); ylabel('幅值'); grid on; %加高斯白噪声wgn SNR=5;%单位为dB snr=10^(SNR/10);%将分贝化成十进制 Am=0.7; N0=((Am^2)/2)/snr;%噪声功率单位为W N0_db=10*log10(N0);%dBw N=wgn(1,num*tnum/2,N0_db); qp1=Qpsk+N; t=1:1:num*tnum/2;%采样频率 subplot(413); plot(t,qp1); axis([1 1000 -2 2]); set(gca,'Xtick',[0:100:1000],'Ytick',-1:1:1);%设置网格线坐标 xlabel('加噪后的QPSK信号'); ylabel('幅值'); grid on; %[B,A] = BUTTER(N,Wn,'high') ---用来设计高通滤波器 %[B,A] = BUTTER(N,Wn,'low') 用来设计低通滤波器 %[B,A] = BUTTER(N,Wn)--带通滤波器 %比如说你的采样频率Fs=1000Hz，设计一个8阶、通带为100-200Hz的带通滤波器： %[b,a]=butter(8,[0.2 0.4])=butter(8,[100/(1000/2) 200/(1000/2) ]) %y=filter(B,A,x) %采样频率为1; [b,a]=butter(1,[0.002 0.04]);%带通滤波器 qp2=filter(b,a,qp1); t=1:1:1000;%采样频率 subplot(414); plot(t,qp2); axis([1 1000 -2 2]); set(gca,'Xtick',[0:100:1000],'Ytick',-1:1:1);%设置网格线坐标 xlabel('带通后QPSK信号'); ylabel('幅值'); grid on; %与载波信号相乘 figure(2); subplot(321); t=1:1:num*tnum/2; I4=qp2.*sin(2*pi*fc*t); %I4=Qpsk.*sin(2*pi*fc*t); t=1:1:1000;%采样频率 plot(t,I4); axis([1 1000 -2 2]); set(gca,'Xtick',[0:100:1000],'Ytick',-1:1:1);%设置网格线坐标 xlabel('I信号与载波相乘'); ylabel('幅值'); grid on; subplot(322); t=1:1:num*tnum/2; Q4=qp2.*cos(2*pi*fc*t); plot(t,Q4); axis([1 1000 -2 2]); set(gca,'Xtick',[0:100:1000],'Ytick',-1:1:1);%设置网格线坐标 xlabel('Q信号与载波信号相乘'); ylabel('幅值'); grid on; %低通滤波 [b1,a1] = butter(1,0.01); I5=filter(b1,a1,I4); t=1:1:1000;%采样频率 subplot(323); plot(t,I5); axis([1 1000 -2 2]); set(gca,'Xtick',[0:100:1000],'Ytick',-1:1:1);%设置网格线坐标 xlabel('低通滤波后的I信号'); ylabel('幅值'); grid on; [b2,a2] = butter(1,0.01); Q5=filter(b2,a2,Q4); t=1:1:1000;%采样频率 subplot(324); plot(t,Q5); axis([1 1000 -2 2]); set(gca,'Xtick',[0:100:1000],'Ytick',-1:1:1);%设置网格线坐标 xlabel('低通滤波后的Q信号'); ylabel('幅值'); grid on; %抽样判决 %JT=zeros(1,num*tnum); i6=[]; I6=[]; for i=1:num/2 if(I5(i*tnum)>0) I6=ones(1,tnum); i6=[i6 I6]; else I6=-1*ones(1,tnum); i6=[i6 I6]; end end t=1:1:1000;%采样频率 %figure(2); subplot(325); plot(t,i6); axis([1 1000 -2 2]); set(gca,'Xtick',[0:100:1000],'Ytick',-1:1:1);%设置网格线坐标 xlabel('解调后的I信号'); ylabel('幅值'); grid on; q6=[]; Q6=[]; for i=1:num/2 if(Q5(i*tnum)>0) Q6=ones(1,tnum); q6=[q6 Q6]; else Q6=-1*ones(1,tnum); q6=[q6 Q6]; end end t=1:1:1000;%采样频率 subplot(326); plot(t,q6); axis([1 1000 -2 2]); set(gca,'Xtick',[0:100:1000],'Ytick',-1:1:1);%设置网格线坐标 xlabel('解调后的Q信号'); ylabel('幅值'); grid on; %解调后的基带信号 JT=[]; for i=1:1:num/2 JT=[JT i6(((i-1)*100+1):i*100)]; JT=[JT q6(((i-1)*100+1):i*100)]; end t=1:1:2000;%采样频率 figure(3); subplot(411); plot(t,JT); axis([1 2000 -2 2]); set(gca,'Xtick',[0:200:2000],'Ytick',-1:1:1);%设置网格线坐标 xlabel('解调后的基带信号'); ylabel('幅值'); grid on;