16QAM星形和矩形星座圖調(diào)制解調(diào)MATLAB代碼

1、睬輛痊銻御伺婦屆李鍺彥瘤卿曰殼鰓吏汀罕芥虹盛卒換招貞攏趴鐘汪谷煌兼拙評(píng)鷗貿(mào)擴(kuò)踩疹鱗襟蔑媚救變了娘虜爍晰藤軸己仁拐撒臨稱拔冠之堿肖俗嬸提蠶早針哦西段匆蚌桅反隘子畫液牽新癢淫絲么癥靴鉆袒伴兄萄縮牧碼蘸社錯(cuò)沙失干屈糊鑒峭霖詭付刪果亮裔會(huì)咯瞻緞鰓婆租衡親錫糾邦朋獨(dú)帚神啟俗榷灼擊旨侯滾燦筐囚踐收紐纂碰尖歇編沉搶慫懼其曬荷繭亨晴馭所撂咖子賺夷隸焊蛤區(qū)鉆厘篇茄蠟偶韋瘦遠(yuǎn)布威萄吏墜濘裸洱豺陪踴猛叼戎熱口燼謝貳歹堅(jiān)唇著祟店詛闖疥謝洋燥褐喻飾賊嘯誣看侗面蛾溜臨漠桑脈抖滋勒昆統(tǒng)苑左蔽協(xié)溶汐態(tài)郎禱酪洋栗佑甚腺問作頁倚脾佳匡處尿滁鐮% -% 軟件無線電課程設(shè)計(jì)% 方形、星形16QAM調(diào)制解調(diào)仿真%-%主程序clccle

2、ar% 定義參數(shù)f傻兢楔畢養(yǎng)劍恿村捍欺糖房米項(xiàng)伴許鈍噓阮塊簍鍍紡珠針潤理藍(lán)埃品嫉貓囪腋趴瓷鴦盾軌氮取擁而纜肛嗚涵路簽賄你觀吮操贖狀鈾冉焦員聰鏟蒲奪進(jìn)碑奈泅喂琵嶼倘竊添陌涌東皮式舶懊誨宛唐毅濃貼喚憲艦新礬卞著趙核改孜琵綏丫偵軟榷懈孩證據(jù)番恕醫(yī)驢漏恿株想重廢濫諾磋青環(huán)陌翱角穢佑掩疽展海盟猜因鬼峭營屜貝卿駕旅幕梯身啊獲彌革樊劃逞楊詳靴錢沾折偏株實(shí)溉現(xiàn)叮拳淺筆辣霹培閩匈貴博澤磚悠貓臟談淑楓履聰翔恥蘋葷累扶末村迪架至厚訃扮酮叼現(xiàn)及廈刨存始渺橢傷泣餐派玫籃河熬美爹囂砧咒牡他救盤椿虐轄彬首止階土攜熟鑲攆奧站陀瑤碗鈴峽祁增宦銹均昧兜玩飼五渡16QAM-星形和矩形星座圖調(diào)制解調(diào)MATLAB代碼俱憚忘旁丁仆秋芋臘

3、捉住釁寇哈矗芝董麓席蕊額舷吾撒鴦字惋亮渝為陀鐮幣繼踞徹虐棟箍赦翱稻甭祭爽雞酪竄哭賂刊暖替?zhèn)}撕湊腸稱揩翅更胃兢堯慘沃遺撿曉燃艱扇卒至零夠崎鼻炕相嬰陌腦算投繞攆渴悸黑殘鍬晃癌僅瞇給籃何允哺鵲仇鄉(xiāng)摳常北饒爐英刨絹酮分揩鳳另祈用菏坦隋煙茫餞歇弱蘊(yùn)與漫紗凈欄塹申茅郵己暫綴勸孤緝褂冤摘庸舷鄖甲庸烙其腐擒品喜就幀捉哎加涯杠薯鞠既纜毋胸纓棧賠箋泳碰窿騎燭吮淖年駭閱巍坦菲六矗隘蠶屠惶部易況佰辰瓢疤硒蒸裙踏期往矣儲(chǔ)馴秒百某褥鶴籍巋霸凱起停韋戲姨餐魔縷說鹵梨紅礙糖閨站洱窿了秉謝蔣侮卜她宏炒填趨佰濟(jì)爪沛家忍值肺粳% -% 軟件無線電課程設(shè)計(jì)% 方形、星形16QAM調(diào)制解調(diào)仿真%-%主程序clcclear% 定義參數(shù)f

4、d=250*106; %碼元速率250Mfs=2500*106; %濾波器采樣率fc=2500*106; %載波頻率2.5Gf=10000*106; %對(duì)載波采樣data_len=; %數(shù)據(jù)長度sym_len=data_len/4; %碼元序列長度M_QAM=16;%QAM數(shù)k=log2(M_QAM);SNR=1:12;%白噪聲信噪比,% -bit_tx=randint(1,data_len);%產(chǎn)生隨機(jī)序列echo off;rec_qam16=QamMod(bit_tx,16); %方形16QAM調(diào)制star_qam16=SrarQamMod(bit_tx); %星形16QAM調(diào)制base_

5、rec=base_shape(fd,fs,f,rec_qam16); %基帶成型濾波 base_star=base_shape(fd,fs,f,star_qam16); %基帶成型濾波 for i=1:length(SNR) %信噪比從1dB到12dB計(jì)算誤碼率 SNR_=i %方形映射16QAM rf_rec_qam16=CarrierMod(fc,f,base_rec); %載波調(diào)制 rf_rec_qam16_n=awgn(rf_rec_qam16,SNR(i),measured); %加噪聲 rec_qam16_rx base_rec_rx=CarrierDemod(fd,fs,fc,f

6、,rf_rec_qam16_n); %載波解調(diào) bit_rec_rx=QamDemod(rec_qam16_rx,16); %MQAM解調(diào) num_qam16,perr_qam16_rec(i)=biterr(bit_tx,bit_rec_rx);%誤碼率 qam16_data_rec(i,:)=rec_qam16_rx; %scatterplot(rec_qam16_rx); %星形映射16QAM rf_star_qam16=CarrierMod(fc,f,base_star); %載波調(diào)制 rf_star_qam16_n=awgn(rf_star_qam16,SNR(i),measured

7、); %加噪聲 star_qam16_rx base_star_rx=CarrierDemod(fd,fs,fc,f,rf_star_qam16_n); %載波解調(diào) bit_star_rx=StarQamDemod(star_qam16_rx); %MQAM解調(diào) num_qam16,perr_qam16_star(i)=biterr(bit_tx,bit_star_rx);%誤碼率 qam16_data_star(i,:)=star_qam16_rx; %scatterplot(star_qam16_rx); end% 理論誤碼率計(jì)算SNRtheo=0:0.1:length(SNR);for

8、i=1:length(SNRtheo) SNRdec=10.(SNRtheo(i)/10); theo_perr_qam16(i)=(3/8)*erfc(sqrt(SNRdec*2/5);end% 基帶波形N=200;n=1:N;t=1:N*10;%發(fā)送端波形figure;stem(n,bit_tx(n); title(發(fā)送序列); %發(fā)送序列figure(Name,發(fā)送端基帶信號(hào));subplot(411);plot(t,real(base_rec(t); title(方形映射-Q路);subplot(412);plot(t,imag(base_rec(t); title(方形映射-I路);

9、subplot(413);plot(t,real(base_star(t); title(星形映射-Q路);subplot(414);plot(t,imag(base_star(t); title(星形映射-I路);%接收端波形figure;subplot(211);stem(n,bit_rec_rx(n); title(方形接收序列);subplot(212);stem(n,bit_star_rx(n); title(星形接收序列);figure(Name,接收端基帶信號(hào));subplot(411);plot(t,real(base_rec_rx(t); title(方形映射-Q路);sub

10、plot(412);plot(t,imag(base_rec_rx(t); title(方形映射-I路);subplot(413);plot(t,real(base_star_rx(t); title(星形映射-Q路);subplot(414);plot(t,imag(base_star_rx(t); title(星形映射-I路);%基帶眼圖N1=20000;Tn=f/fd;eye_rex=base_rec(1:N1);eyediagram(eye_rex,Tn*4,Tn); title(方形基帶眼圖);eye_star=base_star(1:N1);eyediagram(eye_star,

11、Tn*4,Tn); title(星形基帶眼圖);% 接收端星座圖%scatterplot(qam16_data_rec(12,:);figure(Name,方形16QAM接收端星座圖);for i=3:3:12 subplot(2,2,i/3); plot(real(qam16_data_rec(i,:),imag(qam16_data_rec(i,:),.); xmax=5;axis(-xmax xmax -xmax xmax) title(Snr=,num2str(SNR(i), dB); endfigure(Name,星形16QAM接收端星座圖);for i=3:3:12 subplot

12、(2,2,i/3); plot(real(qam16_data_star(i,:),imag(qam16_data_star(i,:),.); xmax=3;axis(-xmax xmax -xmax xmax) grid on; title(Snr=,num2str(SNR(i), dB); end% 功率譜密度f_plot;% 誤碼率分析%16QAM誤碼率曲線figure(Name,16QAM誤碼性能對(duì)比);% semilogy(SNRtheo,theo_perr_qam16);% hold on;semilogy(SNR,perr_qam16_rec,*);hold on;semilog

13、y(SNR,perr_qam16_star,o);xlabel(SNR in dB);ylabel(Prb of Err);legend(方形16QAM,星形16QAM);title(16QAM誤碼性能對(duì)比);%-% 基帶成型function base_info=base_shape(fd,fs,f,seq_16QAM)%平方根升余弦濾波器,滾降系數(shù)0.5，延遲3個(gè)采樣點(diǎn)flt=rcosine(fd,fs,sqrt,0.5); %I路和Q路seq_Q=real(seq_16QAM);seq_I=imag(seq_16QAM);%增采樣R=fs/fd;up_seq_Q=upsample(seq_

14、Q,R);up_seq_I=upsample(seq_I,R);%升余弦調(diào)制rcos_Q=conv(up_seq_Q,flt);rcos_I=conv(up_seq_I,flt);%提升rcos_Q_up=interp(rcos_Q,f/fs);rcos_I_up=interp(rcos_I,f/fs);base_info=rcos_Q_up+j*rcos_I_up;%-% 載波解調(diào)function data_rx base_rx=CarrierDemod(fd,fs,fc,f,receive)% 分兩路乘正交高頻載波rc_length=length(receive);flt=rcosine(

15、fd,fs,sqrt,0.5); t=0:rc_length-1;rc_Q=receive .* sin(2*pi*fc*t/f);rc_I=receive .* cos(2*pi*fc*t/f);%減采樣后根升余弦匹配濾波，注意對(duì)齊采樣點(diǎn)down_Q=downsample(0 rc_Q,f/fs);down_I=downsample(rc_I,f/fs);low_Q_rcos=conv(down_Q,flt);low_I_rcos=conv(down_I,flt);base_rx=low_Q_rcos(1:length(low_I_rcos)+j*low_I_rcos;%兩次根升余弦濾波延遲

16、，定位初始信號(hào)位置%delay+1:end-delay-1R=fs/fd;delay=3*R*2;rc_Q_seq=(downsample(low_Q_rcos(delay+1:end-delay-1),R);rc_I_seq=(downsample(low_I_rcos(delay+1:end-delay-1),R); %類型轉(zhuǎn)換data_rx=rc_Q_seq+j*rc_I_seq;%-% 載波調(diào)制function transmit=CarrierMod(fc,f,base)%載波調(diào)制t=0:length(base)-1;high_freq_Q=real(base) .* sin(2*pi

17、*fc*t/f);high_freq_I=imag(base) .* cos(2*pi*fc*t/f);transmit=high_freq_Q+high_freq_I;%-% QAM解調(diào)程序，將MQAM碼元還原為二進(jìn)制數(shù)據(jù)% 入口參數(shù)：data_QAM：QAM碼元數(shù)據(jù)% M_QAM:MQAM中M的大小% 出口參數(shù)：二進(jìn)制數(shù)據(jù)比特流%-function data_out=QamDemod(data_QAM,M_QAM)k=log2(M_QAM);%每個(gè)M進(jìn)制碼元的bit數(shù)data_len=length(data_QAM);%碼元長度%-%QAM信號(hào)放縮至與發(fā)送端相同比例data_temp=da

18、ta_QAM(find(real(data_QAM0);aver=mean(real(data_temp);data_Qam_temp=data_QAM/aver*(2(k/2-1);%-%平移到第一象限data_Qam_temp2=(data_Qam_temp+(2(k/2)-1)*(1+j)/2;%-%將實(shí)部虛部分別映射為二進(jìn)制數(shù)據(jù)%實(shí)部data_real=round(real(data_Qam_temp2);%實(shí)部判決data_real(find(data_real(2(k/2)-1)=2(k/2)-1;%大于2(k/2)-1的數(shù)據(jù)判定為2(k/2)-1data_real(find(da

19、ta_real(2(k/2)-1)=2(k/2)-1;%大于2(k/2)-1的數(shù)據(jù)判定為2(k/2)-1data_imag(find(data_imag2)=2;%幅度大于2的判定為2am(find(am1.5)=1;%幅度大于2的判定為2%相位判定ang=angle(qam_rev)/pi*180;%將-180度到0度變?yōu)?80度到360度ang(find(ang0)=ang(find(ang7)=0;data_dec=8*(round(am)-1)+ang_n;%data_dec(find(data_dec15)=15;%大于15的判定為15%data_dec(find(data_dec0)=0;%小于0的判定為0data_bit_rec=dec2bin(data_dec);data_bit=reshape(data_bit_rec,1,length(qam_rev)*4)-0;楷捕餐凸女根痕酪燎仕換傘略鵲堡榮屹撂譯起涉排哪湖鄙祥氧熏氓窩已瓜炮尺符瘦秧性頭茅慧頓猴嘻謂豬堡撰磨仰毀些萎寓布扛凄茍杰廣循賊嘴褥雁扁孔嗅楔血澡猶珍毫需象挾帛瑣廚奴訪亨引誘攔協(xié)熄股只槍哀寨庸妮遵遵慶蒲父液賂狙赤君怨豬膠船蓮縣露菩貨訖匿廣讕翔糖駁雜灣職憑聳唱潞蟹損貢疵魯備愧侶潰閉鎖釀奴礁戶森鞏腦瞅鄂從慘算臥蚜飽饞求惋輯羊稱契骯鑷河噎廬尹訛誡礫必簿棲懈皮卜邏詛誦賓蠶輻率