与时间无关的故事-CSDN博客

原创如何基于matlab软件读取sac文件中的地震数据

基于matlab读取标准sac文件的代码

2022-08-26 09:21:31 2483

原创（matlab代码）从sac文件中读取数据然后导出到txt文件里

for i = 1:30ieval([‘a=rsaca(’’’,num2str(i),’.sac’’);’]);b = zeros(13634399-756523,1);b(:,1) = a(756524:13634399,2);eval([‘save(’’’,num2str(i),’.txt’’,’‘b’’,’’-ASCII’’);’]);end

2021-11-25 15:39:43 1600

原创（matlab代码）wigb命令使用方法，地震记录规范出图

load ob_data_1_depth0.mat %读取地震数据dt=0.00025; %地震数据的采样间隔x=[1 21:2:196]; %所有地震道集的偏移距time=[0:dt:100000*dt-dt]; %采样点的时间轴(所有采样点对应的时刻)figurewigb(data1,1,x,time)xlabel(‘偏移距(m)’)ylabel(‘时间(s)’)如上图所示，共89道，第1道和第2道之间间隔20m，其他道间距均为2m...

2021-11-25 09:55:02 6039 1

原创（matlab代码）有关论文出图格式的命令

figureh1=plot(Xc,ga,‘ks-’,Xc1,ga1,‘ks-’,Xc2,ga2,‘ks-’,‘linewidth’,2)hold onh2=plot(Xc,GC(index,:),‘ro-’,Xc1,GC1(index,:),‘ro-’,Xc2,GC2(index,:),‘ro-’)% 设置坐标区的字体格式set(gca,‘fontsize’,10,‘fontname’,‘Times’);% 设置横纵坐标轴标签的字体格式ylabel(‘相速度(m/s)’,‘FontSize’,

2021-04-02 13:29:53 528

原创（matlab代码）绘制地震记录的F-K谱

% 需要读入数据：地震记录DATA_DISP 采样时间dt 道间距dxdatatt = DATA_DISP’;%转置后横轴为时间% dt=0.005;dv = 1;dx = 1; %【改！道间距】vmax=vmax_plot;vmin=5; % 设置速度间隔最大最小扫描速度nt=size(datatt,2); % 采样点数f_samplt = 1/dt; % 采样频率nx=size(datatt,1); % 道数kx=1/dx; % 波数nx_plus = 4096;nt_plus

2021-03-24 19:30:53 5449 2

原创（matlab代码）地震数据带通滤波

DATA_for_bandpass=data1; % 选择要分析的数据f1=0.001; f2=1; f3=15; f4=20;for j=1:size(DATA_for_bandpass,2)data_bandpass(:,j)=bandpass_filter(DATA_for_bandpass(:,j),dt,f1,f2,f3,f4);enddata3=data_bandpass;figureplot(data3,‘k’)% axis([0,size(data3,1),-1,1]);

2020-12-29 21:35:21 1477 1

原创（matlab代码）频谱分析(1)

data_pin=s1(:,2); % 选择要分析的数据dt=0.00025; % 采样间隔fftplot(data_pin,dt,‘b’,1,1);axis([0,100,0,1]); % 设置坐标轴显示范围grid on; % 显示网格线

2020-12-29 21:18:19 879

旅人2870的博客

原创（matlab）地震数据频谱分析-频谱图代码