再舉一例，圖2

代碼在這里，僅僅修改幾個參數而已。頻率為100,500,800的三個正弦波，采樣為10k，基音周期不理想，低于直接FFT

nt=4000;

fs=10000;

t=(0:nt-1)*1/fs;

f1=100;

f2=500;

f3=800;

x=0.6*sin(2*pi*f1*t)+0.8*sin(2*pi*f2*t)+1.2*sin(2*pi*f3*t);

x=x-mean(x);? ?? ?? ?? ?? ?? ?? ?? ?? ?? ???% 消去直流分量

x=x/max(abs(x));? ?? ?? ?? ?? ?? ?? ?? ?? ? % 幅值歸一化

wlen=320; inc=80;? ?? ?? ?? ???% 分幀的幀長和幀移

overlap=wlen-inc;? ?? ?? ?? ???% 幀之間的重疊部分

T1=0.05;? ?? ?? ?? ?? ?? ?? ???% 設置基音端點檢測的參數

y??= enframe(x,wlen,inc)';? ?? ?? ?? ?? ?? ?% 分幀

fn??= size(y,2);? ?? ?? ?? ?? ?? ?? ?? ?? ? % 取得幀數

[voiceseg,vosl,SF,Ef]=pitch_vad1(y,fn,T1);

frameTime = frame2time(fn, wlen, inc, fs);??% 計算各幀對應的時間坐標

lmin=fix(fs/1000);? ?? ?? ?? ?? ?? ?? ?? ?? ?% 基音周期的最小值

lmax=fix(fs/60);? ?? ?? ?? ?? ?? ?? ?? ?? ? % 基音周期的最大值

period=zeros(1,fn);? ?? ?? ?? ?? ?? ?? ?? ? % 基音周期初始化

for k=1:fn

if SF(k)==1? ?? ?? ?? ?? ?? ?? ?? ?? ???% 是否在有話幀中

y1=y(:,k).*hamming(wlen);? ?? ?? ???% 取來一幀數據加窗函數

xx=fft(y1);? ?? ?? ?? ?? ?? ?? ?? ? % FFT

a=2*log(abs(xx)+eps);? ?? ?? ?? ?? ?% 取模值和對數

%B(:,k)=ifft(a);

b=ifft(a);? ?? ?? ?? ?? ?? ?? ?? ???% 求取倒譜

[R(k),Lc(k)]=max(b(lmin:lmax));? ???% 在lmin和lmax區間中尋找最大值

period(k)=Lc(k)+lmin-1;? ?? ?? ?? ? % 給出基音周期

end

end

% 作圖

subplot 311, plot(t,x,'k');??title('語音信號')

axis([0 max(t) -1 1]); ylabel('幅值');

subplot 312; plot(frameTime,period,'k');

xlim([0 max(t)]);ylim([0,150]); title('基音周期');

grid; xlabel('時間/s'); ylabel('樣點數');

subplot 313; plot(frameTime,fs./period,'k');

xlim([0 max(t)]);ylim([0,1000]); title('基音周期');

grid; xlabel('時間/s'); ylabel('頻率');

圖2.jpg

(107.72 KB, 下載次數: 0)

2018-8-29 16:53 上傳

第二個圖

總結

以上是生活随笔為你收集整理的matlab中frame是什么意思,frame2time结果是什么含义？的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。