spmain2.m

自己写的关于lpc分析与合成的Demo程序

% spmain.m
OrderP=10;
FrameLen=20             % 20 ms per Frame;
NPoint=FrameLen*8;       %

[Y0,Fs0,Bits0]=wavread('bj.wav');
DataInP01=Y0(1:17/20*length(Y0));

[Y,Fs,Bits]=wavread('bj2.wav');
DataInP1=Y(1:17/20*length(Y));
DataInP2=DataInP1(1:length(DataInP1)/5);

set (0, 'DefaultAxesFontSize', 8, ... 
        'DefaultAxesFontName', 'Times'); 
figure(1); 
set (gcf, 'Units', 'centimeter'); 
pos = get (gcf, 'Position'); 
pos(3) = 10; 
pos(4) = 8; 
set (gcf, 'Position', pos);   
set (gca, 'Units', 'centimeter'); 
set (gca, 'Position', [1.5 1.2 8 6]); 

subplot(211);
plot((1:length(DataInP1))/Fs,DataInP01(1:length(DataInP1)));
xlabel('Time, s');ylabel('Amplitude');title('Original Waveform');  
subplot(212);
plot((1:length(DataInP1))/Fs,DataInP1);
xlabel('Time, s');ylabel('Amplitude');title('Syn-Waveform');

set (gcf, 'PaperPositionMode', 'auto');   
print -depsc  bjwavform2.eps; 


DataInWin=DataInP2.*window(@hamming,length(DataInP2));
Dfft=fft(DataInP2);
plotlen=floor(length(Dfft)/2);


TotalPoint = length(DataInWin);
NFrame = floor(TotalPoint/NPoint);

for i=4:4:4
    
DataF = DataInWin(NPoint*i+1:NPoint*(i+1));
DataWin = DataF.*window(@hamming,NPoint);
DataFFT = fft(DataWin);

A = lpc(DataWin,OrderP);
[H1]=freqz(1,A,160);


set (0, 'DefaultAxesFontSize', 8, ... 
        'DefaultAxesFontName', 'Times'); 
figure(2); 
set (gcf, 'Units', 'centimeter'); 
pos = get (gcf, 'Position'); 
pos(3) = 10; 
pos(4) = 8; 
set (gcf, 'Position', pos);   
set (gca, 'Units', 'centimeter'); 
set (gca, 'Position', [1.5 1.2 8 6]); 

plotlen=floor(length(DataFFT)/2);
plot(4000*(1:plotlen)/plotlen,10*log10(abs(DataFFT(1:plotlen))/max(abs(DataFFT))));grid;
hold on;
plot((1:160)/160*4000,10*log10(abs(H1)/max(abs(H1))),'-.');
xlabel('Frequency, Hz');ylabel('Normlized Amplitude /dB');
set (gcf, 'PaperPositionMode', 'auto');   
print -depsc  bjSpetrum2.eps; 
end;