[信息與通信]基于matlab的多進(jìn)制正交幅度調(diào)制系統(tǒng)的仿真終稿-資料下載頁(yè)
2025-08-20 17:15本頁(yè)面
【正文】 SNR_in_dB=8:2:24。 %AWGN信道信噪比f(wàn)or j=1:length(SNR_in_dB) y_add_noise=awgn(y2,SNR_in_dB(j))。 %加入不同強(qiáng)度的高斯白噪聲 y_output=qamdet(y_add_noise,fs,fb,fc)。 %對(duì)已調(diào)信號(hào)進(jìn)行解調(diào) numoferr=0。 for i=1:N if (y_output(i)~=info(i)), numoferr=numoferr+1。 end。 end。 Pe(j)=numoferr/N。 %統(tǒng)計(jì)誤碼率end。figure。semilogy(SNR_in_dB,Pe,39。red*39。)。 grid on。xlabel(39。SNR in dB39。)。 ylabel(39。Pe39。)。title(39。16QAM調(diào)制在不同信道噪聲強(qiáng)度下的誤碼率39。)。%產(chǎn)生二進(jìn)制信源隨機(jī)序列function [info]=random_binary(N) if nargin == 0, %如果沒(méi)有輸入?yún)?shù)，則指定信息序列為10000個(gè)碼元 N=10000。end。for i=1:N, temp=rand。 if (temp), info(i)=0。 % 1/2的概率輸出為0 else info(i)=1。 % 1/2的概率輸出為1 endend。function [y,I,Q]=qam(x,Kbase,fs,fb,fc)。%T=length(x)/fb。 m=fs/fb。 nn=length(x)。dt=1/fs。 t=0:dt:Tdt。%串/并變換分離出I分量、Q分量，然后再分別進(jìn)行電平映射I=x(1:2:nn1)。 [I,In]=two2four(I,4*m)。 Q=x(2:2:nn)。 [Q,Qn]=two2four(Q,4*m)。 if Kbase==2。 %基帶成形濾波 I=bshape(I,fs,fb/4)。 Q=bshape(Q,fs,fb/4)。 end。 y=I.*cos(2*pi*fc*t)Q.*sin(2*pi*fc*t)。 %調(diào)制%QAM信號(hào)解調(diào)function [xn,x]=qamdet(y,fs,fb,fc)。dt=1/fs。 t=0:dt:(length(y)1)*dt。I=y.*cos(2*pi*fc*t)。 Q=y.*sin(2*pi*fc*t)。 [b,a]=butter(2,2*fb/fs)。 %設(shè)計(jì)巴特沃斯濾波器I=filtfilt(b,a,I)。 Q=filtfilt(b,a,Q)。m=4*fs/fb。 N=length(y)/m。 n=(.6:1:N)*m。 n=fix(n)。In=I(n)。 Qn=Q(n)。 xn=four2two([In Qn])。 %I分量Q分量并/串轉(zhuǎn)換，最終恢復(fù)成碼元序列xnnn=length(xn)。 xn=[xn(1:nn/2)。xn(nn/2+1:nn)]。 xn=xn(:)。 xn=xn39。%基帶升余弦成形濾波器function y=bshape(x,fs,fb,N,alfa,delay)。%設(shè)置默認(rèn)參數(shù)if nargin6。 delay=8。 end。if nargin5。 alfa=。 end。if nargin4。 N=16。 end。b=firrcos(N,fb,2*alfa*fb,fs)。y=filter(b,1,x)。%二進(jìn)制轉(zhuǎn)換成四進(jìn)制function [y,yn]=two2four(x,m)。T=[0 1。3 2]。 n=length(x)。 ii=1。for i=1:2:n1。 xi=x(i:i+1)+1。 yn(ii)=T(xi(1),xi(2))。 ii=ii+1。end。yn=。 y=yn。 for i=1:m1。 y=[y。yn]。end。 y=y(:)39。 %；；；%四進(jìn)制轉(zhuǎn)換成二進(jìn)制function xn=four2two(yn)。y=yn。 ymin=min(y)。 ymax=max(y)。 ymax=max([ymax abs(ymin)])。ymin=abs(ymax)。 yn=(yymin)*3/(ymaxymin)。 %設(shè)置門(mén)限電平，判決I0=find(yn )。 yn(I0)=zeros(size(I0))。I1=find(yn= amp。 yn)。 yn(I1)=ones(size(I1))。I2=find(yn= amp。 yn)。 yn(I2)=ones(size(I2))*2。I3=find(yn=)。 yn(I3)=ones(size(I3))*3。%一位四進(jìn)制碼元轉(zhuǎn)換為兩位二進(jìn)制碼元T=[0 0。0 1。1 1。1 0]。 n=length(yn)。 for i=1:n。 xn(i,:)=T(yn(i)+1,:)。end。 xn=xn39。 xn=xn(:)。 xn=xn39。%畫(huà)出星座圖function c=constel(x,fs,fb,fc)。N=length(x)。 m=2*fs/fb。 n=fs/fc。 i1=mn。 i=1。 ph0=(i11)*2*pi/n。 while i = N/m。 xi=x(i1:i1+n1)。 y=2*fft(xi)/n。 c(i)=y(2)。 i=i+1。 i1=i1+m。end。 %如果無(wú)輸出，則作圖if nargout1。 cmax=max(abs(c))。 ph=(0:5:360)*pi/180。 plot(*cos(ph),*sin(ph),39。c39。)。 hold on。 for i=1:length(c)。 ph=ph0angle(c(i))。 a=abs(c(i))/cmax*。 plot(a*cos(ph),a*sin(ph),39。r*39。)。 end。 plot([ ],[0 0],39。k:39。,[0 0],[ ],39。k:39。)。 hold off。 axis equal。 axis([ ])。end。指導(dǎo)老師評(píng)語(yǔ)
電大資料相關(guān)推薦
鄂ICP備17016276號(hào)-1