标题: 语音信号的时频认知(噪音为1/20的高斯白噪音)GUI界面和MATLAB代码 [打印本页]
作者: cyguet    时间: 2018-12-22 14:00
标题: 语音信号的时频认知(噪音为1/20的高斯白噪音)GUI界面和MATLAB代码
设计基于MATLAB，添加的噪音为1/20的高斯白噪音


MATLAB源程序如下:
  1. function varargout = untitled(varargin)
  2. gui_Singleton = 1;
  3. gui_State = struct('gui_Name',       mfilename, ...
  4.                    'gui_Singleton',  gui_Singleton, ...
  5.                    'gui_OpeningFcn', @untitled_OpeningFcn, ...
  6.                    'gui_OutputFcn',  @untitled_OutputFcn, ...
  7.                    'gui_LayoutFcn',  [] , ...
  8.                    'gui_Callback',   []);
  9. if nargin && ischar(varargin{1})
  10.     gui_State.gui_Callback = str2func(varargin{1});
  11. end

  13. if nargout
  14.     [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
  15. else
  16.     gui_mainfcn(gui_State, varargin{:});
  17. end
  18. function untitled_OpeningFcn(hObject, ~, handles, varargin)
  19. handles.output = hObject;
  20. guidata(hObject, handles);
  21. function varargout = untitled_OutputFcn(~, ~, handles)
  22. varargout{1} = handles.output;
  23. function pushbutton1_Callback(~, ~, handles)
  24. global bit SingleC Number_of_Wave_Data_Point SingleCn;
  25. [file,path]=uigetfile('*.m4a;*.mp3;*.wav;*.flac');
  26. set(handles.text3,'String',[path,file]);   % 将路径及文件名显示在文本框中
  27. [y,Fs] = audioread([path,file]);
  28. bit=Fs;
  29. SingleC=y(:,1); % 抽取第 1 声道  即抽取其中一维数组
  30. SingleCn=SingleC;
  31. Single_Channel_Max = max(abs(SingleC));  %幅值数据归一化
  32. Single_Channel_Wave = SingleC/Single_Channel_Max;   %幅值重构
  33. Number_of_Wave_Data_Point = length(Single_Channel_Wave);   
  34. Time_Single_Channel_Wave=(0:Number_of_Wave_Data_Point-1)/bit;  %计算数据点时刻 ， 需要将横轴的纯数字点转化为时间（秒）
  35. axes(handles.axes1);                                                                        
  36. plot(Time_Single_Channel_Wave,Single_Channel_Wave);                  %调取数据点，绘制图像
  37. xlabel('时间/s','FontSize',10,'FontWeight','bold');
  38. ylabel('幅值','FontSize',10,'FontWeight','bold');
  39. axes(handles.axes2);
  40. Single_Channel_Wave_Fft = fft(SingleC);                    
  41. Fft_Xabel=(0 : Number_of_Wave_Data_Point/2)*bit / Number_of_Wave_Data_Point;     %提取一般数据，以绘制单边谱波形
  42. Single_Channel_Wave_Fft_Max = max(abs(Single_Channel_Wave_Fft));  %幅值数据归一化
  43. Single_Channel_Wave_Fft = Single_Channel_Wave_Fft/Single_Channel_Wave_Fft_Max;   %幅值重构
  44. plot(Fft_Xabel,abs(Single_Channel_Wave_Fft (1 : Number_of_Wave_Data_Point / 2 +1)));%调用数据，绘制图像，以Fft_Xabel的数据为横轴数据，
  45. xlabel('频率/Hz','FontSize',10,'FontWeight','bold');
  46. ylabel('频率幅值','FontSize',10,'FontWeight','bold');
  47. function nc_Callback(~, ~, handles)
  48. global bit SingleC Number_of_Wave_Data_Point SingleCn;
  49. noise=wgn(Number_of_Wave_Data_Point,1,1);
  50. SingleCn=SingleC+0.05*noise;  %加入高斯白噪声
  51. Single_Channel_Max = max(abs(SingleCn));  %幅值数据归一化
  52. Single_Channel_Wave = SingleCn/Single_Channel_Max;   %幅值重构   
  53. Time_Single_Channel_Wave=(0:Number_of_Wave_Data_Point-1)/bit;  %计算数据点时刻 ， 需要将横轴的纯数字点转化为时间（秒）
  54. axes(handles.axes1);                                                                        
  55. plot(Time_Single_Channel_Wave,Single_Channel_Wave);                  %调取数据点，绘制图像
  56. xlabel('时间/s','FontSize',10,'FontWeight','bold');
  57. ylabel('幅值','FontSize',10,'FontWeight','bold');
  58. axes(handles.axes2);
  59. Single_Channel_Wave_Fft = fft(SingleCn);                    
  60. Fft_Xabel=(0 : Number_of_Wave_Data_Point/2)*bit / Number_of_Wave_Data_Point;     %提取一般数据，以绘制单边谱波形
  61. Single_Channel_Wave_Fft_Max = max(abs(Single_Channel_Wave_Fft));  %幅值数据归一化
  62. Single_Channel_Wave_Fft = Single_Channel_Wave_Fft/Single_Channel_Wave_Fft_Max;   %幅值重构
  63. plot(Fft_Xabel,abs(Single_Channel_Wave_Fft (1 : Number_of_Wave_Data_Point / 2 +1)));%调用数据，绘制图像，以Fft_Xabel的数据为横轴数据，
  64. xlabel('频率/Hz','FontSize',10,'FontWeight','bold');
  65. ylabel('频率幅值','FontSize',10,'FontWeight','bold');
  66. set(handles.uipanel2,'Title','加噪后的时域波形和频谱图');
  67. Time_Noise=(0:Number_of_Wave_Data_Point-1)/bit;  %计算数据点时刻
  68. Noise_Max = max(abs(SingleCn));  %幅值数据归一化
  69. SingleCn = SingleCn/Noise_Max;   %幅值重构
  70. axes(handles.axes3);
  71. plot(Time_Noise,0.05*noise,'r');
  72. xlabel('时间/s');
  73. ylabel('幅值');
  74. function stp_Callback(~, ~, ~)
  75. clear sound;
  76. function pushbutton11_Callback(~, ~, handles)
  77. global bit SingleC Number_of_Wave_Data_Point SingleCn;
  78. SingleCn=SingleC;
  79. Single_Channel_Max = max(abs(SingleCn));  %幅值数据归一化
  80. Single_Channel_Wave = SingleCn/Single_Channel_Max;   %幅值重构   
  81. Time_Single_Channel_Wave=(0:Number_of_Wave_Data_Point-1)/bit;  %计算数据点时刻 ， 需要将横轴的纯数字点转化为时间（秒）
  82. axes(handles.axes1);                                                                        
  83. plot(Time_Single_Channel_Wave,Single_Channel_Wave);                  %调取数据点，绘制图像
  84. xlabel('时间/s','FontSize',10,'FontWeight','bold');
  85. ylabel('幅值','FontSize',10,'FontWeight','bold');
  86. axes(handles.axes2);
  87. Single_Channel_Wave_Fft = fft(SingleCn);                    
  88. Fft_Xabel=(0 : Number_of_Wave_Data_Point/2)*bit / Number_of_Wave_Data_Point;     %提取一般数据，以绘制单边谱波形
  89. Single_Channel_Wave_Fft_Max = max(abs(Single_Channel_Wave_Fft));  %幅值数据归一化
  90. Single_Channel_Wave_Fft = Single_Channel_Wave_Fft/Single_Channel_Wave_Fft_Max;   %幅值重构
  91. plot(Fft_Xabel,abs(Single_Channel_Wave_Fft (1 : Number_of_Wave_Data_Point / 2 +1)));%调用数据，绘制图像，以Fft_Xabel的数据为横轴数据，
  92. xlabel('频率/Hz','FontSize',10,'FontWeight','bold');
  93. ylabel('频率幅值','FontSize',10,'FontWeight','bold');
  94. set(handles.uipanel2,'Title','原始信号的时域波形和频谱图');
  95. axes(handles.axes3);
  96. plot(0);
  97. xlabel('时间/s');
  98. ylabel('幅值');
  99. function pc_Callback(~, ~, handles)
  100. global SingleCn bit
  101. sound(SingleCn,0.5*get(handles.spd,'value')*bit);
复制代码

所有资料51hei提供下载:
信号大作业.rar (62.34 KB, 下载次数: 19)






欢迎光临 (http://www.51hei.com/bbs/) Powered by Discuz! X3.1