【正文】 filter an input at frequency domain by an ideal bandpass filter%Inputs:% f: frequency samples% sf: input data spectrum samples% B1: bandpass39。T = 1/df。bf2 = floor( length(f)/2 ) bf。[t,st]=F2T(f,yf)。 %時(shí)間采樣間隔 fm=1。 mt = sqrt(2)*cos(2*pi*fm*t)。%noise = noise_nb(fc,B,N0,t)。 %畫(huà)出DSB信號(hào)波形plot(t,mt,39。DSB調(diào)制信號(hào)39。)。[t,rt] = lpf(f,rf,2*fm)。r39。)。psf = (abs(sf).^2)/T。DSB信號(hào)功率譜39。)。T = 1/df。yf=hf.*sf。clear all。 %載波中心頻率T=5。 %白噪單邊功率譜密度 %AM modulationA=2。 %窄帶高斯噪聲產(chǎn)生%s_am = s_am + noise。r39。)。%AM demodulationrt = s_am.*cos(2*pi*fc*t)。 %低通濾波subplot(312)plot(t,rt)。)。xlabel(39。plot(f,psf)。)。例 19%顯示模擬調(diào)制的波形及解調(diào)方法SSB，%信源close all。 %信源最高頻率fc=10。 %信源%N0 = 。%s_ssb=s_ssb+noise。r39。)。 %SSB demodulationrt = s_ssb.*cos(2*pi*fc*t)。 subplot(312)plot(t,rt)。)。xlabel(39。plot(f,psf)。)。例 20%顯示模擬調(diào)制的波形及解調(diào)方法VSB，%信源close all。 %信源最高頻率fc=20。 %信源%VSB modulations_vsb = mt.*cos(2*pi*fc*t)。figure(1)subplot(311)plot(t,s_vsb)。)。xlabel(39。 [f,rf] = T2F(t,rt)。plot(t,mt/2,39。相干解調(diào)后的信號(hào)波形與輸入信號(hào)的比較39。)subplot(313)[f,sf]=T2F(t,s_vsb)。title(39。f39。T = 1/df。f1 = bf1 + floor( length(f)/2 ) 。hf(f2)=1。hf(f4)=1。例 21%顯示模擬調(diào)制的波形及解調(diào)方法AM、DSB、SSB, %信源close all。fc=10。 %AM modulationA=2。s_am = s_am + noise。r39。[f,rf] = T2F(t,rt)。)。subplot(322)plot(t,rt)。)。xlabel(39。B=2*fm。hold on。title(39。t39。[f,rf] = T2F(t,rt)。plot(t,mt/2,39。DSB解調(diào)信號(hào)39。)。s_ssb=s_ssb+noise。)。%SSB demodulationrt = s_ssb.*cos(2*pi*fc*t)。subplot(326)plot(t,rt)。)。xlabel(39。Fmx = 1/dt。[f,rf] = T2F(t,rn)。 Kf = 5。t = 0:dt:T。 %信源信號(hào)%FM 調(diào)制A = sqrt(2)。figure(1)subplot(311)。r39。)。plot(f, abs(sf))。ylabel(39。subplot(313)plot(t,rt)。)。ylabel(39。Ts=1。 %碼元數(shù)t = 0:dt:(N*N_sample1)*dt。mt3 = sinc((t5)/Ts)。 %對(duì)序列間隔插入N_sample1個(gè)0st1 = conv(data,gt1)。 st3 = conv(data,gt3)。 figure(1)subplot(321)plot(t,[st1(1:length(t))] )。)。 ylabel(39。axis([0 20 ])。subplot(324)plot(f,10*log10(abs(st2f).^2/T))。)。雙極性sinc波形39。)。sinc波形功率譜密度(dB/Hz)39。)。out = reshape(out,1,M*N)。N_sample=8。 gt = ones(1,N_sample)。 %數(shù)字基帶信號(hào) ht1 = gt。 figure(1)subplot(321)plot( t,st(1:length(t)) )。)。輸入數(shù)字序列39。方波濾波后輸出39。stem( t,ddd(1:length(t))/8 )。)。t/Ts39。)。axis([0 20 ])。 ylabel(39。 close all。 N_data=1000。dd= sigexpand(d,N_sample)。tt = 5*Ts:dt:(N_data+5)*N_sample*dtdt。xlabel(39。部分響應(yīng)基帶信號(hào)39。for k=5:50 ss = st(k*N_sample+1:(k+eye_num)*N_sample)。end%plot(ttt,ss)。ylabel(39。close all。N_sample = 8。 %波形采樣間隔t = 0:dt:N*Tsdt。gt = ones(1,fc*N_sample)。 axis([0 10 0 ])。 subplot(222)。ylabel(39。s_2ask = d_NRZ(1:Lt).*ht。2ASK39。axis([fc4 fc+4 50 10])。 figure(2)%2PSK信號(hào)d_2psk = 2*d_NRZ1。 ylabel(39。plot( f,10*log10(abs(s_2pskf).^2/T) )。)。subplot(223)plot(t,s_2fsk)。)。plot(f,10*log10(abs(s_2fskf).^2/T))。)。例 27%QPSK amp。Ts= 1。 dt = 1/fc/N_sample。 %功率譜密度1py2f = zeros(1,length(t))。 %QPSK調(diào)制 s1 = sigexpand(d1,fc*N_sample)。 s1 = b1(1:length(s1))。 st_oqpsk= s1.*cos(2*pi*fc*t) s2_delay.*sin(2*pi*fc*t)。 [t y2] = bpf(f,y2f,fc1/Ts,fc+1/Ts)。t39。)。)。)。 axis([5 15 ])。 %經(jīng)過(guò)非線性電路 y1 = *tanh(2*y1)。 py1f = py1f + abs(y1f).^2/T。 %OQPSK 100段功率譜密度平均 subplot(223)。)。title(39。 subplot(224)plot(f,10*log10(py2f))。ylabel(39。經(jīng)過(guò)非線性電路后的功率譜密度39。y=*tanh(2*x)。)。A=。 endend figure(1)plot(x,yu,39。u Law39。ylabel(39。yy=[1,7/8,6/8,5/8,4/8,3/8,2/8,1/8,1/8,2/8,3/8,4/8,5/8,6/8,7/8,1]。stem(xx,yy,39。u律壓縮特性39。 figure(2)plot(x,ya,39。A Law39。ylabel(39。yy=[1,7/8,6/8,5/8,4/8,3/8,2/8,1/8,1/8,2/8,3/8,4/8,5/8,6/8,7/8,1]。stem(xx,yy,39。A律壓縮特性39。例29%show the pcm encode and decodeclear allclose allt=0::10。 v=1。 yy = pcm_decode(y,v)。 drawnow subplot(211) plot(t,x)。 subplot(212) plot(t,yy) title(39。endfigure(2)plot(vm1,snrq)。%4096x4096for i=1:n if x(i)0 out(i,1)=1。out(i,3)=0。 elseif 32=abs(x(i)) amp。step=2。out(i,3)=1。 elseif 128=abs(x(i)) amp。step=8。out(i,3)=0。 elseif 512=abs(x(i)) amp。step=32。out(i,3)=1。 elseif 2048=abs(x(i)) amp。step=128。out(i,4)=1。 else tmp = floor( (abs(x(i))st)/step )。,1,8*n)。slot(1) = 0。slot(5) = 256。 step(1) = 2。step(5) = 16。 for i=1:n/8 ss = 2*in(i,1)1。 out(i)=ss*(st+dt)/4096*v。N_sample = 8。 gt1 = ones(1,N_sample)。 % sin(pi*t/Ts)/(pi*t/Ts)波形，截段取10個(gè)碼元gt3 = mt3(1:10*N_sample)。st2 = conv(data,gt2)。 xt = st1。[t1,yt1] = F2T(f,yf1)。%相頻失真、群時(shí)延無(wú)失真信道%hf3 = exp(j*pi*f+j**pi)。 yf3 = xf.*hf3。[t4,yt4]=F2T(f,yf4)。)。grid on。)。ylabel(39。幅頻失真信道39。f39。xlabel(39。plotyy(f,abs(hf3),f,angle(hf3)/pi)。title(39。subplot(222)。)。ylabel(39。相頻失真、群時(shí)延失真信道39。f39。axis([0,20, ])。)。N_sample = 8。dLen = length(t)。 %對(duì)序列間隔插入N_sample1個(gè)0st1 = conv(data,gt1)。hf = m(1)*exp(j*2*pi*f*tao(1)) +m(2)* exp(j*2*pi*f*tao(2))+... m(3)* exp(j*2*pi*f*tao(3))。LineWidth39。輸入信號(hào)39。,2)。)。subplot(222)。axis([2 2 0 60])。subplot(224)。axis([2 2 0 60])。xlabel(39。LineWidth39。信道幅頻特性39。)。)。axis([2 2 1 1])。... 0 1 1 1 0 1 0。0 0 0 1。0 1 0 1。1 0 0 1。1 1 0 1。 E = dec2bin([0:2^71],7)48。for k=1:8 e(k,:)=find(s==k1)。 %對(duì)應(yīng)S＝010的Ee3 = dec2bin( e(4,:)1,7 )48。 %對(duì)應(yīng)S＝110的Ee7 = dec2bin( e(8,:)1,7 )48。EbN0dB = 0:6。 %卷積碼的結(jié)構(gòu)trellis = poly2trellis(3,[7,5])。 d = [d1 zeros(1,2)]。 %誤碼計(jì)數(shù) error(k) = error(k) + sum( abs( dd(1:length(d1)) d1) ) n = n+1。for k=1:length(din) d(k) = ( curState+1,din(k)+1 )。out = reshape(out39。numStates = 。next_metric = ones(1,numStates)。 dout = 2*dout1。 if next_metric(nextState)0 %判斷下一狀態(tài)是否新寫(xiě)入？ next_metric(nextState) = x。 end end end end cur_metric = next_metric。例 34%m序列發(fā)生器及其自相關(guān) clear all。% state=1000L = 10