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

以下是基于瑞利平坦衰落信道建立频率选择性衰落信道模型的 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

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