请写出OFDM调制解调代码，要求包含成型滤波，载波调制

OFDM调制解调是一项比较复杂的工作，代码实现需要考虑多个因素，包括成型滤波、载波调制、IFFT/FFT变换等。以下是一个简单的OFDM调制解调代码示例，供参考：

%% OFDM调制解调代码

%% 参数设置
N = 64; % 子载波数量
M = 16; % 星座点数
K = log2(M); % 每个子载波的比特数
P = 4; % 周期扩展数
SNR = 10; % 信噪比

%% 信源编码
bits = randi([0 1], N*K, 1); % 生成随机比特序列
tx_data = reshape(bits, K, N).'; % 将比特序列重塑为矩阵

%% 加载调制
qam_mod = comm.RectangularQAMModulator('ModulationOrder', M, 'BitInput', true, 'NormalizationMethod', 'Average Power');
qam_demod = comm.RectangularQAMDemodulator('ModulationOrder', M, 'BitOutput', true, 'NormalizationMethod', 'Average Power');

%% OFDM调制
tx_signal = ofdm_mod(tx_data, P, qam_mod);

%% 信道模拟
rx_signal = awgn(tx_signal, SNR, 'measured');

%% OFDM解调
rx_data = ofdm_demod(rx_signal, P, qam_demod);

%% 信道编码
rx_bits = reshape(rx_data.', [], 1); % 将解调后的数据重塑为比特序列

%% BER计算
ber = sum(xor(bits, rx_bits))/length(bits);

%% 函数实现
function tx_signal = ofdm_mod(tx_data, P, qam_mod)
    N = size(tx_data, 2);
    % IFFT变换
    tx_ifft = ifft(tx_data, [], 2);
    % 周期扩展
    tx_cp = [tx_ifft(:, end-P+1:end) tx_ifft];
    % 并列转串行
    tx_serial = tx_cp(:);
    % 成型滤波
    h = ones(1, P+1);
    tx_filter = filter(h, 1, tx_serial);
    % 载波调制
    tx_signal = qam_mod(tx_filter);
end

function rx_data = ofdm_demod(rx_signal, P, qam_demod)
    N = length(rx_signal)/(N+P);
    % 载波解调
    rx_qam = qam_demod(rx_signal);
    % 成型滤波
    h = ones(1, P+1);
    rx_filter = filter(h, 1, rx_qam);
    % 串行转并列
    rx_parallel = reshape(rx_filter, N+P, []);
    % 去掉循环前缀
    rx_nocp = rx_parallel(P+1:end, :);
    % FFT变换
    rx_fft = fft(rx_nocp, [], 2);
    % 提取数据
    rx_data = reshape(rx_fft, [], 1);
end

以上代码中，`ofdm_mod`函数实现了OFDM调制的过程，包括IFFT变换、周期扩展、成型滤波和载波调制；`ofdm_demod`函数实现了OFDM解调的过程，包括载波解调、成型滤波、去掉循环前缀、FFT变换和提取数据。在代码中，我们使用了MATLAB中的`comm`库来实现载波调制和解调的过程，更加方便快捷。