✅博主简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,Matlab项目合作可私信。
🍎个人主页:海神之光
🏆代码获取方式:
海神之光Matlab王者学习之路—代码获取方式
⛳️座右铭:行百里者,半于九十。
更多Matlab仿真内容点击👇
Matlab图像处理(进阶版)
路径规划(Matlab)
神经网络预测与分类(Matlab)
优化求解(Matlab)
语音处理(Matlab)
信号处理(Matlab)
车间调度(Matlab)
⛄一、傅立叶变换语音降噪简介
什么是傅里叶变换?
法国科学家傅里叶提出,任何一条周期曲线,无论多么跳跃或不规则,都能表示成一组光滑正弦曲线叠加之和。傅里叶变换即是把一条不规则的曲线拆解成一组光滑正弦曲线的过程。
傅里叶变换的目的是可将时域(即时间域)上的信号转变为频域(即频率域)上的信号,随着域的不同,对同一个事物的了解角度也就随之改变,因此在时域中某些不好处理的地方,在频域就可以较为简单的处理。这就可以大量减少处理信号存储量。
example:弹钢琴
假设有一时间域函数:y = f(x),根据傅里叶的理论它可以被分解为一系列正弦函数的叠加,他们的振幅A,频率ω或初相位φ不同:
所以傅里叶变换可以把一个比较复杂的函数转换为多个简单函数的叠加,看问题的角度也从时间域转到了频率域,有些的问题处理起来就会比较简单。
⛄二、部分源代码
function varargout = DSP(varargin)
% DSP MATLAB code for DSP.fig
% DSP, by itself, creates a new DSP or raises the existing
% singleton*.
%
% H = DSP returns the handle to a new DSP or the handle to
% the existing singleton*.
%
% DSP(‘CALLBACK’,hObject,eventData,handles,…) calls the local
% function named CALLBACK in DSP.M with the given input arguments.
%
% DSP(‘Property’,‘Value’,…) creates a new DSP or raises the
% existing singleton*. Starting from the left, property value pairs are
% applied to the GUI before DSP_OpeningFcn gets called. An
% unrecognized property name or invalid value makes property application
% stop. All inputs are passed to DSP_OpeningFcn via varargin.
%
% *See GUI Options on GUIDE’s Tools menu. Choose “GUI allows only one
% instance to run (singleton)”.
%
% See also: GUIDE, GUIDATA, GUIHANDLES
% Edit the above text to modify the response to help DSP
% Last Modified by GUIDE v2.5 13-Dec-2015 17:32:39
% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct(‘gui_Name’, mfilename, …
‘gui_Singleton’, gui_Singleton, …
‘gui_OpeningFcn’, @DSP_OpeningFcn, …
‘gui_OutputFcn’, @DSP_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
% End initialization code - DO NOT EDIT
% — Executes just before DSP is made visible.
function DSP_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% varargin command line arguments to DSP (see VARARGIN)
% Choose default command line output for DSP
handles.output = hObject;
% Update handles structure
guidata(hObject, handles);
% UIWAIT makes DSP wait for user response (see UIRESUME)
% uiwait(handles.main);
% — Outputs from this function are returned to the command line.
function varargout = DSP_OutputFcn(hObject, eventdata, handles)
% varargout cell array for returning output args (see VARARGOUT);
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Get default command line output from handles structure
varargout{1} = handles.output;
% --------------------------------------------------------------------
function file_Callback(hObject, eventdata, handles)
% hObject handle to file (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% --------------------------------------------------------------------
function noise_Callback(hObject, eventdata, handles)
% hObject handle to noise (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% --------------------------------------------------------------------
function filter_Callback(hObject, eventdata, handles)
% hObject handle to filter (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% --------------------------------------------------------------------
function help_Callback(hObject, eventdata, handles)
% hObject handle to help (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% --------------------------------------------------------------------
function help_doc_Callback(hObject, eventdata, handles)
% hObject handle to help_doc (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% --------------------------------------------------------------------
function about_Callback(hObject, eventdata, handles)
% hObject handle to about (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
string =str2mat(’ 声音处理示例程序’,[],’ 版本:1.0’,[],’ 开发:项晓强 褚倩云 徐如林’,[],’ 电子邮箱:xrlin@163.com’,[],’ 特别鸣谢:安大航模队大力赞助’);
msgbox(string);
% --------------------------------------------------------------------
function IIR_Callback(hObject, eventdata, handles)
% hObject handle to IIR (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% --------------------------------------------------------------------
function FIR_Callback(hObject, eventdata, handles)
% hObject handle to FIR (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% --------------------------------------------------------------------
function flat_noise_Callback(hObject, eventdata, handles)
% hObject handle to flat_noise (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
set(handles.L1,‘Visible’,‘off’);
set(handles.L2,‘Visible’,‘off’);
set(handles.L3,‘Visible’,‘off’);
set(handles.L4,‘Visible’,‘off’);
set(handles.E1,‘Visible’,‘off’);
set(handles.E2,‘Visible’,‘off’);
set(handles.E3,‘Visible’,‘off’);
set(handles.E4,‘Visible’,‘off’);
set(handles.single_run,‘Enable’,‘off’);
set(handles.multi_run,‘Enable’,‘off’);
set(handles.IIR_run,‘Enable’,‘off’);
set(handles.FIR_run,‘Enable’,‘off’);
wav = evalin(‘base’,‘data’);
fs = evalin(‘base’,‘Fs’);
reswav = awgn(wav,20);%添加信噪比为10dB的高斯噪声
assignin(‘base’,‘resdata’,reswav);
fftwav = abs(fft(reswav));
axes(handles.Time);
x = (0:length(reswav)-1)/fs;
handles.Line1 = plot(x,reswav);
guidata(hObject,handles);%保存值
set(handles.Time,‘XMinorTick’,‘on’);
grid on;
xlabel(‘时间/s’);
ylabel(‘幅度’);
title(‘时域图’);
axes(handles.Freq);
xf = (0:length(reswav)-1)'*fs/length(fftwav);
handles.Line2 = plot(xf,fftwav);
guidata(hObject,handles);%保存值
set(handles.Freq,‘XMinorTick’,‘on’);
grid on;
xlabel(‘频率/Hz’);
ylabel(‘幅度’);
title(‘频域图’);
assignin(‘base’,‘flag’,1);
handles.sou = audioplayer(reswav,fs);
guidata(hObject,handles);%保存值
% --------------------------------------------------------------------
function single_Callback(hObject, eventdata, handles)
% hObject handle to single (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% --------------------------------------------------------------------
function multi_Callback(hObject, eventdata, handles)
% hObject handle to multi (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% --------------------------------------------------------------------
function open_Callback(hObject, eventdata, handles)
% hObject handle to open (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
[filename,pathname] = uigetfile(‘.wav;.mp3’,‘请选择声音文件’, ‘MultiSelect’, ‘off’);
[wav,fs]=audioread(strcat(pathname,filename));
wav = wav(:,1);
handles.sou = audioplayer(wav,fs);
guidata(hObject,handles);%保存值
assignin(‘base’,‘data’,wav);
assignin(‘base’,‘Fs’,fs);
set(handles.L1,‘Visible’,‘off’);
set(handles.L2,‘Visible’,‘off’);
set(handles.L3,‘Visible’,‘off’);
set(handles.L4,‘Visible’,‘off’);
set(handles.E1,‘Visible’,‘off’);
set(handles.E2,‘Visible’,‘off’);
set(handles.E3,‘Visible’,‘off’);
set(handles.E4,‘Visible’,‘off’);
set(handles.button_play,‘Visible’,‘off’);
set(handles.Time,‘Visible’,‘off’);
set(handles.Freq,‘Visible’,‘off’);
set(handles.single_run,‘Enable’,‘off’);
set(handles.multi_run,‘Enable’,‘off’);
set(handles.IIR_run,‘Enable’,‘off’);
set(handles.FIR_run,‘Enable’,‘off’);
set(handles.button_play,‘Visible’,‘on’);
set(handles.Time,‘Visible’,‘on’);
set(handles.Freq,‘Visible’,‘on’);
set(handles.noise,‘Enable’,‘on’);
set(handles.filter,‘Enable’,‘on’);
set(handles.button_play,‘String’,‘播放’);
assignin(‘base’,‘ps’,0);
reswav = wav;
assignin(‘base’,‘resdata’,reswav);
fftwav = abs(fft(reswav));
axes(handles.Time);
x = (0:length(reswav)-1)/fs;
handles.Line1 = plot(x,reswav);
guidata(hObject,handles);%保存值
set(handles.Time,‘XMinorTick’,‘on’);
grid on;
xlabel(‘时间/s’);
ylabel(‘幅度’);
title(‘时域图’);
axes(handles.Freq);
xf = (0:length(reswav)-1)'*fs/length(fftwav);
handles.Line2 = plot(xf,fftwav);
guidata(hObject,handles);%保存值
set(handles.Freq,‘XMinorTick’,‘on’);
grid on;
xlabel(‘频率/Hz’);
ylabel(‘幅度’);
title(‘频域图’);
assignin(‘base’,‘flag’,1);
% --------------------------------------------------------------------
function close_Callback(hObject, eventdata, handles)
% hObject handle to close (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
clc;
clear all;
close;
% — Executes on button press in button_play.
function button_play_Callback(hObject, eventdata, handles)
% hObject handle to button_play (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
ps = evalin(‘base’,‘ps’);
if(ps == 0)
set(handles.button_play,‘String’,‘停止’);
assignin(‘base’,‘ps’,1);
play(handles.sou);
end;
if(ps == 1)
set(handles.button_play,‘String’,‘播放’);
assignin(‘base’,‘ps’,0);
stop(handles.sou);
end;
%sound(evalin(‘base’,‘resdata’),evalin(‘base’,‘Fs’));
% --------------------------------------------------------------------
function record_Callback(hObject, eventdata, handles)
% hObject handle to record (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
set(handles.L1,‘Visible’,‘off’);
set(handles.L2,‘Visible’,‘off’);
set(handles.L3,‘Visible’,‘off’);
set(handles.L4,‘Visible’,‘off’);
set(handles.E1,‘Visible’,‘off’);
set(handles.E2,‘Visible’,‘off’);
set(handles.E3,‘Visible’,‘off’);
set(handles.E4,‘Visible’,‘off’);
set(handles.button_play,‘Visible’,‘off’);
set(handles.Time,‘Visible’,‘off’);
set(handles.Freq,‘Visible’,‘off’);
set(handles.noise,‘Enable’,‘off’);
set(handles.filter,‘Enable’,‘off’);
set(handles.single_run,‘Enable’,‘off’);
set(handles.multi_run,‘Enable’,‘off’);
set(handles.IIR_run,‘Enable’,‘off’);
set(handles.FIR_run,‘Enable’,‘off’);
set(handles.button_play,‘String’,‘播放’);
assignin(‘base’,‘ps’,0);
%--------------------------------%
if (evalin(‘base’,‘flag’)==1)
set(handles.Line1,‘Visible’,‘off’);
set(handles.Line2,‘Visible’,‘off’);
end;
⛄三、运行结果
⛄四、matlab版本及参考文献
1 matlab版本
2014a
2 参考文献
[1]韩纪庆,张磊,郑铁然.语音信号处理(第3版)[M].清华大学出版社,2019.
[2]柳若边.深度学习:语音识别技术实践[M].清华大学出版社,2019.
3 备注
简介此部分摘自互联网,仅供参考,若侵权,联系删除
🍅 仿真咨询
1 各类智能优化算法改进及应用
生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化、公交排班优化、充电桩布局优化、车间布局优化、集装箱船配载优化、水泵组合优化、解医疗资源分配优化、设施布局优化、可视域基站和无人机选址优化
2 机器学习和深度学习方面
卷积神经网络(CNN)、LSTM、支持向量机(SVM)、最小二乘支持向量机(LSSVM)、极限学习机(ELM)、核极限学习机(KELM)、BP、RBF、宽度学习、DBN、RF、RBF、DELM、XGBOOST、TCN实现风电预测、光伏预测、电池寿命预测、辐射源识别、交通流预测、负荷预测、股价预测、PM2.5浓度预测、电池健康状态预测、水体光学参数反演、NLOS信号识别、地铁停车精准预测、变压器故障诊断
3 图像处理方面
图像识别、图像分割、图像检测、图像隐藏、图像配准、图像拼接、图像融合、图像增强、图像压缩感知
4 路径规划方面
旅行商问题(TSP)、车辆路径问题(VRP、MVRP、CVRP、VRPTW等)、无人机三维路径规划、无人机协同、无人机编队、机器人路径规划、栅格地图路径规划、多式联运运输问题、车辆协同无人机路径规划、天线线性阵列分布优化、车间布局优化
5 无人机应用方面
无人机路径规划、无人机控制、无人机编队、无人机协同、无人机任务分配
6 无线传感器定位及布局方面
传感器部署优化、通信协议优化、路由优化、目标定位优化、Dv-Hop定位优化、Leach协议优化、WSN覆盖优化、组播优化、RSSI定位优化
7 信号处理方面
信号识别、信号加密、信号去噪、信号增强、雷达信号处理、信号水印嵌入提取、肌电信号、脑电信号、信号配时优化
8 电力系统方面
微电网优化、无功优化、配电网重构、储能配置
9 元胞自动机方面
交通流 人群疏散 病毒扩散 晶体生长
10 雷达方面
卡尔曼滤波跟踪、航迹关联、航迹融合