Matlab信号处理-基于FastICA的声源分离系统

function varargout = ICA_fast(varargin) % ICA_FAST MATLAB code for ICA_fast.fig % ICA_FAST, by itself, creates a new ICA_FAST or raises the existing % singleton*. % % H = ICA_FAST returns the handle to a new ICA_FAST or the handle to % the existing singleton*. % % ICA_FAST('CALLBACK',hObject,eventData,handles,...) calls the local % function named CALLBACK in ICA_FAST.M with the given input arguments. % % ICA_FAST('Property','Value',...) creates a new ICA_FAST or raises the % existing singleton*. Starting from the left, property value pairs are % applied to the GUI before ICA_fast_OpeningFcn gets called. An % unrecognized property name or invalid value makes property application % stop. All inputs are passed to ICA_fast_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 ICA_fast % Last Modified by GUIDE v2.5 21-May-2021 08:57:10 % Begin initialization code - DO NOT EDIT gui_Singleton = 1; gui_State = struct('gui_Name', mfilename, ... 'gui_Singleton', gui_Singleton, ... 'gui_OpeningFcn', @ICA_fast_OpeningFcn, ... 'gui_OutputFcn', @ICA_fast_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 ICA_fast is made visible. function ICA_fast_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 ICA_fast (see VARARGIN) movegui(gcf,'center') warning off % Choose default command line output for ICA_fast handles.output = hObject; % Update handles structure guidata(hObject, handles); % UIWAIT makes ICA_fast wait for user response (see UIRESUME) % uiwait(handles.figure1); % --- Outputs from this function are returned to the command line. function varargout = ICA_fast_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; % --- Executes on button press in pushbutton2. function pushbutton2_Callback(hObject, eventdata, handles) % hObject handle to pushbutton2 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global MixedS FS FileName PathName [FileName,PathName] = uigetfile('*.wav'); pathall=strcat(PathName,FileName);%获得wav路径 [MixedS,FS]=audioread(pathall);%读取wav文件的数据和频率 t=(0:length(MixedS)-1)/FS;%变换得到播放时间 MixedS1=MixedS'; plot(handles.axes1,t,MixedS1(1,:)); title(handles.axes1,'混合信号1（声道1）'); plot(handles.axes2,t,MixedS1(2,:)); title(handles.axes2,'混合信号2（声道2）'); N=length(t); f=(0:N-1)*FS/N; MixedS_1=abs(fft(MixedS1(1,:))); MixedS_2=abs(fft(MixedS1(2,:))); plot(handles.axes3,f,MixedS_1); title(handles.axes3,'混合信号1频谱'); plot(handles.axes4,f,MixedS_2); title(handles.axes4,'混合信号2频谱'); % --- Executes on button press in pushbutton3. function pushbutton3_Callback(hObject, eventdata, handles) % hObject handle to pushbutton3 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global MixedS FS ICAedS MixedS=MixedS'; MixedS_bak=MixedS; %标准化 MixedS_mean=zeros(2,1); for i=1:2 MixedS_mean(i)=mean(MixedS(i,:)); end % 计算MixedS的均值 for i=1:2 for j=1:size(MixedS,1) MixedS(i,j)=MixedS(i,j)-MixedS_mean(i); end end %白化 MixedS_cov=cov(MixedS'); % cov为求协方差的函数 [E,D]=eig(MixedS_cov); % 对信号矩阵的协方差函数进行特征值分解 Q=inv(sqrt(D))*(E)'; % Q为白化矩阵 MixedS_white=Q*MixedS; % MixedS_white为白化后的信号矩阵 IsI=cov(MixedS_white'); % IsI应为单位阵 %FASTICA算法 X=MixedS_white; % 以下算法将对X进行操作 [VariableNum,SampleNum]=size(X); numofIC=VariableNum; % 在此应用中，独立元个数等于变量个数 B=zeros(numofIC,VariableNum); % 初始化列向量w的寄存矩阵,B=[b1 b2 ... bd] for r=1:numofIC i=1;maxIterationsNum=1000; % 设置最大迭代次数（即对于每个独立分量而言迭代均不超过此次数） IterationsNum=0; b=rand(numofIC,1)-.5; % 随机设置b初值 b=b/norm(b); % 对b标准化 norm(b):向量元素平方和开根号 while i<=maxIterationsNum+1 if i == maxIterationsNum % 循环结束处理 fprintf('\n第%d分量在%d次迭代内并不收敛。', r,maxIterationsNum); break; end bOld=b; a2=1; u=1; t=X'*b; g=t.*exp(-a2*t.^2/2); dg=(1-a2*t.^2).*exp(-a2*t.^2/2); b=((1-u)*t'*g*b+u*X*g)/SampleNum-mean(dg)*b; % 核心公式 b=b-B*B'*b; % 对b正交化 b=b/norm(b); if abs(abs(b'*bOld)-1)<1e-9 % 如果收敛，则 B(:,r)=b; % 保存所得向量b break; end i=i+1; end end for k=1:numofIC W(:,k)=B(:,k)/5^k; % 得到解混矩阵W end %ICA计算的数据复原并构图 ICAedS=W'*Q*MixedS_bak; % 计算ICA后的矩阵 % 将混合矩阵重新排列并输出 t=(0:length(ICAedS)-1)/FS;%变换得到播放时间 plot(handles.axes1,t,ICAedS(1,:)); title(handles.axes1,'分离信号1'); plot(handles.axes2,t,ICAedS(2,:)); title(handles.axes2,'分离信号2'); N=length(t); f=(0:N-1)*FS/N; ICAedS_1=abs(fft(ICAedS(1,:))); ICAedS_2=abs(fft(ICAedS(2,:))); plot(handles.axes3,f,ICAedS_1); title(handles.axes3,'分离信号1频谱'); plot(handles.axes4,f,ICAedS_2); title(handles.axes4,'分离信号2频谱'); % --- Executes on button press in pushbutton4. function pushbutton4_Callback(hObject, eventdata, handles) % hObject handle to pushbutton4 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global FileName PathName hunhe=fullfile(PathName,FileName); set(handles.activex1,'URL',hunhe); % --- Executes on button press in pushbutton5. function pushbutton5_Callback(hObject, eventdata, handles) % hObject handle to pushbutton5 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global ICAedS FS filename = 'ICA1.mp3'; audiowrite(filename, ICAedS(1,:), FS); set(handles.activex1,'URL',filename); % --- Executes on button press in pushbutton6. function pushbutton6_Callback(hObject, eventdata, handles) % hObject handle to pushbutton6 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global ICAedS FS filename = '