%去噪頻譜
close all;
clc;

PV=[235 180 664 463];

m_gui_win_h=figure('Position',PV,...
		'Resize','off',...
		'Numbertitle','off',...%'color',BACK_COLOR,...
		'Name','頻域參數(shù)分析');
    
set(m_gui_win_h,'MenuBar','none');

uicontrol('Style','frame',...
      'Position',[5 5 180 455],...
      'Callback',[]);

     m_gui_pb_time_h=uicontrol('Style','pushbutton',...
      'Position',[10 330 165 30],...
      'FontSize',10.0,...
      'String','載入語音',...
      'Callback','loadsp');

   
   m_gui_pb_vocalfold_h=uicontrol('Style','pushbutton',...
      'Position',[10 270 165 30],...
      'FontSize',10.0,...
      'String','原始信號頻譜',...
      'Callback','yuanshipinpu');

   
   m_gui_pb_formanttrk_h=uicontrol('Style','pushbutton',...
      'Position',[10 210 165 30],...
      'FontSize',10.0,...
      'String','去噪信號頻譜',...
      'Callback','quzaopinpu');
uicontrol('Style','pushbutton',...
      'Position',[10 150 165 30],...
      'FontSize',10.0,...
      'String','返回',...
      'Callback','main_can_shu_fen_xi');
ax1=axes('Position',[235 180 664 463]);
axes('Position',[.32 .5 .65 .35],'Box','on','Visible','on');
ffs=fftshift(fft(A));
plot(abs(ffs));
title('原始信號頻譜');

axes('Position',[.32 .05 .65 .35],'Box','on','Visible','on');
n=5;fc1=60;fc2=900;
[b,a]=butter(n,[fc1/(FS/2) fc2/(FS/2)]);
y=filter(b,a,A);
ffb=fftshift(fft(y));
plot(abs(ffb));
title('去噪信號頻譜');