資源簡(jiǎn)介
實(shí)現(xiàn)利用matlab,輸入GPR數(shù)據(jù)最終輸出GPR波形堆積圖像

代碼片段和文件信息
clear?all;
clc
[num1?num2?num3]=textread(‘F:\liumaomao\mythesisdata\softwaredata\Matlab\boxingduiji\xd007.txt‘‘%f?%f?%f‘‘headerlines‘0);
%讀取導(dǎo)出數(shù)據(jù),num1位置,num2采樣時(shí)間,num3幅值
ts=320;?????%時(shí)窗
fs=1024;????%每道波采樣點(diǎn)數(shù)
tc=ts/fs;???%采樣間隔
nlen=length(num1);??%數(shù)據(jù)長(zhǎng)度
n=nlen/fs;??????????%道數(shù)
n3=save_in(num3nfs);%創(chuàng)建對(duì)應(yīng)于num2的保存數(shù)組并分列保存
[MN]?=?size(n3);?
M=[1:1:M];??????%時(shí)間序列
M=M‘;
t=tc.*M;????????%時(shí)間坐標(biāo)
ds1=0.02;???????%道間距
sx1=(1:n)*ds1;??%測(cè)線位置
figure(1);
imagesc(sx1tn3);%繪制灰度圖
colormap(gray);
xlabel(‘測(cè)線位置/m‘);ylabel(‘t/ns‘);title(‘灰度圖‘)
%生成波形堆積圖,由于采樣時(shí)道間距較小,顯示時(shí)適當(dāng)增大道間距
n4=n3(:1:4:n);??%選取適當(dāng)列數(shù)據(jù)顯示,注意與道間距對(duì)應(yīng)
[BC]?=?size(n4);?
ds2=0.08;???????%道間距
sx2=(1:C)*ds2;??%測(cè)線位置
figure(2);
wigb(n41sx2t)?%?生成波形圖,豎坐標(biāo)改成距離了
xlabel(‘測(cè)線位置/m‘);ylabel(‘t/ns‘);title(‘波形堆積圖‘)
?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----
?????文件????????953??2018-04-25?10:20??波形堆積\boxingduiji?-?副本.m
?????文件????????675??2015-09-22?18:25??波形堆積\duqushuju?-?副本.m
?????目錄??????????0??2018-07-05?19:53??波形堆積
-----------?---------??----------?-----??----
?????????????????1628????????????????????3
評(píng)論
共有 條評(píng)論