dplresolution.m

在MATLAB环境下设计一个加速度计动态校准的数据处理软件附有程序源代码,压缩包中有软件使用说明

  set(findobj('Tag','text41'),'String',[num2str(Horizon_Time),'s']);
  set(findobj('Tag','text42'),'String',CHSTD_Status);
  set(findobj('Tag','edit24'),'String',num2str(countSTD));
  set(findobj('Tag','edit30'),'String',num2str(Filter_Freq));
  set(findobj('Tag','edit17'),'String',num2str(Filter_Freq));
  set(findobj('Tag','edit18'),'String',num2str(Filter_Freq));
  x=1:countSTD;
  plot(x,DPR_DataSTD,'r')
  set(gca,'Color',[0,0,0],'XColor',[0,1/2,0],'YColor',[0,1/2,0]);
  grid on;zoom off;zoom xon;
  hold on;
  load DPR_DataSupport DPRCover_PopMenuTEA_Value DPRCover_PopMenuTEB_Value CHTEA_Status CHTEB_Status;
  if DPRCover_PopMenuTEA_Value==4 & DPRCover_PopMenuTEB_Value~=4
      if any(CHTEB_Status)
          set(findobj('Tag','pushbutton9'),'Enable','on');
      end
  elseif DPRCover_PopMenuTEA_Value~=4 & DPRCover_PopMenuTEB_Value==4
      if any(CHTEA_Status)
          set(findobj('Tag','pushbutton9'),'Enable','on');
      end
  elseif DPRCover_PopMenuTEA_Value~=4 & DPRCover_PopMenuTEB_Value~=4
      if any(CHTEA_Status) & any(CHTEB_Status)
          set(findobj('Tag','pushbutton9'),'Enable','on');
      end
  elseif DPRCover_PopMenuTEA_Value==4 & DPRCover_PopMenuTEB_Value==4
      set(findobj('Tag','pushbutton9'),'Enable','on');    
  end
  handles.countSTD=countSTD;
  handles.DPR_DataSTD_Backup=DPR_DataSTD;
  handles.DPR_DataSTD=DPR_DataSTD;
  handles.Horizon_Interval=Horizon_Interval;
  handles.Vertical_Factor_STD=Vertical_Factor_STD;
  handles.output = hObject;
  guidata(hObject, handles);
end
  
% --- Executes on button press in pushbutton17.
function pushbutton17_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton17 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
load DPR_DataSupport DPRCover_PopMenuCOL_Value
if DPRCover_PopMenuCOL_Value==1 
  load DPR_DataSupport path;
  name='*.wfm';
  filename=[path,name];
  [FileTEA_Name,path]=uigetfile(filename,'打开被校器具数据文件');        
  filename=[path,FileTEA_Name];                                       
  status=copyfile(filename,'FileTEA.wfm');
  !cnvrtwfm -r -8 FileTEA.wfm
  [fid,message]=fopen('FileTEA.wav','rb');
  [DPR_DataTEA,countTEA]=fread(fid,inf,'int8');
  status=fclose(fid);
  save DPR_DataSupport FileTEA_Name path countTEA -append;
  set(hObject,'String',FileTEA_Name);
  !cnvrtwfm -p -d- -b- FileTEA.wfm
  [fid,message]=fopen('FileTEA.asc','rt');
  Inf=fread(fid,inf,'char')';
  InfStr=char(Inf);
  FidPos=findstr(InfStr,'YMULT');
  status=fseek(fid,FidPos+4,-1);
  FlagStr=fread(fid,7,'char')';
  Vertical_Factor_TEA=sscanf(char(FlagStr),'%f');
  FidPos=findstr(InfStr,'WFID');
  status=fseek(fid,FidPos+10,-1);
  FlagStr=fread(fid,26,'char')';
  CHTEA_Status=sscanf(char(FlagStr),'%c');
  save DPR_DataSupport CHTEA_Status -append;
  status=fclose(fid);
  set(findobj('Tag','text44'),'String',CHTEA_Status);
  x=1:countTEA;
  plot(x,DPR_DataTEA,'c')
  set(gca,'Color',[0,0,0],'XColor',[0,1/2,0],'YColor',[0,1/2,0]);
  grid on;zoom off;zoom xon;
  hold on;
  load DPR_DataSupport DPRCover_PopMenuTEA_Value DPRCover_PopMenuTEB_Value CHSTD_Status CHTEB_Status;
  if DPRCover_PopMenuTEB_Value==4
      if any(CHSTD_Status)
          set(findobj('Tag','pushbutton9'),'Enable','on');
      end
  elseif DPRCover_PopMenuTEB_Value~=4
      if any(CHTEB_Status) & any(CHSTD_Status)
          set(findobj('Tag','pushbutton9'),'Enable','on');
      end
  end
  handles.DPR_DataTEA_Backup=DPR_DataTEA;
  handles.DPR_DataTEA=DPR_DataTEA;
  handles.Vertical_Factor_TEA=Vertical_Factor_TEA;
  handles.output = hObject;
  guidata(hObject, handles);
elseif DPRCover_PopMenuCOL_Value==2 
  load DPR_DataSupport path;
  name='*.wfm';
  filename=[path,name];
  [FileTEA_Name,path]=uigetfile(filename,'打开被校器具数据文件');        
  filename=[path,FileTEA_Name];                                       
  status=copyfile(filename,'FileTEA.wfm');
  Vertical_Factor_TEA=1;
  CHTEA_Status='NA';
  save DPR_DataSupport CHTEA_Status -append;
  !ConvertTekWfm FileTEA.wfm /MLAB
  [fid,message]=fopen('FileTEA.dat','r');
  [DPR_DataTEA,countTEA]=fscanf(fid,'%f',inf);
  status=fclose(fid);
  save DPR_DataSupport FileTEA_Name path countTEA -append;
  set(hObject,'String',FileTEA_Name);
  x=1:countTEA;
  plot(x,DPR_DataTEA,'c')
  set(gca,'Color',[0,0,0],'XColor',[0,1/2,0],'YColor',[0,1/2,0]);
  grid on;zoom off;zoom xon;
  hold on;
  load DPR_DataSupport DPRCover_PopMenuTEA_Value DPRCover_PopMenuTEB_Value CHSTD_Status CHTEB_Status;
  if DPRCover_PopMenuTEB_Value==4
      if any(CHSTD_Status)
          set(findobj('Tag','pushbutton9'),'Enable','on');
      end
  elseif DPRCover_PopMenuTEB_Value~=4
      if any(CHTEB_Status) & any(CHSTD_Status)
          set(findobj('Tag','pushbutton9'),'Enable','on');
      end
  end
  handles.DPR_DataTEA_Backup=DPR_DataTEA;
  handles.DPR_DataTEA=DPR_DataTEA;
  handles.Vertical_Factor_TEA=Vertical_Factor_TEA;
  handles.output = hObject;
  guidata(hObject, handles);
elseif DPRCover_PopMenuCOL_Value==3 
  load DPR_DataSupport path;
  name='*.txt';
  filename=[path,name];
  [FileTEA_Name,path]=uigetfile(filename,'打开被校器具数据文件');        
  filename=[path,FileTEA_Name];                                       
  [fid,message]=fopen(filename,'rt');
  flag1=fscanf(fid,'%c',625);
  flagTX1=findstr(flag1,'XInc:');
  status=fseek(fid,flagTX1+20,'bof');
  Horizon_Interval=fscanf(fid,'%c',24);
  Horizon_Interval=sscanf(char(Horizon_Interval),'%f');
  flagTX2=findstr(flag1,'Data:')
  status=fseek(fid,flagTX2+25,'bof');
  [DPR_DataTEA,countTEA]=fscanf(fid,'%f',inf);
  status=fclose(fid);
  save DPR_DataSupport FileTEA_Name path countTEA -append;
  set(hObject,'String',FileTEA_Name);
  Vertical_Factor_TEA=1;
  CHTEA_Status='NA';
  save DPR_DataSupport CHTEA_Status -append;
  x=1:countTEA;
  plot(x,DPR_DataTEA,'c')
  set(gca,'Color',[0,0,0],'XColor',[0,1/2,0],'YColor',[0,1/2,0]);
  grid on;zoom off;zoom xon;
  hold on;
  load DPR_DataSupport DPRCover_PopMenuTEA_Value DPRCover_PopMenuTEB_Value CHSTD_Status CHTEB_Status;
  if DPRCover_PopMenuTEB_Value==4
      if any(CHSTD_Status)
          set(findobj('Tag','pushbutton9'),'Enable','on');
      end
  elseif DPRCover_PopMenuTEB_Value~=4
      if any(CHTEB_Status) & any(CHSTD_Status)
          set(findobj('Tag','pushbutton9'),'Enable','on');
      end
  end
  handles.DPR_DataTEA_Backup=DPR_DataTEA;
  handles.DPR_DataTEA=DPR_DataTEA;
  handles.Vertical_Factor_TEA=Vertical_Factor_TEA;
  handles.output = hObject;
  guidata(hObject, handles);
end
  
% --- Executes on button press in pushbutton18.
function pushbutton18_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton18 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
load DPR_DataSupport DPRCover_PopMenuCOL_Value
if DPRCover_PopMenuCOL_Value==1 
  load DPR_DataSupport path;
  name='*.wfm';
  filename=[path,name];
  [FileTEB_Name,path]=uigetfile(filename,'打开标准器具数据文件');        
  filename=[path,FileTEB_Name];                                       
  status=copyfile(filename,'FileTEB.wfm');
  !cnvrtwfm -r -8 FileTEB.wfm
  [fid,message]=fopen('FileTEB.wav','rb');
  [DPR_DataTEB,countTEB]=fread(fid,inf,'int8');
  status=fclose(fid);
  save DPR_DataSupport FileTEB_Name path countTEB -append;
  if countTEB ~=handles.countSTD
    errordlg('数据长度错误')
    return
  else
    set(hObject,'String',FileTEB_Name);
  end
  !cnvrtwfm -p -d- -b- FileTEB.wfm
  [fid,message]=fopen('FileTEB.asc','rt');
  Inf=fread(fid,inf,'char')';
  InfStr=char(Inf);
  FidPos=findstr(InfStr,'YMULT');
  status=fseek(fid,FidPos+4,-1);
  FlagStr=fread(fid,7,'char')';
  Vertical_Factor_TEB=sscanf(char(FlagStr),'%f');
  FidPos=findstr(InfStr,'WFID');
  status=fseek(fid,FidPos+10,-1);
  FlagStr=fread(fid,26,'char')';
  CHTEB_Status=sscanf(char(FlagStr),'%c');
  save DPR_DataSupport CHTEB_Status -append;
  status=fclose(fid);
  set(findobj('Tag','text44'),'String',CHTEB_Status);
  x=1:countTEB;
  plot(x,DPR_DataTEB,'y')
  set(gca,'Color',[0,0,0],'XColor',[0,1/2,0],'YColor',[0,1/2,0]);
  grid on;zoom off;zoom xon;
  hold on;
  load DPR_DataSupport DPRCover_PopMenuTEA_Value DPRCover_PopMenuTEB_Value CHSTD_Status CHTEA_Status;
  if DPRCover_PopMenuTEA_Value==4
      if any(CHSTD_Status)
          set(findobj('Tag','pushbutton9'),'Enable','on');
      end
  elseif DPRCover_PopMenuTEA_Value~=4
      if any(CHTEA_Status) & any(CHSTD_Status)
          set(findobj('Tag','pushbutton9'),'Enable','on');
      end
  end
  handles.DPR_DataTEB_Backup=DPR_DataTEB;
  handles.DPR_DataTEB=DPR_DataTEB;
  handles.Vertical_Factor_TEB=Vertical_Factor_TEB;
  handles.output = hObject;
  guidata(hObject, handles);
elseif DPRCover_PopMenuCOL_Value==2 
  load DPR_DataSupport path;
  name='*.wfm';
  filename=[path,name];
  [FileTEB_Name,path]=uigetfile(filename,'打开标准器具数据文件');        
  filename=[path,FileTEB_Name];                                       
  status=copyfile(filename,'FileTEB.wfm');
  Vertical_Factor_TEB=1;
  CHTEB_Status='NA';
  save DPR_DataSupport CHTEB_Status -append;
  !ConvertTekWfm FileTEB.wfm /MLAB
  [fid,message]=fopen('FileTEB.dat','r');
  [DPR_DataTEB,countTEB]=fscanf(fid,'%f',inf);
  save DPR_DataSupport FileTEB_Name path countTEB -append;
  if countTEB ~=handles.countSTD
    errordlg('数据长度错误')
    return
  else
    set(hObject,'String',FileTEB_Name);
  end
  x=1:countTEB;
  plot(x,DPR_DataTEB,'y')
  set(gca,'Color',[0,0,0],'XColor',[0,1/2,0],'YColor',[0,1/2,0]);
  grid on;zoom off;zoom xon;
  hold on;
  load DPR_DataSupport DPRCover_PopMenuTEA_Value DPRCover_PopMenuTEB_Value CHSTD_Status CHTEA_Status;
  if DPRCover_PopMenuTEA_Value==4
      if any(CHSTD_Status)
          set(findobj('Tag','pushbutton9'),'Enable','on');
      end
  elseif DPRCover_PopMenuTEA_Value~=4
      if any(CHTEA_Status) & any(CHSTD_Status)
          set(findobj('Tag','pushbutton9'),'Enable','on');
      end
  end
  handles.DPR_DataTEB_Backup=DPR_DataTEB;
  handles.DPR_DataTEB=DPR_DataTEB;
  handles.Vertical_Factor_TEB=Vertical_Factor_TEB;
  handles.output = hObject;
  guidata(hObject, handles);
elseif DPRCover_PopMenuCOL_Value==3 
  load DPR_DataSupport path;
  name='*.txt';
  filename=[path,name];
  [FileTEB_Name,path]=uigetfile(filename,'打开标准器具数据文件');        
  filename=[path,FileTEB_Name];                                       
  [fid,message]=fopen(filename,'rt');
  [DPR_DataTEB,countTEB]=fscanf(fid,'%f',inf);
  status=fclose(fid);
  save DPR_DataSupport FileTEB_Name path countTEB -append;
  if countTEB ~=handles.countSTD
    errordlg('数据长度错误')
    return
  else
    set(hObject,'String',FileTEB_Name);
  end
  Vertical_Factor_TEB=1;
  CHTEB_Status='NA';
  save DPR_DataSupport CHTEB_Status -append;
  x=1:countTEB;
  plot(x,DPR_DataTEB,'y')
  set(gca,'Color',[0,0,0],'XColor',[0,1/2,0],'YColor',[0,1/2,0]);
  grid on;zoom off;zoom xon;
  hold on;
  load DPR_DataSupport DPRCover_PopMenuTEA_Value DPRCover_PopMenuTEB_Value CHSTD_Status CHTEA_Status;
  if DPRCover_PopMenuTEA_Value==4
      if any(CHSTD_Status)
          set(findobj('Tag','pushbutton9'),'Enable','on');
      end
  elseif DPRCover_PopMenuTEA_Value~=4
      if any(CHTEA_Status) & any(CHSTD_Status)
          set(findobj('Tag','pushbutton9'),'Enable','on');
      end
  end
  handles.DPR_DataTEB_Backup=DPR_DataTEB;
  handles.DPR_DataTEB=DPR_DataTEB;
  handles.Vertical_Factor_TEB=Vertical_Factor_TEB;
  handles.output = hObject;
  guidata(hObject, handles);
end


% --- Executes on button press in pushbutton39.
function pushbutton39_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton39 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
cla;
load DPR_DataSupport CHSTD_Status CHTEA_Status CHTEB_Status countSTD;
load DPR_DataSupport DPRCover_PopMenuSTD_Value DPRCover_PopMenuTEA_Value DPRCover_PopMenuTEB_Value;
if DPRCover_PopMenuTEA_Value==4 & DPRCover_PopMenuTEB_Value~=4
    if any(CHSTD_Status) & any(CHTEB_Status)
        x=1:countSTD;
        plot(x,handles.DPR_DataSTD,'r',x,handles.DPR_DataTEB,'y')
        set(gca,'Color',[0,0,0],'XColor',[0,1/2,0],'YColor',[0,1/2,0]);
        grid on;zoom off;zoom xon;
    elseif any(CHSTD_Status) & ~any(CHTEB_Status)
        x=1:countSTD;
        plot(x,handles.DPR_