循環(huán)碼的性能分析課程設(shè)計(jì)任務(wù)書(shū)matlab課程設(shè)計(jì)任務(wù)書(shū)-閱讀頁(yè)
2025-01-31 16:28本頁(yè)面
【正文】 （ 2—17） 并從中得到信息組 ()mx? 。 （ 2）根據(jù)伴隨式 s(x)找出對(duì)應(yīng)的估值錯(cuò)誤圖樣 ()ex? ； （ 3）計(jì)算 )()( ^^ )( xexc xy ?? ，得到估計(jì)碼字 ()cx? 。 由于 g(x) 的次數(shù)為 n k 次， g(x) 除 E(x) 后得余式（即伴隨式）的最高次數(shù)為 nk1 次，故 S(x) 共有 2nk 個(gè)可能的表達(dá)式，每一個(gè)表達(dá)式對(duì)應(yīng)一個(gè)錯(cuò)誤格式。 BCH（ 7， 4）循環(huán)碼錯(cuò)誤圖樣表： 武漢理工大學(xué) matlab《循環(huán)碼的性能分析》課程設(shè)計(jì)任務(wù)書(shū) 5 否 錯(cuò)誤圖樣 錯(cuò)誤圖樣碼字 伴隨式 S(x) 伴隨式 E6(x)=x6 1000000 x2 100 E5(x)=x5 0100000 x2+x 110 E4(x)=x4 0010000 x2+x+1 111 E3(x)=x3 0001000 x+1 011 E2(x)=x2 0000100 x2+1 101 E1(x)=x1 0000010 x 010 E0(x)=x0 0000001 1 001 E(x)=0 0000000 0 000 上式指出了系統(tǒng)循環(huán)碼的譯碼方法：將收到的碼字 R(x) 用 g(x) 去除，如果除盡則無(wú)錯(cuò)；否則有錯(cuò)。 根據(jù)前面的討論，可得（ 7， 4）循環(huán)碼譯碼的程序框圖如圖 所示 ： 表 1 BCH（ 7， 4）循環(huán)碼錯(cuò)誤圖樣表 初始化 由 R(x)確定 S(x)： R(x)/g(x)的余式 S(x)=0，無(wú)誤碼誤碼 由 S(x)確定錯(cuò)誤圖樣 E(x) 糾錯(cuò) )()()( XRXExc ?? 存儲(chǔ) c(x) 圖 譯碼程序框圖 武漢理工大學(xué) matlab《循環(huán)碼的性能分析》課程設(shè)計(jì)任務(wù)書(shū) 6 3 程序與仿真 程序函數(shù)介紹 MATLAB 中提供了循環(huán)碼的編碼和譯碼函數(shù)，本程序直接調(diào)用進(jìn)行編程。其中 msg 代表信息碼元， 是一個(gè) K 列矩陣， N 是編碼后的碼字長(zhǎng)度； K 是信息位的長(zhǎng)度； opt 是有些編碼方式需要的參數(shù)。 說(shuō)明： 這個(gè)函數(shù)對(duì)接收到的碼字進(jìn)行譯碼，恢復(fù)出原始的信息，譯碼參數(shù)和方式必須和編碼時(shí)采用的嚴(yán)格相同。 （ 3） cyclpoly 函數(shù) 功能：生成循環(huán)碼的生成多項(xiàng)式。 p=cyclpoly（ N,K,fd_flag） 。 （ 4） randint 函數(shù) 功能：引起一致地分布的任意整數(shù)矩陣 語(yǔ)法： out = randint(m) out = randint(m,n) out = randint(m,n,rg) 武漢理工大學(xué) matlab《循環(huán)碼的性能分析》課程設(shè)計(jì)任務(wù)書(shū) 7 out = randint(m,n,rg,state) （ 5） awgn 函數(shù) 功能：在某一信號(hào)中加入高斯白噪聲 語(yǔ)法 ： y = awgn(x,SNR)。 x 的強(qiáng)度假定為 0dBW。 （ 6） quantiz 函數(shù) 功能：產(chǎn)生一個(gè)量化序號(hào)和輸出量化值。 [index ,quantiz]=quantiz(sig,partition,codebook)。 說(shuō)明：這個(gè)函數(shù)根據(jù)給定的區(qū)間參數(shù)量化信號(hào)，間隔矢量必須按遞增順序排列。 n=2^m1。 fs=100。%設(shè)置信噪比單位： dB p=cyclpoly(7,4)。%設(shè)置 8個(gè)信息組，每組 k個(gè) code=encode(msg,n,k,39。,p)。measured39。%信道輸出 %將浮點(diǎn)數(shù)轉(zhuǎn)化為 2進(jìn)制 for index=1:length(code) if(code_addnoise(index)) noisecode(index)=0。 end end %譯碼輸出 newmsg=decode(noisecode,n,k,39。)。 fft_new=fftshift(fft(newmsg,N))。 %畫(huà)圖輸出波形 figure(1)。stem(msg)。編碼器輸入信號(hào) 39。 subplot(2,2,2),stem(code),title(39。)。譯碼器輸出信號(hào) 39。 subplot(2,2,4),plot(f,abs(fft_new)/max(abs(fft_new)))。頻率 :Hz39。 ylabel(39。)。譯碼器輸出后頻譜圖 39。 figure(2) plot(code_addnoise),title(39。)。 n = 2^m1。 %信息位長(zhǎng) Fs=40。 %碼速率 N=Fs/Fd。 %進(jìn)制數(shù) for SNR = 1:100。 code = encode(msg,n,k,39。)。measured39。dB39。 %將浮點(diǎn)數(shù)轉(zhuǎn)化為二 進(jìn)制，波形整形過(guò)程 for i=1:100 for a=1:n if noisycode(i,a) noisycode(i,a) = 0。 end end end %譯碼 newmsg = decode(noisycode,n,k,39。)。 result(SNR)=ratio。The bit error rate is39。 end %不同信噪比下循環(huán)碼經(jīng)過(guò)加性高斯白噪聲信道的誤碼率 figure(1) stem(result)。循環(huán)碼在不同信噪比下的誤碼率 39。誤碼率 39。*39。信噪比 39。 ylabel(39。)。此外，對(duì)循環(huán)碼的理論知識(shí)有了進(jìn)一步的認(rèn)識(shí)。 這次 MATLAB 課程設(shè)計(jì)不僅僅加深了對(duì)課程理論知識(shí)的了解，并更加熟悉了計(jì)算機(jī)語(yǔ)言，軟件的應(yīng)用。 武漢理工大學(xué) matlab《循環(huán)碼的性能分析》課程設(shè)計(jì)任務(wù)書(shū) 14 5. 參考文獻(xiàn) [1]樊昌信， 曹麗娜 通信原理（第六版 ).北京：國(guó)防工業(yè)出版社， [2]徐明遠(yuǎn)，邵玉斌 MATLAB 仿真在通信與電子工程中的應(yīng)用 .西安：西安電子科技大學(xué)學(xué)出版社， [3]唐向宏，岳恒立，鄭雪峰 MATLAB 及在電子信息類(lèi)課程中的應(yīng)用 (第二版） .北京：電子工業(yè)出版社，
公司管理相關(guān)推薦
鄂ICP備17016276號(hào)-1