【正文】 最高，且實(shí)現(xiàn)的難度大。這種擴(kuò)展只需在輸入端與調(diào)制器間增加一些數(shù)字基帶處理模塊，如信源編碼、加密、信道編碼等，在解調(diào)后增加相應(yīng)的解碼解密器即可。fc=200。t=[ts:ts:t0]。%生成2PSKfor i=1:10 x=rand。 if x= mi(l)=1。 end endend%加噪聲signal_power=spower(y(1:length(t)))。 % generate noiser=y+noise。f=[0:df1:df1*(length(m)1)]fs/2。C=C*ts。NOISE=NOISE*ts。H(1:n_cutoff)=2*ones(1,n_cutoff)。YY=YY*ts。yout=dem。 else yout(1,i)=0。YOUT=YOUT*ts。subplot(2,2,1)。xlabel(39。源信號波形39。xlabel(39。源信號頻譜39。grid。)。subplot(2,2,4)。)。pause。axis([0 2 2])。title(39。plot(f,abs(fftshift(NOISE)))。title(39。plot(t,y(1:length(t)))。Time39。)。Frequency39。)。axis([0 2 2])。title(39。plot(f,abs(fftshift(R)))。 title(39。figure(3)。axis([0 2 2])。title(39。plot(f,abs(fftshift(YY)))。 title(39。plot(t,dem(1:length(t)))。)。subplot(2,2,4)。)。pause。axis([0 2 2])。title(39。plot(f,abs(fftshift(YOUT)))。 title(39。plot(t,mi(1:length(t)))。Time39。)。Frequency39。)。fc=200。t=[ts:ts:t0]。%生成2DPSKbn=[]。 else m(i)=0。 c1(l)=cos(2*pi*fc*t(l)+pi)。 end if bn(i+1)==1。 y(l)=c1(l)。 % pute noise standard deviationnoise=noise_std*randn(1,length(y))。M=M*ts。[Y,y,df1]=fftseq(y,ts,df)。[R,r,df1]=fftseq(r,ts,df)。%解調(diào)f_cutoff=70。H(length(f)n_cutoff+1:length(f))=2*ones(1,n_cutoff)。DEM=H.*YY。y1=ytemp。 else y1(1,i)=0。%碼反變換 k=y1((i1)*15+1:i*15)。 end endy3(1)=xor(y2(1),1)。 for j=(1:15)。 else yout(l)=0。YOUT=YOUT*ts。subplot(3,2,1)。xlabel(39。源信號波形39。xlabel(39。源信號頻譜39。grid。)。subplot(3,2,4)。)。subplot(3,2,5)。xlabel(39。載波波形39。xlabel(39。載波頻譜39。subplot(3,2,1)。xlabel(39。噪聲波形39。xlabel(39。噪聲頻譜39。grid。)。subplot(3,2,4)。)。subplot(3,2,5)。xlabel(39。加噪聲調(diào)制波形39。xlabel(39。加噪聲調(diào)制頻譜39。subplot(3,2,1)。xlabel(39。相干解調(diào)后波形39。xlabel(39。相干解調(diào)后頻譜39。grid。title(39。plot(f,abs(fftshift(DEM)))。 title(39。plot(t,y1(1:length(t)))。Time39。)。Frequency39。)。plot(t,yout(1:length(t)))。Time39。)。Frequency39。)。axis([0 2 2])。title(39。plot(f,abs(fftshift(M)))。 title(39。t0=。fs=1/ts。nm=length(t)。 s(l)=sin(2*pi*fc*t(l))。 else mi(l)=3。 s(l)=sin(2*pi*fc*t(l))。 else mq(l)=3。ask1=mi.*c。 % power in modulated signalnoise_power=signal_power/snr_lin。 % add noise to the modulated signalyi=r.*c。ml=mq。M2=M2/fs。S=S/fs。ASK2=ASK2/fs。u=y。R=R/fs。Yq=Yq/fs。H(1:n_cutoff)=2*ones(1,n_cutoff)。DEM2=H.*Yq。%抽樣判決for i=(1:10)。 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)。 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)。%仿真pause。xlabel(39。源信號1波形39。Frequncy39。)。xlabel(39。源信號2波形39。Frequncy39。)。axis([0 2 2])。title(39。xlabel(39。載波1頻譜39。axis([0 2 2])。title(39。xlabel(39。載波2頻譜39。grid。)。subplot(2,2,2)plot(f,abs(fftshift(ASK1)))。title(39。grid。)。subplot(2,2,4)plot(f,abs(fftshift(ASK1)))。title(39。figure(4)subplot(3,2,1)plot(t,noise(1:length(t)))。Time39。)。)。subplot(3,2,3)plot(t,u(1:length(t)))。Time39。)。)。subplot(3,2,5)plot(t,r(1:length(t)))。Time39。)。)。pause。Time39。)。)。subplot(2,2,3)plot(t,dem1(1:length(t)))。)。subplot(2,2,4)plot(f,abs(fftshift(DEM1)))。title(39。figure(6)subplot(2,2,1)plot(t,yq(1:length(t)))。)。subplot(2,2,2)plot(f,abs(fftshift(Yq)))。title(39。grid。title(39。xlabel(39。低通后信號2頻譜39。grid。)。subplot(2,2,2)plot(f,abs(fftshift(Y1)))。title(39。grid。)。subplot(2,2,4)plot(f,abs(fftshift(M1)))。title(39。figure(8)subplot(2,2,1)plot(t,y2(1:length(t)))。Time39。)。)。grid。)。subplot(2,2,4)plot(f,abs(fftshift(M2)))。title(39。if nargin == 2 n1=0。M=fft(m,n)。df=fs/n。endn2=length(m)。)。Frequncy39。源信號2波形39。xlabel(39。抽樣判決后信號2頻譜39。xlabel(39。title(39。axis([0 ])。)。Frequncy39。源信號1波形39。xlabel(39。)。Frequency39。抽樣判決后信號1波形39。xlabel(39。pause。)。)。Time39。)。Frequency39。相干解調(diào)后信號2波形39。xlabel(39。)。Frequency39。低通后信號1波形39。xlabel(39。相干解調(diào)后信號1頻譜39。xlabel(39。title(39。grid。加噪聲調(diào)制16QAM頻譜39。xlabel(39。title(39。axis([0 5 5])。未加噪聲調(diào)制16QAM頻譜39。xlabel(39。title(39。axis([0 5 5])。噪聲頻譜39。xlabel(39。title(39。axis([0 2 2])。)。Frequncy39。4ASK信號2波形39。xlabel(39。)。Frequncy39。4ASK信號1波形39。xlabel(39。pause。)。)。Time39。subplot(2,2,3)plot(t,s(1:length(t)))。)。)。Time39。figure(2)subplot(2,2,1)plot(t,c(1:length(t)))。title(39。subplot(2,2,4)plot(f,abs(fftshift(M2)))。)。grid。title(39。subplot(2,2,2)plot(f,abs(fftshift(M1)))。)。grid。[Y2,y2,df1]=fftseq(y2,ts,df)。 x2(i)=mean(k)。 x1(i)=mean(k)。dem1=dem1(1:length(t))。DEM1=H.*Yi。n_cutoff=floor(70/df1)。Yi=Yi/fs。U=U/fs。[NOISE,noise,df1]=fftseq(noise,t