【正文】 hold off。r*39。 plot(*cos(ph),*sin(ph),39。 y=2*fft(xi)/n。 N=length(x)。 n=length(yn)。 yn)。 ymin=abs(ymax)。yn]。 xi=x(i:i+1)+1。 b=firrcos(N,fb,2*alfa*fb,fs)。 delay=8。 xn=four2two([In Qn])。 %設(shè)計巴特沃斯濾波器 I=filtfilt(b,a,I)。 end。 t=0:dt:Tdt。 if (temp), info(i)=0。)。plot(y_add_noise)。 numoferr=0。 %對已調(diào)信號進行解調(diào) %傅里葉變換，求出解調(diào)信號的頻譜 n1=length(y_output)。 %畫出16QAM調(diào)制方式對應(yīng)的星座圖 subplot(224)。 plot(f,y)。 y=fft(y)/n。title(39。 y2=y。N=100。這里采用正交調(diào)幅法。解調(diào)后序列39。axis([0 100 2 2])。subplot(2,1,2)。原序列39。plot(FQPSK_rc)。 endendp=error/1000。 bit_recover=[bit_recover 1]。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%抽樣判決data_recover=[]。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%解調(diào)I_demo=QPSK_rc_n0.*cos(2*pi*f1*t1)。n0=rand(size(t1))。for i=1:1e4data_I2(i)=data_I1(i)。bit_I=bit_in(1:2:1e3)。ylabel(39。 title(39。 xlabel(39。 len1=length(so)。 title(39。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%濾波器系數(shù) a1=[1 ]。 axis([0 100e6 2 2])。 axis([0 100e6 2 2])。 gain=.135。 grid on。)。 pskn=(psk+noise)。 noise=vn*(randn(size(t)))。 %方波幅度phi=0。 %抽樣頻率 fm=20e3。2PSK解調(diào)原理框圖由于PSK信號的功率譜中五載波分量，所以必須采用相干解調(diào)的方式。 %顯示psk調(diào)制信號頻譜plot(f,Fpsk)。 %已調(diào)信號頻譜帶寬subplot(312)。 %坐標范圍title(39。 %方波幅度phi=0。 %抽樣頻率 fm=20e3。plot(t,st1)。at(1,j)=1。title(39。title(39。fsk信號加高斯白噪聲輸出信號39。2FSK信號39。)。 %顯示載波1與基帶信號相乘波形plot(t,F1)。)。plot(t,st2)。 %對信號進行傅里葉變換%%%%%%%%%%%%%%%%%%%%已調(diào)信號加噪聲后的頻譜Ffsk1=abs(Ffsk1(1:length(Ffsk1)/2+1))。 else st2(n)=1。 %載波2帶寬%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%產(chǎn)生基帶信號st1=t。 %基帶信號碼元數(shù)j=5000。 其中，為的反碼，即若，則；若，則。例如，1碼用頻率f1來傳輸，0碼用頻率f2來傳輸，而其振幅和初始相位不變。抽樣判決后的數(shù)字信號波形39。)。 %運用模擬濾波器對信號進行濾波Fx=fft(x)。 %已調(diào)信號頻譜帶寬、%Butterworth 模擬低通濾波器設(shè)計%獲得Butterworth 模擬低通原型濾波器的階數(shù)及3dB 截止頻率[Ord,omega_c]=buttord(4*pi*f*,4*pi*f*,2,30,39。 %長度tau=(0:N1)/(M1)/f。 %一個碼元寬度內(nèi)的信號長度carry2=kron(ones(size(m)),carry1)。 %序列的長度F=150。為確保無失真還原信號，必須在接收端提供一個與調(diào)制載波嚴格同步的本地載波，這是整個解調(diào)過程能否順利完好進行的關(guān)鍵。已調(diào)信號頻譜39。 %數(shù)字信號波形subplot(2,2,1) %子圖分割plot(tau,dsig),title(39。 %對信號進行傅里葉變換Fask=abs(Fask(1:length(Fask)/2+1))。 %一個碼元寬度內(nèi)的正弦載波信號Lcarry1=length(carry1)。 %數(shù)字信號序列Lm=length(m)?？梢灾苯诱{(diào)用,用戶也可以將自己編寫的實用程序?qū)氲組ATLAB函數(shù)庫中方便自己以后調(diào)用，此外許多的MATLAB愛好者都編寫了一些經(jīng)典的程序，用戶可以直接進行下載就可以用。通常使用鍵控法來實現(xiàn)數(shù)字調(diào)制，比如對載波的振幅、頻率和相位進行鍵控。其次，還有4PSK、16QAM等調(diào)制方式。2ASK調(diào)制解調(diào)原理與系統(tǒng)仿真、2ASK調(diào)制原理與仿真數(shù)字幅度調(diào)制又稱幅度鍵控（ASK），二進制幅度鍵控記作2ASK。 %序列的長度F=150。 %一個碼元寬度內(nèi)的信號長度carry2=kron(ones(size(m)),carry1)。 %已調(diào)信號頻譜f1=(0:length(Fask)1)*f/length(Fask)/2。數(shù)字信號波形39。)。解調(diào)過程如圖5所示。 %數(shù)字信號的帶寬f=600。 %載波信號ask=kron(m,carry1)。 %時間Tmin=min(tau)。s39。 %對信號進行傅里葉變換Fx=abs(Fx(1:length(Fx)/2+1))。 %解調(diào)信號頻譜grid on %添加網(wǎng)格axis([Tmin Tmax ]) %設(shè)置坐標范圍th=。)%畫出抽樣判決%后%%的數(shù)字%%信號波形axis([min(n_tau) max(n_tau) ]) %設(shè)置坐標范圍grid on subplot(3,2,5) %子圖分割plot(tau,dsig),title(39。故其表示式為： 式中，假設(shè)碼元的初始相位分別為和；和為兩個不同的碼元的角頻率；幅度為A為一常數(shù)，表示碼元的包絡(luò)為矩形脈沖。2FSK信號的頻譜特性： 由于相位離散的2FSK信號可看成是兩個2ASK信號之和，所以，這里可以直接應(yīng)用2ASK信號的頻譜分析結(jié)果，比較方便，即 2FSK信號帶寬為 式中，是基帶信號的帶寬。a=round(rand(1,i))。for n=1:10 if a(n)1。 endend。%%%%%%%%%%%%%已調(diào)信號加噪聲信號后的頻譜帶寬f3=(0:length(Ffsk1)1)*fm/length(Ffsk1)/2 subplot(421)。title(39。 %顯示載波波形subplot(424),plot(s2)。title(39。fsk=F1+F2。)subplot(428),plot(f,Ffsk),title(39。)subplot(322),plot(f3,Ffsk1),title(39。fsk與載波1相乘后波形39。fsk與載波2相乘后波形39。endendend。title(39。 %基帶頻率 n=2*(6*fs/fm)。 x = A * cos(twopi_fc_t + phi)?；鶐盘?9。 %顯示psk調(diào)制信號 plot(t,psk)。 title(39。在相干解調(diào)中，如何得到同頻同相的本地載波是個關(guān)鍵問題。 %基帶頻率 n=2*(6*fs/fm)。 x = A * cos(twopi_fc_t + phi)。 %產(chǎn)生噪音 subplot(211)。 %調(diào)制后加噪 Fpskn=fft(pskn)。 grid on。 %畫格子%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%帶通濾波 fBW=40e3。 Hz=gain*(z+1).*(z1)./(z.^2(p))。 title(39。 grid on。%(z(p)) b1=[ ]。通過低通濾波器后輸出39。 for ii=1:len1 if so(ii) = vt Vs(ii)=High。時間 (s)39。原始基帶信號39。幅度(V)39。bit_Q=bit_in(2:2:1e3)。da