MATLAB音频波形调制时出现矩阵维度错误求助

嘿，我太懂这种调试到头疼的感觉了——矩阵维度不匹配简直是MATLAB处理信号时的“常客”。结合你描述的情况，我帮你梳理几个最可能的原因和对应的解决办法：

最常见的元凶：立体声音频的维度问题 #

MATLAB的audioread默认会返回立体声数据（N×2的矩阵，N是采样点数量，2代表左右声道），但你的调制逻辑大概率是针对单声道（N×1的列向量）设计的，直接用就会触发维度错误。

排查&解决： #

先在函数开头加一行代码确认音频维度：

[ss_raw, Fs] = audioread('C:\Users\furka\Documents\MATLAB\sample.wav');
disp(size(ss_raw)); % 看看输出是[N 1]还是[N 2]

如果输出是[N 2]，把它转成单声道就行，两种方法选一个：

% 方法1：取左右声道的平均值（最常用）
ss = mean(ss_raw, 2);
% 方法2：直接取左声道或者右声道
ss = ss_raw(:, 1); % 左声道

容易忽略的细节：时间轴t的长度必须完全匹配 #

你说已经调整了时间轴，但要确保t的长度和音频采样点数量完全一致。比如：

N = length(ss); % 现在ss是N×1的单声道信号
t = linspace(0, (N-1)/Fs, N)'; % 生成和ss长度相同的列向量时间轴

这里要注意用(N-1)/Fs而不是N/Fs，因为采样点是从0到第N-1个，这样时间轴的每个点都能对应到一个音频采样点，避免出现长度差1的情况。

载波信号的维度也要对齐 #

假设你生成载波的代码是carrier = cos(2*pi*fc*t);，如果t是行向量，生成的carrier也会是行向量，和列向量的ss相乘就会报错。解决方法很简单：

% 要么把t转成列向量
t = linspace(0, (N-1)/Fs, N)';
% 要么把载波转成列向量
carrier = cos(2*pi*fc*t)';

给你一个可运行的调试示例 #

我把这些要点整合到一个AM调制的示例里，你可以参考着改你的代码：

function [sm] = modulation(ss, fc, mtype)
    % 加载音频并处理维度
    [ss_raw, Fs] = audioread('C:\Users\furka\Documents\MATLAB\sample.wav');
    % 转单声道
    ss = mean(ss_raw, 2);
    N = length(ss);
    % 生成匹配长度的时间轴（列向量）
    t = linspace(0, (N-1)/Fs, N)';
    
    % 根据调制类型生成已调信号
    switch mtype
        case 'AM'
            % 生成载波（自动是列向量）
            carrier = cos(2*pi*fc*t);
            % 归一化原始信号避免过调制
            ss_norm = ss / max(abs(ss));
            % 常规AM调制
            sm = (1 + ss_norm) .* carrier;
        case 'FM'
            % FM调制示例，你可以根据需求调整参数
            kf = 50; % 调频系数
            integral_ss = cumsum(ss) / Fs;
            sm = cos(2*pi*fc*t + 2*pi*kf*integral_ss);
        otherwise
            error('抱歉，暂不支持该调制类型，请输入''AM''或''FM''');
    end
    
    % 绘制波形对比
    figure('Color','w');
    subplot(2,1,1); plot(t, ss); 
    title('原始音频信号'); xlabel('时间(s)'); ylabel('幅度');
    subplot(2,1,2); plot(t, sm); 
    title('已调信号'); xlabel('时间(s)'); ylabel('幅度');
end

调用的时候记得传入合理的载波频率，比如你的音频采样率是44100Hz，载波频率可以设为10000Hz（远大于音频的最高频率，一般音频最高到20kHz，但载波要选合适的数值避免混叠）。

内容的提问来源于stack exchange，提问作者Furkan Ülger