Tworzenie fali sinusoidalnej lub kwadratowej w C#

Question

Jak wygenerować sinusoidę lub falę prostokątną o danej częstotliwości?

Mam nadzieję, że zrobię to, aby skalibrować sprzęt, więc jak precyzyjne byłyby te fale?

Answer 1

Można użyć NAudio i utworzyć wyprowadzony WaveStream, który wysyła fale sinusoidalne lub prostokątne, które można wypisać na karcie dźwiękowej lub zapisać w pliku WAV. Jeśli używałeś 32-bitowych próbek zmiennoprzecinkowych, możesz zapisać wartości bezpośrednio z funkcji sin, bez konieczności skalowania, ponieważ jest to już od -1 do 1.

Co do dokładności, masz na myśli dokładnie odpowiednią częstotliwość, lub dokładnie właściwy kształt fali? Nie ma czegoś takiego jak prawdziwa fala prostokątna, a nawet fala sinusoidalna będzie prawdopodobnie miała kilka bardzo cichych artefaktów na innych częstotliwościach. Jeśli liczy się dokładność częstotliwości, jesteś zależny od stabilności i dokładności zegara w twojej karcie dźwiękowej. Powiedziawszy to, wyobrażam sobie, że dokładność byłaby wystarczająca dla większości zastosowań.

Oto niektóre przykładowy kod, który sprawia, że ​​próbkę 1   kHz z częstotliwością próbkowania 8   kHz i 16 próbek bitowych (to znaczy nie zmiennoprzecinkowe):

int sampleRate = 8000; 
short[] buffer = new short[8000]; 
double amplitude = 0.25 * short.MaxValue; 
double frequency = 1000; 
for (int n = 0; n < buffer.Length; n++) 
{ 
    buffer[n] = (short)(amplitude * Math.Sin((2 * Math.PI * n * frequency)/sampleRate)); 
} 

Answer 2

To pozwala dać częstotliwość, czas trwania i amplituda, i jest to 100% .NET CLR kodu. Brak zewnętrznych bibliotek DLL. Działa poprzez utworzenie formatu MemoryStream w formacie WAV, który przypomina tworzenie pliku tylko w pamięci, bez przechowywania go na dysku. Następnie odtwarza to MemoryStream z System.Media.SoundPlayer.

using System; 
using System.Collections.Generic; 
using System.IO; 
using System.Linq; 
using System.Windows.Forms; 

public static void PlayBeep(UInt16 frequency, int msDuration, UInt16 volume = 16383) 
{ 
    var mStrm = new MemoryStream(); 
    BinaryWriter writer = new BinaryWriter(mStrm); 

    const double TAU = 2 * Math.PI; 
    int formatChunkSize = 16; 
    int headerSize = 8; 
    short formatType = 1; 
    short tracks = 1; 
    int samplesPerSecond = 44100; 
    short bitsPerSample = 16; 
    short frameSize = (short)(tracks * ((bitsPerSample + 7)/8)); 
    int bytesPerSecond = samplesPerSecond * frameSize; 
    int waveSize = 4; 
    int samples = (int)((decimal)samplesPerSecond * msDuration/1000); 
    int dataChunkSize = samples * frameSize; 
    int fileSize = waveSize + headerSize + formatChunkSize + headerSize + dataChunkSize; 
    // var encoding = new System.Text.UTF8Encoding(); 
    writer.Write(0x46464952); // = encoding.GetBytes("RIFF") 
    writer.Write(fileSize); 
    writer.Write(0x45564157); // = encoding.GetBytes("WAVE") 
    writer.Write(0x20746D66); // = encoding.GetBytes("fmt ") 
    writer.Write(formatChunkSize); 
    writer.Write(formatType); 
    writer.Write(tracks); 
    writer.Write(samplesPerSecond); 
    writer.Write(bytesPerSecond); 
    writer.Write(frameSize); 
    writer.Write(bitsPerSample); 
    writer.Write(0x61746164); // = encoding.GetBytes("data") 
    writer.Write(dataChunkSize); 
    { 
     double theta = frequency * TAU/(double)samplesPerSecond; 
     // 'volume' is UInt16 with range 0 thru Uint16.MaxValue (= 65 535) 
     // we need 'amp' to have the range of 0 thru Int16.MaxValue (= 32 767) 
     double amp = volume >> 2; // so we simply set amp = volume/2 
     for (int step = 0; step < samples; step++) 
     { 
      short s = (short)(amp * Math.Sin(theta * (double)step)); 
      writer.Write(s); 
     } 
    } 

    mStrm.Seek(0, SeekOrigin.Begin); 
    new System.Media.SoundPlayer(mStrm).Play(); 
    writer.Close(); 
    mStrm.Close(); 
} // public static void PlayBeep(UInt16 frequency, int msDuration, UInt16 volume = 16383) 

Answer 3

Try od Creating sine and save to wave file in C#

private void TestSine() 
{ 
    IntPtr format; 
    byte[] data; 
    GetSineWave(1000, 100, 44100, -1, out format, out data); 
    WaveWriter ww = new WaveWriter(File.Create(@"d:\work\sine.wav"), 
     AudioCompressionManager.FormatBytes(format)); 
    ww.WriteData(data); 
    ww.Close(); 
} 

private void GetSineWave(double freq, int durationMs, int sampleRate, short decibel, out IntPtr format, out byte[] data) 
{ 
    short max = dB2Short(decibel);//short.MaxValue 
    double fs = sampleRate; // sample freq 
    int len = sampleRate * durationMs/1000; 
    short[] data16Bit = new short[len]; 
    for (int i = 0; i < len; i++) 
    { 
     double t = (double)i/fs; // current time 
     data16Bit[i] = (short)(Math.Sin(2 * Math.PI * t * freq) * max); 
    } 
    IntPtr format1 = AudioCompressionManager.GetPcmFormat(1, 16, (int)fs); 
    byte[] data1 = new byte[data16Bit.Length * 2]; 
    Buffer.BlockCopy(data16Bit, 0, data1, 0, data1.Length); 
    format = format1; 
    data = data1; 
} 

private static short dB2Short(double dB) 
{ 
    double times = Math.Pow(10, dB/10); 
    return (short)(short.MaxValue * times); 
}