【正文】 axis([ ])。k:39。 end。 a=abs(c(i))/cmax*。)。 cmax=max(abs(c))。 i=i+1。 while i = N/m。 n=fs/fc。 xn=xn39。 xn(i,:)=T(yn(i)+1,:)。1 1。 I3=find(yn=)。 yn(I1)=ones(size(I1))。 %設(shè)置門限電平，判決 I0=find(yn )。 ymax=max(y)。 y=y(:)39。 for i=1:m1。 ii=ii+1。 ii=1。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%二進制轉(zhuǎn)換成四進制 function [y,yn]=two2four(x,m)。 N=16。 if nargin5。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%基帶升余弦成形濾波器 function y=bshape(x,fs,fb,N,alfa,delay)。 xn=[xn(1:nn/2)。 In=I(n)。 m=4*fs/fb。 Q=y.*sin(2*pi*fc*t)。 %調(diào)制 %%%%%%%%%%%%%%%%%%%%%%%%%%%%% %QAM%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%信號解調(diào) function [xn,x]=qamdet(y,fs,fb,fc)。 %基帶成形濾波 I=bshape(I,fs,fb/4)。 [I,In]=two2four(I,4*m)。 nn=length(x)。 % 1/2的概率輸出為1 end end。 end。plot(Fy_output)。title(39。已調(diào)信號加噪后輸出39。 figure。 end。 f1=m/n1。 Fy_output=abs(Fy_output(1:fix(n1/2)))*2。 %AWGN信道信噪比 for j=1:length(SNR_in_dB) y_add_noise=awgn(y2,SNR_in_dB(j))。 title(39。f/fb39。 title(39。 f=0:f1:(length(y)1)*f1。 q=find(y1e04)。)。)subplot(222)。 subplot(221)。 m=fs/fb。 %產(chǎn)生二進制信號序列 [y,I,Q]=qam(info,Kbase,fs,fb,fc)。 %基帶信號頻率 fs=32。 圖2 16QAM正交相干解調(diào) 程序代碼：clc。 圖中串/并變換器將速率為Rb的二進制碼元序列分為兩路,速率為Rb/,4電平信號與正交載波相乘,完成正交調(diào)制,經(jīng)24電平變換器輸出為4電平信號,即M=,輸出16個信號狀態(tài),即16QAM，RS=Rb/log216=RB/4。16進制的正交振幅調(diào)制是一種振幅相位聯(lián)合鍵控信號。運行結(jié)果： 分析可知：解調(diào)后恢復原調(diào)制信號。axis([0 100 2 2])。原序列39。figure(3)subplot(2,1,1)。axis([0 100 2 2])。高斯白噪聲39。figure(2)subplot(2,1,1)。axis([0 100 2 2])。4psk信號頻譜39。)。plot(t1,QPSK_rc)。for i=1:2e4 ddd2(i)=ddd1(i)。dd=2*bit_in+1。 bit_recover=[bit_recover 1]。 end。t2=0::。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%低通濾波I_recover=conv(I_demo,xrc)。FQPSK_rc=fft(QPSK_rc)。I_rc=data_I2_rc.*cos(2*pi*f1*t1)。data_Q2_rc=conv(data_Q2,xrc)/。end。data_Q1=repmat(data_Q39。data_I=2*bit_I+1。close all。)。時間 (s)39。)。subplot(326)。)。解調(diào)后輸出信號39。 else Vs(ii)=Low。 vt=0。)。 plot(t,so)。%gain*[1 1] so=filter(b1,a1,cm)。 gain1=。 title(39。 %解調(diào) subplot(323)。通過帶通濾波后輸出39。 subplot(322)。 %Hz(Hz==0)=10^(8)。 a=.8547。 f=[0:3e3:4e5]。psk加噪后信號頻譜39。 figure(2)。 title(39。 %對信號進行傅里葉變換Fpskn=abs(Fpskn(1:length(Fpskn)/2+1))。)。 plot(t,noise)。 %載波 psk=x.*car。 % 方波 am=1。 %耐奎斯特頻率 %用正弦波產(chǎn)生方波 twopi_fc_t=2*pi*fm*t。final=(1/fs)*(n1)。 close all。只有對PSK信號進行非線性變換，才能產(chǎn)生載波分量。而就鍵控法來說，用數(shù)字基帶信號s(t)控制開關(guān)電路，選擇不同相位的載波輸出，這時s(t)為單極性NRZ或雙極性NRZ脈沖序列信號均可。PSK信號頻譜39。 grid on。 axis([0 200e6 2 2])。 %對信號進行傅里葉變換Fpsk=abs(Fpsk(1:length(Fpsk)/2+1))。)。 %顯示基帶信號圖像plot(t,x)。 % 方波 am=1。 %耐奎斯特頻率 %用正弦波產(chǎn)生方波 twopi_fc_t=2*pi*fm*t。 final=(1/fs)*(n1)。 %清除close all。原基帶信號39。抽樣判決后波形39。subplot(325)。at(1,j)=0。)。subplot(324)。)。subplot(323)。調(diào)制信號加高斯噪聲后頻%譜%39。plot(t,fsk1)。已調(diào)信號頻譜39。plot(t,fsk)。Ffsk=fft(fsk)。title(39。s1*st139。 %加入載波1F2=st2.*s2。title(39。title(39?；鶐盘柗创a39。%顯示基帶波形axis([0,5,1,2])。plot(t,st1)。 %高斯噪聲fsk1=fsk+nosie。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%載波信號s1=cos(2*pi*f1*t)s2=cos(2*pi*f2*t)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%調(diào)制F1=st1.*s1。 if st1(n)=1。 for m=j/i*(n1)+1:j/i*n st1(m)=0。 %基帶信號頻率B1=2*f1。 %產(chǎn)生隨機序列t=linspace(0,5,j)。假設(shè)信道噪聲為高斯白噪聲。、2FSK解調(diào)原理2FSK的解調(diào)方法有包絡檢波法、相干檢波法、過零檢測法、差分檢測法，本實驗用到相干檢測法。這兩種方法產(chǎn)生的2FSK信號的波形基本相同，只有一點差異，即由調(diào)頻器產(chǎn)生的2FSK信號在相鄰碼元之間的相位是連續(xù)的，而鍵控法產(chǎn)生的2FSK信號，則分別有兩個獨立的頻率源產(chǎn)生兩個不同頻率的信號，故相鄰碼元的相位不一定是連續(xù)的。2FSK信號的產(chǎn)生方法有兩種：（1） 模擬法，即用數(shù)字基帶信號作為調(diào)制信號進行調(diào)頻。解調(diào)抽樣后恢復原二進制序列。原始信號波形與解調(diào)后的信號作對比39。 %抽樣判決后的數(shù)字信號波形n_tau=tau+。 %抽樣判決的閾值設(shè)置t_judge=(0:Lm1)*Lcarry1+Lcarry1/2。plot(f1,Fx),title(39。 %解調(diào)信號頻譜f1=(0:length(Fx)1)*f/length(Fx)/2。 h=tf(num,den)。)。 %對信號進行傅里葉變換Fask=abs(Fask(1:length(Fask)/2+1))。 %最小時刻Tmax=max(tau)。 %加高斯白噪聲后的已調(diào)信號sig_mul=ask.*carry2。 %調(diào)制后的信號 V=。 %一個正弦信號周期內(nèi)的時間carry1=repmat(A*sin(2*pi*f*t),1,Q)。 %正弦載波信號的頻率A=1。 %清除空間變量m=[1 0 1 1 0 1 1 0]。S( t )定時脈沖 帶通 相乘 低通 抽樣判決eo( t )cosωc t2fs圖5 相干解調(diào)帶通 整流 低通 抽樣判決判決定時脈沖 S( t ) eo( t ) abcd包絡檢波器通常由整流器和低通濾波器組成，可以直接從已調(diào)波中提取原始基帶信號，結(jié)構(gòu)簡單，如圖6所示。信號加噪后波形：、2ASK解調(diào)原理與