How to insert a segment of noise to music file

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

为了研究噪声和音乐对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. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  12. %打开一个名为“请选择背景音乐”的对话框，显示指定路径下所有的.wav文件
  13. [file_name,file_path]=uigetfile('F:\MatlabProgram\*.wav','请选择背景音乐');  
  15. sound_path=strcat(file_path,file_name);
  17. %将选中的音乐文件的数据和采样率读到sound_data和Fs 
  18. [sound_data,Fs,bits]=wavread(sound_path);    
  20. sound_length=length(sound_data);
  22. %求出音乐的播放时间长度
  23. sound_time=sound_length/Fs;       
  25. fprintf(1,'背景音乐播放长度    %6.2f(秒)\n',sound_time);
  27. %检查是否定义噪声时间， 否则输入噪声时间
  28. if ~exist('noise_time')                                   
  29.    noise_time=input('请输入噪声长度(秒)：');             
  30. end
  32. noise_length=noise_time*Fs;
  34. %输出一个噪声强度为1dBw，noise_length x 1的高斯白噪声矩阵数据
  35. noise_data=wgn(noise_length,1,1);             
  37. %检查是否输入插入噪声的起始时间 
  38. if ~exist('noise_start_time')                            
  39.    noise_start_time=input('请输入您要设置的噪声起始时间(秒)：');               
  40. end
  42. %计算噪声数据插入位置
  43. noise_position=noise_start_time*Fs;     
  45. %插入噪声数据                    
  46. sound_data(noise_position:noise_position+noise_length-1)=noise_data;   
  48. %dlmwrite('F:\notify.txt',y);                               %将插入噪声后的音乐数据写入文本文件
  49. %soundData=csvread('F:\notify.txt');                        %读出音乐数据至soundData
  50. object= analogoutput('winsound');                           %建立硬件对象
  51. addchannel(object,[1 2]);                                   %创建声音输出通道
  52. set(object,'SampleRate',Fs)                                 %设置采样率
  53. data1 =sound_data(:,1);                                     %双声道
  54. data2 =sound_data(:,1);
  55. putdata(object,[data1 data2]);                              %往声卡输送音乐数据
  56. %start(object);                                             %开始输出音乐
  57. %stop(object)                                               %停止输出音乐




posted @ 2015-08-23 16:33  AndyBlog  阅读(262)  评论(0)    收藏  举报