【文章內(nèi)容簡(jiǎn)介】 源信號(hào)相比都發(fā)生了明顯變化，說(shuō)明都產(chǎn)生了誤碼，如圖３－１４所示?！「鞣N調(diào)制比較 （１）2PSK、2DPSK檢測(cè)這兩種信號(hào)時(shí)判決器均可工作在最佳門(mén)限電平（零電平）。2DPSK抗噪聲性能不及2PSK。2PSK系統(tǒng)存在“反向工作”問(wèn)題，而2DPSK系統(tǒng)不存在“反向工作”問(wèn)題，因此在實(shí)際應(yīng)用中，真正作為傳輸用的數(shù)字調(diào)相信號(hào)幾乎都是DPSK信號(hào)。2PSK信號(hào)的帶寬、頻帶利用率和2DPSK信號(hào)的相同。（２）16QAM與2PSK、2DPSK相比，在相同的碼元速率情況下其信息速率較高，頻帶利用率高，但誤碼率最高，且實(shí)現(xiàn)的難度大。第四章 結(jié)束語(yǔ)數(shù)字調(diào)制技術(shù)的發(fā)展日新月異，如今在現(xiàn)實(shí)中應(yīng)用的數(shù)字調(diào)制系統(tǒng)大部分是經(jīng)過(guò)改進(jìn)的，性能較好的系統(tǒng)，但是，作為理論發(fā)展最成熟的調(diào)制方式，PSK，DPSK，QAM等的研究仍然具有很重大的意義，因此，我們選擇了這幾種調(diào)制方式做仿真研究。仿真這幾種理論已經(jīng)很成熟的數(shù)字調(diào)制方式，一方面，可以更容易將仿真結(jié)果與成熟的理論進(jìn)行比較，從而驗(yàn)證仿真的合理性；另一方面，也可以以此為基礎(chǔ)將仿真系統(tǒng)進(jìn)行改進(jìn)擴(kuò)展，使其成為仿真更多的數(shù)字調(diào)制方式的模板。另外，數(shù)字調(diào)制解調(diào)系統(tǒng)只是通信系統(tǒng)的一個(gè)重要組成部分，因此，我們所設(shè)計(jì)的數(shù)字調(diào)制仿真系統(tǒng)也可以擴(kuò)展成通信系統(tǒng)的仿真。這種擴(kuò)展只需在輸入端與調(diào)制器間增加一些數(shù)字基帶處理模塊，如信源編碼、加密、信道編碼等，在解調(diào)后增加相應(yīng)的解碼解密器即可。參考文獻(xiàn)[1] 樊昌信等.　通信原理（第6版）　2008　國(guó)防工業(yè)出版社[2]　鄧華等.　MATLAB通信仿真及應(yīng)用實(shí)例詳解　2003　人民郵電出版社[3] 王立寧,樂(lè)光新,詹菲.　MATLAB與通信仿真　2000　人民郵電出版社致 謝附 錄2PSK調(diào)制解調(diào)程序clear。t0=。ts=。fc=200。snr=10。fs=1/ts。df=。t=[ts:ts:t0]。snr_lin=10^(snr/10)。y=[]。nm=length(t)。%生成2PSKfor i=1:10 x=rand。 for j=1:15 l=(i1)*15+j。 c(l)=cos(2*pi*fc*t(l))。 c1(l)=cos(2*pi*fc*t(l)+pi)。 if x= mi(l)=1。 y(l)=c(l)。 elseif x mi(l)=0。 y(l)=c1(l)。 end endend%加噪聲signal_power=spower(y(1:length(t)))。 % power in modulated signalnoise_power=signal_power/snr_lin。 % pute noise powernoise_std=sqrt(noise_power)。 % pute noise standard deviationnoise=noise_std*randn(1,length(y))。 % generate noiser=y+noise。 % add noise to the modulated signal%時(shí)域頻域轉(zhuǎn)換m=mi。[M,m,df1]=fftseq(m,ts,df)。M=M*ts。f=[0:df1:df1*(length(m)1)]fs/2。[Y,y,df1]=fftseq(y,ts,df)。Y=Y*ts。[C,c,df1]=fftseq(c,ts,df)。C=C*ts。[R,r,df1]=fftseq(r,ts,df)。R=R*ts。[NOISE,noise,df1]=fftseq(noise,ts,df)。NOISE=NOISE*ts。%解調(diào)f_cutoff=70。%70Hz低通濾波器n_cutoff=floor(70/df1)。H=zeros(size(f))。H(1:n_cutoff)=2*ones(1,n_cutoff)。H(length(f)n_cutoff+1:length(f))=2*ones(1,n_cutoff)。yy=r.*c。[YY,yy,df1]=fftseq(yy,ts,df)。YY=YY*ts。DEM=H.*YY。%濾波dem=real(ifft(DEM))/fs。dem=dem(1:length(t))。yout=dem。for i=(1:length(t))。% 判決，得到解調(diào)結(jié)果 if yout(1,i)0。 yout(1,i)=1。 else yout(1,i)=0。 end。end。 [YOUT,yout,df1]=fftseq(yout,ts,df)。YOUT=YOUT*ts。pause。%仿真pause。figure(1)。subplot(2,2,1)。plot(t,mi(1:length(t)))。grid。axis([0 2 2])。xlabel(39。Time39。)。title(39。源信號(hào)波形39。)。subplot(2,2,2)。plot(f,abs(fftshift(M)))。xlabel(39。Frequency39。)。 title(39。源信號(hào)頻譜39。)。subplot(2,2,3)。plot(t,c(1:length(t)))。grid。axis([0 2 2])。xlabel(39。Time39。)。title(39。載波波形39。)。subplot(2,2,4)。plot(f,abs(fftshift(C)))。xlabel(39。Frequency39。)。title(39。載波頻譜39。)。pause。figure(2)subplot(3,2,1)。plot(t,noise(1:length(t)))。grid。axis([0 2 2])。xlabel(39。Time39。)。title(39。噪聲波形39。)。subplot(3,2,2)。plot(f,abs(fftshift(NOISE)))。xlabel(39。Frequency39。)。title(39。噪聲頻譜39。)。subplot(3,2,3)。plot(t,y(1:length(t)))。grid。axis([0 2 2])。xlabel(39。Time39。)。title(39。未加噪聲調(diào)制波形39。)。subplot(3,2,4)。plot(f,abs(fftshift(Y)))。xlabel(39。Frequency39。)。 title(39。未加噪聲調(diào)制頻譜39。)。subplot(3,2,5)。plot(t,r(1:length(t)))。grid。axis([0 2 2])。xlabel(39。Time39。)。title(39。加噪聲調(diào)制波形39。)。subplot(3,2,6)。plot(f,abs(fftshift(R)))。xlabel(39。Frequency39。)。 title(39。加噪聲調(diào)制頻譜39。)。pause。figure(3)。subplot(2,2,1)。plot(t,yy(1:length(t)))。grid。axis([0 2 2])。xlabel(39。Time39。)。title(39。相干解調(diào)后波形39。)。subplot(2,2,2)。plot(f,abs(fftshift(YY)))。xlabel(39。Frequency39。)。 title(39。相干解調(diào)后頻譜39。)。subplot(2,2,3)。plot(t,dem(1:length(t)))。grid。xlabel(39。Time39。)。title(39。低通后波形39。)。subplot(2,2,4)。plot(f,abs(fftshift(DEM)))。xlabel(39。Frequency39。)。 title(39。低通后頻譜39。)。pause。figure(4)subplot(2,2,1)。plot(t,yout(1:length(t)))。grid。axis([0 2 2])。xlabel(39。Time39。)。title(39。抽樣判決后波形39。)。subplot(2,2,2)。plot(f,abs(fftshift(YOUT)))。xlabel(39。Frequency39。)。 title(39。抽樣判決后頻譜39。)。subplot(2,2,3)。plot(t,mi(1:length(t)))。grid。axis([0 2 2])。xlabel(39。Time39。)。title(39。源信號(hào)波形39。)。subplot(2,2,4)。plot(f,abs(fftshift(M)))。xlabel(39。Frequency39。)。 title(39。源信號(hào)頻譜39。)。2DPSK調(diào)制解調(diào)程序clear。t0=。ts=。fc=200。snr=10。fs=1/ts。df=。t=[ts:ts:t0]。snr_lin=10^(snr/10)。y=[]。nm=length(t)。%生成2DPSKbn=[]。bn(1)=1。%bn是差分碼，設(shè)bn的第一個(gè)符號(hào)為1for i=1:10 x=rand。 if x= m(i)=1。 else m(i)=0。 end bn(i+1)=xor(m(i),bn(i))。 for j=1:15 l=(i1)*15+j。 c(l)=cos(2*pi*fc*t(l))。 c1(l)=cos(2*pi*fc*t(l)+pi)。 if m(i)==1。 mi(l)=1。 else mi(l)=0。 end if bn(i+1)==1。 bn1(l)=1。 y(l)=c(l)。 else bn1(l)=0。 y(l)=c1(l)。 end endend%加噪聲signal_power=spower(y(1:length(t)))。 % power in modulated signalnoise_power=signal_power/snr_lin。 % pute noise powernoise_std=sqrt(noise_power)。 % pute noise standard deviationnoise=noise_std*randn(1,length(y))。 % generate noiser=y+noise。 % add noise to the modulated signal%時(shí)域頻域轉(zhuǎn)換m=mi。[M,m,df1]=fftseq(m,ts,df)。M=M*ts。[BN1,bn1,df1]=fftse