Sir i have to find the lpc coefficient for an audio. I tried a code now i am getting 20x1 coefficients. But i need 209x20 lpc coefficients. Because i have to find the coefficients for each frame.
3 views (last 30 days)
Show older comments
How i do that ? I tried one code and it is given below,
startIndex=round(56000);
[y, fs] = audioread('5pulsating.wav');
time=(1:length(y))/fs;
frameSize=floor(44*fs/1000);
endIndex=startIndex+frameSize-1;
frame = y(startIndex:endIndex);
frameSize=length(frame);
v=hamming(length(frame));
frame2=frame.*v;
spectral=20*log10(abs(fft(frame2)));
len=length(spectral);
ncoeff=20;
a=lpc(frame2,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
figure;
plot(a);
0 Comments
Answers (0)
See Also
Categories
Find more on Signal Modeling in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!