dtmf_fsk音频数据测试

最新推荐文章于 2024-06-09 13:40:41 发布
最新推荐文章于 2024-06-09 13:40:41 发布
阅读量510 收藏 3
点赞数
分类专栏： test 文章标签： 音视频 c语言 c++
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_45698138/article/details/129876706
版权

文章目录

音频测试

1.dtmf调制成音频数据

  1. 输入任意按键值，输出是该按键对应的dtmf调制的音频信号，输出结果保存到文件中。
  2. 要求dtmf使用的基波响度为-12dB，采样率为16000Hz，生成的音频长度为100ms。

test.c

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

#define FILE_NAME_SIZE_MAX  50
#define SAMPLE_RATE         16000 
#define AMPLITUDE           32767 
#define PI                  3.14159265358979323846
#define AUDIO_TIME_SEC_LEN  (0.1)
#define AUDIO_LOUDNESS_DBFS (-12) 

typedef struct {
    char key;
    int freq1;
    int freq2;
} dtmf_st;

static const dtmf_st dtmf_table[] = {
    {'1', 697, 1209}, {'2', 697, 1336}, {'3', 697, 1477}, {'A', 697, 1633},
    {'4', 770, 1209}, {'5', 770, 1336}, {'6', 770, 1477}, {'B', 770, 1633},
    {'7', 852, 1209}, {'8', 852, 1336}, {'9', 852, 1477}, {'C', 852, 1633},
    {'#', 941, 1209}, {'0', 941, 1336}, {'D', 941, 1477}
};

static const int dtmf_table_size = sizeof(dtmf_table) / sizeof(dtmf_st);

int main(void) {
    char key = 0;
    short *audio_data = NULL;
    int num_samples = 0;
    int freq1 = -1, freq2 = -1;
    double dtmf_amp = 0;
    double time_in_sec = 0;
    double sin_wave1 = 0, sin_wave2 = 0;
    char filename[FILE_NAME_SIZE_MAX] = {0};
    FILE *fp = NULL;
    int i = 0;

    while(1) {
        printf("\nPlease input key:");
        if(scanf("%c", &key) != 1) {
            while(getchar() != '\n');
            printf("Input error! \nPlease input (0-9,A-D,#,*)\n");
            continue;
        }

        for(i=0;i<dtmf_table_size; i++) {
            if(key == dtmf_table[i].key) {
                break;
            }
        }
        if(i < dtmf_table_size) {
            freq1 = dtmf_table[i].freq1;
            freq2 = dtmf_table[i].freq2;
            // printf("%d %d\n",freq1,freq2);
            break;
        } 
        else {
            printf("Invalid key!\n");
        }
    }

    num_samples = SAMPLE_RATE * AUDIO_TIME_SEC_LEN;

    audio_data = malloc(sizeof(short) * num_samples);
    if(audio_data == NULL) {
        goto ERROR;
    }

    dtmf_amp = pow(10, AUDIO_LOUDNESS_DBFS/20.0) * AMPLITUDE;
    for (int i = 0; i < num_samples; i++) {
        time_in_sec = i / (double) SAMPLE_RATE;
        sin_wave1 = sin(2 * PI * freq1 * time_in_sec);
        sin_wave2 = sin(2 * PI * freq2 * time_in_sec);
        audio_data[i] = dtmf_amp * (sin_wave1 + sin_wave2);
    }

    sprintf(filename, "key_%c.raw", key);
    fp = fopen(filename, "wb");
    if (!fp) {
        perror("fopen");
        goto ERROR;
    }

    if(fwrite(audio_data, sizeof(short), num_samples, fp) < 0)
    {
        perror("fwrite");
        goto ERROR;
    }

ERROR:
    if (fp) {
        fclose(fp);
        fp = NULL;
    }
    if (audio_data) {
        free(audio_data);
        audio_data = NULL;
    }
    return 0;
}

sound forge打开

在这里插入图片描述

2.使用fsk调成音频数据

使用C语言编写一段程序，程序可以实现将输入的16进制数据使用fsk调制成音频数据。

1.输入为一个16进制数，输出为一个音频文件。例如输入0x55，对应的二进制为01010101，fsk调制后的频率组合为 低频-高频-低频-高频-低频-高频-低频-高频。

2.音频采样率为16000Hz。

3.fsk调制时使用的两个信号分别为1209Hz和697Hz。

4.fsk使用的调制信号的幅值为-12dB。

5.fsk调制时每个频率的信号长度为10ms。

6.两个不同频率的信号拼接时要注意保持相位上的连续性。

在linux系统中，数据保存为raw格式

test.c

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <string.h>
#include <ctype.h>
#include <stdint.h>
#include <stdbool.h>

#define AMPLITUDE             32767.0 
#define SAMPLE_RATE           (16000.0)                    
#define LOW_FREQ              697                         
#define HIGH_FREQ             1209                       
#define AUDIO_LOUDNESS_DBFS   (-12.0) 
#define SIGNAL_LENGTH         0.01                   
#define BIN_STR_SIZE_MAX      64
#define PI                    3.14159265358979323846            
#define OUTPUT_FILE           "output.raw"

struct wave {
    int freq;               
    double amplitude;       
};
static const struct wave waves[2] = {
        {LOW_FREQ, pow(10, AUDIO_LOUDNESS_DBFS / 20.0)*AMPLITUDE}, 
        {HIGH_FREQ, pow(10, AUDIO_LOUDNESS_DBFS / 20.0)*AMPLITUDE}
    };

static const char *hex_table[] = {"0000", "0001", "0010", "0011", "0100",
                                  "0101", "0110", "0111", "1000", "1001",
                                  "1010", "1011", "1100", "1101", "1110",
                                  "1111"};

void hex_to_bin(char *hex_string, char *bin_string) 
{
    int i = 0;
    char hex_digit = 0;
    for (i = 0; hex_string[i] != '\0'; i++) {
        hex_digit = hex_string[i];
        if (hex_digit >= '0' && hex_digit <= '9') {
            strcat(bin_string, hex_table[hex_digit - '0']);
        } else if (hex_digit >= 'A' && hex_digit <= 'F') {
            strcat(bin_string, hex_table[hex_digit - 'A' + 10]);
        } else if (hex_digit >= 'a' && hex_digit <= 'f') {
            strcat(bin_string, hex_table[hex_digit - 'a' + 10]);
        } else {
            printf("error invalid digit '%c'\n", hex_digit);
            goto ERROR;
        }
    }

    if (strlen(bin_string) == 0) {
        printf("error input is null\n");
        goto ERROR;
    }

    if (strlen(bin_string) % 4 != 0) {
        printf("error：invalid len %d\n", (int)strlen(bin_string));
        goto ERROR;
    }
ERROR:
    return ;
}


double generate_fsk_signal(int fi,int freq,int amplitude,int num_samples,short *audio_data) 
{
    static int pos = 0;
    double t = 0;
    for (int i = 0; i < num_samples; i++)
    {
        audio_data[pos++] = (short)(amplitude * sin(2 * PI * freq * t + fi));
        t += 1/SAMPLE_RATE;
    }
    return (2 * PI * freq * (t+1/SAMPLE_RATE)  + fi); //相位对齐
}

void generate_sin_wave(int n,char *bin_string,short *audio_data,int num_samples) 
{
    int len = strlen(bin_string);
    int i=0 ,j=0;
    double fi = 0;
    int index = 0;
    for(i=0;i<len;i++)
    {
        index = bin_string[i] - '0';
        if(index == 0) {
            fi = generate_fsk_signal(fi,waves[index].freq,waves[index].amplitude,num_samples,audio_data);
        } 
        else if(index == 1) {
            fi = generate_fsk_signal(fi,waves[index].freq,waves[index].amplitude,num_samples,audio_data);
        }
    }  
}



char *get_hex_str(char *str)
{
    int len = strlen(str);

    if(str == NULL) {
        goto ERROR;
    }

    if (len == 0 || len > BIN_STR_SIZE_MAX / 4) {
        goto ERROR;
    }

    if(strstr(str,"0x") || strstr(str,"0X") ) {
        str += 2;
    }

    for(int i=0; i<strlen(str); i++){
        if(!isxdigit(str[i])){
            printf("error: invalid hex_string\n");
            goto ERROR;
        }
    }

    return str;
ERROR:
    return NULL;
}

int main(int argc, char **argv) 
{
    int i = 0, j = 0;
    int index = 0;
    int num_samples = 0;
    FILE *fp = NULL;
    int signal_length = 0;
    char *hex_str = NULL;
    double last_fi = 0;
    char bin_string[BIN_STR_SIZE_MAX] = {0};
    short *audio_data = NULL;
    
    if (argc < 2) {
        printf("Usage: %s hex_string\n", argv[0]);
        goto ERROR;
    }

    hex_str = get_hex_str(argv[1]);
    if(hex_str == NULL){
        goto ERROR;
    }

    hex_to_bin(hex_str, bin_string);
    num_samples = SIGNAL_LENGTH * SAMPLE_RATE;
    signal_length = strlen(bin_string) * num_samples;
    audio_data = malloc(sizeof(short) * signal_length);
    memset(audio_data,0,sizeof(short) * signal_length);
  
    generate_sin_wave(signal_length,bin_string,audio_data,num_samples);

    fp = fopen(OUTPUT_FILE, "wb");
    if (fp == NULL) {
        perror("fopen");
        goto ERROR;
    }

    if(fwrite(audio_data, sizeof(short), signal_length, fp) < 0) {
        perror("fwrite");
        goto ERROR;
    }
ERROR:
    if (fp) {
        fclose(fp);
        fp = NULL;
    }
    if (audio_data) {
        free(audio_data);
        audio_data = NULL;
    }

    return 0;
}

在这里插入图片描述

相关知识：

dtmf:

在这里插入图片描述

响度:

响度常用的单位是dBSPL、dBFS、dBpa、dBm0。

dB计算：对于功率，dB = 10 * lg(A/B)。对于电压或电流，dB= 20 * lg(A/B)。此处A，B代表参与比较的功率值或者电流、电压值。

以dBfs举例说明：

dBfs = 20 * lg(x/32767)。

其中32767是16bit采样精度时最大的值，x取值0 ~ 32767。所以dBfs的最大值为当x取值为32767时的值，即0dBfs。

采样率：

采样率是指对模拟音频信号进行取样的每秒次数，通常以赫兹 (Hz) 为单位表示。当我们要将模拟音频信号转换成数字形式时（例如录制音频或将其从 CD 转换成数字格式），需要对其进行采样，即在时间轴上不断地对音频信号进行抽样，获取一系列数字值。采样率越高，每秒中采集到的样本数量就越多，这可以提供更精确的音频表示，但同时也会增加所需的存储空间和处理能力。

振幅：

振幅是指波形的最大偏移量，通常用单位分贝 (dB) 或电平（以电压为单位）来表示。在音频信号中，振幅描述了声音的响度或强度。例如，在一段声音信号中，振幅越高，则声音会听起来更响亮、更强烈，振幅越低，则声音会听起来更柔和、更安静。与采样率不同，振幅是代表音频信号能量的单一数值，其大小与声音响度直接相关。

正弦波：

y = A * sin(ωt + ϕ)

其中，A代表振幅，即波形的最大偏移量；

ω代表角频率，即单位时间内的变化角度；

t代表时间；

ϕ代表初相位，即初始状态下正弦波图像的左右平移位置

φ = arcsin(y/A) - wt

每秒中采集到的样本数量就越多，这可以提供更精确的音频表示，但同时也会增加所需的存储空间和处理能力。

振幅：

振幅是指波形的最大偏移量，通常用单位分贝 (dB) 或电平（以电压为单位）来表示。在音频信号中，振幅描述了声音的响度或强度。例如，在一段声音信号中，振幅越高，则声音会听起来更响亮、更强烈，振幅越低，则声音会听起来更柔和、更安静。与采样率不同，振幅是代表音频信号能量的单一数值，其大小与声音响度直接相关。

正弦波：

y = A * sin(ωt + ϕ)

其中，A代表振幅，即波形的最大偏移量；

ω代表角频率，即单位时间内的变化角度；

t代表时间；

ϕ代表初相位，即初始状态下正弦波图像的左右平移位置

φ = arcsin(y/A) - wt

Yengi
关注
  • 0
    点赞
  • 踩
  • 3
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
dtmf音频解析
12-27
dtmf 基于matable可视化编程实现音频解析
dtmf（双音多频）音频分析软件
07-07
DTMF2NUM是一个小工具，可以帮助你从PCM解码波形文件和MF的DTMF音调。它支持任何文件类型的波（频率，频道和8，16，24和32位），自动优化（直流偏置调整和正常化）和两个WAV和原始的PCM数据。
DTMF双音频信令检测，基于Goertzel（戈泽尔算法），FPGA实现
Fairchild_1947的博客
05-10 6691
本文阐述了在DTMF双音频信令检测中使用的Goertzel（戈泽尔）算法的理论知识，并从公式计算的角度入手，一步步深入，从理论走向实现，再从实现走向切实可用。 本文在文末分享了本文提到的所有代码文件，感兴趣等小伙伴可以下载玩玩哦。
DTMF /FSK
让我们的手指敲起锋利的代码勇敢前进！
07-12 2287
DTMF     DTMF（Dual Tone Multi Frequency），双音多频，由高频群和低频群组成，高低频群各包含 4个频率。一个高频信号和一个低频信号叠加组成一个组合信号，代表一个数字。 DTMF信令有 16个编码。利用 DTMF信令可选择呼叫相应的对讲机双音多频信号（ DTMF），电话系统中电话机与交换机之间的一种用户信令，通常用于发送被叫号码。     在使用双音多频信号
DTMF -- 双音多频编码，即电话过程种的按键控制事件
weixin_34167819的博客
11-16 1166
d
DTMF.zip_DTMF c++_DTMF音频_dtmf_双音频
09-22
关于dtmf应用，DTMF双音频发生器，DTMF双音频接收器
FSK解码原理及实现方法.rar_FSK Simulink_dtmf_fsk c++_fsk解码_stft
09-24
FSK解码原理及实现方法
DTMF.rar_DTMF c语言_dtmf_dtmf 声卡_modem 音频_sunydc
09-24
C语言源代码，用计算机声卡发出电话双音频信号，可用于无modem的计算机电话拨号控制智能电器
DTMF.rar_FSK-DTMF解码_HT95_dtmf_dtmf解码_emc78
09-23
通过单片机（例如：HT95系列、EMC78系列、MCS51系列等）和外围电路的设计应用三级管的3个特性区，结合程序来实现来电显示的解码。...在设计应用中，免去了因采用FSK或DTMF专有解码芯片的高额成本。
DTMF.rar_DTMF TO_dtmf_necessaryyu5_信号能量检测_双音频信号
09-24
数组x为产生的双音频信号，长度102点，32位浮点型； 4、数组z为DTMF检测中能量模板 ，以完成按键的行位置及列位置检测； 6、变量outkeycol 及变量outkeyrow分别为按键行位置及列位置的检测结果，均为整 型变量；...
利用MATLAB仿真DTMF音频的产生及其解码.pdf DTMF音频的产生及其解码.zip
05-20
利用MATLAB仿真DTMF音频的产生及其解码.pdf DTMF音频的产生及其解码。
音频接口Android FSK通信
05-02
音频接口与Android实现通信，利用2fsk编码。可以设计在便携式医疗电子、可穿戴式设备。
中国来电显示双制式（FSK和DTMF）标准
12-24
中国来电显示双制式（FSK和DTMF）标准，最新权威，数据格式，利于参考与编程实现指导。
信号调制疑问_FSK调制中Frequency Deviation、DataRate及RxBW参数设置
weixin_42216454的博客
01-11 3233
　　很多朋友在使用FSK调制方式的射频通信模块时，经常会有一些疑问：射频配置参数中的FrequencyDev、DataRate及RxBW具体概念是什么意思？这几个参数之间有什么关系？在不同的应用需求中，该如何合理正确的设置它们？在本文中将对这些常见问题一一给出答案。　　1， 概念的理解　　FrequencyDeviation（频率偏差）　　在FSK信号调制中，信号‘0’和信号‘1’是用不同的频率来...
FSK数据格式
artmcu的专栏
10-30 2603
来电显示  FSK  有效的串行数据格式1 格式 字节数 说明 举例 消 息 消息类型 1 80H，解码方式 80H 消息长度 1 所有参数的数据长度 13H（十进制时为19） 参 数 参数1类型 1
音视频笔记-----三种数字调制的形式之FSK
qq_36187285的博客
02-09 7475
数字调制技术 数字调制是指用数字数据调制模拟信号，主要有三种形式：移频键控法FSK(Frequency-shift keying)、移幅键控法ASK()、移相键控法PSK()。计算机在处理信息时，只能处理数字信号，所以光信号，声信号，电信号等都需要转换成数字信号，即计算机能够识别的0，1二进制数字格式，但当数字信号传输到终端时，我们又需要对数据进行再次调制，使数字信号恢复成模拟信号，最显著的即...
【gnuradio 仿真音频数据ASK调制和解调】
qq_38529833的博客
03-12 1625
然后使用threshold模块进行判决，threshold模块的判决参数设置为0.1，这里设为0.1的原因是因为上图中，经过rms计算后，调制信号部分整体上移约0.2，这里设为0.1，把高于0.1的判决为1，低于0.1的判决为0，这样就可以检波出原始信号了。由于我们的音频数据是处理过的，因此要进行复原，kepp 1 in N的参数要设置为100，因为刚刚repeat中参数设置的为100，每个bit重复了100次。通过对比原始音频波形，和解调后波形，会发现是完全一致的，这是由于理想环境没有噪声。
话机按键音效——DTMF生成
qq_41281869的博客
11-18 2772
拨号音效——DTMF前言干货 前言 通常制作话机功能都需要用到0-9*#拨号音效，网上的资源少之又少，有的只是短暂时间的按键音效，文件大小大，音效时间无法自控，因博主最近也在寻找拨号音效，途中频繁碰壁不是要钱也是音效不对劲，最后找到了一个外国网址，可自定义制作DTMF音效，果然好资源都得找外网，话不多说献上干货。 干货 点击进入自定义DTMF网站 按键对应的双音频频率图标 下面示例按键1的DTMF音效怎么生成(写入两个频率，调整持续时间，采样率越高，占内存越大，其他随缘)。 对应内容 频率1：697 频
UniAnimate：华科提出人类跳舞视频生成新框架，支持合成一分钟高清视频
最新发布
2401_84033492的博客
06-09 832
其次，研究团队还提出了一种统一的噪声输入，其支持随机噪声输入和基于第一帧的条件噪声输入，随机噪声输入可以配合参考图像和姿态序列生成一段视频，而基于第一帧的条件噪声输入（First Frame Conditioning）则以视频第一帧作为条件输入延续生成后续的视频。随着视频生成技术的快速发展，特别是生成模型的迭代演化，跳舞视频生成任务取得了前所未有的进展，并展示了广泛的应用潜力。：UniAnimate 框架通过统一的视频扩散模型，消除了对额外参考网络的依赖，降低了训练难度和模型参数的数量。
multisim中DTMF_KEYPAD怎么用
05-28
在Multisim中，DTMF_KEYPAD模块用于产生DTMF信号，可以模拟电话键盘的输入。下面是DTMF_KEYPAD的使用步骤： 1. 打开Multisim软件，创建一个新的电路图。 2. 从左侧的“元件”库中选择“DTMF_KEYPAD”模块，将其拖动到电路图中。 3. 连接DTMF_KEYPAD模块的输入端和输出端，如下图所示： ![DTMF_KEYPAD](https://img-blog.csdnimg.cn/20211015104510321.png) 4. 右键单击DTMF_KEYPAD模块，选择“属性”，在弹出的对话框中可以设置DTMF信号的频率和振幅等参数。 5. 点击“运行”按钮或按下F5键，就可以开始模拟DTMF信号的产生了。 如果需要将DTMF信号输出到外部设备，可以使用Multisim中的其他模块来完成，比如声音放大器或者数字信号处理器等模块。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Yengi

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值