測量數(shù)據(jù)處理實驗四 利用傅立葉分析進(jìn)行測量數(shù)據(jù)處理

《測量數(shù)據(jù)處理實驗四 利用傅立葉分析進(jìn)行測量數(shù)據(jù)處理》由會員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。

1、《測量數(shù)據(jù)處理》實驗報告（測繪工程專業(yè)使用）班級：測繪112班學(xué)號：2011122712姓名：許艷超淮海工學(xué)院測繪工程學(xué)院實驗四利用傅立葉分析進(jìn)行測量數(shù)據(jù)處理一、實驗?zāi)康暮鸵?.理解離散傅立葉分析；2.了解快速傅立葉變換；3.能運(yùn)用該方法進(jìn)行簡單的測量數(shù)據(jù)處理，并能夠進(jìn)行總結(jié)。二、實驗時間、地點測量數(shù)據(jù)處理實驗室三、實驗內(nèi)容1.對某點進(jìn)行位移觀測，得到8個不同時間的位移值，見下表，試進(jìn)行傅立葉變換，并作出頻譜圖。表5.1觀測數(shù)據(jù)時間k01234567位移f(k)2.33.03.84.25.14.53.21.82.對某變形點高程進(jìn)行觀測，得數(shù)據(jù)如下時間12345678910

2、1112131415高程1621751621561741571541771591711661621711831631617181920152165175167170試進(jìn)行傅立葉變換，并作出頻譜圖。答：1、設(shè)計如下代碼并輸入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、對已知數(shù)據(jù)進(jìn)行傅立葉變換，利用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('期數(shù)')ylabel('累計沉降量');pval=p-fval;subplot(312);plot(t,pval);title('去除逐式項')xlabel('期數(shù)')ylabel('殘差值')y=f

4、ft(pval);m=abs(y);t=1:length(pval)subplot(313);plot(t,m);title('頻譜分析');xlabel('頻率值');ylabel('能量值');最后得到結(jié)果為下圖：四、實驗心得這次實驗可以發(fā)現(xiàn)利用傅立葉變換在進(jìn)行處理數(shù)據(jù)時能帶來極大的方便，理解離散傅立葉分析，了解快速傅立葉變換，能運(yùn)用該方法進(jìn)行簡單的測量數(shù)據(jù)處理，并能夠進(jìn)行總結(jié)所以學(xué)會利用傅立葉變換進(jìn)行簡單的數(shù)據(jù)處理非常重要。成績：日期：指導(dǎo)教師：