现在我有一段编码后的2进制0和1序列，现在要对编码后的序列进行ppm调制，然后通过符合对数正态分布的衰落信道和高斯白噪声，然后再对处理后的接收信号进行ppm解调得到开始的编码后的2进制序列，其中调制后传输需要脉冲成形，注意接收时候为了过滤噪声要进行采用和阈值判断，根据以上描述给我一段matlab代码

下面是一个简单的MATLAB代码，演示了如何对编码后的二进制序列进行PPM调制，通过符合对数正态分布的衰落信道和高斯白噪声进行传输，然后进行PPM解调，得到原始编码后的二进制序列：

% 编码后的二进制序列
bits = [1 0 1 1 0 1 0 0 0 1 0 1 0 1 1 0 1 0 1 0];

% PPM调制
fs = 100000;         % 采样频率
T = 1/fs;            % 采样周期
fc = 10000;          % 载波频率
Tc = 1/fc;           % 载波周期
n = 4;               % 每个样本值对应的二进制位数
Ts = n*T;            % PPM信号周期
A = 1;               % 调制信号幅度
B = 0.5;             % PPM信号幅度
t = 0:T:Ts-T;        % 时间序列
f = zeros(1,length(bits)); % 待调制的调制信号
for i = 1:length(bits)
    if bits(i) == 0
        f(i) = -A;
    else
        f(i) = A;
    end
end
L = length(f);       % 采样点数
bits = 2^n-1;        % 每个样本值对应的二进制位数
t_sample = 0:T:n*T-T;% 样本时刻序列
f_sample = f(1:n:end);% 在样本时刻下采样到的信号
sample_values = round((f_sample+B)/B*bits/2); % 对样本值进行量化
sample_values_bin = dec2bin(sample_values,n); % 将样本值转换为二进制
sample_values_bin = reshape(sample_values_bin',1,[]); % 将二进制串展开
ppm_signal_bin = reshape(sample_values_bin,bits,[]); % 将二进制串分割成每个样本值对应的二进制串
ppm_signal_bin = ppm_signal_bin'; % 将二进制串转置
ppm_signal_bin = bin2dec(ppm_signal_bin); % 将二进制串转换为整数
ppm_signal = B*(1+ppm_signal_bin/bits).*sin(2*pi*fc*t(1:bits:end)); % PPM信号
ppm_signal = ppm_signal(:)'; % 转换为行向量
figure;
plot(t,f,'b',t(1:bits:end),ppm_signal,'r');
xlabel('Time (s)');
ylabel('Amplitude');
title('PPM Modulation');
legend('Modulating Signal','PPM Signal');

% 符合对数正态分布的衰落信道和高斯白噪声
SNR = 10;            % 信噪比
Lp = 10;             % 衰落信道的平均损耗
sigma = sqrt(A^2/2*10^(-Lp/10)/SNR); % 噪声标准差
ppm_signal_fading = ppm_signal.*lognrnd(0,1,[1,length(ppm_signal)]); % 加入对数正态分布衰落
ppm_signal_noisy = ppm_signal_fading + sigma*randn(1,length(ppm_signal_fading)); % 加入高斯白噪声
figure;
plot(t(1:bits:end),ppm_signal_noisy(1:bits:end),'b',t(1:bits:end),ppm_signal_fading(1:bits:end),'r');
xlabel('Time (s)');
ylabel('Amplitude');
title('PPM Signal with Fading and Noise');
legend('Noisy PPM Signal','Fading PPM Signal');

% PPM信号解调
ppm_demod = zeros(1,length(ppm_signal_noisy)/bits);
for i = 1:length(ppm_demod)
    ppm_bin = round((ppm_signal_noisy((i-1)*bits+1:i*bits)+B)/B*bits/2); % 将每个样本值对应的脉冲位置解调出来
    ppm_demod(i) = sum(ppm_bin.*(2.^(bits-1:-1:0))); % 将二进制数转换成十进制数，得到对应的样本值
end
demod_bits = zeros(1,length(ppm_demod)*n); % 解调出的二进制序列
for i = 1:length(ppm_demod)
    demod_bits((i-1)*n+1:i*n) = dec2bin(ppm_demod(i),n)-'0'; % 将每个样本值转换成对应的二进制序列
end
figure;
plot(1:length(bits),bits,'b',1:length(demod_bits),demod_bits,'r');
xlabel('Bit Index');
ylabel('Value');
title('PPM Demodulation');
legend('Original Bits','Demodulated Bits');

需要注意的是，代码中的信噪比SNR和衰落信道的平均损耗Lp需要根据具体情况进行设置。