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