测量数据处理实验四 利用傅立叶分析进行测量数据处理

《测量数据处理实验四 利用傅立叶分析进行测量数据处理》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、《测量数据处理》实验报告（测绘工程专业使用）班级：测绘112班学号：2011122712姓名：许艳超淮海工学院测绘工程学院实验四利用傅立叶分析进行测量数据处理一、实验目的和要求1.理解离散傅立叶分析；2.了解快速傅立叶变换；3.能运用该方法进行简单的测量数据处理，并能够进行总结。二、实验时间、地点测量数据处理实验室三、实验内容1.对某点进行位移观测，得到8个不同时间的位移值，见下表，试进行傅立叶变换，并作出频谱图。表5.1观测数据时间k01234567位移f(k)2.33.03.84.25.14.53.21.82.对某变形点高程进行观测，得数据如下时间12345678910

2、1112131415高程1621751621561741571541771591711661621711831631617181920152165175167170试进行傅立叶变换，并作出频谱图。答：1、设计如下代码并输入matlab，得到频谱图：clc;clear;f=[2.33.03.84.25.14.53.21.8]';%e:expF=fft(f)FM=abs(F)t1=(0:1:7)fori=1:8t=0:0.001:FM(i);plot(t1(i),t,'b-')holdonEnd1、对已知数据进行傅立叶变换，利用matlab作出频谱图：p=[162,175,16

3、2,156,174,157,154,177,159,171,166,162,171,183,163,152,165,175,167,170];t=[1:length(p)];pfit=polyfit(t,p,5);fval=polyval(pfit,t);subplot(311);plot(t,p,'*',t,fval,'-');title('采用多项式拟合')xlabel('期数')ylabel('累计沉降量');pval=p-fval;subplot(312);plot(t,pval);title('去除逐式项')xlabel('期数')ylabel('残差值')y=f

4、ft(pval);m=abs(y);t=1:length(pval)subplot(313);plot(t,m);title('频谱分析');xlabel('频率值');ylabel('能量值');最后得到结果为下图：四、实验心得这次实验可以发现利用傅立叶变换在进行处理数据时能带来极大的方便，理解离散傅立叶分析，了解快速傅立叶变换，能运用该方法进行简单的测量数据处理，并能够进行总结所以学会利用傅立叶变换进行简单的数据处理非常重要。成绩：日期：指导教师：