Matlab 常用读写命令#

1
%% read the 2-binary data
2
function data=read_bin_file(filename, num)
3
     fp = fopen(filename,'rb');
4
     data = fread(fp,num,"float");
5
    fclose(fp);
6
end
7

8
%% Ricker wavelet
9
function [wave]=ricker(fm, t, st,amp)
10
    wave = amp.*(1-2.*(pi.*fm.*(t-st)).^2).*exp(-(pi.*fm.*(t-st)).^2);
11
end
12

13
function [winer]= wiener(wave1f, wave2f)
14
    amp = real(wave2f).*real(wave2f)+imag(wave2f).*imag(wave2f);
15
%     amp = amp+ max(amp)*0.01;
16
    winerreal =  real(wave1f).*real(wave2f)+imag(wave1f).*imag(wave2f);
17
    winerimag =  imag(wave1f).*real(wave2f)-real(wave1f).*imag(wave2f);
18
    winerimag = winerimag./amp;
19
    winerreal = winerreal./amp;
20
    winer = winerreal + 1i*winerimag;
21
end
22

23
%% convert the .mat to .rsf
24
clear;clc;
25
maindir = uigetdir("choose a folder!");
26
matfiles = dir(fullfile(maindir));
27
LengthFiles = length(matfiles);
28

29
for i=3:LengthFiles   % note : change with the depth of folder
30
    name=matfiles(i).name;
31
    folder=matfiles(i).folder;
32
    data = importdata([folder,'/',name]);
33
    data(data==data(1,1)) = 1500.0;
34
    rsfname = [folder,'/rsf/'];
35
    rsfnamedata = [rsfname,'/data/'];
36
    mkdir(rsfname); mkdir(rsfnamedata);
37
    fp = fopen([rsfnamedata,name(1:end-4),'.rsf@'],"wb");
38
    fwrite(fp,data,"float");
39
    fclose(fp);
40
    strings = ['n1=201 n2=301 d1=10 d2=10 o1=0 o2=0 data_format=native_float in=./data/',name(1:end-4),'.rsf@']; % change with the paramaters
41
    fp2 = fopen([rsfname,name(1:end-4),'.rsf'],"w");
42
    fprintf(fp2,strings);
43
    fclose(fp2);
44
end