【正文】 itle(39。)。 figure(2)。 plot(f1,Fpskn)。psk加噪后信號(hào)頻譜39。 grid on。 f=[0:3e3:4e5]。 z=exp(w*j)。 a=.8547。 gain=.135。 %Hz(Hz==0)=10^(8)。%[1 0 p] b=[ 0 ]。 subplot(322)。 axis([0 100e6 2 2])。通過帶通濾波后輸出39。 grid on。 %解調(diào) subplot(323)。 axis([0 100e6 2 2])。 title(39。)。 gain1=。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%濾波器系數(shù) a1=[1 ]。%gain*[1 1] so=filter(b1,a1,cm)。%add gain so=somean(so)。 plot(t,so)。 title(39。)。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%比較 High=。 vt=0。 len1=length(so)。 else Vs(ii)=Low。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%一下為顯示波形subplot(325)。解調(diào)后輸出信號(hào)39。 xlabel(39。)。幅度(V)39。subplot(326)。 title(39。)。 %坐標(biāo)范圍grid on。時(shí)間 (s)39。ylabel(39。)。4PSK調(diào)制解調(diào)原理與系統(tǒng)仿真、4PSK調(diào)制原理、4PSK解調(diào)原理信號(hào)在解調(diào)之前，經(jīng)過信道會(huì)受到信道噪聲的干擾，假設(shè)信道噪聲為高斯白噪聲。close all。bit_I=bit_in(1:2:1e3)。data_I=2*bit_I+1。data_I1=repmat(data_I39。data_Q1=repmat(data_Q39。for i=1:1e4data_I2(i)=data_I1(i)。end。xrc=+*cos(pi*f)。data_Q2_rc=conv(data_Q2,xrc)/。n0=rand(size(t1))。I_rc=data_I2_rc.*cos(2*pi*f1*t1)。QPSK_rc=(sqrt(1/2).*I_rc+sqrt(1/2).*Q_rc)。FQPSK_rc=fft(QPSK_rc)。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%解調(diào)I_demo=QPSK_rc_n0.*cos(2*pi*f1*t1)。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%低通濾波I_recover=conv(I_demo,xrc)。I=I_recover(11:10010)。t2=0::。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%抽樣判決data_recover=[]。 end。for i=1:20:20000 if sum(data_recover(i:i+19))0。 bit_recover=[bit_recover 1]。 bit_recover=[bit_recover 1]。dd=2*bit_in+1。]。for i=1:2e4 ddd2(i)=ddd1(i)。 endendp=error/1000。plot(t1,QPSK_rc)。title(39。)。plot(FQPSK_rc)。4psk信號(hào)頻譜39。subplot(3,1,1)。axis([0 100 2 2])。原序列39。figure(2)subplot(2,1,1)。axis([0 100 2 2])。高斯白噪聲39。subplot(2,1,2)。axis([0 100 2 2])。4psk加白噪聲波形39。figure(3)subplot(2,1,1)。axis([0 100 2 2])。原序列39。subplot(2,1,2)。axis([0 100 2 2])。解調(diào)后序列39。運(yùn)行結(jié)果： 分析可知：解調(diào)后恢復(fù)原調(diào)制信號(hào)。它是2ASK體制的推廣，和2ASK相比，這種體制的優(yōu)點(diǎn)在于信息傳輸速率高。16進(jìn)制的正交振幅調(diào)制是一種振幅相位聯(lián)合鍵控信號(hào)。這里采用正交調(diào)幅法。 圖中串/并變換器將速率為Rb的二進(jìn)制碼元序列分為兩路,速率為Rb/,4電平信號(hào)與正交載波相乘,完成正交調(diào)制,經(jīng)24電平變換器輸出為4電平信號(hào),即M=,輸出16個(gè)信號(hào)狀態(tài),即16QAM，RS=Rb/log216=RB/4。然后經(jīng)過低通濾波器，低通濾波器LPF濾除乘法器產(chǎn)生的高頻分量，獲得有用信號(hào)，低通濾波器LPF輸出經(jīng)抽樣判決可恢復(fù)出電平信號(hào)。 圖2 16QAM正交相干解調(diào) 程序代碼：clc。N=100。 %基帶信號(hào)頻率 fs=32。 %載波頻率,為便于觀察已調(diào)信號(hào)，我們把載波頻率設(shè)的較低 Kbase=2。 %產(chǎn)生二進(jìn)制信號(hào)序列 [y,I,Q]=qam(info,Kbase,fs,fb,fc)。 y2=y。 m=fs/fb。 dt=1/fs。 subplot(221)。title(39。)subplot(222)。 title(39。)。 y=fft(y)/n。 q=find(y1e04)。 y=20*log10(y)。 f=0:f1:(length(y)1)*f1。 plot(f,y)。 title(39。)。f/fb39。 %畫出16QAM調(diào)制方式對(duì)應(yīng)的星座圖 subplot(224)。 title(39。)。 %AWGN信道信噪比 for j=1:length(SNR_in_dB) y_add_noise=awgn(y2,SNR_in_dB(j))。 %對(duì)已調(diào)信號(hào)進(jìn)行解調(diào) %傅里葉變換，求出解調(diào)信號(hào)的頻譜 n1=length(y_output)。 Fy_output=abs(Fy_output(1:fix(n1/2)))*2。 Fy_output(q)=1e04。 f1=m/n1。 numoferr=0。 end。 Pe(j)=numoferr/N。 figure。plot(y_add_noise)。已調(diào)信號(hào)加噪后輸出39。subplot(312)。title(39。)。plot(Fy_output)。解調(diào)信號(hào)頻譜39。 end。 if (temp), info(i)=0。 % 1/2的概率輸出為1 end end。 T=length(x)/fb。 nn=length(x)。 t=0:dt:Tdt。 [I,In]=two2four(I,4*m)。 [Q,Qn]=two2four(Q,4*m)。 %基帶成形濾波 I=bshape(I,fs,fb/4)。 end。 %調(diào)制 %%%%%%%%%%%%%%%%%%%%%%%%%%%%% %QAM%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%信號(hào)解調(diào) function [xn,x]=qamdet(y,fs,fb,fc)。 t=0:dt:(length(y)1)*dt。 Q=y.*sin(2*pi*fc*t)。 %設(shè)計(jì)巴特沃斯濾波器 I=filtfilt(b,a,I)。 m=4*fs/fb。 n=(.6:1:N)*m。 In=I(n)。 xn=four2two([In Qn])。 xn=[xn(1:nn/2)。 xn=xn(:)。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%基帶升余弦成形濾波器 function y=bshape(x,fs,fb,N,alfa,delay)。 delay=8。 if nargin5。 end。 N=16。 b=firrcos(N,fb,2*alfa*fb,fs)。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%二進(jìn)制轉(zhuǎn)換成四進(jìn)制 function [y,yn]=two2four(x,m)。3 2]。 ii=1。 xi=x(i:i+1)+1。 ii=ii+1。 yn=。 for i=1:m1。yn]。 y=y(:)39。 y=yn。 ymax=max(y)。 ymin=abs(ymax)。 %設(shè)置門限電平，判決 I0=find(yn )。 I1=find(yn= amp。 yn(I1)=ones(size(I1))。 yn)。 I3=find(yn=)。 %一位四進(jìn)制碼元轉(zhuǎn)換為兩位二進(jìn)制碼元 T=[0 0。1 1。 n=length(yn)。 xn(i,:)=T(yn(i)+1,:)。 xn=xn39。 xn=xn39。 N=length(x)。 n=fs/fc。 i=1。 while i = N/m。 y=2*fft(xi)/n。 i=i+1。 end。 cmax=max(abs(c))。 plot(*cos(ph),*sin(ph),39。)。 for i=1:length(c)。 a=abs(c(i))/cmax*。r*39。 end。k:39。k:39。 hold off。 axis([ ])。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%運(yùn)行結(jié)果： 課程設(shè)計(jì)心得體會(huì)參考文獻(xiàn)《通信電子電路》 謝沅清 鄧鋼編著 電子工業(yè)出版社出版；《通信原理matlab仿真教程》 趙鴻圖 茅艷主編 人民郵電出版社出版；《matlab程序設(shè)計(jì)與應(yīng)用》 劉衛(wèi)國主編 高等教育出版社出版；《matlab教程（2012a）》 張志泳 楊祖櫻編著 北京航空航天大學(xué)出版社；百度資料：40 /