从R驱动声卡

Question

我正在做一个利用声卡播放声音并从多通道声卡输出电信号的实验。从MatLab，我可以直接驱动每个通道，但是我更喜欢使用R。我已经能够制作一个多通道文件并命名通道，但是我无法使其正常播放。我已将声卡设置为默认值，并且在“声音”窗口的配置窗口中播放测试时，但是在运行代码时却无法播放。该代码将打开WMP，并播放一个文件，该文件的长度与我创建的信号的正确长度相同，但是我没有从扬声器中听到声音，也无法从已设置的其他通道中测量电信号。 Windows 7的 声卡：Lynx E44 R 3.4.4

library("seewave")
library("tuneR")
s1 <- 10*sin(2*pi*440*seq(0,1,length.out=8000))
s1_multi <- WaveMC(data = cbind(s1, s1, s1, s1),
               samp.rate=8000, 
               bit=16)
 colnames(s1_multi) <- c("FL", "FR", "FC", "LF")
setWavPlayer(shQuote("C:/Program Files (x86)/Windows Media Player/wmplayer.exe"))

listen(s1_multi)

我做错了什么吗？有没有办法直接从R驱动声卡的通道？

---

我在Seewave开发者留言板上发布了相同的问题，并被告知为什么我的确切代码不起作用（WaveMC的listen（）仅播放第一个频道）。但是我更大的问题是，是否有一种方法可以直接从R驱动声卡，或者是否必须通过第三方播放器。到目前为止，我已经尝试过tuneR和seewave。

谢谢。

Answer 1

您可以尝试使用audio软件包。

带有R的星球大战帝国三月

让我们写下一些笔记和持续时间：

p <- "A3 
      A3 A3 A3 F3 C3# A3 F3 C3# A3
      F4 F4 F4 F4# C3# G3 F3 C3# A3"

d <- c(2, 
       1, 1, 1, 0.75, 0.25, 1, 0.75, 0.25, 2,
       1, 1, 1, 0.75, 0.25, 1, 0.75, 0.25, 2)

制作音频波数据帧：

w <- make_wave(pitch = p, duration = d, tempo = 120, sample_rate= 44100)

这使用的是Nick Kennedy's phantastic function的稍微改编的版本。

...并直接从控制台播放它：

audio::play(w)

---

说明

您基本上需要的是带有频率和持续时间的表格。这也可以通过查找表来实现。获得每个音符的频率和持续时间后，您可以使用make_sine()将其轻松转换为音频波并进行播放。

注意：tuneR还具有一个称为play()的函数，该函数的工作原理不同，并且可能由于屏蔽而导致错误。

---

功能

make_wave <- function(pitch, duration, tempo, sample_rate){
    # Credit to Nick Kennedy
    require(dplyr)

    notes <- c(A = 0, B = 2, C = 3, D = 5, E = 7, F = 8, G = 10)

    # remove leading, trailing, and consecutive spaces (and some line breaks)
    pitch <- gsub("\n", " ", pitch)
    pitch <- gsub("(?<=[\\s])\\s*|^\\s+|\\s+$", "", pitch, perl=TRUE)

    # split notes, add duration
    x <- data_frame(pitch = strsplit(pitch, " ")[[1]],
                    duration = duration)

    # calculate frequency
    x <-
        x %>%
        mutate(octave = substring(pitch, nchar(pitch)) %>%
        {suppressWarnings(as.numeric(.))} %>%
            ifelse(is.na(.), 4, .),
        note = notes[substr(pitch, 1, 1)],
        note = note + grepl("#", pitch) -
            grepl("b", pitch) + octave * 12 +
            12 * (note < 3),
        freq = 2 ^ ((note - 60) / 12) * 440)


    # create fade to remove clipping
    fade <- seq(0, 1, 50 / sample_rate)


    # turn frequency and duration into sound
    make_sine <- function(freq, duration) {
        # create audio wave
        wave <- sin(seq(0, duration / tempo * 60, 1 / sample_rate) *
                        freq * 2 * pi)

        # apply fade to audio wave
        wave * c(fade, rep(1, length(wave) - 2 * length(fade)), rev(fade))
    }

    # create waveform
    wave_df <-
        mapply(make_sine, x$freq, x$duration) %>%
        do.call("c", .)

    return(wave_df) 
}