【正文】 end。 axis equal。)。,[0 0],[ ],39。 plot([ ],[0 0],39。)。 plot(a*cos(ph),a*sin(ph),39。 ph=ph0angle(c(i))。 hold on。c39。 ph=(0:5:360)*pi/180。%如果無輸出，則作圖 if nargout1。 i1=i1+m。 c(i)=y(2)。 xi=x(i1:i1+n1)。 ph0=(i11)*2*pi/n。 i1=mn。 m=2*fs/fb。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%畫出星座圖 function c=constel(x,fs,fb,fc)。 xn=xn(:)。 end。 for i=1:n。1 0]。0 1。 yn(I3)=ones(size(I3))*3。 yn(I2)=ones(size(I2))*2。 I2=find(yn= amp。 yn)。 yn(I0)=zeros(size(I0))。 yn=(yymin)*3/(ymaxymin)。 ymax=max([ymax abs(ymin)])。 ymin=min(y)。 %；；； %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%四進(jìn)制轉(zhuǎn)換成二進(jìn)制 function xn=four2two(yn)。 end。 y=[y。 y=yn。 end。 yn(ii)=T(xi(1),xi(2))。 for i=1:2:n1。 n=length(x)。 T=[0 1。 y=filter(b,1,x)。 end。 if nargin4。 alfa=。 end。 %設(shè)置默認(rèn)參數(shù) if nargin6。 xn=xn39。xn(nn/2+1:nn)]。 %I分量Q分量并/串轉(zhuǎn)換，最終恢復(fù)成碼元序列xn nn=length(xn)。 Qn=Q(n)。 n=fix(n)。 N=length(y)/m。 Q=filtfilt(b,a,Q)。 [b,a]=butter(2,2*fb/fs)。 I=y.*cos(2*pi*fc*t)。 dt=1/fs。 y=I.*cos(2*pi*fc*t)Q.*sin(2*pi*fc*t)。 Q=bshape(Q,fs,fb/4)。 if Kbase==2。 Q=x(2:2:nn)。 %串/并變換分離出I分量、Q分量，然后再分別進(jìn)行電平映射 I=x(1:2:nn1)。 dt=1/fs。 m=fs/fb。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%function [y,I,Q]=qam(x,Kbase,fs,fb,fc)。 % 1/2的概率輸出為0 else info(i)=1。 for i=1:N, temp=rand。)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%一下為需要用到的M文件%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%產(chǎn)生二進(jìn)制信源隨機(jī)序列 function [info]=random_binary(N) if nargin == 0, %如果沒有輸入?yún)?shù)，則指定信息序列為10000個(gè)碼元 N=10000。title(39。subplot(313)。解調(diào)輸出信號39。plot(y_output)。)。title(39。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%一下為顯示波形代碼subplot(311)。 %統(tǒng)計(jì)誤碼率 end。 end。 for i=1:N if (y_output(i)~=info(i)), numoferr=numoferr+1。 f2=0:f1:(length(Fy_output)1)*f1。 Fy_output=20*log10(Fy_output)。 q=find(Fy_output1e04)。 Fy_output=fft(y_output)/n1。 %加入不同強(qiáng)度的高斯白噪聲 y_output=qamdet(y_add_noise,fs,fb,fc)。 SNR_in_dB=8:2:24。16QAM調(diào)制方式對應(yīng)的星座圖39。 constel(y1,fs,fb,fc)。)。 xlabel(39。已調(diào)信號頻譜39。 grid on。 subplot(223)。 f1=m/n。 y(q)=1e04。 y=abs(y(1:fix(n/2)))*2。 %%%%%%%%%%%%%%%%%%傅里葉變換，求出已調(diào)信號的頻譜 n=length(y)。已調(diào)信號輸出波形39。 %顯示的已調(diào)信號 %由于測試信號碼元數(shù)量為10000個(gè)，在這里我們只顯示其總數(shù)的1/10 plot(t(1:1000),y(1:1000),t(1:1000),I(1:1000),t(1:1000),Q(1:1000)[0 20],[0 0])。原序列39。plot(info)。 t=0:dt:Tdt。nn=length(info)。 %備份信號，供后續(xù)仿真用 T=length(info)/fb。 %對基帶信號進(jìn)行16QAM調(diào)制 y1=y。 % Kbase=2,基帶經(jīng)成形濾波器濾波后，再進(jìn)行調(diào)制 info=random_binary(N)。 %抽樣頻率 fc=4。 %設(shè)定碼元數(shù)量 fb=1。close all。16QAM正交相干解調(diào)如圖2所示。 、16QAM解調(diào)原理16QAM信號采取正交相干解調(diào)的方法解調(diào)，解調(diào)器首先對收到的16QAM信號進(jìn)行正交相干解調(diào)，一路與tcωcos相乘，一路與tcωsin相乘。16QAM正交調(diào)制的原理如下圖1所示。16QAM的產(chǎn)生有2種方法：（1）正交調(diào)幅法，它是有2路正交的四電平振幅鍵控信號疊加而成；（2）復(fù)合相移法：它是用2路獨(dú)立的四相位移相鍵控信號疊加而成。 正交幅度調(diào)制是利用多進(jìn)制振幅鍵控（MASK）和正交載波調(diào)制相結(jié)合產(chǎn)生的。16QAM調(diào)制解調(diào)原理與系統(tǒng)仿真、16QAM調(diào)制原理 16QAM是用兩路獨(dú)立的正交4ASK信號疊加而成，4ASK是用多電平信號去鍵控載波而得到的信號。)。title(39。plot(t2,data_recover_a)。)。title(39。plot(t2,ddd2)。)。title(39。plot(t1,QPSK_rc_n0)。)。title(39。plot(n0)。)。title(39。plot(t2,ddd2)。)。title(39。subplot(3,1,3)。4psk信號波形39。axis([0 100 2 2])。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%一下均為波形顯示部分figure(1)subplot(3,1,2)。endfor i=1:1e3 if bit_recover(i)~=ddd(i) error=error+1。ddd1=repmat(ddd,20,1)。ddd=[dd39。 endenderror=0。 else data_recover_a(i:i+19)=1。 data_recover_a(i:i+19)=1。bit_recover=[]。for i=1:20:10000 data_recover=[data_recover I(i:1:i+19) Q(i:1:i+19)]。t3=0::。Q=Q_recover(11:10010)。Q_recover=conv(Q_demo,xrc)。Q_demo=QPSK_rc_n0.*sin(2*pi*f1*t1)。FQPSK_rc=abs(FQPSK_rc(1:length(FQPSK_rc)/2+1))。QPSK_rc_n0=QPSK_rc+n0。Q_rc=data_Q2_rc.*sin(2*pi*f1*t1)。 %產(chǎn)生高斯噪聲%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%調(diào)制f1=1。t1=0::1e3+。data_I2_rc=conv(data_I2,xrc)/。f=0::1。data_Q2(i)=data_Q1(i)。,20,1)。,20,1)。data_Q=2*bit_Q+1。bit_Q=bit_in(2:2:1e3)。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%產(chǎn)生雙極性碼bit_in=randint(1e3,1,[0 1])。仿真代碼：clear all。運(yùn)行結(jié)果： 分析可知：解調(diào)后恢復(fù)原始信號。幅度(V)39。)。xlabel(39。 axis([0 2e4 2 2])。原始基帶信號39。 plot (t,x)。)。 %坐標(biāo)軸標(biāo)記 ylabel(39。時(shí)間 (s)39。), axis([0 2e4 5 5]) grid on。 plot (t,Vo), title(39。 end end Vo=Vs。 for ii=1:len1 if so(ii) = vt Vs(ii)=High。 %設(shè)立比較標(biāo)準(zhǔn) error=0。 Low=。 grid on。通過低通濾波器后輸出39。 axis([0 6e4 ])。%removes DC ponent subplot(324)。 so=so*10。%(z(p)) b1=[ ]。%gain=(1p)/2 Hz1=gain1*(z+1)./(z(p))。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%低通濾波器 p=。通過相乘器后輸出39。 grid