帮我编写一段matlab程序，实现音频数字水印的嵌入，并且画出嵌入前、后的频谱

以下是一个简单的MATLAB程序，用于嵌入音频数字水印并绘制其频谱。它使用了Discrete Wavelet Transformation（DWT）来嵌入水印。

% 读入音频文件
audio = audioread('original_audio.wav');

% 设定水印（长度为64）
watermark = [0 1 0 1 1 0 1 0 1 1 0 1 1 0 1 0 1 1 0 1 0 1 1 0 1 1 0 1 0 1 1 0 1 1 0 1 0 1 1 0 1 1 0 1 0 1 1 0 1 1 0 1 0 1 1 0 1 1 0 1 0 1 1 0 1];

% 设定DWT变换级数
level = 2;

% 对音频文件进行DWT变换
[coeffs,~,~] = wavedec(audio,level,'db4');

% 将水印嵌入到第一层DWT系数中
coeffs(1:length(watermark)) = embed_watermark(coeffs(1:length(watermark)), watermark);

% 重构嵌入了水印的音频文件
watermarked_audio = waverec(coeffs, 'db4');

% 绘制原始音频和嵌入了水印的音频的频谱
subplot(2,1,1);
spectrogram(audio);
title('Original Audio');
subplot(2,1,2);
spectrogram(watermarked_audio);
title('Watermarked Audio');

function [new_coeffs] = embed_watermark(coeffs, watermark)
% 将水印嵌入到DWT系数中
% coeffs: DWT系数
% watermark: 要嵌入的水印

    % 将系数进行量化
    quantized_coeffs = quantize_coeffs(coeffs);
    
    % 将水印嵌入到第一层DWT系数中
    quantized_coeffs(1:length(watermark)) = embed_watermark_bits(quantized_coeffs(1:length(watermark)), watermark);
    
    % 反量化系数
    new_coeffs = dequantize_coeffs(quantized_coeffs);
end

function [new_coeffs] = quantize_coeffs(coeffs)
% 对DWT系数进行量化

    % 设定量化步长
    step_size = 0.1;
    
    % 进行线性量化
    new_coeffs = round(coeffs/step_size);
end

function [new_coeffs] = dequantize_coeffs(quantized_coeffs)
% 对量化的DWT系数进行反量化

    % 设定量化步长
    step_size = 0.1;
    
    % 进行线性反量化
    new_coeffs = quantized_coeffs*step_size;
end

function [new_coeffs] = embed_watermark_bits(coeffs, watermark)
% 将水印嵌入到DWT系数中
% coeffs: DWT系数
% watermark: 要嵌入的水印

    % 获取要嵌入的系数
    bits_to_embed = watermark(:)';
    
    % 将系数转换为二进制
    bits = dec2bin(abs(coeffs));
    
    % 将要嵌入的系数替换到二进制编码中
    bits(1:length(bits_to_embed),end) = num2str(bits_to_embed);
    
    % 将二进制编码转换回十进制系数
    new_coeffs = bin2dec(bits);
    
    % 保留原始系数的符号
    new_coeffs = sign(coeffs).*new_coeffs;
end

请注意，此代码仅用于演示目的，可能需要进行更多的优化和改进，以实现更好的水印嵌入和提取性能。