请帮我用matlab代码实现turbo码在AWGN、瑞利衰落信道的误码率的分析对比。请注意代码的准确性，要求能在matlab2022a上准确运行。

以下是一个基本的Turbo码在AWGN信道和瑞利衰落信道下误码率分析的Matlab代码实现，你可以根据实际需求进行修改和调整：

% 生成Turbo码
% 定义Turbo码参数
K = 1024;   % 数据块长度
rate = 1/3; % 编码速率
M = 4;      % 星座大小
decoding_iterations = 4; % 解码迭代次数

% 生成随机数据
data = randi([0 M-1], K, 1);

% 编码
tb = comm.TurboEncoder('InterleaverIndicesSource', 'Input port');
tb = comm.TurboEncoder('TrellisStructure', poly2trellis(4, [13 15], 13), 'InterleaverIndicesSource', 'Input port');
tb = comm.TurboEncoder('TrellisStructure', poly2trellis(4, [13 15], 13), 'InterleaverIndices', randperm(K), 'PuncturePatternSource', 'Property', 'PuncturePattern', [1;1;0]);
turbo_code = step(tb, data);

% 仿真参数
EbN0 = 0:2:12;          % 信噪比范围
num_errors_awgn = zeros(size(EbN0));    % AWGN信道误码个数
num_errors_rayleigh = zeros(size(EbN0)); % 瑞利衰落信道误码个数
num_bits = length(turbo_code);   % 位数

% 创建AWGN信道和瑞利衰落信道
awgn_channel = comm.AWGNChannel('NoiseMethod', 'Signal to noise ratio (SNR)');
rayleigh_channel = comm.RayleighChannel('PathDelays', [0 1], 'AveragePathGains', [0 -3], 'MaximumDopplerShift', 100, 'SampleRate', 1e6, 'PathGainsOutputPort', true);

% 模拟传输和解码
for i = 1:length(EbN0)
    % 传输并添加噪声
    channel_output_awgn = step(awgn_channel, turbo_code, 10^(EbN0(i)/10));
    channel_output_rayleigh = step(rayleigh_channel, turbo_code);
    
    % 解码
    td = comm.TurboDecoder('TrellisStructure', poly2trellis(4, [13 15], 13), 'InterleaverIndices', randperm(K), 'NumIterations', decoding_iterations);
    decoded_bits_awgn = step(td, channel_output_awgn);
    decoded_bits_rayleigh = step(td, channel_output_rayleigh);
    
    % 统计误码率
    num_errors_awgn(i) = biterr(decoded_bits_awgn, data);
    num_errors_rayleigh(i) = biterr(decoded_bits_rayleigh, data);
end

% 绘制误码率曲线
ber_awgn = num_errors_awgn / num_bits;
ber_rayleigh = num_errors_rayleigh / num_bits;
semilogy(EbN0, ber_awgn, 'r', EbN0, ber_rayleigh, 'b');
xlabel('Eb/N0 (dB)');
ylabel('BER');
legend('AWGN channel', 'Rayleigh fading channel');

这段代码首先生成一个Turbo码，然后模拟了在AWGN信道和瑞利衰落信道下的传输和解码，并统计了误码率。最后，通过绘制误码率曲线进行对比分析。你需要根据实际情况，调整代码中的参数和配置，以便正确地执行和分析。