資源簡介
Matlab關(guān)于BPSK仿真-BPSK仿真.rar
該程序是老師給的,內(nèi)容很詳細(xì),每行都有解釋,很適合新手看
所含文件:
Figure47.jpg
該程序是老師給的,內(nèi)容很詳細(xì),每行都有解釋,很適合新手看
所含文件:
Figure47.jpg
代碼片段和文件信息
function?y=ask_2(snr_in);
snr=10^(snr_in/10);
x=rand(1001);%產(chǎn)生隨機(jī)序列
h=length(x);%計(jì)算序列長度
fc=10000;fs=100000;%設(shè)定采樣頻率以及高頻和低頻(調(diào)制頻率)?
ts=10/fc;%設(shè)定碼元寬度;每碼元100個(gè)點(diǎn)?
t=0:1/fs:ts*h-1/fs;%產(chǎn)生全部碼元總的時(shí)間序列起始時(shí)刻為0步長為采樣頻率的倒數(shù)即采樣時(shí)間結(jié)束時(shí)刻為序列結(jié)束時(shí)刻
l=length(t)/h;%計(jì)算出這個(gè)序列的長度
t=reshape(t‘lh);%將序列處理成矩陣形式,每行元素為一個(gè)碼元內(nèi)的所有采樣點(diǎn)l,行數(shù)為碼元數(shù)n??:reshape(lh)
t=t‘;%因矩陣是按列存儲(chǔ),故須作轉(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;%加噪以后的信號(hào)
t=reshape(t‘l*h1);
subplot(311);plot(ty);
grid?on
title(‘2ASK信號(hào)的時(shí)域圖形‘)
?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件???????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
評(píng)論
共有 條評(píng)論