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.
2 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
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 commentaires
Réponses (0)
Voir également
Catégories
En savoir plus sur Signal Modeling dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!