【正文】 ta_Q2(i)=data_Q1(i)。 %產(chǎn)生高斯噪聲%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%調(diào)制f1=1。Q_demo=QPSK_rc_n0.*sin(2*pi*f1*t1)。for i=1:20:10000 data_recover=[data_recover I(i:1:i+19) Q(i:1:i+19)]。 endenderror=0。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%一下均為波形顯示部分figure(1)subplot(3,1,2)。title(39。)。plot(t1,QPSK_rc_n0)。title(39。)。16QAM正交調(diào)制的原理如下圖1所示。 %設(shè)定碼元數(shù)量 fb=1。 %備份信號，供后續(xù)仿真用 T=length(info)/fb。原序列39。 y=abs(y(1:fix(n/2)))*2。 grid on。 constel(y1,fs,fb,fc)。 Fy_output=fft(y_output)/n1。 for i=1:N if (y_output(i)~=info(i)), numoferr=numoferr+1。title(39。subplot(313)。 % 1/2的概率輸出為0 else info(i)=1。 %串/并變換分離出I分量、Q分量，然后再分別進(jìn)行電平映射 I=x(1:2:nn1)。 y=I.*cos(2*pi*fc*t)Q.*sin(2*pi*fc*t)。 Q=filtfilt(b,a,Q)。 %I分量Q分量并/串轉(zhuǎn)換，最終恢復(fù)成碼元序列xn nn=length(xn)。 end。 y=filter(b,1,x)。 yn(ii)=T(xi(1),xi(2))。 end。 yn=(yymin)*3/(ymaxymin)。 yn(I2)=ones(size(I2))*2。 for i=1:n。 m=2*fs/fb。 c(i)=y(2)。c39。)。 axis equal。)。 plot(a*cos(ph),a*sin(ph),39。 ph=(0:5:360)*pi/180。 xi=x(i1:i1+n1)。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%畫出星座圖 function c=constel(x,fs,fb,fc)。1 0]。 I2=find(yn= amp。 ymax=max([ymax abs(ymin)])。 y=[y。 for i=1:2:n1。 end。 %設(shè)置默認(rèn)參數(shù) if nargin6。 Qn=Q(n)。 [b,a]=butter(2,2*fb/fs)。 Q=bshape(Q,fs,fb/4)。 dt=1/fs。 for i=1:N, temp=rand。解調(diào)輸出信號39。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%一下為顯示波形代碼subplot(311)。 f2=0:f1:(length(Fy_output)1)*f1。 %加入不同強度的高斯白噪聲 y_output=qamdet(y_add_noise,fs,fb,fc)。)。 subplot(223)。 %%%%%%%%%%%%%%%%%%傅里葉變換，求出已調(diào)信號的頻譜 n=length(y)。plot(info)。 %對基帶信號進(jìn)行16QAM調(diào)制 y1=y。close all。16QAM的產(chǎn)生有2種方法：（1）正交調(diào)幅法，它是有2路正交的四電平振幅鍵控信號疊加而成；（2）復(fù)合相移法：它是用2路獨立的四相位移相鍵控信號疊加而成。title(39。plot(t2,ddd2)。)。title(39。subplot(3,1,3)。endfor i=1:1e3 if bit_recover(i)~=ddd(i) error=error+1。 else data_recover_a(i:i+19)=1。t3=0::。FQPSK_rc=abs(FQPSK_rc(1:length(FQPSK_rc)/2+1))。t1=0::1e3+。,20,1)。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%產(chǎn)生雙極性碼bit_in=randint(1e3,1,[0 1])。)。 plot (t,x)。), axis([0 2e4 5 5]) grid on。 %設(shè)立比較標(biāo)準(zhǔn) error=0。 axis([0 6e4 ])。%gain=(1p)/2 Hz1=gain1*(z+1)./(z(p))。 plot(t,cm)。 plot(t,fpskn)。%BW=2(1a)/sqrt(a) p=(j^2*a^2)。)。加噪后信號39。 axis([0 .2e3 1 1])。 %載波調(diào)制vn=。 A=1。 fs=8e5。 2PSK信號屬于DSB信號，它的解調(diào)，不再能采用包絡(luò)檢測的方法，只能進(jìn)行相干解調(diào)。 subplot(313)。 %已調(diào)信號頻譜f=(0:length(Fpsk)1)*fs/length(Fpsk)/2。 axis([0 2e4 2 2])。 A=1。 %關(guān)閉窗口fs=8e5。)subplot(326)。endelsefor j=m*500+1:(m+1)*500。plot(t,st22)。plot(t,st11)。title(39。title(39。s2*st239。 %加入載波2subplot(425)。載波信號139。subplot(422)。 %fsk信號加高斯白噪聲Ffsk1=fft(fsk1)。 st2(n)=0。 %載波1帶寬B2=2*f2。2FSK調(diào)制解調(diào)代碼：clear allclose alli=10。 (a) (b) 2FSK信號產(chǎn)生原理圖 由鍵控法產(chǎn)生原理可知，一位相位離散的2FSK信號可看成不同頻率交替發(fā)送的兩個2ASK信號之和，即 其中是脈寬為的矩形脈沖表示的NRZ數(shù)字基帶信號。2FSK調(diào)制解調(diào)原理與系統(tǒng)仿真、2FSK調(diào)制原理 2FSK是利用數(shù)字基帶信號控制在波的頻率來傳送信息。 %抽樣判決后的信號對應(yīng)的時間subplot(3,2,4) %子圖分割plot(n_tau,y_value), title(39。解調(diào)信號頻譜39。 %獲得濾波器傳遞函數(shù)%濾波x=lsim(h,sig_mul,tau)。 %已調(diào)信號頻譜f1=(0:length(Fask)1)*f/length(Fask)/2。 %已調(diào)信號與載波信號相乘N=length(ask)。 %一個碼元寬度內(nèi)的正弦載波信號Lcarry1=length(carry1)。 %數(shù)字信號序列Lm=length(m)。相干解調(diào)需要將載頻位置的已調(diào)信號頻譜重新搬回原始基帶位置，因此用相乘器與載波相乘來實現(xiàn)。plot(f1,Fask),title(39。 %一個數(shù)字信號‘1’dsig=kron(m,T)。 %加高斯白噪聲后的已調(diào)信號Fask=fft(ask)。 %一個正弦信號周期內(nèi)的時間carry1=repmat(A*sin(2*pi*f*t),1,Q)。 %清除空間變量m=[1 0 1 1 0 1 1 0]。在新的版本中也加入了對C，F(xiàn)ORTRAN，C++，JAVA的支持。這種用數(shù)字基帶信號控制載波，把數(shù)字基帶信號變換為數(shù)字帶通信號的過程稱為數(shù)字調(diào)制。任務(wù)和要求、課程設(shè)計目的?通過課程設(shè)計加深對《通信原理》課程基本內(nèi)容的理解?學(xué)習(xí)與掌握通信系統(tǒng)仿真的基本方法、課程設(shè)計要求1）完成2ASK、2FSK、2PSK、4PSK、16QAM五種調(diào)制系統(tǒng)的仿真，要求有調(diào)制、信道、解調(diào)三個部分；2）信道為高斯信道；3）系統(tǒng)參數(shù)可自行設(shè)置；4）仿真結(jié)果要求給出已調(diào)信號時域波形和頻譜，解調(diào)后信號的時域波形與頻譜。2ASK是利用代表數(shù)字信息“0”或“1”的基帶矩形脈沖去鍵控一個連續(xù)的載波，使載波時斷時續(xù)地輸出。 %數(shù)字信號的帶寬f=600。 %載波信號ask=kron(m,carry1)。 %已調(diào)信號頻譜帶寬N=length(ask)。) %畫出數(shù)字信號的波形圖grid on axis([Tmin Tmax ]) %設(shè)置坐標(biāo)范圍subplot(2,2,2) %子圖分割plot(tau,carry2),title (39。figure(2)subplot(111),plot(tau,ask1) %畫出調(diào)制加噪波形title(39。S( t )定時脈沖 帶通 相乘 低通 抽樣判決eo( t )cosωc t2fs圖5 相干解調(diào)帶通 整流 低通 抽樣判決判決定時脈沖 S( t ) eo( t ) abcd