【導(dǎo)讀】數(shù)字調(diào)制是通信系統(tǒng)中最為重要的環(huán)節(jié)之一，數(shù)字調(diào)制技術(shù)的改進(jìn)也是通信系統(tǒng)性能提高的重要途徑。本文首先分析了數(shù)字調(diào)制系統(tǒng)的幾種基本調(diào)制解調(diào)方法，然后，運(yùn)用Matlab設(shè)計(jì)了這幾種數(shù)字調(diào)制解調(diào)方法的仿真程序，主要包括PSK，DPSK和16QAM。通過仿真，分析了這三種調(diào)制解調(diào)過程中各環(huán)節(jié)時(shí)域和頻域的波形，并考慮了信道噪聲的影響。通過仿真更深刻地理解了數(shù)字調(diào)制解調(diào)系統(tǒng)基本原理。Matlab是一種交互式的、以矩陣為基礎(chǔ)的軟件開發(fā)環(huán)境，它用于科學(xué)和工程的計(jì)算與可視化。應(yīng)用Matlab可方便地解決復(fù)雜數(shù)值計(jì)算問題。數(shù)字調(diào)制是指用數(shù)字基帶信號(hào)對(duì)載波的某些參量進(jìn)行控制，使載波的這些參量隨基帶信號(hào)的變化而變化。根據(jù)控制的載波參量的不同，數(shù)字調(diào)制有調(diào)幅、調(diào)相和調(diào)頻三種基本形式，并可以派生出多種其他形式。為了進(jìn)行長途傳輸，必須對(duì)數(shù)字信號(hào)進(jìn)行載波調(diào)制，將信號(hào)頻譜搬移到高頻處才能在信道中傳輸。另外，由于數(shù)字通信具有建網(wǎng)靈活，容易采用數(shù)字差錯(cuò)控制技術(shù)和數(shù)

　　

【正文】 ,n_cutoff)。DEM1=H.*Yi。dem1=real(ifft(DEM1))/fs。DEM2=H.*Yq。dem2=real(ifft(DEM2))/fs。dem1=dem1(1:length(t))。dem2=dem2(1:length(t))。%抽樣判決for i=(1:10)。 k=dem1((i1)*15+1:i*15)。 x1(i)=mean(k)。 for j=(1:15) l=(i1)*15+j。 if x1(i)=0 x1(i)=max(k) if x1(i)3*10^6 y1(l)=3 else y1(l)=1 end end if x1(i)0 x1(i)=min(k) if x1(i)3*10^6 y1(l)=3 else y1(l)=1 end end end endfor i=(1:10)。 k=dem2((i1)*15+1:i*15)。 x2(i)=mean(k)。 for j=(1:15) l=(i1)*15+j。 if x2(i)=0 x2(i)=max(k) if x2(i)3*10^6 y2(l)=3 else y2(l)=1 end end if x2(i)0 x2(i)=min(k) if x2(i)3*10^6 y2(l)=3 else y2(l)=1 end end end end [Y1,y1,df1]=fftseq(y1,ts,df)。Y1=Y1/fs。[Y2,y2,df1]=fftseq(y2,ts,df)。Y2=Y2/fs。%仿真pause。figure(1)subplot(2,2,1)plot(t,mi(1:length(t)))。grid。axis([0 ])。xlabel(39。Time39。)。title(39。源信號(hào)1波形39。)。subplot(2,2,2)plot(f,abs(fftshift(M1)))。xlabel(39。Frequncy39。)。title(39。源信號(hào)1頻譜39。)。subplot(2,2,3)plot(t,mq(1:length(t)))。grid。axis([0 ])。xlabel(39。Time39。)。title(39。源信號(hào)2波形39。)。subplot(2,2,4)plot(f,abs(fftshift(M2)))。xlabel(39。Frequncy39。)。title(39。源信號(hào)2頻譜39。)。pause。figure(2)subplot(2,2,1)plot(t,c(1:length(t)))。grid。axis([0 2 2])。xlabel(39。Time39。)。title(39。載波1波形39。)。subplot(2,2,2)plot(f,abs(fftshift(C)))。xlabel(39。Frequncy39。)。title(39。載波1頻譜39。)。subplot(2,2,3)plot(t,s(1:length(t)))。grid。axis([0 2 2])。xlabel(39。Time39。)。title(39。載波2波形39。)。subplot(2,2,4)plot(f,abs(fftshift(S)))。xlabel(39。Frequncy39。)。title(39。載波2頻譜39。)。pause。figure(3)subplot(2,2,1)plot(t,ask1(1:length(t)))。grid。axis([0 5 5])。xlabel(39。Time39。)。title(39。4ASK信號(hào)1波形39。)。subplot(2,2,2)plot(f,abs(fftshift(ASK1)))。xlabel(39。Frequncy39。)。title(39。4ASK信號(hào)1頻譜39。)。subplot(2,2,3)plot(t,ask2(1:length(t)))。grid。axis([0 5 5])。xlabel(39。Time39。)。title(39。4ASK信號(hào)2波形39。)。subplot(2,2,4)plot(f,abs(fftshift(ASK1)))。xlabel(39。Frequncy39。)。title(39。4ASK信號(hào)2頻譜39。)。pause。figure(4)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,u(1:length(t)))。grid。axis([0 5 5])。xlabel(39。Time39。)。title(39。未加噪聲調(diào)制16QAM波形39。)。subplot(3,2,4)plot(f,abs(fftshift(U)))。xlabel(39。Frequncy39。)。title(39。未加噪聲調(diào)制16QAM頻譜39。)。subplot(3,2,5)plot(t,r(1:length(t)))。grid。axis([0 5 5])。xlabel(39。Time39。)。title(39。加噪聲調(diào)制16QAM波形39。)。subplot(3,2,6)plot(f,abs(fftshift(R)))。xlabel(39。Frequncy39。)。title(39。加噪聲調(diào)制16QAM頻譜39。)。pause。figure(5)subplot(2,2,1)plot(t,yi(1:length(t)))。grid。xlabel(39。Time39。)。title(39。相干解調(diào)后信號(hào)1波形39。)。subplot(2,2,2)plot(f,abs(fftshift(Yi)))。xlabel(39。Frequncy39。)。title(39。相干解調(diào)后信號(hào)1頻譜39。)。subplot(2,2,3)plot(t,dem1(1:length(t)))。grid。xlabel(39。Time39。)。title(39。低通后信號(hào)1波形39。)。subplot(2,2,4)plot(f,abs(fftshift(DEM1)))。xlabel(39。Frequency39。)。title(39。低通后信號(hào)1頻譜39。)。pause。figure(6)subplot(2,2,1)plot(t,yq(1:length(t)))。grid。xlabel(39。Time39。)。title(39。相干解調(diào)后信號(hào)2波形39。)。subplot(2,2,2)plot(f,abs(fftshift(Yq)))。xlabel(39。Frequency39。)。title(39。相干解調(diào)后信號(hào)2頻譜39。)。subplot(2,2,3)plot(t,dem2(1:length(t)))。grid。xlabel(39。Time39。)。title(39。低通后信號(hào)2波形39。)。subplot(2,2,4)plot(f,abs(fftshift(DEM2)))。xlabel(39。Frequency39。)。title(39。低通后信號(hào)2頻譜39。)。pause。figure(7)subplot(2,2,1)plot(t,y1(1:length(t)))。grid。axis([0 ])。xlabel(39。Time39。)。title(39。抽樣判決后信號(hào)1波形39。)。subplot(2,2,2)plot(f,abs(fftshift(Y1)))。xlabel(39。Frequency39。)。title(39。抽樣判決后信號(hào)1頻譜39。)。subplot(2,2,3)plot(t,mi(1:length(t)))。grid。axis([0 ])。xlabel(39。Time39。)。title(39。源信號(hào)1波形39。)。subplot(2,2,4)plot(f,abs(fftshift(M1)))。xlabel(39。Frequncy39。)。title(39。源信號(hào)1頻譜39。)。pause。figure(8)subplot(2,2,1)plot(t,y2(1:length(t)))。grid。axis([0 ])。xlabel(39。Time39。)。title(39。抽樣判決后信號(hào)2波形39。)。subplot(2,2,2)plot(f,abs(fftshift(Y2)))。xlabel(39。Frequency39。)。title(39。抽樣判決后信號(hào)2頻譜39。)subplot(2,2,3)plot(t,mq(1:length(t)))。grid。axis([0 ])。xlabel(39。Time39。)。title(39。源信號(hào)2波形39。)。subplot(2,2,4)plot(f,abs(fftshift(M2)))。xlabel(39。Frequncy39。)。title(39。源信號(hào)2頻譜39。)。function [M,m,df]=fftseq(m,ts,df) % [M,m,df]=fftseq(m,ts,df)% [M,m,df]=fftseq(m,ts)%FFTSEQ generates M, the FFT of the sequence m.% The sequence is zero padded to meet the required frequency resolution df.% ts is the sampling interval. The output df is the final frequency resolution.% Output m is the zero padded version of input m. M is the FFT.fs=1/ts。if nargin == 2 n1=0。else n1=fs/df。endn2=length(m)。n=2^(max(nextpow2(n1),nextpow2(n2)))。M=fft(m,n)。m=[m,zeros(1,nn2)]。df=fs/n。function p=spower(x)% p=spower(x)%SPOWER returns the power in signal xp=(norm(x)^2)/length(x)。