地磁低点位移.ppt

新建台站增大密度不现实。可把青海、宁夏、陕西的地磁台观测数据整合起来，一是可以形成一个完整的闭合区域(矩形或圆形)利于绘制突变界线，二是可以凑出足够多的台站数量保证突变界线的精确度。 加强台站的维护管理，使得数据不要间断。 改进地磁低点位移地震预报的程序(软件)，使其功能更完善。 clc; clear; vv=load(D:\MapSIS\H01结果_数据\62001_2_31230_分.TXT); vv=num2str(vv(end,1)); v1=str2double(vv(1:4)); v2=str2double(vv(5:6)); v3=str2double(vv(7:8)); vv=[v1 v2 v3]; %******地磁数据使用的是格林尼治时间，加上7小时之后是兰州当地时间 path=D:\MapSIS\H01结果_数据\; %输入路径** ww=[2014,07,16]; %输入绘图日日期** 绘图日期必须在工作日期之前 %*********************************** rr=datenum(ww)-datenum([2014,1,1])+1; [XX,nn]=inp(path,vv); gg=draw(XX,ww,nn,rr); %% 地图模块--------------------- nl={兰州1,兰州2,兰州3,山丹4,天水5,天水6,天水7,嘉峪关8,嘉峪关9,... 嘉峪关10,嘉峪关11,合作12,合作13,玛曲14,玛曲15,舟曲16,舟曲17,英鸽18,... 寺滩19,芦阳20,松山21,古丰22,黄羊23,横梁24,临夏25,肃北26,肃北27};%27个 [ax,tl]=drap(XX,gg,nl,rr); %% ------------------------ noo=14; tt=dras(XX,noo,nl,ww,rr); function [XX,nn,yy]=inp(path,vv) %UNTITLED2 此处显示有关此函数的摘要 % 此处显示详细说明 FileName=dir(fullfile(path,*.TXT)); FileName=struct2cell(FileName); Q=FileName(1,:); %提取文件名矩阵的第一行信息 nn=length(Q); pp=datenum(vv)-datenum([2014,1,1])+1;%工作天数 XX=zeros(nn,pp); %文件夹数(台站数)*工作天数 for k=1:nn; %num(2)=27个文件夹=台站数 QQ=strcat(path,Q{k}); SS=load(QQ); DD=SS(:,2:end); DD(DD==999999.0)=NaN; %把矩阵DAE中所有999999的数值用NaN来代替 [yy,XX(k,:)]=min(DD.); %输出单个台站从年初到工作日每日的地磁地点位移时间序列 end XX(XX==1)=NaN; XX=XX+420; %格林尼治时间转化为北京时间7*60=420 子程序inp function gg=draw(XX,ww,nn,rr) %UNTITLED2 此处显示有关此函数的摘要 % 此处显示详细说明 figure; plot(1:nn,XX(:,rr),*-b); set(gca,xtick,1:1:nn,xlim,[1 nn],ytick,(4:9)*120,ylim,[1 3]*360,yticklabel,... {8点,10点,12点,14点,16点,18点},FontName,Times New Roman,FontSize,15); xlabel(仪器序号, FontName,Times New Roman,FontSize,16,FontWeight,Bold); ylabel(地磁低点位移时间,FontName,Times New Roman,FontSize,16,FontWeight,Bold); hold on; grid on; plot([1 nn],[1 1]*max(XX(:,rr)),r,[1 nn],[1 1]*min(XX(:,rr)),r); md=max(XX(:,rr))-min(XX(:,rr)); uu=strcat(md=,num2str(md),分钟); text(nn/2,max(XX(:,rr)),uu,FontName,Times New Roman,FontSize,16,FontWeight,Bold); if md120 er=\color{