从最简单的入手,1.脚本命令M文件,一般没有输入输出,2.函数M文件,一般有输入输出.
脚本命令M文件简单,下面给出两种M文件的形式,先不看内容,内容还需要修行才能看懂。
如:%一个简单的从excel中读取经纬度来计算震中距的程序
clear
clc
fileDir='F:\MATLAB WORK\111\';
property='wenchuan_para.xlsx';
loglat=xlsread([fileDir,property],'para','I2:J1215');
du=pi/180;
R=6371.393;
alfa1=103.4*du;
beta1=31*du;
for i=1:length(loglat(:,1))
alfa2=loglat(i,1)*du;
beta2=loglat(i,2)*du;
Garc(i,1)=R*acos(cos(beta1)*cos(beta2)*cos(alfa1-alfa2)+sin(beta1)*sin(beta2));
end
函数M文件:
%进行基线校正和滤波处理,并输出为新的地震波
%对不同文件夹和不同地震要在第5、6行设置文件位置和文件通配符
function correct
clear;clc;
fdir='E:\records and programs for
phd\2.NIED(KIK-NETK-NET)\s