语音信号的时频认知(噪音为1/20的高斯白噪音)GUI界面和MATLAB代码

ID:452077 · 发表于 2018-12-22 14:00

设计基于MATLAB，添加的噪音为1/20的高斯白噪音

MATLAB源程序如下:

function varargout = untitled(varargin)
gui_Singleton = 1;
gui_State = struct('gui_Name', mfilename, ...
'gui_Singleton', gui_Singleton, ...
'gui_OpeningFcn', @untitled_OpeningFcn, ...
'gui_OutputFcn', @untitled_OutputFcn, ...
'gui_LayoutFcn', [] , ...
'gui_Callback', []);
if nargin && ischar(varargin{1})
gui_State.gui_Callback = str2func(varargin{1});
end
if nargout
[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
gui_mainfcn(gui_State, varargin{:});
end
function untitled_OpeningFcn(hObject, ~, handles, varargin)
handles.output = hObject;
guidata(hObject, handles);
function varargout = untitled_OutputFcn(~, ~, handles)
varargout{1} = handles.output;
function pushbutton1_Callback(~, ~, handles)
global bit SingleC Number_of_Wave_Data_Point SingleCn;
[file,path]=uigetfile('*.m4a;*.mp3;*.wav;*.flac');
set(handles.text3,'String',[path,file]); % 将路径及文件名显示在文本框中
[y,Fs] = audioread([path,file]);
bit=Fs;
SingleC=y(:,1); % 抽取第 1 声道即抽取其中一维数组
SingleCn=SingleC;
Single_Channel_Max = max(abs(SingleC)); %幅值数据归一化
Single_Channel_Wave = SingleC/Single_Channel_Max; %幅值重构
Number_of_Wave_Data_Point = length(Single_Channel_Wave);
Time_Single_Channel_Wave=(0:Number_of_Wave_Data_Point-1)/bit; %计算数据点时刻，需要将横轴的纯数字点转化为时间（秒）
axes(handles.axes1);
plot(Time_Single_Channel_Wave,Single_Channel_Wave); %调取数据点，绘制图像
xlabel('时间/s','FontSize',10,'FontWeight','bold');
ylabel('幅值','FontSize',10,'FontWeight','bold');
axes(handles.axes2);
Single_Channel_Wave_Fft = fft(SingleC);
Fft_Xabel=(0 : Number_of_Wave_Data_Point/2)*bit / Number_of_Wave_Data_Point; %提取一般数据，以绘制单边谱波形
Single_Channel_Wave_Fft_Max = max(abs(Single_Channel_Wave_Fft)); %幅值数据归一化
Single_Channel_Wave_Fft = Single_Channel_Wave_Fft/Single_Channel_Wave_Fft_Max; %幅值重构
plot(Fft_Xabel,abs(Single_Channel_Wave_Fft (1 : Number_of_Wave_Data_Point / 2 +1)));%调用数据，绘制图像，以Fft_Xabel的数据为横轴数据，
xlabel('频率/Hz','FontSize',10,'FontWeight','bold');
ylabel('频率幅值','FontSize',10,'FontWeight','bold');
function nc_Callback(~, ~, handles)
global bit SingleC Number_of_Wave_Data_Point SingleCn;
noise=wgn(Number_of_Wave_Data_Point,1,1);
SingleCn=SingleC+0.05*noise; %加入高斯白噪声
Single_Channel_Max = max(abs(SingleCn)); %幅值数据归一化
Single_Channel_Wave = SingleCn/Single_Channel_Max; %幅值重构
Time_Single_Channel_Wave=(0:Number_of_Wave_Data_Point-1)/bit; %计算数据点时刻，需要将横轴的纯数字点转化为时间（秒）
axes(handles.axes1);
plot(Time_Single_Channel_Wave,Single_Channel_Wave); %调取数据点，绘制图像
xlabel('时间/s','FontSize',10,'FontWeight','bold');
ylabel('幅值','FontSize',10,'FontWeight','bold');
axes(handles.axes2);
Single_Channel_Wave_Fft = fft(SingleCn);
Fft_Xabel=(0 : Number_of_Wave_Data_Point/2)*bit / Number_of_Wave_Data_Point; %提取一般数据，以绘制单边谱波形
Single_Channel_Wave_Fft_Max = max(abs(Single_Channel_Wave_Fft)); %幅值数据归一化
Single_Channel_Wave_Fft = Single_Channel_Wave_Fft/Single_Channel_Wave_Fft_Max; %幅值重构
plot(Fft_Xabel,abs(Single_Channel_Wave_Fft (1 : Number_of_Wave_Data_Point / 2 +1)));%调用数据，绘制图像，以Fft_Xabel的数据为横轴数据，
xlabel('频率/Hz','FontSize',10,'FontWeight','bold');
ylabel('频率幅值','FontSize',10,'FontWeight','bold');
set(handles.uipanel2,'Title','加噪后的时域波形和频谱图');
Time_Noise=(0:Number_of_Wave_Data_Point-1)/bit; %计算数据点时刻
Noise_Max = max(abs(SingleCn)); %幅值数据归一化
SingleCn = SingleCn/Noise_Max; %幅值重构
axes(handles.axes3);
plot(Time_Noise,0.05*noise,'r');
xlabel('时间/s');
ylabel('幅值');
function stp_Callback(~, ~, ~)
clear sound;
function pushbutton11_Callback(~, ~, handles)
global bit SingleC Number_of_Wave_Data_Point SingleCn;
SingleCn=SingleC;
Single_Channel_Max = max(abs(SingleCn)); %幅值数据归一化
Single_Channel_Wave = SingleCn/Single_Channel_Max; %幅值重构
Time_Single_Channel_Wave=(0:Number_of_Wave_Data_Point-1)/bit; %计算数据点时刻，需要将横轴的纯数字点转化为时间（秒）
axes(handles.axes1);
plot(Time_Single_Channel_Wave,Single_Channel_Wave); %调取数据点，绘制图像
xlabel('时间/s','FontSize',10,'FontWeight','bold');
ylabel('幅值','FontSize',10,'FontWeight','bold');
axes(handles.axes2);
Single_Channel_Wave_Fft = fft(SingleCn);
Fft_Xabel=(0 : Number_of_Wave_Data_Point/2)*bit / Number_of_Wave_Data_Point; %提取一般数据，以绘制单边谱波形
Single_Channel_Wave_Fft_Max = max(abs(Single_Channel_Wave_Fft)); %幅值数据归一化
Single_Channel_Wave_Fft = Single_Channel_Wave_Fft/Single_Channel_Wave_Fft_Max; %幅值重构
plot(Fft_Xabel,abs(Single_Channel_Wave_Fft (1 : Number_of_Wave_Data_Point / 2 +1)));%调用数据，绘制图像，以Fft_Xabel的数据为横轴数据，
xlabel('频率/Hz','FontSize',10,'FontWeight','bold');
ylabel('频率幅值','FontSize',10,'FontWeight','bold');
set(handles.uipanel2,'Title','原始信号的时域波形和频谱图');
axes(handles.axes3);
plot(0);
xlabel('时间/s');
ylabel('幅值');
function pc_Callback(~, ~, handles)
global SingleCn bit
sound(SingleCn,0.5*get(handles.spd,'value')*bit);

复制代码

所有资料51hei提供下载:

信号大作业.rar (62.34 KB, 下载次数: 19)

帐号		自动登录	找回密码
密码			立即注册

语音信号的时频认知(噪音为1/20的高斯白噪音)GUI界面和MATLAB代码

评分