fork download
  1. > #ライブラリの読み込み。wavファイルのプレイソフト名の登録
  2. > library(tuneR)
  3. tuneR >= 1.0 has changed its Wave class definition.
  4. Use updateWave(object) to convert Wave objects saved with previous versions of tuneR.
  5. > library(seewave)
  6.  
  7. ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  8. Welcome to seewave !
  9. The package is regularly updated, please check for new version [http://r...content-available-to-author-only...n.fr/seewave]
  10. Thanks to use the right reference when citing seewave in publications
  11. See citation('seewave')
  12. ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  13.  
  14. > setWavPlayer("play")
  15. >
  16. > #波形作成
  17. > t <- seq(0,2,1/44100) # 時間 44100サンプル/秒 のデータを2秒分
  18. > s1 <- (2^15-1)*sin(2*pi*t*440)*(1:length(t))/length(t) # 440Hzのsin波を時間とともに振幅線形的に増大
  19. > s2 <- (2^15-1)*sin(2*pi*t*880)*(length(t):1)/length(t) # 880Hzのsin波を時間とともに振幅線形的に減少
  20. >
  21. > # 音データを作成
  22. > w <- Wave(s1+s2,bit=16)
  23. 警告メッセージ:
  24. In .local(left, ...) : 'samp.rate' not specified, assuming 44100Hz
  25. > writeWave(w,"test002.wav") # wavファイルとして書きだし。
  26. 警告メッセージ:
  27. In writeWave(w, "test002.wav") :
  28. channels' data will be rounded to integers for writing the wave file
  29. >
  30. > # 演奏
  31. > play(w)
  32.  
  33. /tmp/RtmpkVHVFa/tuneRtemp.wav:
  34.  
  35. File Size: 176k Bit Rate: 706k
  36. Encoding: Signed PCM
  37. Channels: 1 @ 16-bit
  38. Samplerate: 44100Hz
  39. Replaygain: off
  40. Duration: 00:00:02.00
  41.  
  42.  
  43. In:0.00% 00:00:00.00 [00:00:02.00] Out:0 [ | ] Clip:0
  44. In:46.4% 00:00:00.93 [00:00:01.07] Out:41.0k [!=====|=====!] Clip:0
  45. In:55.7% 00:00:01.11 [00:00:00.89] Out:49.2k [======|======] Hd:0.0 Clip:0
  46. In:65.0% 00:00:01.30 [00:00:00.70] Out:57.3k [======|======] Hd:0.0 Clip:0
  47. In:74.3% 00:00:01.49 [00:00:00.51] Out:65.5k [======|======] Hd:0.0 Clip:0
  48. In:83.6% 00:00:01.67 [00:00:00.33] Out:73.7k [!=====|=====!] Hd:0.0 Clip:0
  49. In:92.9% 00:00:01.86 [00:00:00.14] Out:81.9k [!=====|=====!] Hd:0.0 Clip:0
  50. In:100% 00:00:02.00 [00:00:00.00] Out:88.2k [!=====|=====!] Hd:0.0 Clip:0
  51. In:100% 00:00:02.00 [00:00:00.00] Out:88.2k [ | ] Hd:0.0 Clip:0
  52. Done.
  53. 警告メッセージ:
  54. In writeWave(object, filename) :
  55. channels' data will be rounded to integers for writing the wave file
  56. >
  57. > # スペクトル解析
  58. > png('test002.png');
  59. > d <- spectro(w,osc=T);
  60. > title('test002');
  61. > dev.off()
  62. null device
  63. 1
  64. >
  65. > # spectroのかえりち。このなかのampを使いたい。
  66. > str(d)
  67. List of 3
  68. $ time: num [1:172] 0 0.0117 0.0234 0.0351 0.0468 ...
  69. $ freq: num [1:256] 0 0.0861 0.1723 0.2584 0.3445 ...
  70. $ amp : num [1:256, 1:172] -74.2 -72.9 -69.9 -65.9 -53.1 ...
  71. > a<-d["amp"] # or a<-d[3]
  72.  
Not running #stdin #stdout 0s 0KB
stdin
Standard input is empty
stdout
Standard output is empty