标题: 从乐谱中提取基音频率matlab代码 [打印本页]

---

作者: mingyang1996    时间: 2017-12-12 21:54
标题: 从乐谱中提取基音频率matlab代码
大家好，这是我们之前做的音乐识别dsp项目用到的乐谱F0提取matlab代码，先在matlab上面仿真后再迁移到CPU上面。看到有的人做dsp用到，发在这里共享。

源程序如下:

1. 
   
2. % [x,Fs] = audioread('DJF0.wav');
   
3. % % y=abs(hilbert(x));
   
4. % % plot(y);title('解调信号');
   
5. %
   
6. % % fs=30;
   
7. % t=0:1/Fs:13.306;
   
8. % x6=x;
   
9. % x66 = hilbert(x6);
   
10. % xx = abs(x66+j*x6);
    
11. % figure(1)
    
12. % hold on
    
13. % plot(t,x6);
    
14. % plot(t,xx,'r')
    
15. % xlim([0 5])
    
16. % hold off
    
17. % % 包络算法，未考虑边界条件
    
18. % d = diff(x6);
    
19. % n = length(d);
    
20. % d1 = d(1:n-1);
    
21. % d2 = d(2:n);
    
22. % indmin = find(d1.*d2<0 & d1<0)+1;
    
23. % indmax = find(d1.*d2<0 & d1>0)+1;
    
24. % envmin = spline(t(indmin),x6(indmin),t);
    
25. % envmax = spline(t(indmax),x6(indmax),t);
    
26. % figure
    
27. % hold on
    
28. % plot(t,x6);
    
29. % plot(t,envmin,'r');
    
30. % plot(t,envmax,'m');
    
31. % hold off
    
32. % xlim([0 5])
    
33. clc,clear;
    
34. % c=readmidi('潇洒走一回MIDI.mid');
    
35. c=readmidi('res.mid');
    
36. b=melcontour(c,0.25,'abs'); %% 0.25拍每个点，即每0.25beat取一个点
    
37. plot(b);%%画出旋律轮廓（melodic contour）,画出F0；
    
38. d=midi2hz(b);%把从midi格式代表的音高转换为频率，基频；
    
39. gettempo(c);%得到文件中音乐的bpm 即每分钟多少拍；从而可以计算出每秒多少拍，即音乐频率
    
40. %在“潇洒走一回”音乐中，得到1beat对应2Hz，
    
41. %则在midi文件中，因为每0.25beat取1个点，则计算得到对应8Hz（合成时需要的fs）
    
42. 
    
43. save mydata b;
    
44. 
    
45. 
    

复制代码

所有资料51hei提供下载:
tiqu.rar (778 Bytes, 下载次数: 27)

2017-12-12 21:53 上传

点击文件名下载附件
下载积分: 黑币 -5

---

欢迎光临 (http://www.51hei.com/bbs/)

Powered by Discuz! X3.1