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

为了研究噪声和音乐对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. bzoj2301-Problem b

    题意 \(T\le 5\times 10^4\) 次询问,每次询问 \(a,b,c,d,k\le 5\times 10^4\),求 \[ \sum _{i=a}^b\sum _{j=c}^d[gcd( ...

  2. canvas - 炫酷的3D星空

    1.国际惯例,先上效果 (⊙o⊙)… 效果图看上去效果并不很炫酷啊,直接戳 这里 看效果吧! 2代码部分 html: <canvas id="canvas" width=&q ...

  3. BZOJ4998 星球联盟(LCT+双连通分量+并查集)

    即要求动态维护边双.出现环时将路径上的点合并即可.LCT维护.具体地,加边成环时makeroot+access+splay一套把这段路径提出来,暴力dfs修改并查集祖先,并将这部分与根断开,视为删除这 ...

  4. P4169 [Violet]天使玩偶/SJY摆棋子

    题目背景 感谢@浮尘ii 提供的一组hack数据 题目描述 Ayu 在七年前曾经收到过一个天使玩偶,当时她把它当作时间囊埋在了地下.而七年后 的今天,Ayu 却忘了她把天使玩偶埋在了哪里,所以她决定仅 ...

  5. NOIP2016愤怒的小鸟 题解报告 【状压DP】

    题目什么大家都清楚 题解 我们知道,三点确定一条抛物线,现在这条抛物线过原点,所以任意两只猪确定一条抛物线.通过运算的出对于两头猪(x1,y1),(x2,y2),他们所在抛物线a=(y1*x2-y2* ...

  6. java规范(三)-----判空----方法内的为空判断

    一般我们判空或者有判断条件时 都是使用 if(条件成立){ 执行代码 } 这样的逻辑. 但是如果对象的字段很深层次时或者条件很多时就容易形成很多个{}的情况,这样就容易分不出哪个花括号属于哪部分.如下 ...

  7. 51nod 1766 树上的最远点对(线段树)

    像树的直径一样,两个集合的最长路也是由两个集合内部的最长路的两个端点组成的,于是我们知道了两个集合的最长路,枚举一下两两端点算出答案就可以合并了,所以就可以用线段树维护一个区间里的最长路了. #inc ...

  8. [知识点]C++中STL容器之vector

    零.STL目录 1.容器之map 2.容器之vector 3.容器之set 一.前言 关于STL和STL容器的概念参见STL系列第一篇——map(见上).今天介绍第二个成员——vector. 二.用途 ...

  9. vue添加属性绑定

    html <div id="app-2"> <span v-bind:title="message"> 鼠标悬停几秒钟查看此处动态绑定的 ...

  10. Java中的Class.forName

    在做JAVA EE开发的过程中,更多的是使用框架来提高开发效率.越来越发现,之前很基础的一些东西,都忘记的差不多了.从今天开始慢慢的复习一下基础.今天在看JDBC的时候,就有一个有趣的地方,之前学的时 ...