实现发送全1序列的双极性不归零功率谱密度的Matlab代码
时间: 2024-02-13 12:05:53 浏览: 35
发送全1序列的双极性不归零码的功率谱密度是一定的,可以通过直接计算得到。具体的MATLAB代码如下:
```
seq = ones(1,1000); % 发送全1序列
t = 0:0.01:length(seq)-0.01; % 采样时间
s = zeros(1,length(t)); % 初始化输出序列
for i = 1:length(seq)
if seq(i) == 0
s((i-1)*100+1:i*100) = -1;
else
s((i-1)*100+1:i*100) = 1;
end
end
s_diff = diff(s); % 对输出序列进行差分
for i = 1:length(s_diff)
if s_diff(i) == -2
s(i+1:end) = s(i+1:end) + 2;
elseif s_diff(i) == 2
s(i+1:end) = s(i+1:end) - 2;
end
s(i*100) = 0; % 在每个码元的中间加入归零位
end
NFFT = 1024; % FFT点数
Pxx = abs(fft(s,NFFT)).^2/length(s)/100; % 计算功率谱密度
f = linspace(-0.5,0.5,NFFT)*100; % 频率坐标
plot(f,10*log10(fftshift(Pxx)),'LineWidth',2);
grid on;
xlabel('Frequency');
ylabel('Power Spectral Density (dB/Hz)');
title('Bipolar NRZI PSD of All "1"s');
```
在这段代码中,我们首先定义了一个长度为1000的全1序列,然后根据该序列生成了双极性不归零码,并计算了其功率谱密度。为了得到更好的频谱图,我们使用了1024点FFT,并使用fftshift将频谱移动到中心位置,最后画出功率谱密度随频率变化的曲线。
相关推荐
![doc](https://img-home.csdnimg.cn/images/20210720083327.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)