【正文】 的基本數(shù)據(jù)單位是矩陣，它的指令表達式與數(shù)學、工程中常用的形式十分相似，故用MATLAB來解算問題要比用C，F(xiàn)ORTRAN等語言完成相同的事情簡捷得多，并且MATLAB也吸收了像Maple等軟件的優(yōu)點，使MATLAB成為一個強大的數(shù)學軟件。 %一個正弦周期內(nèi)的采樣點數(shù)t=(0:M1)/M/f。 %最大時刻T=ones(size(carry1))。信號加噪后波形：、2ASK解調(diào)原理與仿真2ASK有兩種基本解調(diào)方法：相干解調(diào)法（同步檢測法）和非相fc干解調(diào)法（包絡檢波法）。 %一個正弦信號周期內(nèi)的時間carry1=repmat(A*sin(2*pi*f*t),1,Q)。 %對信號進行傅里葉變換Fask=abs(Fask(1:length(Fask)/2+1))。plot(f1,Fx),title(39。解調(diào)抽樣后恢復原二進制序列。假設信道噪聲為高斯白噪聲。 if st1(n)=1。%顯示基帶波形axis([0,5,1,2])。 %加入載波1F2=st2.*s2。plot(t,fsk)。subplot(323)。at(1,j)=0。 %清除close all。 %顯示基帶信號圖像plot(t,x)。 grid on。 close all。 %載波 psk=x.*car。 title(39。 a=.8547。 %解調(diào) subplot(323)。 plot(t,so)。解調(diào)后輸出信號39。時間 (s)39。data_Q1=repmat(data_Q39。FQPSK_rc=fft(QPSK_rc)。 bit_recover=[bit_recover 1]。)。高斯白噪聲39。axis([0 100 2 2])。 圖2 16QAM正交相干解調(diào) 程序代碼：clc。 subplot(221)。 f=0:f1:(length(y)1)*f1。 %AWGN信道信噪比 for j=1:length(SNR_in_dB) y_add_noise=awgn(y2,SNR_in_dB(j))。 figure。 end。 %基帶成形濾波 I=bshape(I,fs,fb/4)。 In=I(n)。 N=16。 for i=1:m1。 yn(I1)=ones(size(I1))。 xn=xn39。 cmax=max(abs(c))。k:39。 end。 i=i+1。 xn(i,:)=T(yn(i)+1,:)。 %設置門限電平，判決 I0=find(yn )。 ii=ii+1。 if nargin5。 m=4*fs/fb。 [I,In]=two2four(I,4*m)。plot(Fy_output)。 end。 title(39。 q=find(y1e04)。 m=fs/fb。 圖中串/并變換器將速率為Rb的二進制碼元序列分為兩路,速率為Rb/,4電平信號與正交載波相乘,完成正交調(diào)制,經(jīng)24電平變換器輸出為4電平信號,即M=,輸出16個信號狀態(tài),即16QAM，RS=Rb/log216=RB/4。原序列39。figure(2)subplot(2,1,1)。plot(t1,QPSK_rc)。 end。I_rc=data_I2_rc.*cos(2*pi*f1*t1)。data_I=2*bit_I+1。)。 else Vs(ii)=Low。%gain*[1 1] so=filter(b1,a1,cm)。通過帶通濾波后輸出39。 f=[0:3e3:4e5]。 %對信號進行傅里葉變換Fpskn=abs(Fpskn(1:length(Fpskn)/2+1))。 % 方波 am=1。只有對PSK信號進行非線性變換，才能產(chǎn)生載波分量。 axis([0 200e6 2 2])。 % 方波 am=1。原基帶信號39。)。調(diào)制信號加高斯噪聲后頻%譜%39。Ffsk=fft(fsk)。title(39。plot(t,st1)。 for m=j/i*(n1)+1:j/i*n st1(m)=0。、2FSK解調(diào)原理2FSK的解調(diào)方法有包絡檢波法、相干檢波法、過零檢測法、差分檢測法，本實驗用到相干檢測法。原始信號波形與解調(diào)后的信號作對比39。 %解調(diào)信號頻譜f1=(0:length(Fx)1)*f/length(Fx)/2。 %最小時刻Tmax=max(tau)。 %正弦載波信號的頻率A=1。figure(2)subplot(111),plot(tau,ask1) %畫出調(diào)制加噪波形title(39。 %已調(diào)信號頻譜帶寬N=length(ask)。 %數(shù)字信號的帶寬f=600。任務和要求、課程設計目的?通過課程設計加深對《通信原理》課程基本內(nèi)容的理解?學習與掌握通信系統(tǒng)仿真的基本方法、課程設計要求1）完成2ASK、2FSK、2PSK、4PSK、16QAM五種調(diào)制系統(tǒng)的仿真，要求有調(diào)制、信道、解調(diào)三個部分；2）信道為高斯信道；3）系統(tǒng)參數(shù)可自行設置；4）仿真結(jié)果要求給出已調(diào)信號時域波形和頻譜，解調(diào)后信號的時域波形與頻譜。在新的版本中也加入了對C，F(xiàn)ORTRAN，C++，JAVA的支持。 %一個正弦信號周期內(nèi)的時間carry1=repmat(A*sin(2*pi*f*t),1,Q)。 %一個數(shù)字信號‘1’dsig=kron(m,T)。相干解調(diào)需要將載頻位置的已調(diào)信號頻譜重新搬回原始基帶位置，因此用相乘器與載波相乘來實現(xiàn)。 %一個碼元寬度內(nèi)的正弦載波信號Lcarry1=length(carry1)。 %已調(diào)信號頻譜f1=(0:length(Fask)1)*f/length(Fask)/2。解調(diào)信號頻譜39。2FSK調(diào)制解調(diào)原理與系統(tǒng)仿真、2FSK調(diào)制原理 2FSK是利用數(shù)字基帶信號控制在波的頻率來傳送信息。2FSK調(diào)制解調(diào)代碼：clear allclose alli=10。 st2(n)=0。subplot(422)。 %加入載波2subplot(425)。title(39。plot(t,st11)。endelsefor j=m*500+1:(m+1)*500。 %關閉窗口fs=8e5。 axis([0 2e4 2 2])。 subplot(313)。 fs=8e5。 %載波調(diào)制vn=。加噪后信號39。%BW=2(1a)/sqrt(a) p=(j^2*a^2)。 plot(t,cm)。 axis([0 6e4 ])。), axis([0 2e4 5 5]) grid on。)。,20,1)。FQPSK_rc=abs(FQPSK_rc(1:length(FQPSK_rc)/2+1))。 else data_recover_a(i:i+19)=1。subplot(3,1,3)。)。title(39。close all。plot(info)。 subplot(223)。 %加入不同強度的高斯白噪聲 y_output=qamdet(y_add_noise,fs,fb,fc)。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%一下為顯示波形代碼subplot(311)。 for i=1:N, temp=rand。 Q=bshape(Q,fs,fb/4)。 Qn=Q(n)。 end。 y=[y。 I2=find(yn= amp。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%畫出星座圖 function c=constel(x,fs,fb,fc)。 ph=(0:5:360)*pi/180。)。)。 c(i)=y(2)。 for i=1:n。 yn=(yymin)*3/(ymaxymin)。 yn(ii)=T(xi(1),xi(2))。 end。 Q=filtfilt(b,a,Q)。 %串/并變換分離出I分量、Q分量，然后再分別進行電平映射 I=x(1:2:nn1)。subplot(313)。 for i=1:N if (y_output(i)~=info(i)), numoferr=numoferr+1。 constel(y1,fs,fb,fc)。 y=abs(y(1:fix(n/2)))*2。 %備份信號，供后續(xù)仿真用 T=length(info)/fb。16QAM正交調(diào)制的原理如下圖1所示。title(39。)。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%一下均為波形顯示部分figure(1)subplot(3,1,2)。for i=1:20:10000 data_recover=[data_recover I(i:1:i+19) Q(i:1:i+19)]。 %產(chǎn)生高斯噪聲%%%%%%%%%%%