baobiao007

木蟲 (職業(yè)作家)

中國特色

應(yīng)助: 201 (大學(xué)生)
金幣: 6482.7
散金: 557
紅花: 40
帖子: 3050
在線: 1009.9小時(shí)
蟲號(hào): 505962
注冊: 2008-02-18
專業(yè): 應(yīng)用地球物理學(xué)

[交流] 分享用matlab顯示地震記錄的波形變面積圖已有9人參與

wigb函數(shù)是網(wǎng)上找的，自己添了個(gè)主程序方便使用。

主程序
function PlotRecord_Wigb(filename,nx,nt,dx,dt,zy)
%波形變面積顯示地震記錄
% clear; clf;clc;
% nx=200; %道數(shù)
% nt=2001; %采樣點(diǎn)數(shù)
% dx=10.;    %道間距
% dt=0.001; %時(shí)間采樣間隔
% zy=1.0;    %增益
% filename='r200-2001-1ms.dat';
clf;
fid=fopen(filename,'r');
[A,COUNT]=fread(fid,[nt,nx],'float32');
fclose(fid);

x=0:dx: (nx-1)*dx;
z=0:dt: (nt-1)*dt;
wigb(A,zy,x,z);

end

子函數(shù)
function wigb (a,scal,x,z,amx)
%WIGB: Plot seismic data using wiggles.
%
%  WIGB(a,scal,x,z,amx)
%
%  IN a:    seismic data
%       scale: multiple data by scale
%       x:    x-axis;
%       z:    vertical axis (time or depth)
%
%    If only 'a' is enter, 'scal,x,z,amn,amx' are decided automatically;
%    otherwise, 'scal' is a scalar; 'x, z' are vectors for annotation in
%    offset and time, amx are the amplitude range.
%
%
%  Author(s): Xingong Li (Exxon-Mobil)
%  Copyright 1998-2003 Xingong
%  Revision: 1.2  Date: Dec/2002
%
%  Signal Analysis and Imaging Group (SAIG)
%  Department of Physics, UofA
%

if nargin == 0, nx=10;nz=10; a = rand(nz,nx)-0.5; end;

[nz,nx]=size(a);

trmx= max(abs(a));
if (nargin <= 4); amx=mean(trmx);  end;
if (nargin <= 2); x=[1:nx]; z=[1:nz]; end;
if (nargin <= 1); scal =1; end;

if nx <= 1 ; disp(' ERR: PlotWig: nx has to be more than 1');return;end;

% take the average as dx
dx1 = abs(x(2:nx)-x(1:nx-1));
dx = median(dx1);

dz=z(2)-z(1);
xmx=max(max(a)); xmn=min(min(a));

if scal == 0; scal=1; end;
a = a * dx /amx;
a = a * scal;

fprintf(' PlotWig: data range [%f, %f], plotted max %f \n',xmn,xmx,amx);

% set display range
x1=min(x)-2.0*dx; x2=max(x)+2.0*dx;
z1=min(z)-dz; z2=max(z)+dz;

set(gca,'NextPlot','add','Box','on', ...
  'XLim', [x1 x2], ...
  'YDir','reverse', ...
  'YLim',[z1 z2]);

fillcolor = [0 0 0];
linecolor = [0 0 0];
linewidth = 0.1;

z=z'; % input as row vector
zstart=z(1);
zend  =z(nz);

for i=1:nx,

  if trmx(i) ~= 0; % skip the zero traces
tr=a(:,i); % --- one scale for all section
    s = sign(tr) ;
    i1= find( s(1:nz-1) ~= s(2:nz) ); % zero crossing points
npos = length(i1);

%12/7/97
zadd = i1 + tr(i1) ./ (tr(i1) - tr(i1+1)); %locations with 0 amplitudes
aadd = zeros(size(zadd));

[zpos,vpos] = find(tr >0);
[zz,iz] = sort([zpos; zadd]); % indices of zero point plus positives
aa = [tr(zpos); aadd];
aa = aa(iz);

% be careful at the ends
if tr(1)>0, a0=0; z0=1.00;
else, a0=0; z0=zadd(1);
end;
if tr(nz)>0, a1=0; z1=nz;
else, a1=0; z1=max(zadd);
end;

zz = [z0; zz; z1; z0];
aa = [a0; aa; a1; a0];

zzz = zstart + zz*dz -dz;

patch( aa+x(i) , zzz,  fillcolor);

line( 'Color',[1 1 1],'EraseMode','background',  ...
      'Xdata', x(i)+[0 0], 'Ydata',[zstart zend]); % remove zero line

%'LineWidth',linewidth, ...
%12/7/97 'Xdata', x(i)+[0 0], 'Ydata',[z0 z1]*dz); % remove zero line

line( 'Color',linecolor,'EraseMode','background',  ...
   'LineWidth',linewidth, ...
   'Xdata', tr+x(i), 'Ydata',z); % negatives line

else % zeros trace
line( 'Color',linecolor,'EraseMode','background',  ...
   'LineWidth',linewidth, ...
      'Xdata', [x(i) x(i)], 'Ydata',[zstart zend]);
end;
end;

[ Last edited by baobiao007 on 2012-4-9 at 10:41 ]

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

地球物理資源集錦

小胖的專輯

» 本帖已獲得的紅花（最新10朵）

逄碩

» 猜你喜歡

調(diào)劑西南科技大學(xué) 環(huán)資學(xué)院已經(jīng)有4人回復(fù)
無法從這份源碼判斷項(xiàng)目是否“獲資助批準(zhǔn)” 已經(jīng)有11人回復(fù)
地球化學(xué)論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有222人回復(fù)
西南科技大學(xué) 調(diào)劑綿陽已經(jīng)有10人回復(fù)
散金幣，求基金中標(biāo)！！已經(jīng)有39人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

如何用MATLAB做出地震正演模型已經(jīng)有5人回復(fù)
地震面波英文已經(jīng)有7人回復(fù)
如何用matlab產(chǎn)生markov波形已經(jīng)有3人回復(fù)
【原創(chuàng)】從皇極數(shù)經(jīng)談地震預(yù)報(bào) 已經(jīng)有6人回復(fù)

我同意叔本華的觀點(diǎn)，人們投身藝術(shù)和科學(xué)領(lǐng)域的強(qiáng)烈愿望之一就是逃離痛苦、殘酷和枯燥無味的現(xiàn)實(shí)生活，逃離自己飄忽不定的七情六欲的桎梏。--愛因斯坦

1樓 2012-04-09 10:35:48

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

smithjohn

金蟲 (正式寫手)

應(yīng)助: 4 (幼兒園)
金幣: 1334.9
散金: 372
紅花: 4
帖子: 673
在線: 148.8小時(shí)
蟲號(hào): 1030430
注冊: 2010-05-28

0.8

板塊做地震的不多
樓主試過實(shí)際地震資料么？

贊一下(1人)

回復(fù)此樓

2樓2012-04-09 18:13:25

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

baobiao007

木蟲 (職業(yè)作家)

中國特色

應(yīng)助: 201 (大學(xué)生)
金幣: 6482.7
散金: 557
紅花: 40
帖子: 3050
在線: 1009.9小時(shí)
蟲號(hào): 505962
注冊: 2008-02-18
專業(yè): 應(yīng)用地球物理學(xué)

沒有，實(shí)際地震資料要根據(jù)具體的數(shù)據(jù)格式來弄，勘探地震有自己行業(yè)標(biāo)準(zhǔn)的數(shù)據(jù)格式。天然地震的就不清楚了

贊一下(1人)

回復(fù)此樓

3樓2012-04-09 18:21:52

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

baobiao007

木蟲 (職業(yè)作家)

中國特色

應(yīng)助: 201 (大學(xué)生)
金幣: 6482.7
散金: 557
紅花: 40
帖子: 3050
在線: 1009.9小時(shí)
蟲號(hào): 505962
注冊: 2008-02-18
專業(yè): 應(yīng)用地球物理學(xué)

引用回帖:

9樓: Originally posted by 逄碩 at 2013-05-05 19:25:49
這個(gè)導(dǎo)入的數(shù)據(jù)應(yīng)該是什么參數(shù)啊

普通二進(jìn)制地震記錄數(shù)據(jù)

贊一下(1人)

回復(fù)此樓

10樓2013-05-05 19:29:27

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

cj000211

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 4
帖子: 13
在線: 5.7小時(shí)
蟲號(hào): 804526
注冊: 2009-07-06
專業(yè): 應(yīng)用地球物理學(xué)

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

多謝了~~~剛下載完，回去試試看

贊一下(1人)

回復(fù)此樓

11樓2013-05-07 11:20:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

lzk55882

謝謝分享

4樓2012-04-11 10:29:05

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

i星辰

銀蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1811.8
帖子: 403
在線: 230小時(shí)
蟲號(hào): 1142640
注冊: 2010-11-09
性別: GG
專業(yè): 應(yīng)用地球物理學(xué)

犀利啊，拜讀了

回復(fù)此樓

5樓2012-04-11 14:59:32

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

sapyren

銀蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 341.5
散金: 427
紅花: 3
帖子: 110
在線: 31.9小時(shí)
蟲號(hào): 1405389
注冊: 2011-09-17
性別: GG
專業(yè): 大地測量學(xué)

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

fname附件貼出來啊，這樣看好沒勁！

贊一下

回復(fù)此樓

做純粹的自己，享受責(zé)任與自由！�。�

6樓2012-04-26 20:07:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

應(yīng)助: 0 (幼兒園)
金幣: 5
帖子: 1
在線: 28分鐘
蟲號(hào): 2135987
注冊: 2012-11-19
專業(yè): 應(yīng)用地球物理學(xué)

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

送紅花一朵

這個(gè)導(dǎo)入的數(shù)據(jù)應(yīng)該是什么參數(shù)啊

贊一下

回復(fù)此樓

9樓2013-05-05 19:25:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 268求調(diào)劑 +7	簡單點(diǎn)0 2026-03-17	7/350	2026-03-21 08:32 by JourneyLucky
[考研] 070300化學(xué)319求調(diào)劑 +7	錦鯉0909 2026-03-17	7/350	2026-03-21 03:46 by JourneyLucky
[考研] 一志愿中國石油大學(xué)（華東）本科齊魯工業(yè)大學(xué) +3	石能偉 2026-03-17	3/150	2026-03-21 02:22 by JourneyLucky
[考研] 279分求調(diào)劑一志愿211 +11	chaojifeixia 2026-03-19	12/600	2026-03-21 01:49 by 星空星月
[考研] 354求調(diào)劑 +5	Tyoumou 2026-03-18	8/400	2026-03-21 00:35 by JourneyLucky
[考研] 296求調(diào)劑 +6	www_q 2026-03-18	10/500	2026-03-20 23:56 by JourneyLucky
[考研] 材料專碩英一數(shù)二306 +7	z1z2z3879 2026-03-18	7/350	2026-03-20 23:48 by JourneyLucky
[考研] 294求調(diào)劑材料與化工專碩 +15	陌の森林 2026-03-18	15/750	2026-03-20 23:28 by JourneyLucky
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8	小材化本科 2026-03-18	8/400	2026-03-20 23:16 by JourneyLucky
[考研] 287求調(diào)劑 +7	晨昏線與星海 2026-03-19	8/400	2026-03-20 22:19 by JourneyLucky
[考研] 085600材料與化工 +8	安全上岸！ 2026-03-16	8/400	2026-03-20 22:13 by luoyongfeng
[考研] A區(qū)線材料學(xué)調(diào)劑 +5	周周無極 2026-03-20	5/250	2026-03-20 21:33 by laoshidan
[考研] 求調(diào)劑 +3	eation27 2026-03-20	3/150	2026-03-20 19:32 by JourneyLucky
[考研] 081700化工學(xué)碩調(diào)劑 +3	【1】 2026-03-16	3/150	2026-03-19 23:40 by edmund7
[考研] 085600材料與化工調(diào)劑 324分 +10	llllkkkhh 2026-03-18	12/600	2026-03-19 14:33 by llllkkkhh
[考研] 311求調(diào)劑 +11	冬十三 2026-03-15	12/600	2026-03-18 14:36 by 星空星月
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國級(jí)證書40+數(shù)一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考博] 26博士申請(qǐng) +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111