【正文】 for n=1:a if rn(n)==1 。 l=l。 %每檢測到一個(gè) 1變一次符號(hào) rn(n)=l。 end end AMI=rn。 stairs(rn)。 grid。 axis([0 length(rn) min(rn) max(rn)+]) title(39。AMI 碼波形 39。)。 %計(jì)算功率譜 % Fs=1024。 %采樣頻率 n=0:1/Fs:1。 nfft=1024。 南華大學(xué)電氣工程學(xué)院課程設(shè)計(jì) 第 22 頁 共 27 頁 cxn=xcorr(rn,39。unbiased39。)。 %計(jì)算序列的自相關(guān)函數(shù) CXk=fft(cxn,nfft)。 Pxx=abs(CXk)。 index=0:round(nfft/21)。 k=index*Fs/nfft。 plot_Pxx=10*log10(Pxx(index+1))。 figure(2) plot(k,plot_Pxx)。 title(39。AMI 碼功率譜 39。)。 %HDB3 碼編碼及功率譜 clear clc a=20。 %a 表示元素個(gè)數(shù)，可以修改碼元為任意個(gè)數(shù)如 1024 N=1024。 x=unidrnd(2,1,a1)。 %隨機(jī)產(chǎn)生離散均勻二值分布（ 1， 2），共 a1 個(gè) n=length(x)。 y=zeros(1,n)。 last_y=1。 %設(shè)定初始參考值 last_v=1。 zero_flag=0。 for i=1:n if(x(i)==1) y(i)=last_y。 last_y=y(i)。 zero_flag=0。 else y(i)=0。 南華大學(xué)電氣工程學(xué)院課程設(shè)計(jì) 第 23 頁 共 27 頁 zero_flag=zero_flag+1。 end if(zero_flag==4) if(last_y==last_v) y(i3)=last_y。 y(i)=y(i3)。 last_y=y(i)。 last_v=y(i)。 else y(i)=last_y。 last_y=y(i)。 last_v=y(i)。 end zero_flag=0。 end end %%%%%%%%%%%%%%%%% 以下繪制波形 %%%%%%%%%%%%%%%%% t0=100。 t=0:1/t0:length(x)。 for i=1:length(x) for j=1:t0 h((i1)*t0+j)=y(i)。 %把等式右邊改為 (y(i)==1)可以繪制 HDB3正極性編碼+H %把等式右邊改為 (y(i)==1)可以繪制 HDB3負(fù)極性編碼 H end end h=[h,h(i*t0)]。 M=max(h)。 南華大學(xué)電氣工程學(xué)院課程設(shè)計(jì) 第 24 頁 共 27 頁 m=min(h)。 figure(1)。 plot(t,zeros(1,length(t)),39。LineWidth39。,1)。hold on plot(t,h,39。LineWidth39。,2)。 set(gca,39。xtick39。,[0:1:length(x)],39。Ytick39。,[0,1]) grid on box axis([0 i M+.1])。 title(39。HDB3 碼波形 39。)。 %計(jì)算功率譜 % Fs=1024。 %采樣頻率 n=0:1/Fs:1。 nfft=1024。 cxn=xcorr(h,39。unbiased39。)。 %計(jì)算序列的自相關(guān)函數(shù) Cxk=fft(cxn,nfft)。 Pxx=abs(Cxk)。 index=0:round(nfft/21)。 k=index*Fs/nfft。 plot_Pxx=10*log10(Pxx(index+1))。 figure(2) plot(k,plot_Pxx)。 title(39。HDB3 碼功率譜 39。)。 %CMI 碼編碼及功率譜 clear clc a=20。 %a 表示元素個(gè)數(shù)，可以修改碼元為任意個(gè)數(shù)如 1024 N=1024。 r=unidrnd(2,1,a1)。 %隨機(jī)產(chǎn)生離散均勻二值分布（ 1， 2），共 a1 個(gè) 南華大學(xué)電氣工程學(xué)院課程設(shè)計(jì) 第 25 頁 共 27 頁 rn=r1。 %每個(gè)元素減 1變?yōu)椋?0， 1）二值隨機(jī)數(shù)列 rn=[1,rn]。 %跟書上一致，從 1 開始隨機(jī)產(chǎn)生離散均勻二值分布（ 1， 2），共 a個(gè) ori=rn。 %將原序列保存起來，便于后面解碼后比較 ori x=rn。 m=1。 %1 碼的 cmi 編碼的判決碼 。 n=1。 %1 碼的 cmi 編碼判決碼的輔助碼 。 q=[0,。]。 %cmi 編碼后的前一位序列數(shù)組 。 h=[0,。]。 %cmi 編碼后的后一位序列數(shù)組 。 for k=1:20 %for 語句 NRZ 碼序列的第一個(gè)到第十個(gè)循環(huán) 。 m=mod(n,2)。 %cmi 編碼結(jié)果 00 和 11 交替的判決語句 。 if (x(k)==0) %如果 NRZ 碼為 0,則進(jìn)行下一步 。 q(k)=0。 %cmi 編碼輸出的前一位為 0。 h(k)=1。 %cmi 編碼輸出的后一位為 1。 end %本 if 語句循環(huán)結(jié)束 。 if (x(k)==1 amp。 m==1) %如果 NRZ碼為 1且編碼判決碼為 1,則進(jìn)行下一步 。 q(k)=1。 %cmi 編碼輸出的前一位為 1。 h(k)=1。 %cmi 編碼輸出的后一位為 1。 n=n+1。 %判決輔助碼加 1。 end %本 if 語句循環(huán)結(jié)束 。 if (x(k)==1 amp。 m==0) %如果 NRZ碼為 1且編碼判決碼為 0,則進(jìn)行下一步 。 q(k)=0。 %cmi 編碼輸出的前一位為 0。 h(k)=0。 %cmi 編碼輸出的后一位為 0。 n=n+1。 %判決輔助碼加 1。 end %本 if語句循環(huán)結(jié)束 。 end %本 for 語句循環(huán)結(jié)束 。 fprintf(39。Display encode result:\n39。)。%顯示 輸出編碼結(jié)果為 :。 for k=1:20 %for 語句 cmi 編碼后序列的第一個(gè)到第十個(gè)的循環(huán) 。 if (q(k)==1amp。h(k)==1)%如果 cmi 編碼后的前一位和后一位都為 1,則進(jìn)行下一南華大學(xué)電氣工程學(xué)院課程設(shè)計(jì) 第 26 頁 共 27 頁 步 。 fprintf(39。 11 39。) %cmi 編碼輸出為 11。 else if (q(k)==0 amp。h(k)==0) %如果 cmi 編碼后的前一位和后一位都為 0,則進(jìn)行下一步 。 fprintf(39。 00 39。) %cmi 編碼輸出為 00。 else %上述情況以外 。 fprintf(39。 01 39。)%cmi 編碼輸出為 01。 end %結(jié)束本 if 語句 。 end %結(jié)束上層 if 語句 。 end for i=1:20 if(i==1) kk(1)=q(1)。 end if(i20) kk(i*2)=h(i)。 kk(i*2+1)=q(i+1)。 end if(i==20) kk(i*2)=h(i)。 end end stairs(kk)。 axis([1 41 ])。 title(39。CMI 編碼輸出波形 39。)。 grid on。 Fs=1024。 %采樣頻率 n=0:1/Fs:1。 nfft=1024。 cxn=xcorr(kk,39。unbiased39。)。 %計(jì)算序列的自相關(guān)函數(shù) Cxk=fft(cxn,nfft)。 Pxx=abs(Cxk)。 南華大學(xué)電氣工程學(xué)院課程設(shè)計(jì) 第 27 頁 共 27 頁 index=0:round(nfft/21)。 k=index*Fs/nfft。 plot_Pxx=10*log10(Pxx(index+1))。 figure(2) plot(k,plot_Pxx)。 title(39。CMI 碼功率譜 39。)