資源簡介
關(guān)于matlab bpsk qpsk 16qam 源代碼 在高斯白噪聲信道中加解調(diào)
代碼片段和文件信息
function?y=ask_2(snr_in);
snr=10^(snr_in/10);
x=rand(1001);%產(chǎn)生隨機(jī)序列
h=length(x);%計算序列長度
fc=10000;fs=100000;%設(shè)定采樣頻率以及高頻和低頻(調(diào)制頻率)?
ts=10/fc;%設(shè)定碼元寬度;每碼元100個點(diǎn)?
t=0:1/fs:ts*h-1/fs;%產(chǎn)生全部碼元總的時間序列起始時刻為0步長為采樣頻率的倒數(shù)即采樣時間結(jié)束時刻為序列結(jié)束時刻
l=length(t)/h;%計算出這個序列的長度
t=reshape(t‘lh);%將序列處理成矩陣形式,每行元素為一個碼元內(nèi)的所有采樣點(diǎn)l,行數(shù)為碼元數(shù)n??:reshape(lh)
t=t‘;%因矩陣是按列存儲,故須作轉(zhuǎn)置處理
y=zeros(hl);
for?i=1:h
????if?x(i)>=0.5
?????????y(i:)=3.5*sin(2*pi*fc.*t(i:));
????else?y(i:)=sin(2*pi*fc.*t(i:));%賦值
????end
end
y0=reshape(y‘l*h1);
power=y0‘*y0;
powerzao=power/snr;%以下是產(chǎn)生高斯白噪聲序列
a=sqrt(powerzao/(l*h));
zao=a*randn(l*h1);
powerzao2=zao‘*zao;
y=y0+zao;%加噪以后的信號
t=reshape(t‘l*h1);
subplot(311);plot(ty);
grid?on
title(‘2ASK信號的時域圖形‘)
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????2154??2008-11-03?16:57??QAM_16.m
?????文件????????683??2008-11-13?14:00??psk_2.m
?????文件????????607??2008-12-22?09:56??main.m
?????文件????????937??2008-12-20?16:03??fsk_4.m
?????文件????????957??2008-12-20?16:03??fsk_2.m
?????文件????????862??2008-11-05?16:13??ask_4.m
?????文件????????904??2008-12-20?16:02??ask_2.m
?????文件????????865??2008-12-04?10:05??TZSB.m
?????文件????????779??2008-12-20?15:54??TX.m
?????文件????????867??2008-12-03?15:53??qpsk.m
-----------?---------??----------?-----??----
?????????????????9615????????????????????10
評論
共有 條評論