MATLAB设计异常延拓源程序.doc

磁异常转换计算作业 设计思想：设计一个实测磁异常延拓用户界面，用以实现对水平圆柱体产生的各个磁异常分量进行延拓的人机交互可视化。界面介绍：1菜单栏中model菜单，包含两个子菜单单模型和组合模型；用户可以选择模型，并读出该模型的磁异常数据，赋给相应变量（T,Za,Ha）。2界面右上方异常选择弹出式菜单，包含三个选项：delT | Za | Ha ；用户可以选择磁异常分量（用鼠标点击），左边的绘图界面将会根据该磁异常分量的实测数据，绘制其剖面图。3界面右下方延拓高度编辑框，用户可以在此输入所需的已选磁异常分量的延拓高度（m）（以Enter键确认），左边的绘图界面将会根据所输入的延拓高度，绘制已选磁异常分量延拓后的剖面图。 此用户界面只适用于向上延拓，即输入的数值应=0; 当输入为0时，绘图界面中的剖面图为实测剖面图；用户在使用该界面时，应先选择模型和磁异常分量，然后才能对其进行延拓；用户界面：程序代码：clc;clear;%设置用户界面f=figure(714);set(f,resize,on,menubar,none,name,my demo,color,0.3 1 0.8);%菜单%模型选择 ( 单模型 | 组合模型 ) ui_model=uimenu(f,Label,model);%单模型 sin_singleui_sin=uimenu(ui_model,Label,单模型,callback,. fid=fopen(D:新建文件夹延拓和分量转换数据_单模型产生delT.dat,r);,. T=fscanf(fid,%g %g %g,3 157);,. fclose(fid);,. fid=fopen(D:新建文件夹延拓和分量转换数据_单模型产生Za.dat,r);,. Za=fscanf(fid,%g %g %g,3 157);,. fclose(fid);,. fid=fopen(D:新建文件夹延拓和分量转换数据_单模型产生Ha.dat,r);,. Ha=fscanf(fid,%g %g %g,3 157);,. fclose(fid); );%组合模型 com_combinationui_com=uimenu(ui_model,Label,组合模型,callback,. fid=fopen(D:新建文件夹延拓和分量转换数据_组合模型产生delT_组合.dat,r);,. T=fscanf(fid,%g %g %g,3 157);,. fclose(fid);,. fid=fopen(D:新建文件夹延拓和分量转换数据_组合模型产生Za_组合.dat,r);,. Za=fscanf(fid,%g %g %g,3 157);,. fclose(fid);,. fid=fopen(D:新建文件夹延拓和分量转换数据_组合模型产生Ha_组合.dat,r);,. Ha=fscanf(fid,%g %g %g,3 157);,. fclose(fid); );%绘图界面subplot(position,0.05,0.2,0.7,0.6);title(实测磁异常延拓用户界面,fontsize,15);%异常选择弹出式菜单 t2=uicontrol(f,. style,text,. string,异常选择:,. backgroundcolor,1 0.5 1,. fontsize,12,. Units,Normalized,. position,0.8,0.7,0.15,0.06,. HorizontalAlignment,left); pop=uicontrol(f,style,popup,. string,delT|Za|Ha,. Units,Normalized,. position,0.8,0.62,0.15,0.06,. HorizontalAlignment,left,. backgroundcolor,1 1 0,. fontsize,12,. enable,on,. callback, val = get(pop,Value);,. if (val = 1);,. W=T;,. elseif (val = 2);,. W=Za;,. elseif (val = 3);,. W=Ha;,. end;); %上延高度编辑框t3=uicontrol(f,. style,text,. string,上延高度(m):,. backgroundcolor,1 0.5 1,. fontsize,12,. Units,Normalized,. position,0.8,0.49,0.15,0.06,. HorizontalAlignment,left); high=uicontrol(f, . Style,edit, . HorizontalAlignment,left, . backgroundcolor,1 1 0,. Units,normalized, . Position,0.8,0.41,0.15,0.06,. fontsize,12,. enable,on,. callback,. plot(W(1,:),W(3,:),r);,. hold on;,. h=str2num(get(gcbo,string);,. A=fft(W(3,:);,. B=A.*exp(-h/102*W(1,:);,. b=ifft(B);,. plot(W(1,:),b);,. hold off;); 界面操作：图1：单个模型Za上延2m和实测数据的剖面对比图图2：单个模型Za上延30m和实测数据的剖面对比图图3：组合模型delT上延4m和实测数据的剖面对比图图4：组合模型delT上延15m和实测数据的剖面对比图实验小结： 通过对实测磁异常延拓用户界面设计，学习了如何通过编程实现用户界面的设计，对MATLAB的知识有了进一步的掌握。使用实测磁异常延拓用户界面，输入不同的参