帮我写一个基于OFDM的水下多径信道的MATLAB代码，可以分别绘制不同多经数下幅度随时延变化的图像

OFDM是一种常用的多载波调制技术，可以有效地对抗多径衰落。在水下通信中，多径衰落是不可避免的，因此OFDM技术也被广泛应用。下面是一个基于OFDM的水下多径信道的MATLAB代码，可以绘制不同多经数下幅度随时延变化的图像。

% OFDM水下多径信道模拟
clear all;
close all;

% 参数设置
N = 64;                     % 子载波数
cp_len = 16;                % 循环前缀长度
EbNo = 10;                  % 信噪比
n_carrier = 32;             % 载波数
n_path = 5;                 % 多径数
delay_spread = [0, 1, 2, 3, 4];     % 不同多径下的时延扩展
path_gains = [1, 0.8, 0.6, 0.4, 0.2];    % 不同多径下的多径增益

% 生成QPSK调制信号
tx_data = randi([0, 1], 1, N*n_carrier);
tx_data = reshape(tx_data, N, n_carrier);
tx_data = 2*tx_data - 1;
tx_data = tx_data/sqrt(2);

% IFFT变换
tx_signal = ifft(tx_data, N);

% 添加循环前缀
tx_signal_cp = [tx_signal(N-cp_len+1:N, :); tx_signal];

% 信道
for i = 1:n_path
    delay = delay_spread(i);
    path_gain = path_gains(i);
    channel = zeros(N+cp_len, n_carrier);
    for j = 1:n_carrier
        h = zeros(N+cp_len, 1);
        h(delay+1) = path_gain;
        h = ifft(h);
        channel(:, j) = [h(N-cp_len+2:N); h];
    end
    rx_signal_cp(:, :, i) = filter(channel, 1, tx_signal_cp);
end

% 加噪声
rx_signal_cp = awgn(rx_signal_cp, EbNo, 'measured');

% 去除循环前缀并进行FFT变换
for i = 1:n_path
    rx_signal(:, :, i) = rx_signal_cp(cp_len+1:end, :, i);
    rx_data(:, :, i) = fft(rx_signal(:, :, i), N);
end

% 绘制幅度随时延变化的图像
figure;
for i=1:n_path
    subplot(n_path, 1, i);
    plot(abs(rx_data(1, :, i)), 'r');
    xlabel('Carrier Index'); ylabel('Magnitude');
    title(sprintf('Delay Spread = %d', delay_spread(i)));
    grid on;
end

这个代码可以模拟OFDM水下多径信道，并绘制不同多径数下幅度随时延变化的图像。你可以根据需要修改参数，比如子载波数、循环前缀长度、信噪比、载波数、多径数等等。