# MatLab - PowerPoint PPT Presentation

1 / 11

MatLab. -et krasjkurs. www.stud.ntnu.no/~gjerding/data. >> load_well; >> well well = depth: [18163x1 double] Vp: [18163x1 double] Vs: [18163x1 double] rho: [18163x1 double] name: {1x9 cell} top: [1x9 double] bottom: [1x9 double].

I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.

MatLab

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

## MatLab

-et krasjkurs

### www.stud.ntnu.no/~gjerding/data

>> well

well =

depth: [18163x1 double]

Vp: [18163x1 double]

Vs: [18163x1 double]

rho: [18163x1 double]

name: {1x9 cell}

top: [1x9 double]

bottom: [1x9 double]

• Plott Vp mot rho for ’Ty’ formasjonen og trekk en regresjonslinje

Top Ty

Bunn Ty

rho

Vp

Gardner et al. (1974):

a?

Brønn:

Top ty

id=

Bunn ty

close all; clear

frm = 4;

id = find( well.depth >= well.top(frm) & well.depth <= well.bottom(frm) ) ;

% ...the regression coefficients

a = sum( well.Vp(id).^.25.* well.rho(id) ) ...

/sum(( well.Vp(id) .^.25).^2);

x = [ min(well.Vp(id)) : max(well.Vp(id)) ];

subplot(2,2,1)

hold on

plot( well.Vp(id) , well.rho(id) ,'r.' )

legend('well',4)

plot(x,a*x.^.25,'b')

xlabel('Vp [m/s]'); ylabel('rho [g/cc]');

title(well.name(frm))

axis tight

grid

str = num2str(a);

text(2380,2.20,strcat('Gardner: rho = ',str,'Vp^0^.^2^5'))

Seiswiggle(data)

mesh(data)

view(60,50)

View(0,90)

### Andre ting

• SEGY => MatLab

• www.stud.ntnu.no/~gjerding

• MatLab – Fortran / C++

• MatLab – unix

• Perl, awk, sed etc.

• unix