Matlab插值计算各时刻磁法勘探日变观测值

Matlab插值计算各时刻磁法勘探日变观测值

在磁法勘探中，消日变影响的改正称为日变改正。进行日变改正时必须设立日变站，观测日变情况。根据日变数据和测点观测时间，对观测数据进行改正。

在本次磁法实习中，日变观测间隔为300s，观测的时间段早于最早的磁力测量，晚于最晚的磁力测量。在整理测量数据进行日变改正时就可以根据进行磁力测量的时刻的日变观测值作日变改正。

假设日变观测数据的记录时刻分别为1000、1300、1600……3700、4000，而进行磁测时的时刻可能为1000-3000之间的任意整数，因此就可以根据日变观测数据进行插值，进而找到任意时刻对应的日变观测场值。Matlab可以很方便地进行这样的计算：

function mT=getAllRibian(ribianTime,ribianT,mTime)
%根据日变站数据和磁测时间数据，计算磁测各时刻对应的日变观测值。
%ribianTime:日变观测的各个时间向量
%ribianT:日变观测的T值
%mTime:磁测各个的时间向量
%mT:磁测各时刻对应的日变观测值