【正文】 v(k)= yh(k)。 end sign=yh(k)。%當(dāng)前輸入為 0，則計(jì)數(shù)器清零 end end re=[xn39。,yh39。,B39。 %HDB3解碼 input=yh。 sign=0。 end sign=input(k)。 subplot(3,1,1)。axis([1 length(xn) 2 2])。 ylabel(39。)。HDB3 碼的編碼前圖形 39。 subplot(3,1,2)。axis([1 length(xn) 2 2])。 ylabel(39。)。HDB3 碼的編碼后（解碼前）圖形 39。 subplot(3,1,3)。axis([1 length(xn) 2 2])。 ylabel(39。)。HDB3 碼的解碼后圖形 39。 code1=yh。 ASK 的調(diào)制與解調(diào) source 函數(shù) function sendSignal=source(n,N) sendSignal=randint(1,n) bit=[]。 else bit1=ones(1,N)。 end figure(1) plot(1:length(bit),bit),title(39。),grid on。 end askModu 函數(shù)源程序 function transmittedSignal=askModu(signal,bitRate,fc,N) %signal=[1 0 1 0 1 0 0 1]。 % fc=1000000。 t=linspace(0,1/bitRate,N)。 transmittedSignal=[]。 end figure(2) plot(1:length(transmittedSignal),transmittedSignal)。ASK調(diào)制波形 39。grid on。 F=fft(transmittedSignal)。ASK仿真頻譜分析 39。 grid on。title(39。)。 End gussian 函數(shù)源程序 function signal=gussian(transmittedSignal,noise) signal=sqrt(2)*transmittedSignal。 figure(5) plot(1:length(signal),signal)。包含噪聲的波形 39。 end CheckRatepe 函數(shù)源程序 function PeWrong=CheckRatePe(signal1,signal2,s) rights=0。 for ki=1:s2 if(signal1(ki)==signal2(ki)) rights=rights+1。 end end PeWrong=wrongs/(wrongs+rights)。 signal2=abs(signal1)。 %LPF,包絡(luò)檢波 168。 %延遲時(shí)間 bitstream=[]。 i=IN+LL/2。 i=i+LL。 plot(1:length(signal1),signal1)。接收端波形 (包含噪聲 )39。grid on。 plot(1:length(signal2),signal2)。整流之后的波形 39。grid on。 plot(1:length(signal3),signal3)。LPF濾波后的包絡(luò)波形 39。grid on。 for i=1:length(bitstream) if bitstream(i)==0 bit1=zeros(1,N)。 end bit=[bit,bit1]。接收端二進(jìn)制波形 39。 axis([0,N*length(bitstream),])。 fpeask=[]。 endn=18。 %} n=16。 %fc=bitRate fc/bitRate為每個(gè)包含 sin周期個(gè)數(shù) bitRate=1000000。 %noise=ti。 signal=source(n,N)。 signal1=gussian(transmittedSignal,noise)。 %{ configueSignal。 end figure(8)。 grid on。ASK— 誤碼率 39。 xlabel(39。)。PeASK39。 save fpeask %} FSK的調(diào)制與解調(diào) source 函數(shù) function sendSignal=source(n,N) sendSignal=randint(1,n) bit=[]。 else bit1=ones(1,N)。 end figure(1) plot(bit),title(39。),grid on。 end FskModu 函數(shù)源程序 function transmittedSignal=fskModu(signal,bitRate,f1,f2,N) t=linspace(0,1/bitRate,N)。 c2=sin(2*pi*t*f2)。 for i=1:length(signal) if signal(i)==1 transmittedSignal=[transmittedSignal,c1]。 end end figure(2) plot(1:length(transmittedSignal),transmittedSignal)。FSK調(diào)制波形 39。grid on。 F=fft(transmittedSignal)。FSK仿真頻譜分析 39。 grid on。 signal=awgn(signal,noise)。包含噪聲的波形 39。 grid on。 wrongs=0。 else wrongs=wrongs+1。 end demoFSK 函數(shù)源程序 function bitstream=demoFSK(receivedSignal,bitRate,f1,f2,N) load FSKnum signal1=receivedSignal。 %通過 HPF，得到高頻分量 signal3=abs(signal2)。 %通過 LPF，形成包絡(luò) bitstream=[]。 %延遲時(shí)間 bitstream1=[]。 %每個(gè) bit的抽樣點(diǎn)數(shù) i=IN1 +LL/2。 i=i+LL。 plot(1:length(signal1),signal1)。接收端波形 (包含噪聲 )39。grid on。 plot(1:length(signal2),signal2)。通過 HPF得到的高頻分量波形 39。grid on。 plot(1:length(signal3),signal3)。通過 LPF后的包絡(luò)波形 39。grid on。 %通過 BPF，得到低頻分量 signal5=abs(signal4)。 %通過 LPF，形成包絡(luò) IN2=fix(length(lowpass)/2)+fix(length(daitong)/2)。 LL=N。 while (i=length(signal5)) %判決 bitstream2=[bitstream2,signal5(i)=]。 end bitstream2 figure(6) subplot(3,1,1)。title(39。)。 subplot(3,1,2)。title(39。)。 subplot(3,1,3)。title(39。)。 for i=1:min(length(bitstream1),length(bitstream2)) %判決 if(bitstream1(i)bitstream2(i)) bitstream(i)=1。 end end bitstream bit=[]。 else bit1=ones(1,N)。 end figure(7) plot(bit),title(39。),grid on。 end FSK_main 函數(shù)源程序 close all clear all %{ ti=0。 startn=6。 for ti=startn:endn n=1000。 f1=18000000。 bitRate=1000000。 %noise=ti。 signal=FSKsource(n,N)。 signal1=FSKgussian(transmittedSignal,noise)。 %{ configueSignal。 end figure(8)。 grid on。Bit Error Rate Of FSK39。 xlabel(39。)。PeFSK39。 load PeRate save fpefsk fpeask %} PSK的調(diào)制與解調(diào) source 函數(shù) function sendSignal=source(n,N) sendSignal=randint(1,n) bit=[]。 else bit1=ones(1,N)。 end figure(1) plot(bit),title(39。),grid on。 end PSKModu 函數(shù)源程序 function transmittedSignal=bpskModu(signal,bitRate,fc,N) t=linspace(0,1/bitRate,N)。 c2=sin(2*pi*t*fc + pi)。 for i=1:length(signal)