根据以下描述生成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(:); % 转换为列向量

% 绘制输入信号和输出信号的波形和功率谱对比图
figure;
subplot(2,2,1);
plot(real(x));
xlabel('样本数');
ylabel('幅度');
title('输入信号波形');
subplot(2,2,2);
plot(real(r));
xlabel('样本数');
ylabel('幅度');
title('输出信号波形');
subplot(2,2,3);
periodogram(x,[],[],N);
title('输入信号功率谱密度');
subplot(2,2,4);
periodogram(r,[],[],N);
title('输出信号功率谱密度');