找回密碼
 立即注冊(cè)

QQ登錄

只需一步，快速開始

搜索
查看: 2960|回復(fù): 0
收起左側(cè)

語(yǔ)音信號(hào)的時(shí)頻認(rèn)知(噪音為1/20的高斯白噪音)GUI界面和MATLAB代碼

[復(fù)制鏈接]
ID:452077 發(fā)表于 2018-12-22 14:00 | 顯示全部樓層 |閱讀模式
設(shè)計(jì)基于MATLAB，添加的噪音為1/20的高斯白噪音
0.png

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 聲道  即抽取其中一維數(shù)組
  30. SingleCn=SingleC;
  31. Single_Channel_Max = max(abs(SingleC));  %幅值數(shù)據(jù)歸一化
  32. Single_Channel_Wave = SingleC/Single_Channel_Max;   %幅值重構(gòu)
  33. Number_of_Wave_Data_Point = length(Single_Channel_Wave);   
  34. Time_Single_Channel_Wave=(0:Number_of_Wave_Data_Point-1)/bit;  %計(jì)算數(shù)據(jù)點(diǎn)時(shí)刻 ， 需要將橫軸的純數(shù)字點(diǎn)轉(zhuǎn)化為時(shí)間（秒）
  35. axes(handles.axes1);                                                                        
  36. plot(Time_Single_Channel_Wave,Single_Channel_Wave);                  %調(diào)取數(shù)據(jù)點(diǎn)，繪制圖像
  37. xlabel('時(shí)間/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;     %提取一般數(shù)據(jù)，以繪制單邊譜波形
  42. Single_Channel_Wave_Fft_Max = max(abs(Single_Channel_Wave_Fft));  %幅值數(shù)據(jù)歸一化
  43. Single_Channel_Wave_Fft = Single_Channel_Wave_Fft/Single_Channel_Wave_Fft_Max;   %幅值重構(gòu)
  44. plot(Fft_Xabel,abs(Single_Channel_Wave_Fft (1 : Number_of_Wave_Data_Point / 2 +1)));%調(diào)用數(shù)據(jù)，繪制圖像，以Fft_Xabel的數(shù)據(jù)為橫軸數(shù)據(jù)，
  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));  %幅值數(shù)據(jù)歸一化
  52. Single_Channel_Wave = SingleCn/Single_Channel_Max;   %幅值重構(gòu)   
  53. Time_Single_Channel_Wave=(0:Number_of_Wave_Data_Point-1)/bit;  %計(jì)算數(shù)據(jù)點(diǎn)時(shí)刻 ， 需要將橫軸的純數(shù)字點(diǎn)轉(zhuǎn)化為時(shí)間（秒）
  54. axes(handles.axes1);                                                                        
  55. plot(Time_Single_Channel_Wave,Single_Channel_Wave);                  %調(diào)取數(shù)據(jù)點(diǎn)，繪制圖像
  56. xlabel('時(shí)間/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;     %提取一般數(shù)據(jù)，以繪制單邊譜波形
  61. Single_Channel_Wave_Fft_Max = max(abs(Single_Channel_Wave_Fft));  %幅值數(shù)據(jù)歸一化
  62. Single_Channel_Wave_Fft = Single_Channel_Wave_Fft/Single_Channel_Wave_Fft_Max;   %幅值重構(gòu)
  63. plot(Fft_Xabel,abs(Single_Channel_Wave_Fft (1 : Number_of_Wave_Data_Point / 2 +1)));%調(diào)用數(shù)據(jù)，繪制圖像，以Fft_Xabel的數(shù)據(jù)為橫軸數(shù)據(jù)，
  64. xlabel('頻率/Hz','FontSize',10,'FontWeight','bold');
  65. ylabel('頻率幅值','FontSize',10,'FontWeight','bold');
  66. set(handles.uipanel2,'Title','加噪后的時(shí)域波形和頻譜圖');
  67. Time_Noise=(0:Number_of_Wave_Data_Point-1)/bit;  %計(jì)算數(shù)據(jù)點(diǎn)時(shí)刻
  68. Noise_Max = max(abs(SingleCn));  %幅值數(shù)據(jù)歸一化
  69. SingleCn = SingleCn/Noise_Max;   %幅值重構(gòu)
  70. axes(handles.axes3);
  71. plot(Time_Noise,0.05*noise,'r');
  72. xlabel('時(shí)間/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));  %幅值數(shù)據(jù)歸一化
  80. Single_Channel_Wave = SingleCn/Single_Channel_Max;   %幅值重構(gòu)   
  81. Time_Single_Channel_Wave=(0:Number_of_Wave_Data_Point-1)/bit;  %計(jì)算數(shù)據(jù)點(diǎn)時(shí)刻 ， 需要將橫軸的純數(shù)字點(diǎn)轉(zhuǎn)化為時(shí)間（秒）
  82. axes(handles.axes1);                                                                        
  83. plot(Time_Single_Channel_Wave,Single_Channel_Wave);                  %調(diào)取數(shù)據(jù)點(diǎn)，繪制圖像
  84. xlabel('時(shí)間/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;     %提取一般數(shù)據(jù)，以繪制單邊譜波形
  89. Single_Channel_Wave_Fft_Max = max(abs(Single_Channel_Wave_Fft));  %幅值數(shù)據(jù)歸一化
  90. Single_Channel_Wave_Fft = Single_Channel_Wave_Fft/Single_Channel_Wave_Fft_Max;   %幅值重構(gòu)
  91. plot(Fft_Xabel,abs(Single_Channel_Wave_Fft (1 : Number_of_Wave_Data_Point / 2 +1)));%調(diào)用數(shù)據(jù)，繪制圖像，以Fft_Xabel的數(shù)據(jù)為橫軸數(shù)據(jù)，
  92. xlabel('頻率/Hz','FontSize',10,'FontWeight','bold');
  93. ylabel('頻率幅值','FontSize',10,'FontWeight','bold');
  94. set(handles.uipanel2,'Title','原始信號(hào)的時(shí)域波形和頻譜圖');
  95. axes(handles.axes3);
  96. plot(0);
  97. xlabel('時(shí)間/s');
  98. ylabel('幅值');
  99. function pc_Callback(~, ~, handles)
  100. global SingleCn bit
  101. sound(SingleCn,0.5*get(handles.spd,'value')*bit);
復(fù)制代碼

所有資料51hei提供下載:
信號(hào)大作業(yè).rar (62.34 KB, 下載次數(shù): 19)


評(píng)分

參與人數(shù) 1黑幣 +50 收起 理由
admin + 50 共享資料的黑幣獎(jiǎng)勵(lì)！

查看全部評(píng)分

回復(fù)

使用道具 舉報(bào)

本版積分規(guī)則

手機(jī)版|小黑屋|51黑電子論壇 |51黑電子論壇6群 QQ 管理員QQ:125739409;技術(shù)交流QQ群281945664

Powered by 單片機(jī)教程網(wǎng)

快速回復(fù) 返回頂部 返回列表