0

我有一个 wav 文件,我将其导入 matlab,并使用以下函数从采样数据中恢复原始信号:

sub1=wavread('jamming.wav');
 magnSub1=abs(fft(sub1));
 phaseSub1=angle(fft(sub1));
 sub1_L=magnSub1.*exp(i*phaseSub1);
 sub1_L=ifft(sub1_L);
 q=(2^15-1)*(abs(sub1_L).*cos(angle(sub1_L)));
 plot(q);
 axis([44000,45000,-2^15+1,+2^15-1]);

但现在我想在现有图表上添加一些散点图来比较采样数据和连续波。这是我到目前为止所拥有的:

[wave,fs]=wavread('jamming.wav'); 

sound(wave,fs);

t=0:1/fs:(length(wave)-1)/fs;
z = t([1],:);
wave1=wave(:,[1]);
scatter(z,wave1);

问题是我无法将 x 轴缩放到第一张图。有人可以帮我解决这个问题吗?我的最终图表应该与本页中的图表相同:http: //www.cheers4all.com/2012/07/analog-to-digital-convertor-matlab-source-code/

非常感谢您的帮助,在此先感谢您!

4

1 回答 1

0

这里的问题是您的第一个图是在该范围内绘制的,1:length(q)而第二个图是在该范围内绘制的0 : 1/fs : (length(wave)-1)/fs

假设您的值qwave1长度相同(看起来可能是),请尝试:

figure; hold all;
plot(z, q);
plot(z, wave1);

然后,您可以根据自己的喜好调整每个绘图命令的线条属性。

如果q并且wave1长度不完全相同,则快速而肮脏的修复将是:

figure; hold all;
plot(1:length(q), q);
plot(1:length(wave1), wave1);

这可能不会为您提供您正在寻找的确切内容,但它应该是朝着正确方向迈出的一步。

于 2013-04-28T22:20:23.057 回答