【數(shù)字基帶信號(hào)】基于matlab實(shí)現(xiàn)單雙極歸零碼和功率譜仿真
【數(shù)字基帶信號(hào)】基于matlab實(shí)現(xiàn)單雙極歸零碼和功率譜仿真
TT_Matlab
博主簡(jiǎn)介：擅長(zhǎng)智能優(yōu)化算法、神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)、信號(hào)處理、元胞自動(dòng)機(jī)、圖像處理、路徑規(guī)劃、無(wú)人機(jī)等多種領(lǐng)域的Matlab仿真，完整matlab代碼或者程序定制加qq1575304183。
1 簡(jiǎn)介
數(shù)字基帶信號(hào)是通信系統(tǒng)中最基本的信號(hào)之一,為了更好地學(xué)習(xí)和掌握該類信號(hào),利用MATLAB軟件從單雙極性歸零編碼和不歸零編碼兩個(gè)角度對(duì)其展開(kāi)研究，并且求各自的功率譜密度。
2 部分代碼
close all
clear all
%采樣點(diǎn)數(shù)的設(shè)置
k=
14
;
%每碼元采樣數(shù)的設(shè)置
L=
32
;
N=
2
^k;
M=N/L;%M為碼元個(gè)數(shù)
dt=
1
/L;%時(shí)域采樣間隔
T=N*dt;%時(shí)域截?cái)鄥^(qū)間
df=
1.0
/T;%頻域采樣間隔
Bs=N*df/
2
;%頻域截?cái)鄥^(qū)間
t=linspace(-T/
2
,T/
2
,N);%產(chǎn)生時(shí)域采樣點(diǎn)
f=linspace(-Bs,Bs,N);%產(chǎn)生頻域采樣點(diǎn)
EP1=zeros(size(f));
EP2=zeros(size(f));
EP3=zeros(size(f));
%程序第
1
部分：隨機(jī)產(chǎn)生
1000
列
0
、
1
信號(hào)序列，分別對(duì)其進(jìn)行雙極性歸零編碼和不歸零%編碼，并且求各自的功率譜密度，求功率譜密度的均值
for
x=
1
:
1000
%取
1000
次樣值
a=round(rand(
1
,M));%產(chǎn)生一個(gè)長(zhǎng)度為M的隨機(jī)序列a，
0
和
1
等概出現(xiàn)
nrz=zeros(L,M);%產(chǎn)生一個(gè)L行M列的nrz矩陣，初始化為全
0
矩陣
rz=zeros(L,M);%產(chǎn)生一個(gè)L行M列的rz矩陣，初始化為全
0
矩陣
for
i=
1
:
M
if
a
(
i
)
==
1
nrz(:,i)=
1
;%使nrz矩陣第i列全部元素都為
1
rz(
1
:L/
2
,i)=
1
;%使rz矩陣第i列前L/
2
個(gè)元素為
1
else
nrz(:,i)=
0
;%使nrz矩陣第i列全部元素都為
-1
rz(
1
:L/
2
,i)=
0
;%使rz矩陣第i列前L/
2
個(gè)元素為
-1
end
end
%分別重排nrz、rz矩陣為
1
行N列的矩陣
nrz=reshape(nrz,
1
,N);
rz=reshape(rz,
1
,N);
%做傅里葉變換并算出功率譜密度
NRZ=t2f(nrz,dt);
P1=NRZ.*conj(NRZ)/T;
RZ=t2f(rz,dt);
P2=RZ.*conj(RZ)/T;
%求功率譜密度的均值
EP1=(EP1*(x
-1
)+P1)/x;
EP2=(EP2*(x
-1
)+P2)/x;
end
%程序第
2
部分：畫(huà)波形圖和功率譜密度曲線
figure(
1
)%開(kāi)啟一個(gè)編號(hào)為
1
的繪圖窗口
subplot(
2
,
2
,
1
);plot(t,nrz)%畫(huà)單極性不歸零碼的時(shí)域圖
axis([
-5
,
5
,min(nrz)
-0.1
,max(nrz)+
0.1
])
title(
’單極性不歸零碼’
,
’fontsize’
,
12
)
xlabel(
’t(ms)’
,
’fontsize’
,
12
)
ylabel(
’nrz(t)’
,
’fontsize’
,
12
)
grid
on
subplot
(
2
,
2
,
2
)
;plot(t,rz)%畫(huà)單極性歸零碼的時(shí)域圖
axis([
-5
,
5
,min(rz)
-0.1
,max(rz)+
0.1
])
title(
’單極性歸零碼’
,
’fontsize’
,
12
)
xlabel(
’t(ms)’
,
’fontsize’
,
12
)
ylabel(
’rz(t)’
,
’fontsize’
,
12
)
grid
on
subplot
(
2
,
2
,
3
)
;plot(f,EP1)%畫(huà)單極性不歸零碼的功率譜密度圖
axis([
-5
,
5
,
0
,
1.2
])
title(
’單極性不歸零碼功率譜密度圖’
,
’fontsize’
,
12
)
xlabel(
’f(kHz)’
,
’fontsize’
,
12
)
ylabel(
’P1(f)’
,
’fontsize’
,
12
)
grid
on
subplot
(
2
,
2
,
4
)
;plot(f,EP2)%畫(huà)單極性歸零碼的功率譜密度圖
axis([
-5
,
5
,
0
,
0.3
])
title(
’單極性歸零碼功率譜密度圖’
,
’fontsize’
,
12
)
xlabel(
’f(kHz)’
,
’fontsize’
,
12
)
ylabel(
’P2(f)’
,
’fontsize’
,
12
)
grid
on
3 仿真結(jié)果
4 參考文獻(xiàn)
[1]宋建華, 朱勇, 王爾馥,等. 基于MATLAB的數(shù)字基帶信號(hào)及其功率譜研究[J]. 信息通信, 2016.
博主簡(jiǎn)介：擅長(zhǎng)智能優(yōu)化算法、神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)、信號(hào)處理、元胞自動(dòng)機(jī)、圖像處理、路徑規(guī)劃、無(wú)人機(jī)等多種領(lǐng)域的Matlab仿真，相關(guān)matlab代碼問(wèn)題可私信交流。
部分理論引用網(wǎng)絡(luò)文獻(xiàn)，若有侵權(quán)聯(lián)系博主刪除。