0

サブバンドコーディングについてコーディングし、信号を遅らせます My code:

N = 18; %Setting the filter length
[h0,h1,g0,g1] = firpr2chfb(N-1,0.4); % MATLAB function for the analysis/synthesis       
num=20000;
[x,fs,nbits] = wavread('sub1.wav',num);
% Analysis part
% Level 1
x0 = filter(h0,1,x); % Lowpass filtering
x1 = filter(h1,1,x); % Highpass filtering 
v0 = downsample(x0,2); % Down-sampling, signal component v_0[n]
v1= downsample(x1,2); % Down-sampling, signal component v_1[n]
% Level 2
x2= v0; % Selecting the lowpass output from Level for the input to Level 
x02 = filter(h0,1,x2); % Lowpass filtering
x12 = filter(h1,1,x2); % Lowpass filtering
v02= downsample(x02,2); % Down-sampling
v12= downsample(x12,2); % Down-sampling
v2= v12; % Signal component v_2[n]
% Level 3
x3= v02; % Selecting the lowpass output from Level for the input to Level 
x03 = filter(h0,1,x3); % Lowpass filtering
x13 = filter(h1,1,x3); % Highpass filtering
v03= downsample(x03,2); % Down-sampling
v13= downsample(x13,2); % Down-sampling
v3= v13; % Signal component v_3[n]
w13= [zeros(size(1:N-1)) v3(1:length(v3)-(N-1))]

しかし、私はエラーを受け取ります

==> octavesubband のエラー w13= [zeros(size(1:N-1))v3(1:length(v3)-(N-1))]; %遅延の挿入 z^(-(N-))

私はそれを修正することを知りません。私を助けてください。ありがとう

4

1 に答える 1

1

連結しようとしている配列のサイズを確認することから始めます。

>> size(v3(1:length(v3)-(N-1)))

ans =

        2483           1

>> size(zeros(size(1:N-1)))

ans =

     1    17

それらの1つを転置する必要があることを示唆しています:

w13= [zeros(size(1:N-1))'; v3(1:length(v3)-(N-1))];

また

w13= [zeros(size(1:N-1)) v3(1:length(v3)-(N-1))'];
于 2013-05-10T04:21:47.737 に答える