【正文】 。df=fs/n。M=fft(m,n)。endn2=length(m)。if nargin == 2 n1=0。)。title(39。Frequncy39。subplot(2,2,4)plot(f,abs(fftshift(M2)))。源信號2波形39。)。xlabel(39。grid。抽樣判決后信號2頻譜39。)。xlabel(39。)。title(39。Time39。axis([0 ])。figure(8)subplot(2,2,1)plot(t,y2(1:length(t)))。)。title(39。Frequncy39。subplot(2,2,4)plot(f,abs(fftshift(M1)))。源信號1波形39。)。xlabel(39。grid。)。title(39。Frequency39。subplot(2,2,2)plot(f,abs(fftshift(Y1)))。抽樣判決后信號1波形39。)。xlabel(39。grid。pause。低通后信號2頻譜39。)。xlabel(39。)。title(39。Time39。grid。)。title(39。Frequency39。subplot(2,2,2)plot(f,abs(fftshift(Yq)))。相干解調(diào)后信號2波形39。)。xlabel(39。figure(6)subplot(2,2,1)plot(t,yq(1:length(t)))。)。title(39。Frequency39。subplot(2,2,4)plot(f,abs(fftshift(DEM1)))。低通后信號1波形39。)。xlabel(39。subplot(2,2,3)plot(t,dem1(1:length(t)))。相干解調(diào)后信號1頻譜39。)。xlabel(39。)。title(39。Time39。grid。pause。加噪聲調(diào)制16QAM頻譜39。)。xlabel(39。)。title(39。Time39。axis([0 5 5])。subplot(3,2,5)plot(t,r(1:length(t)))。未加噪聲調(diào)制16QAM頻譜39。)。xlabel(39。)。title(39。Time39。axis([0 5 5])。subplot(3,2,3)plot(t,u(1:length(t)))。噪聲頻譜39。)。xlabel(39。)。title(39。Time39。axis([0 2 2])。figure(4)subplot(3,2,1)plot(t,noise(1:length(t)))。)。title(39。Frequncy39。subplot(2,2,4)plot(f,abs(fftshift(ASK1)))。4ASK信號2波形39。)。xlabel(39。grid。)。title(39。Frequncy39。subplot(2,2,2)plot(f,abs(fftshift(ASK1)))。4ASK信號1波形39。)。xlabel(39。grid。pause。載波2頻譜39。)。xlabel(39。)。title(39。Time39。axis([0 2 2])。subplot(2,2,3)plot(t,s(1:length(t)))。載波1頻譜39。)。xlabel(39。)。title(39。Time39。axis([0 2 2])。figure(2)subplot(2,2,1)plot(t,c(1:length(t)))。)。title(39。Frequncy39。subplot(2,2,4)plot(f,abs(fftshift(M2)))。源信號2波形39。)。xlabel(39。grid。)。title(39。Frequncy39。subplot(2,2,2)plot(f,abs(fftshift(M1)))。源信號1波形39。)。xlabel(39。grid。%仿真pause。[Y2,y2,df1]=fftseq(y2,ts,df)。 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)。 x2(i)=mean(k)。 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)。 x1(i)=mean(k)。%抽樣判決for i=(1:10)。dem1=dem1(1:length(t))。DEM2=H.*Yq。DEM1=H.*Yi。H(1:n_cutoff)=2*ones(1,n_cutoff)。n_cutoff=floor(70/df1)。Yq=Yq/fs。Yi=Yi/fs。R=R/fs。U=U/fs。u=y。[NOISE,noise,df1]=fftseq(noise,ts,df)。ASK2=ASK2/fs。ASK1=ASK1/fs。S=S/fs。C=C/fs。M2=M2/fs。M1=M1/fs。ml=mq。%%時(shí)域頻域轉(zhuǎn)換。 % add noise to the modulated signalyi=r.*c。 % pute noise standard deviationnoise=noise_std*randn(1,length(y))。 % power in modulated signalnoise_power=signal_power/snr_lin。y=mi.*c+mq.*s。ask1=mi.*c。mi(nm)=mi(nm1)。 else mq(l)=3。 elseif x= mq(l)=1。 s(l)=sin(2*pi*fc*t(l))。 for j=1:15 l=(i1)*15+j。 else mi(l)=3。 elseif x= mi(l)=1。 s(l)=sin(2*pi*fc*t(l))。 for j=1:15 l=(i1)*15+j。nm=length(t)。t=[ts:ts:t0]。fs=1/ts。fc=200。t0=。)。 title(39。Frequency39。plot(f,abs(fftshift(M)))。)。title(39。Time39。axis([0 2 2])。plot(t,mi(1:length(t)))。)。 title(39。Frequency39。plot(f,abs(fftshift(YOUT)))。)。title(39。Time39。axis([0 2 2])。plot(t,yout(1:length(t)))。figure(4)。)。 title(39。Frequency39。plot(f,abs(fftshift(Y1)))。)。title(39。Time39。axis([0 2 2])。plot(t,y1(1:length(t)))。)。 title(39。Frequency39。plot(f,abs(fftshift(DEM)))。)。title(39。Time39。grid。subplot(3,2,3)。相干解調(diào)后頻譜39。)。xlabel(39。subplot(3,2,2)。相干解調(diào)后波形39。)。xlabel(39。grid。subplot(3,2,1)。pause。加噪聲調(diào)制頻譜39。)。xlabel(39。subplot(3,2,6)。加噪聲調(diào)制波形39。)。xlabel(39。grid。subplot(3,2,5)。未加噪聲調(diào)制頻譜39。)。xlabel(39。subplot(3,2,4)。未加噪聲調(diào)制波形39。)。xlabel(39。grid。subplot(3,2,3)。噪聲頻譜39。)。xlabel(39。subplot(3,2,2)。噪聲波形39。)。xlabel(39。grid。subplot(3,2,1)。pause。載波頻譜39。)。xlabel(39。subplot(3,2,6)。載波波形39。)。xlabel(39。grid。subplot(3,2,5)。碼變換后頻譜39。)。xlabel(39。subplot(3,2,4)。碼變換后波形39。)。xlabel(39。grid。subplot(3,2,3)。源信號頻譜39。)。xlabel(39。subplot(3,2,2)。源信號波形39。)。xlabel(39。grid。subplot(3,2,1)。%仿真pause。YOUT=YOUT*ts。Y1=Y1*ts。 else yout(l)=0。 if y3(i)==1。 for j=(1:15)。 y3(i)=xor(y2(i1),y2(i))。 end endy3(1)=xor(y2(1),1)。 y2(i)=1。%碼反變換 k=y1((i1)*15+1:i*15)。end。 else y1(1,i)=0。% 判決，得到解調(diào)結(jié)果 if y1(1,i)0。y1=ytemp。dem=dem(1:length(t))。DEM=H.*YY。[YY,yy,df1]=fftseq(yy,ts,df)。H(length(f)n_cutoff+1:length(f))=2*ones(1,n_cutoff)。H=zeros(size(f))。%解調(diào)f_cutoff=70。[NOISE,noise,df1]=fftseq(noise,ts,df)。[R,r,df1]=fftseq(r,ts,df)。[C,c,df1]=fftseq(c,ts,df)。[Y,y,df1]=fftseq(y,ts,df)。BN1=BN1*ts。M=M*ts。 % add noise to the modula