MATLAB绘制AMI编码的功率谱

好的，下面是 MATLAB 实现 AMI 编码的功率谱绘制：

首先，生成随机二进制序列，代码如下：

seq = (rand(1, 1000) > 0.7);

然后，实现 AMI 编码，代码如下：

last = 0;
ami_seq = zeros(1, length(seq));
for i = 1:length(seq)
    if seq(i) == 1
        ami_seq(i) = -last;
        last = -last;
    end
end

接下来，根据 AMI 编码生成波形，代码如下：

symbol_width = 1;  % 脉冲宽度
symbol_rate = 1;  % 符号率
sampling_rate = 8 * symbol_rate;  % 采样率

signal = zeros(1, length(ami_seq) * symbol_width * sampling_rate);
last = 0;
for i = 1:length(ami_seq)
    if ami_seq(i) == 0
        signal((i-1)*symbol_width*sampling_rate+1:i*symbol_width*sampling_rate) = 0;
    else
        signal((i-1)*symbol_width*sampling_rate+1:i*symbol_width*sampling_rate/2) = last;
        signal(i*symbol_width*sampling_rate/2+1:i*symbol_width*sampling_rate) = -last;
        last = -last;
    end
end

最后，通过对波形做 FFT 并计算功率谱，绘制出功率谱图，代码如下：

[Pxx, F] = periodogram(signal, [], [], sampling_rate);
plot(F, Pxx);
title('Power Spectrum');
xlabel('Frequency (Hz)');
ylabel('Power');

完整代码如下：