我需要

• 读取.wav文件
• 使用 STFT 从中制作 2D 矩阵
• 生成其奇异值分解
• 做一些水印

但我无法通过矩阵乘法再次构建原始矩阵：

 x = audioread('0.wav');
 stft = spectrogram(x); 
 [U,S,V] = svd(stft);
 stftb=U*S*V.';

问题是 STFT 和 STFTb 不一样。