設(shè)計(jì)或其他-語(yǔ)言信號(hào)處理chap10作業(yè)

Chap10作(**)切音符的函數(shù)myNoteSegment.m，其note=myNoteSegment(pitchVec,timeUnit,pitchTh,pitchVec:輸入的音高向量timeUnitpitchTh：note:切出的音符，以「音高,音長(zhǎng),音高,音長(zhǎng)...」的格式來，其中音高的單位是Semitone，音長(zhǎng)的單位是1/64秒。，的長(zhǎng)度可以小于minNoteDuration。如何測(cè)試此函數(shù)：wave檔案以及midi。請(qǐng)使用pvPlay.m來原始音高向量，并用notePlay來播(**) 的函數(shù):upperRatio,resolution,distanceType):pitchVec2數(shù)據(jù)庫(kù)中的音高向量lowerRatio:最低伸縮比例upperRatio:最高伸縮比例resolution:伸縮總次數(shù)distanceType:所用的距離函數(shù)，1L1norm，2L2norm。minDist:線性伸縮的最短距離bestVec1pitchVec1allDist:線性伸縮的所有距離vec1vec2，則不再計(jì)算距離。linScalingMexmyLinearScaling，看看是否得到同提示：可用interp1指令來進(jìn)行內(nèi)差，median指令來計(jì)算中位數(shù)，mean指令來計(jì)算平均值，normLpnorminterp1.m的使用有疑問，Example1:Originalvector原始向量%Thevectorinterpolationshouldhavealengthequalto1.8theoriginal %Lengthofthescaledvectorx2=linspace(1,n,m);y2=interp1(x,y,x2);%Vectorafterinterpolationsubplot(2,1,1),plot(x,y,'o-',x2,y2,'o:');legend('Originalvector','Interpolatedvector','location','northOutside','orientation',subplot(2,1,2),plot(1:length(y),y,'o-',1:length(y2),y2,legend('Originallegend('Originalvector','Scaledvector','location','northOutside','orientation',fprintf('Desiredlengthratio=%f\n',fprintf('Actuallengthratio=%f\n',Desiredlengthratio=1.750000Actuallengthratio=1.800000(***)DTW3程序代碼:請(qǐng)寫一個(gè) 函數(shù)myDtw3.m來實(shí)做type-3[minDist, th,dtwTable]=myDtw3(vec1, %dtw3:Dynamictime withlocalpathsof0and45 endCorner,plotOpt)%%vec1:testingvector,whichshouldbeapitchvec2:referencevector,whichshouldbeavectorofnote minDistance:minimundistanceof th:optimalpathofdynamicalprogrammingthroughtheDTW (Itssizeis2xk,wherekisthepath ifnargin<4,endCorner=1;end%Ifinputisvector,makeitrowifsize(vec1,1)==1|size(vec1,2)==1,vec1=vec1(:)';endifsize(vec2,1)==1|size(vec2,2)==1,vec2=vec2(:)';%======ConstructDTWtable%======ConstructthefirstelementoftheDTWtabledtwTable(1,1)=vecDist(vec1(:,1),vec2(:,1));%======ConstructthefirstrowoftheDTWtable(xyview)fori=2:size1%======ConstructthefirstcolumnoftheDTWtable(xyview)ifbeginCorner==0fordtwTable(1,j)=vecDist(vec1(:,1),%======ConstructalltheotherrowsofDTWtablefori=2:size1forpointDist=vecDist(vec1(:,i),%======Check45-degree%======Check0-degree%======Findtheoveralloptimum%==========Subfunctionfunctiondistance=vecDist(x,y)(***)DTW3程序代碼：改進(jìn)方案:請(qǐng)寫一個(gè)函數(shù)來重復(fù)上一題，(***)旋律辨識(shí)的效能測(cè)試之一:請(qǐng)先閱讀本節(jié)的最后一個(gè)習(xí)題，充分了解如LSresolution為可變參數(shù)，使其從5增加到12。畫出的圖應(yīng)該類似下圖：調(diào)移位次數(shù)為可變參數(shù)（此參數(shù)也是每次比對(duì)一首wav時(shí)，所必須執(zhí)行DTW的個(gè)數(shù)，請(qǐng)使值為1,3,5,7,9。畫出的圖應(yīng)該類似下圖：LSDTW在辨識(shí)效能（含辨識(shí)率及辨 (***)旋律辨識(shí)的