我正在尝试绘制陷波滤波器的图像,即陷波滤波之前和之后的信号的频率响应。但是,我无法显示两个倒数图,只能显示第一个(凹口)。因为该错误:'plot:输入参数#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)播放经过过滤的音频;但是实际上我只有:一个图(陷波图)
答案 0 :(得分:0)
我认为您只需更改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);
错误消失了,但是我不确定代码的正确性(无法帮助您)...
我使用了立体声采样,所以我不得不添加另一个小的更改。
这是我的代码测试示例:
{{1}}