【導讀】提高的重要途徑。本文首先分析了數(shù)字調制系統(tǒng)的幾種基本調制解調方法，然后，運用。Matlab設計了這幾種數(shù)字調制解調方法的仿真程序，主要包括PSK，DPSK和16QAM。仿真，分析了這三種調制解調過程中各環(huán)節(jié)時域和頻域的波形，并考慮了信道噪聲的影響。通過仿真更深刻地理解了數(shù)字調制解調系統(tǒng)基本原理。Matlab的編程功能簡單，并且很容易擴展和創(chuàng)造新的命令與函數(shù)。Matlab具有強大的Simulink動態(tài)仿真環(huán)境，可以實現(xiàn)可。Simulink支持連續(xù)、離散及兩者混合的線。計和分析的功能，可以在Matlab環(huán)境下獨立使用，也可以配合Simulink使用。因此，Matlab在通信系統(tǒng)仿真中得到了廣泛應用，本文也選用該工具。例如，2021年我國的移動電話用戶首次超過了固定電話用戶。統(tǒng)計數(shù)據(jù)，到2021年底移動電話用戶近4億。信號的變化而變化。根據(jù)控制的載波參量的不同，數(shù)字調制有調幅、調相和調頻三種基本。形式，并可以派生出多種其他形式。由于傳輸失真、傳輸損耗以及保證帶內特性的原因，

　　

【正文】 K2=ASK2/fs。 f=[0:df1:df1*(length(mk)1)]fs/2。 [NOISE,noise,df1]=fftseq(noise,ts,df)。 35 NOISE=NOISE*ts。 u=y。 [U,u,df1]=fftseq(u,ts,df)。 U=U/fs。 [R,r,df1]=fftseq(r,ts,df)。 R=R/fs。 [Yi,yi,df1]=fftseq(yi,ts,df)。 Yi=Yi/fs。 [Yq,yq,df1]=fftseq(yq,ts,df)。 Yq=Yq/fs。 %濾波 f_cutoff=70。 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)。 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 end for 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 36 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。源信號 1波形 39。)。 subplot(2,2,2) plot(f,abs(fftshift(M1)))。xlabel(39。Frequncy39。)。title(39。源信號 1頻譜 39。)。 subplot(2,2,3) plot(t,mq(1:length(t)))。grid。axis([0 ])。xlabel(39。Time39。)。title(39。源信號 2波形 39。)。 subplot(2,2,4) plot(f,abs(fftshift(M2)))。xlabel(39。Frequncy39。)。title(39。源信號 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信號 1波形 39。)。 subplot(2,2,2) plot(f,abs(fftshift(ASK1)))。xlabel(39。Frequncy39。)。title(39。4ASK信號 1頻譜 39。)。 subplot(2,2,3) plot(t,ask2(1:length(t)))。grid。axis([0 5 5])。xlabel(39。Time39。)。title(39。4ASK信號 2波形 39。)。 subplot(2,2,4) 37 plot(f,abs(fftshift(ASK1)))。xlabel(39。Frequncy39。)。title(39。4ASK信號 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。未加噪聲調制 16QAM波形 39。)。 subplot(3,2,4) plot(f,abs(fftshift(U)))。xlabel(39。Frequncy39。)。title(39。未加噪聲調制 16QAM頻譜 39。)。 subplot(3,2,5) plot(t,r(1:length(t)))。grid。axis([0 5 5])。xlabel(39。Time39。)。title(39。加噪聲調制 16QAM波形 39。)。 subplot(3,2,6) plot(f,abs(fftshift(R)))。xlabel(39。Frequncy39。)。title(39。加噪聲調制 16QAM頻譜 39。)。 pause。 figure(5) subplot(2,2,1) plot(t,yi(1:length(t)))。grid。xlabel(39。Time39。)。title(39。相干解調后信號 1波形 39。)。 subplot(2,2,2) plot(f,abs(fftshift(Yi)))。xlabel(39。Frequncy39。)。title(39。相干解調后信號 1頻譜 39。)。 subplot(2,2,3) plot(t,dem1(1:length(t)))。grid。xlabel(39。Time39。)。title(39。低通后信號 1波形 39。)。 subplot(2,2,4) plot(f,abs(fftshift(DEM1)))。xlabel(39。Frequency39。)。title(39。低通后信號 1頻譜 39。)。 pause。 figure(6) subplot(2,2,1) plot(t,yq(1:length(t)))。grid。xlabel(39。Time39。)。title(39。相干解調后信號 2波形 39。)。 subplot(2,2,2) plot(f,abs(fftshift(Yq)))。xlabel(39。Frequency39。)。title(39。相干解調后信號 2頻譜 39。)。 subplot(2,2,3) plot(t,dem2(1:length(t)))。grid。xlabel(39。Time39。)。title(39。低通后信號 2波形 39。)。 subplot(2,2,4) plot(f,abs(fftshift(DEM2)))。xlabel(39。Frequency39。)。title(39。低通后信號 2頻譜 39。)。 pause。 figure(7) subplot(2,2,1) plot(t,y1(1:length(t)))。grid。axis([0 ])。xlabel(39。Time39。)。title(39。抽樣判決后信號 1波形 39。)。 subplot(2,2,2) plot(f,abs(fftshift(Y1)))。xlabel(39。Frequency39。)。title(39。抽樣判決后信號 1頻譜 39。)。 subplot(2,2,3) plot(t,mi(1:length(t)))。grid。axis([0 ])。xlabel(39。Time39。)。title(39。源信號 1波形 39。)。 subplot(2,2,4) plot(f,abs(fftshift(M1)))。xlabel(39。Frequncy39。)。title(39。源信號 1頻譜 39。)。 pause。 figure(8) subplot(2,2,1) plot(t,y2(1:length(t)))。grid。axis([0 ])。xlabel(39。Time39。)。title(39。抽樣判決后信號 2波形 39。)。 subplot(2,2,2) 38 plot(f,abs(fftshift(Y2)))。xlabel(39。Frequency39。)。title(39。抽樣判決后信號 2頻譜 39。) subplot(2,2,3) plot(t,mq(1:length(t)))。grid。axis([0 ])。xlabel(39。Time39。)。title(39。源信號 2波形 39。)。 subplot(2,2,4) plot(f,abs(fftshift(M2)))。xlabel(39。Frequncy39。)。title(39。源信號 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。 end n2=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 x p=(norm(x)^2)/length(x)。 載波