iPhone音頻相移

是否可以在iPhone上實時移相音頻？iPhone音頻相移

我有一個簡單的應用程序設置，實時輸出麥克風的音頻。我試圖做的是處理這個音頻，即做相移。

我正確的唯一方法就是採取我的樣本空間，做一個FFT，相移，然後逆fft？

我知道vDSP庫，但它看起來像一個簡單的任務太多的開銷。

更新：我有一個來自Elec Eng的DSP基礎，是的，我真的想做一個相位轉換。我不需要進行頻率轉換或過濾，他們是我將在稍後實施的單獨過程。

來源

2012-07-31 Colin

你是指頻移？相移是毫無意義和瑣碎的。這只是原版的延遲版本。 – TJD 2012-07-31 21:38:04

是的，我想要相移。 – Colin 2012-08-01 19:12:37

如果您只想執行相移，則根本不需要任何DSP，因爲這只是輸入的時間延遲。只需讀取音頻數據並輸出即可 - 這會爲您帶來相位轉換。 – duskwuff 2012-08-08 17:33:40

是的，這是一種方法。

此示例C代碼表示它：

#include <stdio.h> 
#include <math.h> 

#ifndef M_PI 
#define M_PI 3.14159265358979324 
#endif 

typedef struct 
{ 
    double x, y; 
} tComplex; 

tComplex complexAdd(const tComplex* a, const tComplex* b) 
{ 
    tComplex c; 
    c.x = a->x + b->x; 
    c.y = a->y + b->y; 
    return c; 
} 

tComplex complexMul(const tComplex* a, const tComplex* b) 
{ 
    tComplex c; 
    c.x = a->x * b->x - a->y * b->y; 
    c.y = a->x * b->y + a->y * b->x; 
    return c; 
} 

void dft(tComplex out[], const tComplex in[], size_t n, int direction) 
{ 
    size_t k, i; 
    for (k = 0; k < n; k++) 
    { 
    tComplex r = { 0, 0 }, e; 
    for (i = 0; i < n; i++) 
    { 
     e.x = cos(-2 * direction * M_PI/n * ((double)k - n/2) * ((double)i - n/2)); 
     e.y = sin(-2 * direction * M_PI/n * ((double)k - n/2) * ((double)i - n/2)); 
     e = complexMul(&e, &in[i]); 
     r = complexAdd(&r, &e); 
    } 
    out[k] = r; 
    } 
} 

double maxAbs(const tComplex in[], size_t n) 
{ 
    double m = 0; 
    while (n--) 
    { 
    double a = hypot(in->x, in->y); 
    if (m < a) m = a; 
    in++; 
    } 
    return m; 
} 

#define SAMPLE_CNT 32 
#define SAMPLE_SHIFT 3 

int main(void) 
{ 
    tComplex signalIn[SAMPLE_CNT]; 
    tComplex signalOut[SAMPLE_CNT]; 
    tComplex tmp[SAMPLE_CNT]; 
    int i; 

    // signalIn[] = square pulse 
    for (i = 0; i < SAMPLE_CNT; i++) 
    { 
    signalIn[i].x = ((i - SAMPLE_CNT/2) >= 0) * ((i - SAMPLE_CNT/2) < SAMPLE_CNT/4); 
    signalIn[i].y = 0; 
    } 

    // tmp[] = DFT{signalIn[]} 
    dft(tmp, signalIn, SAMPLE_CNT, 1); 

    // tmp[] = DFT{signalIn[]} * exp(j * 2 * PI * f * TimeShift) 
    for (i = 0; i < SAMPLE_CNT; i++) 
    { 
    tComplex e; 
    e.x = cos(2 * M_PI * (i - SAMPLE_CNT/2) * SAMPLE_SHIFT/SAMPLE_CNT); 
    e.y = -sin(2 * M_PI * (i - SAMPLE_CNT/2) * SAMPLE_SHIFT/SAMPLE_CNT); 
    tmp[i] = complexMul(&tmp[i], &e); 
    } 

    // signalOut[] = IDFT{tmp[]} 
    dft(signalOut, tmp, SAMPLE_CNT, -1); 

    printf("  Re{In[]} .  Im{In[]} |" 
     " |DFT{In[]}| |" 
     " Re{Out[]} . Im{Out[]}\n"); 

    for (i = 0; i < SAMPLE_CNT; i++) 
    { 
    int j, s; 

    s = signalIn[i].x/maxAbs(signalIn, SAMPLE_CNT) * 8 + .5; 
    for (j = -8; j <= 8; j++) printf("%c", " *"[s == j]); 
    printf("."); 
    s = signalIn[i].y/maxAbs(signalIn, SAMPLE_CNT) * 8 + .5; 
    for (j = -8; j <= 8; j++) printf("%c", " *"[s == j]); 

    printf("|"); 

    s = hypot(tmp[i].x, tmp[i].y)/maxAbs(tmp, SAMPLE_CNT) * 8; 
    for (j = -8; j <= 8; j++) printf("%c", " *"[s == j]); 

    printf("|"); 

    s = signalOut[i].x/maxAbs(signalOut, SAMPLE_CNT) * 8 + .5; 
    for (j = -8; j <= 8; j++) printf("%c", " *"[s == j]); 
    printf("."); 
    s = signalOut[i].y/maxAbs(signalOut, SAMPLE_CNT) * 8 + .5; 
    for (j = -8; j <= 8; j++) printf("%c", " *"[s == j]); 

    printf("\n"); 
    } 

    return 0; 
}

輸出：

 Re{In[]} .  Im{In[]} | |DFT{In[]}| | Re{Out[]} . Im{Out[]} 
     *  .  *  |  *  |  *  .  *   
     *  .  *  |  *  |  *  .  *   
     *  .  *  |   *  |  *  .  *   
     *  .  *  |  *  |  *  .  *   
     *  .  *  |  *  |  *  .  *   
     *  .  *  |  *  |  *  .  *   
     *  .  *  |   *  |  *  .  *   
     *  .  *  |  *  |  *  .  *   
     *  .  *  |  *  |  *  .  *   
     *  .  *  |   *  |  *  .  *   
     *  .  *  |   *  |  *  .  *   
     *  .  *  |   *  |  *  .  *   
     *  .  *  |  *  |  *  .  *   
     *  .  *  |   *  |  *  .  *   
     *  .  *  |    * |  *  .  *   
     *  .  *  |    * |  *  .  *   
       *.  *  |    *|  *  .  *   
       *.  *  |    * |  *  .  *   
       *.  *  |    * |  *  .  *   
       *.  *  |   *  |    *.  *   
       *.  *  |  *  |    *.  *   
       *.  *  |   *  |    *.  *   
       *.  *  |   *  |    *.  *   
       *.  *  |   *  |    *.  *   
     *  .  *  |  *  |    *.  *   
     *  .  *  |  *  |    *.  *   
     *  .  *  |   *  |    *.  *   
     *  .  *  |  *  |  *  .  *   
     *  .  *  |  *  |  *  .  *   
     *  .  *  |  *  |  *  .  *   
     *  .  *  |   *  |  *  .  *   
     *  .  *  |  *  |  *  .  *

可以看到，輸出是由SAMPLE_SHIFT樣品（3，在這裏）在時間上偏移的輸入的複製品。

您也可以將SAMPLE_SHIFT更改爲非整數值，從而得到一個小數的樣本。

輸出，用於SAMPLE_SHIFT = 2.5：

 Re{In[]} .  Im{In[]} | |DFT{In[]}| | Re{Out[]} . Im{Out[]} 
     *  .  *  |  *  |  *  .  *   
     *  .  *  |  *  |  *  .  *   
     *  .  *  |   *  |  *  .  *   
     *  .  *  |  *  |  *  .  *   
     *  .  *  |  *  |  *  .  *   
     *  .  *  |  *  |  *  .  *   
     *  .  *  |   *  |  *  .  *   
     *  .  *  |  *  |  *  .  *   
     *  .  *  |  *  |  *  .  *   
     *  .  *  |   *  |  *  .  *   
     *  .  *  |   *  |  *  .  *   
     *  .  *  |   *  |  *  .  *   
     *  .  *  |  *  |  *  .  *   
     *  .  *  |   *  |  *  .  *   
     *  .  *  |    * |  *  .  *   
     *  .  *  |    * |  *  .  *   
       *.  *  |    *|  *  .  *   
       *.  *  |    * |  *  .  *   
       *.  *  |    * |   *  .  *   
       *.  *  |   *  |    *.  *   
       *.  *  |  *  |    * .  *   
       *.  *  |   *  |    * .  *   
       *.  *  |   *  |    * .  *   
       *.  *  |   *  |    * .  *   
     *  .  *  |  *  |    * .  *   
     *  .  *  |  *  |    *.  *   
     *  .  *  |   *  |   *  .  *   
     *  .  *  |  *  |  *  .  *   
     *  .  *  |  *  |  *  .  *   
     *  .  *  |  *  |  *  .  *   
     *  .  *  |   *  |  *  .  *   
     *  .  *  |  *  |  *  .  *

來源

2012-08-09 08:59:55