matlab - 如何解决此错误：“绘图：输入参数 #2 和 #3 的尺寸错误：尺寸不兼容。”？

Question

我正在尝试绘制陷波滤波器的图像，即滤波前和陷波滤波后信号的频率响应。但是我不能显示最后两个图表，只有第一个（Notch）。因为那个错误：'情节：输入参数＃2和＃3的尺寸错误：尺寸不兼容。'

我试图改变频率的间隔，但不起作用

wavread("C:\Users\kaline\Desktop\Disciplinas\2019.1\PDS\2° etapa\Trabalho final\sirene.wav","size")
[y,Fs,bits]=wavread("C:\Users\kaline\Desktop\Disciplinas\2019.1\PDS\2° etapa\Trabalho final\sirene.wav");Fs,bits


y=loadwave("C:\Users\kaline\Desktop\Disciplinas\2019.1\PDS\2° etapa\Trabalho final\sirene.wav");
//playsnd(y)

fs=12000;
n=[0:5999];
a=[0.940809 -1.105987 0.940809];
b=[1 -1.105987 0.881618];
[xm,fr]=frmag(a,b,512);
subplot(3,1,1)
plot(fr,abs(xm))
Y=fft(y);
f=n*fs/6000;
subplot(3,1,2)
plot(f,abs(Y)),
z=filter(a,b,y);
z=filter(a,b,z);
z=filter(a,b,z);
w=fft(z);
subplot(3,1,3)
plot(f,abs(w),'r')

playsnd(w);

我希望树图在同一个窗口中，并使用 playnsd(w) 播放过滤后的音频；但实际上我所拥有的是：只有一个图（Notch 的图）

标签： matlabaudioscilab

Answer 1

我认为您只需要更改 的大小n以匹配Y.

Y=fft(y);在代替之后添加以下代码行f=n*fs/6000;：

n_samp = size(Y,2);
n = [0:n_samp-1]
f=n*fs/n_samp;

错误消失了，但我不确定代码的正确性（并且无法帮助您）...

---

我使用了立体声样本，所以我不得不添加另一个小改动。
这是我的代码测试示例：

//wavread("C:\Users\kaline\Desktop\Disciplinas\2019.1\PDS\2° etapa\Trabalho final\sirene.wav","size")
wavread("c:\Program Files\MATLAB\R2016a\toolbox\audio\samples\Turbine-16-44p1-mono-22secs.wav","size")
//[y,Fs,bits]=wavread("C:\Users\kaline\Desktop\Disciplinas\2019.1\PDS\2° etapa\Trabalho final\sirene.wav");Fs,bits
[y,Fs,bits]=wavread("c:\Program Files\MATLAB\R2016a\toolbox\audio\samples\WashingMachine-16-44p1-stereo-10secs.wav");Fs,bits

//y=loadwave("C:\Users\kaline\Desktop\Disciplinas\2019.1\PDS\2° etapa\Trabalho final\sirene.wav");
y=loadwave("c:\Program Files\MATLAB\R2016a\toolbox\audio\samples\Engine-16-44p1-stereo-20sec.wav");
//playsnd(y)

fs=12000;
n=[0:5999];
a=[0.940809 -1.105987 0.940809];
b=[1 -1.105987 0.881618];
[xm,fr]=frmag(a,b,512);
subplot(3,1,1)
plot(fr,abs(xm))
Y=fft(y);
n_samp = size(Y,2);
n = [0:n_samp-1]
f=n*fs/n_samp;
subplot(3,1,2)
plot(f,abs(Y)),
z=filter(a,b,y(1,:));
z=filter(a,b,z(1,:));
z=filter(a,b,z(1,:));
w=fft(z);
subplot(3,1,3)
plot(f,abs(w),'r')

playsnd(w);