sir i tried one code for framing and windowing but the code showing error like this Error in lpcoefficient (line 21) wind(i,:)=fr_ws(i,:).*r(i,:); how i rectify the error

Question

Suchithra K S 2019-5-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/461105-sir-i-tried-one-code-for-framing-and-windowing-but-the-code-showing-error-like-this-error-in-lpcoef

评论： Suchithra K S 2019-5-14

the code is given below,

[x,fs]=audioread('1calf.wav');
s = x(1:93098,:);
N=length(s);
frameduration=0.00475;
framelength=frameduration*fs;
num_frames=floor(N/framelength);
count=0;
for k= 1: num_frames
frames(k,:)=s((k-1)*framelength +1 : framelength*k);
    count=count +1;
    
end
m_amp = abs(max(frames,[],2)); % find maximum of each frame
id = find(m_amp > 0.06); % finding ID of frames with max amp > 0.03
fr_ws = frames(id,:); % frames without silence
for i=1:num_frames
    v=length(fr_ws);
    r=hamming(v);
    wind(i,:)=fr_ws(i,:).*r(i,:);
   
end

2 个评论
显示无隐藏无

KSSV 2019-5-9

What error? Without specifying error and without input data, how you think you will get help?

Suchithra K S 2019-5-9

the error is like this "Index exceeds matrix dimensions.

Error in lpcoefficient (line 21)

wind(i,:)=fr_ws(i,:).*r(i,:);

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Jan 2019-5-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/461105-sir-i-tried-one-code-for-framing-and-windowing-but-the-code-showing-error-like-this-error-in-lpcoef#answer_374264

编辑：Jan 2019-5-9

在 MATLAB Online 中打开

One hint: In

for i=1:num_frames
    v=length(fr_ws);
    r=hamming(v);
    wind(i,:)=fr_ws(i,:).*r(i,:);
end

v and r are constant. Then move the definition out of the loop to avoid a repeated computation. Maybe wind has the wrong size, because it is not pre-allocated - postinmg the error message would reveal this detail, so share this important message with the readers. Then:

v    = length(fr_ws);
r    = hamming(v);
wind = zeros(num_frames, v);
for i = 1:num_frames
    wind(i,:) = fr_ws(i,:) .* r(i,:);
end

Now r(i, :) is a scalar. Is this wanted, or do you mean:

for i = 1:num_frames
    wind(i,:) = fr_ws(i,:) .* r.';
end

I guess, you can omit the loop completely using the modern auto-expanding:

wind = fr_ws .* r.';

7 个评论
显示 5更早的评论隐藏 5更早的评论

Suchithra K S 2019-5-9

在 MATLAB Online 中打开

Actually I tried a matlab code for LPC coefficient for an audio. I want 20 lpc coeffiecients for 209 frames. ie the matrix 209x20. the code is given below,

[x,fs]=audioread('1calf.wav');
s = x(1:93098,:);
N=length(s);
frameduration=0.00475;
framelength=frameduration*fs;
num_frames=floor(N/framelength);
count=0;
for k= 1: num_frames
frames(k,:)=s((k-1)*framelength +1 : framelength*k);
    count=count +1;
    
end
m_amp = abs(max(frames,[],2)); % find maximum of each frame
id = find(m_amp > 0.06); % finding ID of frames with max amp > 0.03
fr_ws = frames(id,:); % frames without silence
    v=length(fr_ws);
    r=hamming(v);
   
%     for i = 1:num_frames
    wind = fr_ws .* r.';
    spectral=20*log10(abs(fft(wind))); 
len=length(spectral);
ncoeff=19;
a=lpc(wind,ncoeff);
[h,f]=freqz(1,a,512,fs);
subplot(3,1,1);
plot(time,y); title(y); axis tight
%ylim=get(gca, 'ylim');
%line([time(startIndex), time(startIndex)], ylim, 'color', 'r');
%line([time(endIndex), time(endIndex)], ylim, 'color', 'r');
axis xy
% xlabel('s');
% ylabel('
subplot(3,1,2);
xj=(1:len/2)*fs/len;
plot(xj,spectral(1:len/2));
axis xy
% xlabel('frequency/Hz');
% ylabel('magnitude/dB');
subplot(3,1,3);
plot(f,20*log10(abs(h)+eps));
legend('LPC');
% xlabel('frequency(Hz)');
% ylabel(magnitude(dB)');
r=roots(a);     
r=r(imag(r)>0.001);
ffreq=sort(atan2(imag(r),real(r))*fs/(2*pi));
ff=ffreq';
for i=1:length(ffreq)
    fprintf('Formant %dFrequency%. 1f\n',i,ffreq(i));
end
    

The error is showing like this:

Matrix dimensions must agree.

Error in lpcoefficient (line 23)

wind = fr_ws .* r.';

How to rectify the error. Can you please help me to complete the code.

Jan 2019-5-10

在 MATLAB Online 中打开

@Suchithra K S: As soon as you use the debugger and post, what

size(wind)
size(fr_ws)
size(r)

replies, I could post a working solution. Maybe you need:

% R2016a - no auto-expanding:
wind = bsxfun(@times, fr_ws, r.');

Suchithra K S 2019-5-14

I tried the above to my code at that time also it showing error

Error using bsxfun

Non-singleton dimensions of the two input arrays must match each other.

Error in lpcoefficient (line 24)

wind = bsxfun(@times, fr_ws, r.');

How to rectify this error?

请先登录，再进行评论。

sir i tried one code for framing and windowing but the code showing error like this Error in lpcoefficient (line 21) wind(i,:)=fr_ws(i,:).*r(i,:); how i rectify the error

2 个评论
显示无隐藏无

回答（1 个）

7 个评论
显示 5更早的评论隐藏 5更早的评论

另请参阅

类别

标签

Community Treasure Hunt

sir i tried one code for framing and windowing but the code showing error like this Error in lpcoefficient (line 21) wind(i,:)=​fr_ws(i,:)​.*r(i,:); how i rectify the error

2 个评论 显示 无隐藏 无

回答（1 个）

7 个评论 显示 5更早的评论隐藏 5更早的评论

另请参阅

类别

标签

Community Treasure Hunt

sir i tried one code for framing and windowing but the code showing error like this Error in lpcoefficient (line 21) wind(i,:)=fr_ws(i,:).*r(i,:); how i rectify the error

2 个评论
显示无隐藏无

7 个评论
显示 5更早的评论隐藏 5更早的评论