基于matlab的qpsk與bpsk信號(hào)性能比較仿真-資料下載頁(yè)
2025-08-10 16:23本頁(yè)面
【正文】 pesystem(m)=pesystem(m)+pe2。 endend peideal=(1/2)*peideal./nx。pesystem=(1/2)*pesystem./nx。 returnfunction QPSKThreeray(p0,p1,p2,delay)NN=256。 %發(fā)送的信號(hào)數(shù)tb=。 %信號(hào)周期fs=16。 %每信號(hào)的采樣數(shù)ebn0db=[1:2:14]。 %信噪比矢量%產(chǎn)生QPSK信號(hào)x=random_binary(NN,fs)+i*random_binary(NN,fs)。delay0=0。delay1=0。delay2=delay。%設(shè)置瑞利信道gain1=sqrt(p1)*abs(randn(1,NN)+i*randn(1,NN))。gain2=sqrt(p2)*abs(randn(1,NN)+i*randn(1,NN))。 for k=1:NN for kk=1:fs index=(k1)*fs+kk。 ggain1(1,index)=gain1(1,k)。 ggain2(1,index)=gain2(1,k)。 endend y1=x。for k=1:delay2 y2(1,k)=y1(1,k)*sqrt(p0)。endfor k=(delay2+1):(NN*fs) y2(1,k)=y1(1,k)*sqrt(p0)+y1(1,kdelay1)*ggain1(1,k)+y1(1,kdelay2)*ggain2(1,k)。end%匹配濾波器b=ones(1,fs)。b=b/fs。a=1。y=filter(b,a,y2)。%應(yīng)用半解析方法，進(jìn)行誤碼率估計(jì)。[cor lags]=vxcorr(x,y)。cmax=max(max(abs(cor)))。nmax=find(abs(cor)==cmax)。timelag=lags(nmax)。corrmag=cmax。theta=angle(cor(nmax))。y=y*exp(i*theta)。 hh=impz(b,a)。ts=1/fs。nbw=(fs/2)*sum(hh.^2)。 index=(10*fs+8:fs:(NN10)*fs+8)。xx=x(index)。yy=y(indextimelag+1)。[n1 n2]=size(y2)。ny2=n1*n2。eb=tb*sum(sum(abs(y2).^2))/ny2。eb=eb/2。[peideal,pesystem]=qpsk_berest(xx,yy,ebn0db,eb,tb,nbw)。figuresemilogy(ebn0db,peideal,39。b*39。,ebn0db,pesystem,39。r+39。)xlabel(39。E_b/N_039。)。ylabel(39。Probability of Error39。)。gridaxis([0 14 10^(10) 1]) return代碼4function bpskThreeray(p0,p1,p2,delay) NN=256。tb=。fs=16。ebn0db=[1:2:14]。 x=random_binary(NN,fs)。 delay0=0。delay1=0。delay2=delay。 gain1=sqrt(p1)*abs(randn(1,NN)+i*randn(1,NN))。gain2=sqrt(p2)*abs(randn(1,NN)+i*randn(1,NN))。 for k=1:NN for kk=1:fs index=(k1)*fs+kk。 ggain1(1,index)=gain1(1,k)。 ggain2(1,index)=gain2(1,k)。 endend y1=x。for k=1:delay2 y2(1,k)=y1(1,k)*sqrt(p0)。endfor k=(delay2+1):(NN*fs) y2(1,k)=y1(1,k)*sqrt(p0)+y1(1,kdelay1)*ggain1(1,k)+y1(1,kdelay2)*ggain2(1,k)。end b=ones(1,fs)。b=b/fs。a=1。y=filter(b,a,y2)。 [cor lags]=vxcorr(x,y)。cmax=max(max(abs(cor)))。nmax=find(abs(cor)==cmax)。timelag=lags(nmax)。corrmag=cmax。theta=angle(cor(nmax))。y=y*exp(i*theta)。 hh=impz(b,a)。ts=1/fs。nbw=(fs/2)*sum(hh.^2)。 index=(10*fs+8:fs:(NN10)*fs+8)。xx=x(index)。yy=y(indextimelag+1)。[n1 n2]=size(y2)。ny2=n1*n2。eb=tb*sum(sum(abs(y2).^2))/ny2。eb=eb/2。[peideal,pesystem]=qpsk_berest(xx,yy,ebn0db,eb,tb,nbw)。figuresemilogy(ebn0db,peideal,39。b*39。,ebn0db,pesystem,39。r+39。)xlabel(39。E_b/N_039。)。ylabel(39。Probability of Error39。)。gridaxis([0 14 10^(10) 1]) return
職業(yè)教育相關(guān)推薦
鄂ICP備17016276號(hào)-1