qpsk瑞利matlab.doc

%%%%qpsk simulation实现QPSK瑞丽信道的各部分图 %***********************初始设置 ******************************** sr=2560; %码符号率，单位symbol/s ml=2; %调制等级，对于qpsk是2 br=sr.*ml; %比特率，单位bit/s nd=20; %码符号数目 ebn0=20; %比特能量，单位db IPIONT=8; %过采样数目 bit_time=1/sr/ml;%比特持续时间 %***********************升余弦滤波器及信道初始化设置***************************** irfn=21;%抽点数目 alfas=0.5; %滚降因子 [xh]=hrollfcoef(irfn,IPIONT,sr,alfas,1);%发送端根升余弦系统 [xh2]=hrollfcoef(irfn,IPIONT,sr,alfas,0);%接受端根升余弦系统 tstp=1/sr/IPIONT; %最小时间间隔 itau = [0]; dlvl = [0]; no=6;%产生单条瑞利信道的需要的不可分辨多径数 th1=[0.0]; itnd0=nd*IPIONT*100; itnd1=[1000]; now1=1; fd=160;%最大多普勒频移,单位HZ flat =1;%1表示仅考虑幅度，0表示幅度和相位都受到衰落 %***********************产生二进制比特流 ********************************* %% 产生二进制比特流 data=rand(1,nd*ml)0.5; %% qpsk调制 [ich,qch]=qpskmod(data,1,nd,ml); figure(1); subplot 311;stairs (data) ,title(码元);axis([0,nd*ml,-2,2]); subplot 312;stairs (ich) ,title(同向支路I);axis([0,nd,-2,2]); subplot 313;stairs (qch) ,title(正交支路Q);axis([0,nd,-2,2]); %**********************脉冲成型滤波器********************************* %% *脉冲成型滤波器 fs=sr*10; % 采样频率 [ich1,qch1]=comoversamp(ich,qch,length(ich),IPIONT);%过采样 data1=ich1+qch1; %% 根升余弦滤波 [ich2,qch2]=comconv(ich1,qch1,xh); data2=ich2+qch2; figure(2); subplot 311; plot(data2); title(发送端采用升余弦波形传输的原始数据);grid on; subplot 312; plot(ich2); title(发送端采用升余弦波形传输的同相支路);grid on; subplot 313; plot(qch2); title(发送端采用升余弦波形传输的正交支路);grid on; h = spectrum.welch; figure; psd(h,data2,fs,fs); title(基带升余弦波形传输信号功率谱密度); xlabel(频率 (MHz) ); %***********************qpsk 载波调制********************************** fc=40*sr; % 载波频率Hz 是符号速率的40倍 ich3=ich2.*cos(2*pi*fc*tstp*(1:length(ich2)));%乘以载波 qch3=qch2.*(-sin(2*pi*fc*tstp*(1:length(qch2)))); data3=ich3+qch3; figure(3); subplot 311;plot (data3) ,title(调制后波形); subplot 312;plot (ich3) ,title(调制后I路波形); subplot 313;plot (qch3) ,title(调制后Q路波形); h = spectrum.welch; figure(4); psd(h,data3,fs,fs); title(经过载波调制后接收端信号功率谱密度); xlabel(频率 (MHz) ); %***********************加性白噪声计算*****************