高程拟合程序.docVIP

function varargout = Open_Callback(h, eventdata, handles, varargin) [filename,pathname]=uigetfile(,打开数据文件：); file=[pathname,filename]; % ----------------------------------------------------------------- function varargout = Exit_Callback(h, eventdata, handles, varargin) close(gcbf); % ----------------------------------------------------------------- function varargout = About_Callback(h, eventdata, handles, varargin) h1=findobj(0,tag,text1); h2=findobj(0,tag,slider2); h3=findobj(0,tag,edit1); h4=findobj(0,tag,axes1); h5=findobj(0,tag,pushbutton1); h6=findobj(0,tag,edit2); h7=findobj(0,tag,text4); set(h1,visible,off); set(h2,visible,off); set(h3,visible,on); set(h3,horizontalalignment,center); set(h4,visible,off); set(h5,visible,off); set(h6,visible,off); set(h7,visible,on); set(h7,string,程序使用说明); % ----------------------------------------------------------------- function varargout = Dxs_Callback(h, eventdata, handles, varargin)%多项式 h1=findobj(0,tag,text1); h2=findobj(0,tag,slider2); h3=findobj(0,tag,edit1); h4=findobj(0,tag,axes1); h5=findobj(0,tag,pushbutton1); h6=findobj(0,tag,edit2); h7=findobj(0,tag,text4); set(h1,visible,off); set(h2,visible,off); set(h3,visible,off); set(h4,visible,on); set(h5,visible,off); set(h6,visible,off); global X Y Z %全局变量 [filename,pathname]=uigetfile({*.txt,文本文件(*.txt),},打开已知点的坐标数据文件(yz.txt)：); file=[pathname,filename]; load(file) b(1:6,1)=1; b(1:6,2)=yz(1:6,1); b(1:6,3)=yz(1:6,2); b(1:6,4)=yz(1:6,1).^2; b(1:6,4)=yz(1:6,1).^2; b(1:6,5)=yz(1:6,2).^2; b(1:6,6)=yz(1:6,1).*yz(1:6,2);%b是系数阵 X c=inv(b);%求逆 [filename,pathname]=uigetfile({*.txt,文本文件(*.txt),},打开已知点的高程异常数据文件(gcyc.txt)：); file=[pathname,filename]; load(file) %load d:\gcyc.txt; d=c*gcyc;%d是所需的系数阵 [filename,pathname]=uigetfile({*.txt,文本文件(*.txt),},打开待求点的坐标数据文件(dq.txt)：); file=[pathname,filename]; load(file) %load d:\dq.txt;%加载待求点的坐标 for i=1:21 x=dq(i,1); y=dq(i,2); f(i,1)=d(1,1)+d(2,1)*x+d(3,1)*y+d(4,1)*x^2+d(5,1)*y^2+d(6,1)*x*y; end [filename,pathname]=uigetfile({*.txt,文本文件(*.t