傅立叶变换函数FFT 是最重要的也是非常典型的,其
调用格式为fft (x )。
程序代码及主要说明如下:
y =fft (x ); %对观测序列x 进行傅立叶变换m =abs (y ); %求能量值t =0:length (x )-1;
p l ot (t,m ); %画出频谱分析结果图,如图2所示
title (‘频谱分析’);xlabel (‘频率值’);ylabel (‘能量值’
)
;图2 能量-频率图3 应用于变形可视化
MAT LAB 具有强大的科学数据可视化能力,能够
将计算的结果以二维、三维图形表现出来。
311 绘制变形过程曲线图
MAT LAB 中的绘图命令p l ot (X,Y,S ),可以方便的
绘制出任意形状的线状图。以沉降监测为例,下列代码
可以绘制出某监测点的下沉过程曲线图(见图3)
。
图3 下沉过程曲线图
p l ot ([t,p,‘-o ’
) %t 为监测时间,p 为累计沉降序列legend (‘累计沉降量’,1)
title (‘下沉过程曲线图’)xlabel (‘时间序列’)ylabel (‘累计沉降值’
)312 绘制等值线图
MAT LAB 中,实现等高线的绘制非常简单,直接用
函数cont our (x,y,z )绘出等高线,函数cont our3(x,y,
z
)生成立体等高线。结果如图4
、图5。
图4 二维等值线图
图5 三维等值线图
313 变形3D 可视化
如果我们对所观测的区域性变形数据,绘制出不同时期的三维等值线图,这样就可以很直观的了解这个区域的整体变形趋势。
MABLAB 不仅提供了卓越的二维图形功能,而且
其三维功能也很强大,MAT LAB 函数可以对不规则
DE M 采样数据点(x,y,z ).的数据进行等距化构造格
网DE M 数据,采用x =lins pace (x m in,x max,n )函数沿x 方向在最小值x m in 和最大值x max 之间均匀设点n