根据以下描述生成matlab代码：基于瑞利平坦衰落信道建立频率选择性衰落信道模型，输入功率时延谱和多普勒频率扩展，生成频率选择性衰落信道，并输入OFDM信号

以下是基于瑞利平坦衰落信道建立频率选择性衰落信道模型并输入OFDM信号的 MATLAB 代码：

% 输入功率时延谱和多普勒频率扩展
power_delay_spectrum = [1, 0.8, 0.6, 0.4, 0.2];
doppler_frequency_spread = 10;

% 生成频率选择性衰落信道
t = 0:0.01:1; % 时间序列
fd = doppler_frequency_spread * randn(); % 多普勒频移
h = zeros(size(t)); % 初始化信道响应
for i = 1:length(power_delay_spectrum)
    tau = i - 1; % 时延
    phi = 2 * pi * fd * tau; % 相位
    h = h + sqrt(power_delay_spectrum(i)) * exp(1j * phi) * sinc(fd*(t - tau));
end

% 生成OFDM信号
N = 64; % 子载波数
M = 16; % 调制阶数
data = randi([0 M-1], N, 1); % 随机生成调制数据
x = sqrt(N)*ifft(data); % IFFT变换
x = x(:); % 转换为列向量

% 传输信号
y = conv(h, x); % 信道传输
y = awgn(y, 20, 'measured'); % 加入高斯白噪声

% 接收信号
r = conv(h, y); % 信道卷积
r = r(1:length(x)); % 去除卷积结果中多余的部分
r = r(:); % 转换为列向量

% OFDM解调
data_hat = (1/sqrt(N))*fft(r); % FFT变换
data_hat = data_hat(:); % 转换为列向量

% 绘制频率选择性衰落信道的幅度响应
f = linspace(-0.5, 0.5, length(t)); % 频率序列
H = fftshift(fft(h)); % 频率域信道响应
figure;
plot(f, abs(H));
xlabel('归一化频率');
ylabel('幅度');
title('频率选择性衰落信道幅度响应');