 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 10154 | 回復(fù): 12 | |||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||||
baobiao007木蟲 (職業(yè)作家)
中國特色
|
[交流]
分享用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朵）
逄碩» 猜你喜歡
地球物理學(xué)和空間物理學(xué)論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有222人回復(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ù)
逄碩
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 5
- 帖子: 1
- 在線: 28分鐘
- 蟲號(hào): 2135987
- 注冊(cè): 2012-11-19
- 專業(yè): 應(yīng)用地球物理學(xué)
9樓2013-05-05 19:25:49
smithjohn
金蟲 (正式寫手)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 1334.9
- 散金: 372
- 紅花: 4
- 帖子: 673
- 在線: 148.8小時(shí)
- 蟲號(hào): 1030430
- 注冊(cè): 2010-05-28
2樓2012-04-09 18:13:25
baobiao007
木蟲 (職業(yè)作家)
中國特色
- 應(yīng)助: 201 (大學(xué)生)
- 金幣: 6482.7
- 散金: 557
- 紅花: 40
- 帖子: 3050
- 在線: 1009.9小時(shí)
- 蟲號(hào): 505962
- 注冊(cè): 2008-02-18
- 專業(yè): 應(yīng)用地球物理學(xué)
3樓2012-04-09 18:21:52
謝謝分享 |
4樓2012-04-11 10:29:05
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿中海洋材料357 +3 | 麥恩莉. 2026-03-30 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +5 | min3 2026-03-24 | 6/300 |
|
|
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +7 | nhj_ 2026-03-25 | 7/350 |
|
|
|
[考研] 085404求調(diào)劑，總分309，本科經(jīng)歷較為豐富 +6 | 來財(cái)aa 2026-03-25 | 6/300 |
|
|
|
[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
|
[考研] 調(diào)劑310 +12 | 溫柔的晚安 2026-03-25 | 13/650 |
|
|
|
[考研] 299求調(diào)劑 +10 | 15188958825 2026-03-25 | 10/500 |
|
|
|
[考研] 2026年華南師范大學(xué)歡迎化學(xué)，化工，生物，生醫(yī)工等專業(yè)優(yōu)秀學(xué)子加入！ +3 | llss0711 2026-03-28 | 6/300 |
|
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
|
[考研] 071000生物學(xué)求調(diào)劑，初試成績343 +7 | 小小甜面團(tuán) 2026-03-25 | 7/350 |
|
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
|
[考研] 331環(huán)境科學(xué)與工程求調(diào)劑 +3 | 熠然好運(yùn)氣 2026-03-27 | 3/150 |
|
|
|
[考研] 330一志愿中國海洋大學(xué) 化學(xué)工程 085602 有讀博意愿 求調(diào)劑 +3 | wywy.. 2026-03-27 | 4/200 |
|
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|