语言信号处理综合实习匹配算法

However,ifthelengthofxisdifferentfromy,thenwecannotusetheaboveformulatocomputethedistance.Instead,weneedamoreflexiblemethodthatcanfindthebestmap fromelementsinxtothoseinyinordertocomputethedistance. 须靠着一套有效的运算方法，才可以DTW是DynamicTimeWar的简称，中文（DynamicProgrammingDP）的方法，Lettandrbetwovectorsoflengthsmandn,respectively.ThegoalofDTWistofinda path{(p1,q1),(p2,q2),...,(pk,qk)}suchthatthedistanceonthismap pathSi=1k∣t(pi)-r(qi)∣isminimized,withthefollowingtrm和n，那么DTW的目标，就是要找到一组路径(p1,q1),(p2,q2),...,(pk,qk)},使得经由上述路径的「点对点」对应距离和Si=1k∣t(pi)-r(qi)∣为Boundaryconditions:(p1,q1)=(1,1),(pk,qk)=(m,n).Thisisatypicalexampleof"anchoredbeginning"and"anchoredend".端点关系：(p1,q1)=(1,1),(pk,qkm,n)。此Localconstraint:Foranygivennode(i,j)inthepath,thepossiblefan-innodesarerestrictedto(i-1,j),(i,j-1),(i-1,j-1).Thislocalconstraintguaranteesthatthemap pathismonotonicallynon-decreasinginitsfirstandsecondarguments.Moreover,foranygivenelementint,weshouldbeabletofindatleastonecorrespondingelementinr,andviceversa.j)，那么其前一点路径只有三种可能：(i-1,j),(i,而且也规定了t的任一个元素至少对应一个r的寻求最佳路径三个步目标函数之定义：定义D(i,j)是t(1:i)和r(1:j)之间的DTW距离，对应的最佳路径是由(11(i,j)。D(i,D(i,j)t(i)r(j) D(i1,jD(i1,j2)最后答案：D(m,D(ij该表示成为D(t(1:ir(1:j))，才能准确地描述D(‧,‧)和t、r的关系在实际运算时，我们通常事先建立一个矩阵，其维度为m×n，先根据端点条件来填入D(1,D(i,j)的值，最后就可以得到我们所要的答案D(m,n)。有上述方法可得知，DTW的计算复杂度大约是O(mn)，比用 知道DTW距离之外，也希望把相关最佳路径出D(m,n)，在反复回推前一个最佳路径的位这个步骤在DP里面称为BackTracking。另一个常用到的localpathconstraint27°-D(D(i,j)t(i)r(j) D(i1,jD(i1,j2)DD(i,j)w(j)t(i)r(j) D(i1,jD(i1,j2)观察››在音符后半部，用户的音高较稳定且WeightedDTW›在音符开始时，权重函数w(j)›在音符后半部，权重函数w(j)较此种localpathconstraint，会DTW路径将会跳掉此2763应Iftherequiredmap pathis"anchoredbeginning,anchoredend",thenthelocalpathconstraintscaninducetheglobalpathconstraints,asshownnext:若要使DTW比对有意义，那么两者长度的比值必须介0.52.0之间，否则我们就不够定义globalpathconstraint，那么在进行DTW计算时，就可以省掉很多不需要计算的时，那么对应的globalpathconstraintTypicalDTWtablesize=128x对应的globalpathconstraint如下：TypicalDTWtablesize128xType-1DTWuses27°-45°-63°localpathType-2DTWuses0°-45°-90°localpathvec1=[7173758080807876757371717175767668767675737170706968687274798080vec2=[69697375798079787673727170696969717375767676767675737170707375808080[minDist,dt th,dtwTable]=dtw(vec1,vec2, thPlot(vec1,vec2, th1,dtwTable1]=dtw(vec1,vec2,[minDist2,dt th2,dtwTable2]=