基于matlab的擴(kuò)頻通信仿真(編輯修改稿)
2025-07-24 18:11 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 de1) %變量含義說(shuō)明 % idata 輸入序列實(shí)部 % qdata 輸入序列虛部 % iout 輸出序列實(shí)部 % qout 輸出序列虛部% code1 擴(kuò)頻碼序列switch nargin case { 0 , 1 } error(39。lack of input argument39。)。 case 2 code1 = qdata。 qdata = idata。 end [hn,vn] = size(idata)。 [hc,vc] = size(code1)。 if hn hc error(39。lack of spread code sequences39。)。 end iout = zeros(hn,vn*vc)。 qout = zeros(hn,vn*vc)。 for ii=1:hn iout(ii,:) = reshape(rot90(code1(ii,:),3)*idata(ii,:),1,vn*vc)。 qout(ii,:) = reshape(rot90(code1(ii,:),3)*qdata(ii,:),1,vn*vc)。 end 信道仿真部分實(shí)現(xiàn)高斯白噪聲代碼函數(shù)function [iout, qout] = b2(idata, qdata, attn) %各變量含義說(shuō)明% idata 輸入序列實(shí)部% qdata 輸入序列虛部 % iout 輸出序列實(shí)部 % qout 輸出序列虛部% attn 根據(jù)信噪比得到的信號(hào)衰減水平 v = length(idata)。 h = length(attn)。 iout = zeros(h,v)。 qout = zeros(h,v)。 for ii=1:h iout(ii,:) = idata + randn(1,v) * attn(ii)。 qout(ii,:) = qdata + randn(1,v) * attn(ii)。 end 接收機(jī)部分 接收機(jī)的解擴(kuò)部分函數(shù)function [iout, qout] = despread(idata, qdata, code1) %各變量含義說(shuō)明% idata 輸入序列實(shí)部 % qdata 輸入序列虛部 % iout 輸出序列實(shí)部 % qout 輸出序列虛部 % code1 擴(kuò)頻碼序列 switch nargin case { 0 , 1 } error(39。lack of input argument39。)。 case 2 code1 = qdata。 qdata = idata。 end [hn,vn] = size(idata)。 [hc,vc] = size(code1)。 vn = fix(vn/vc)。 iout = zeros(hc,vn)。 qout = zeros(hc,vn)。 for ii=1:hc iout(ii,:) = rot90(flipud(rot90(reshape(idata(ii,:),vc,vn)))*rot90(code1(ii,:),3))。 qout(ii,:) = rot90(flipud(rot90(reshape(qdata(ii,:),vc,vn)))*rot90(code1(ii,:),3))。 end 接收機(jī)解調(diào)部分的函數(shù)function [demodata]=qpskdemod(idata,qdata,para,nd,ml) %各變量含義說(shuō)明% idata 輸入數(shù)據(jù)的實(shí)部 % qdata 輸入數(shù)據(jù)的虛部% demodata 解調(diào)后的數(shù)據(jù)% para 并行的信道數(shù)% nd 輸入數(shù)據(jù)個(gè)數(shù)% ml 調(diào)制階數(shù)demodata=zeros(para,ml*nd)。 demodata((1:para),(1:ml:ml*nd1))=idata((1:para),(1:nd))=0。 demodata((1:para),(2:ml:ml*nd))=qdata((1:para),(1:nd))=0。 誤碼計(jì)算和繪圖部分%誤碼率分析 noe2 = sum(sum(abs(datademodata)))。 nod2 = user * n
規(guī)章制度相關(guān)推薦
鄂ICP備17016276號(hào)-1