如何向音频文件中插入噪声

为了研究噪声和音乐对EEG的的影响,实验前需要准备一段夹杂噪声的音乐。

  1. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  2. % Program to insert the white Gaussian noise into muisc
  3. %y 音乐数据
  4. %Fs 采样率
  5. %bits 采样位,默认为16位
  6. %wgn(m,n,p) 获取m * n、噪声强度为p的高斯白噪声
  7. % author: Andy Wu
  8. % date: 6/10/2009
  9. % modify: 9/10/2009
  10. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  11. %打开一个名为“请选择背景音乐”的对话框,显示指定路径下所有的.wav文件
  12. [file_name,file_path]=uigetfile('F:\MatlabProgram\*.wav','请选择背景音乐');
  13. sound_path=strcat(file_path,file_name);
  14. %将选中的音乐文件的数据和采样率读到sound_data和Fs
  15. [sound_data,Fs,bits]=wavread(sound_path);
  16. sound_length=length(sound_data);
  17. %求出音乐的播放时间长度
  18. sound_time=sound_length/Fs;
  19. fprintf(1,'背景音乐播放长度 %6.2f(秒)\n',sound_time);
  20. %检查是否定义噪声时间, 否则输入噪声时间
  21. if ~exist('noise_time')
  22. noise_time=input('请输入噪声长度(秒):');
  23. end
  24. noise_length=noise_time*Fs;
  25. %输出一个噪声强度为1dBw,noise_length x 1的高斯白噪声矩阵数据
  26. noise_data=wgn(noise_length,1,1);
  27. %检查是否输入插入噪声的起始时间
  28. if ~exist('noise_start_time')
  29. noise_start_time=input('请输入您要设置的噪声起始时间(秒):');
  30. end
  31. %计算噪声数据插入位置
  32. noise_position=noise_start_time*Fs;
  33. %插入噪声数据
  34. sound_data(noise_position:noise_position+noise_length-1)=noise_data;
  35. %dlmwrite('F:\notify.txt',y); %将插入噪声后的音乐数据写入文本文件
  36. %soundData=csvread('F:\notify.txt'); %读出音乐数据至soundData
  37. object= analogoutput('winsound'); %建立硬件对象
  38. addchannel(object,[1 2]); %创建声音输出通道
  39. set(object,'SampleRate',Fs) %设置采样率
  40. data1 =sound_data(:,1); %双声道
  41. data2 =sound_data(:,1);
  42. putdata(object,[data1 data2]); %往声卡输送音乐数据
  43. %start(object); %开始输出音乐
  44. %stop(object) %停止输出音乐

How to insert a segment of noise to music file的更多相关文章

  1. 源码浅析:MySQL一条insert操作,会写哪些文件?包括UNDO相关的文件吗?

    DML操作的大致流程 在解答上述疑惑之前,我们来梳理一下DML操作的大致流程: 1.语法解析.语义解析 2.生成执行计划 3.事务修改阶段 1) 激活事务,事务状态由not_active变为activ ...

  2. 【UISegmentedControl】-  分段控件

    一.初始化 二.常见的属性 1.segmentedControlStyle属性:设置基本的样式 2.momentary属性:设置在点击后是否恢复原样 . 3.numberOfSegments属性:只读 ...

  3. UI--普通控件总结1--控件使用

    本文目录 0.UIView常用的属性和操作 0_1.UIView常见的属性 0_2.UIView状态 0_3.UIView常用的方法 1.文本框UITextField和文本视图UITextView 1 ...

  4. 『GreenPlum系列』GreenPlum 4节点集群安装(图文教程)

      目标架构如上图   一.硬件评估 cpu主频,核数推荐CPU核数与磁盘数的比例在12:12以上Instance上执行时只能利用一个CPU核资源进行计算,推荐高主频 内存容量 网络带宽重分布操作 R ...

  5. 【mysql】关于IO/内存方面的一些优化

    这里使用的是mysql  Ver 14.14 Distrib 5.6.19, for Linux (i686) using  EditLine wrapper 一.mysql目录文件 ibdata1: ...

  6. TOYS - POJ 2318(计算几何,叉积判断)

    题目大意:给你一个矩形的左上角和右下角的坐标,然后这个矩形有 N 个隔板分割成 N+1 个区域,下面有 M 组坐标,求出来每个区域包含的坐标数.   分析:做的第一道计算几何题目....使用叉积判断方 ...

  7. lucene原理及源码解析--核心类

    马云说:大家还没搞清PC时代的时候,移动互联网来了,还没搞清移动互联网的时候,大数据时代来了. 然而,我看到的是:在PC时代搞PC的,移动互联网时代搞移动互联网的,大数据时代搞大数据的,都是同一伙儿人 ...

  8. InnoDB存储引擎介绍-(7) Innodb数据页结构

    数据页结构 File Header 总共38 Bytes,记录页的头信息 名称 大小(Bytes) 描述 FIL_PAGE_SPACE 4 该页的checksum值 FIL_PAGE_OFFSET 4 ...

  9. 学习笔记:The Best of MySQL Forum

    http://mysql.rjweb.org/bestof.html I have tagged many of the better forum threads. 'Better' is based ...

随机推荐

  1. hdu3507 Print Article(斜率优化入门)(pascal)

    Problem Description Zero has an old printer that doesn't work well sometimes. As it is antique, he s ...

  2. AtCoder Grand Contest 019 A: Ice Tea Store

    tourist出的题诶!想想就很高明,老年选手可能做不太动.不过A题还是按照惯例放水的. AtCoder Grand Contest 019 A: Ice Tea Store 题意:买0.25L,0. ...

  3. 洛谷 P1352 没有上司的舞会

    树形动规入门题 先放题面 题目描述 某大学有N个职员,编号为1~N.他们之间有从属关系,也就是说他们的关系就像一棵以校长为根的树,父结点就是子结点的直接上司.现在有个周年庆宴会,宴会每邀请来一个职员都 ...

  4. 深入理解JVM一性能监控工具

    一.前言 工欲善其事必先利其器,性能优化和故障排查在我们大都数人眼里是件比较棘手的事情,一是需要具备一定的原理知识作为基础,二是需要掌握排查问题和解决问题的流程.方法.本文就将介绍利用性能监控工具,帮 ...

  5. [HNOI2002]跳蚤 【容斥】

    题目描述 Z城市居住着很多只跳蚤.在Z城市周六生活频道有一个娱乐节目.一只跳蚤将被请上一个高空钢丝的正中央.钢丝很长,可以看作是无限长.节目主持人会给该跳蚤发一张卡片.卡片上写有N+1个自然数.其中最 ...

  6. acid(数据库事务正确执行的四个基本要素的缩写)

    ACID,指数据库事务正确执行的四个基本要素的缩写.包含:原子性(Atomicity).一致性(Consistency).隔离性(Isolation).持久性(Durability).一个支持事务(T ...

  7. LOJ分块⑨题汇总

    从零开始的分块学习系列(感谢hzwer) 题目顺序是我建议的做题顺序 先说一句:分块的核心思想(其实本身分块就可以说是一种思想)是:均摊(或者说平衡/权衡?)复杂度,同时这种思想本身不只局限于序列分块 ...

  8. 【CF962E】Byteland, Berland and Disputed Cities

    Portal! ---> 几句话题意 数轴上面有三种点(B点,R点,P点),现在要将其中的某些点连起来,满足将所有B点去掉之后,所有P点和R点都连通&将所有R点去掉之后,所有B点和P点都 ...

  9. 玲珑杯”ACM比赛 Round #19 B 维护单调栈

    1149 - Buildings Time Limit:2s Memory Limit:128MByte Submissions:588Solved:151 DESCRIPTION There are ...

  10. 04-树4. Root of AVL Tree-平衡查找树AVL树的实现

    对于一棵普通的二叉查找树而言,在进行多次的插入或删除后,容易让树失去平衡,导致树的深度不是O(logN),而接近O(N),这样将大大减少对树的查找效率.一种解决办法就是要有一个称为平衡的附加的结构条件 ...