【文章內(nèi)容簡(jiǎn)介】 os(2*pi*fc*t)。rt = rtmean(rt)。[f,rf] = T2F(t,rt)。[t,rt] = lpf(f,rf,2*fm)。subplot(326)plot(t,rt)。 hold on。plot(t,mt/2,39。r39。)。title(39。SSB解調(diào)信號(hào)39。)。xlabel(39。t39。)。function [out] = noise_nb(fc,B,N0,t)%output the narrow band gaussian noise sample with singlesided power spectrum N0%at carrier frequency equals fc and bandwidth euqals Bdt = t(2)t(1)。Fmx = 1/dt。 n_len = length(t)。p = N0*Fmx。rn = sqrt(p)*randn(1,n_len)。[f,rf] = T2F(t,rn)。 [t,out] = bpf(f,rf,fcB/2,fc+B/2)。例 22%FM modulation and demodulation，clear all。close all。 Kf = 5。fc = 10。T=5。dt=。t = 0:dt:T。 %信源fm= 1。%mt = cos(2*pi*fm*t) + *sin(2*pi**fm*t)。 %信源信號(hào)mt = cos(2*pi*fm*t)。 %信源信號(hào)%FM 調(diào)制A = sqrt(2)。%mti = 1/2/pi/fm*sin(2*pi*fm*t) 3/4/pi/*cos(2*pi**fm*t)。 %mt的積分函數(shù)mti = 1/2/pi/fm*sin(2*pi*fm*t) 。 %mt的積分函數(shù)st = A*cos(2*pi*fc*t + 2*pi*Kf*mti)。figure(1)subplot(311)。plot(t,st)。 hold on。plot(t,mt,39。r39。)。xlabel(39。t39。)。ylabel(39。調(diào)頻信號(hào)39。) subplot(312)[f sf] = T2F(t,st)。plot(f, abs(sf))。axis([25 25 0 3])xlabel(39。f39。)。ylabel(39。調(diào)頻信號(hào)幅度譜39。) %FM 解調(diào)for k=1:length(st)1 rt(k) = (st(k+1)st(k))/dt。endrt(length(st))=0。subplot(313)plot(t,rt)。 hold on。plot(t,A*2*pi*Kf*mt+A*2*pi*fc,39。r39。)。xlabel(39。t39。)。ylabel(39。調(diào)頻信號(hào)微分后包絡(luò)39。)例 23%數(shù)字基帶信號(hào)的功率譜密度 clear all。 close all。Ts=1。N_sample = 8。 %每個(gè)碼元的抽樣點(diǎn)數(shù)dt = Ts/N_sample。 %抽樣時(shí)間間隔N = 1000。 %碼元數(shù)t = 0:dt:(N*N_sample1)*dt。 gt1 = ones(1,N_sample)。 %NRZ非歸零波形gt2 = ones(1,N_sample/2)。 %RZ歸零波形gt2 = [gt2 zeros(1,N_sample/2)]。mt3 = sinc((t5)/Ts)。 % sin(pi*t/Ts)/(pi*t/Ts)波形，截段取10個(gè)碼元gt3 = mt3(1:10*N_sample)。d = ( sign( randn(1,N) ) +1 )/2。data = sigexpand(d,N_sample)。 %對(duì)序列間隔插入N_sample1個(gè)0st1 = conv(data,gt1)。 %Matlab自帶卷積函數(shù)st2 = conv(data,gt2)。d = 2*d1。 %變成雙極性序列data= sigexpand(d,N_sample)。 st3 = conv(data,gt3)。 [f,st1f] = T2F(t,[st1(1:length(t))])。[f,st2f] = T2F(t,[st2(1:length(t))])。[f,st3f] = T2F(t,[st3(1:length(t))])。 figure(1)subplot(321)plot(t,[st1(1:length(t))] )。gridaxis([0 20 ])。ylabel(39。單極性NRZ波形39。)。subplot(322)。plot(f,10*log10(abs(st1f).^2/T) )。gridaxis([5 5 40 10])。 ylabel(39。單極性NRZ功率譜密度(dB/Hz)39。)。 subplot(323)plot(t,[st2(1:length(t))] )。axis([0 20 ])。gridylabel(39。單極性RZ波形39。)。subplot(324)plot(f,10*log10(abs(st2f).^2/T))。axis([5 5 40 10])。gridylabel(39。單極性RZ功率譜密度(dB/Hz)39。)。 subplot(325)plot(t5,[st3(1:length(t))] )。axis([0 20 2 2])。gridylabel(39。雙極性sinc波形39。)。xlabel(39。t/Ts39。)。subplot(326)plot(f,10*log10(abs(st3f).^2/T))。axis([5 5 40 10])。gridylabel(39。sinc波形功率譜密度(dB/Hz)39。)。xlabel(39。f*Ts39。)。 function [out]=sigexpand(d,M)%將輸入的序列擴(kuò)展成間隔為N1個(gè)0的序列N = length(d)。out = zeros(M,N)。out(1,:) = d。out = reshape(out,1,M*N)。例 24%數(shù)字基帶信號(hào)接收示意 clear all。close all。 N =100。N_sample=8。 %每碼元抽樣點(diǎn)數(shù)Ts=1。dt = Ts/N_sample。t=0:dt:(N*N_sample1)*dt。 gt = ones(1,N_sample)。 %數(shù)字基帶波形d = sign(randn(1,N))。 %輸入數(shù)字序列a = sigexpand(d,N_sample)。 st = conv(a,gt)。 %數(shù)字基帶信號(hào) ht1 = gt。rt1 = conv(st,ht1)。 ht2 = 5*sinc(5*(t5)/Ts)。rt2 = conv(st,ht2)。 figure(1)subplot(321)plot( t,st(1:length(t)) )。axis([0 20 ])。 ylabel(39。輸入雙極性NRZ數(shù)字基帶波形39。)。subplot(322)stem( t,a)。axis([0 20 ])。 ylabel(39。輸入數(shù)字序列39。) subplot(323)plot( t,[0 rt1(1:length(t)1)]/8 )。axis([0 20 ])。ylabel(39。方波濾波后輸出39。)。subplot(324)dd = rt1(N_sample:N_sample:end)。ddd= sigexpand(dd,N_sample)。stem( t,ddd(1:length(t))/8 )。axis([0 20 ])。ylabel(39。方波濾波后抽樣輸出39。)。 subplot(325)plot(t5, [0 rt2(1:length(t)1)]/8 )。axis([0 20 ])。xlabel(39。t/Ts39。)。 ylabel(39。理想低通濾波后輸出39。)。subplot(326)dd = rt2(N_sample1:N_sample:end)。ddd=sigexpand(dd,N_sample)。stem( t5,ddd(1:length(t))/8 )。axis([0 20 ])。 xlabel(39。t/Ts39。)。 ylabel(39。理想低通濾波后抽樣輸出39。)。例25%部分響應(yīng)信號(hào)眼圖示意,clear all。 close all。Ts=1。N_sample=16。eye_num = 11。 N_data=1000。 dt = Ts/N_sample。t = 5*Ts:dt:5*Ts。 %產(chǎn)生雙極性數(shù)字信號(hào)d = sign(randn(1,N_data))。dd= sigexpand(d,N_sample)。%部分響應(yīng)系統(tǒng)沖擊響應(yīng)ht = sinc((t+eps)/Ts)./(1 (t+eps)./Ts)。ht( 6*N_sample+1 ) = 1。st = conv(dd,ht)。tt = 5*Ts:dt:(N_data+5)*N_sample*dtdt。 figure(1)subplot(211)。plot(tt,st)。axis([0 20 3 3])。xlabel(39。t/Ts39。)。ylabel(39。部分響應(yīng)基帶信號(hào)39。)。subplot(212)%畫眼圖ss=zeros(1,eye_num*N_sample)。ttt = 0:dt:eye_num*N_sample*dtdt。for k=5:50 ss = st(k*N_sample+1:(k+eye_num)*N_sample)。 drawnow。 plot(ttt,ss)。 hold on。end%plot(ttt,ss)。xlabel(39。t/Ts39。)。ylabel(39。部分響應(yīng)信號(hào)眼圖39。)。例26%2ASK,2PSK，clear all。close all。 A=1。fc = 2。 %2Hz。N_sample = 8。 N = 500。 %碼元數(shù)Ts = 1。 %1 baud/s dt = Ts/fc/N_sample。 %波形采樣間隔t = 0:dt:N*Tsdt。Lt = length(t)。 %產(chǎn)生二進(jìn)制信源d = sign(randn(1,N))。dd = sigexpand((d+1)/2,fc*N_sample)。gt = ones(1,fc*N_sample)。 %NRZ波形 figure(1)subplot(221)。 %輸入NRZ信號(hào)波形(單極性）d_NRZ = conv(dd,gt)。plot(t,d_NRZ(1:length(t)))。 axis([0 10 0 ])。 ylabel(39。輸入信號(hào)39。)。 subplot(222)。 %輸入NRZ頻譜[f,d_NRZf]=T2F( t,d_NRZ(1:length(t)) )。plot(f,10*log10(abs(d_NRZf).^2/T))。axis([2 2 50 10])。ylabel(39。輸入信號(hào)功率譜密度(dB/Hz)39。)。 %2ASK信號(hào)ht = A*cos(2*pi*fc*t)。s_2ask = d_NRZ(1:Lt).*ht。 subplot(223)plot(t,s_2ask)。axis([0 10 ])。 ylabel(39。2ASK39。)。 [f,s_2askf]=T2F(t,s_2ask )。subplot(224)plot(f,10*log10(abs(s_2askf).^2/T))。axis([fc4 fc+4 50 10])。ylabel(39。2ASK功率譜密度(dB/Hz)39。)。 figure(2)%2PSK信號(hào)d_2psk = 2*d_NRZ1。s_2psk = d_2psk(1:Lt).*ht。subplot(221)plot(t,s_2psk)。axis([0 10 ])。 ylabel(39。2PSK39。)。 subplot(222)[f,s_2pskf] = T2F(t,s_2psk)。plot( f,10*log10(abs(s_2pskf).^2/T) )。axis([fc4 fc+4 50 10])。ylabel(39。2PSK功率譜密度(dB/Hz)39。)。 % 2FSK% s_2fsk = Acos(2*pi*fc*t + int(2*d_NRZ1) )。 sd_2fsk = 2*d_NRZ1。 s_2fsk = A*cos(2*pi*fc*t + 2*pi*sd_2fsk(1:length(t)).*t )。subplot(223)plot(t,s_2fsk)。axis([0 10 ])。xlabel(39。t39。)。 ylabel(39。2FSK39。)subplot(224)[f,s_2fskf] = T2F(t,s_2fsk)。plot(f,10*log10(abs(s_2