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

% 需要读入数据：地震记录DATA_DISP 采样时间dt 道间距dx
datatt = 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=5000;
datatt(nx_plus,nt_plus)=0;%道数方向，时间方向充零
fk = abs(fft2(datatt));%t-x转化到f-k域
fk = fftshift(fk);%零频部分移动至中心
% 将双侧频谱fk转化到单侧频谱fk_side 其频率波数间隔df dk,横纵坐标范围:f k 坐标点个数fn kn
f_max=20; % 选择频散图中的最大频率
df = f_samplt/nt_plus;
dk = kx/nx_plus;
fk_side = fk(1:end/2,end/2:end/2+floor(f_max/df));
kn=size(fk_side,1);fn=size(fk_side,2);
k = linspace(0,kx/2,kn);
f2 = linspace(0,f_max,fn);
% figure
% resolution=500;
% [~,ch]=contourf(f2,fliplr(k),fk_side,resolution);
% set(ch,‘linestyle’,‘none’);
% title(‘F-K域能量图’);
% xlabel(‘频率 HZ’);ylabel(‘波数 m-1’);
% axis([0 100 0 0.25])
% colormap(jet)

（matlab代码）绘制地震记录的F-K谱相关推荐

非均匀三次B样条曲线插值实现及MATLAB代码
这篇博客跟我上一篇博客<均匀三次B样条曲线插值实现及MATLAB代码>的内容有点像,只是在基函数的计算上不同,造成均匀/非均匀的区别. 参考资料: [1](这个PPT讲得很通俗,但对于多插 ...
均匀三次B样条曲线插值实现及MATLAB代码
参考资料: [1](这个PPT讲得很通俗,但对于多插值点分段曲线的内容漏讲了一个知识点)三次周期B样条曲线的算法 - 百度文库 (baidu.com) [2](这个介绍只有两个插值点的三次B样条曲线, ...
VMD分解，matlab代码，包络线，包络谱，中心频率，峭度值，能量熵，近似熵，包络熵，频谱图，希尔伯特变换，包含所有程序MATLAB代码，-西储大学数据集为例
目录 1.选取数据 2.VMD函数-matlab代码 3.采用matlab脚本导入数据并做VMD分解 4. VMD分解图 5.计算中心频率 6.画包络线 7. 画包络谱 8. 计算峭度值 9.计算能量 ...
使用Matlab绘制语音信号的语谱图
本文绘制语音信号的语谱图主要使用了spectrogram函数,spectrogram是一个MATLAB函数,使用短时傅里叶变换得到信号的频谱图.当使用时无输出参数,会自动绘制频谱图:有输出参数,则会返 ...
单频信号的相位谱计算与误差修正-附Matlab代码
一.问题描述我们在实际处理时经常遇到只有一个正弦信号的情况,其频率为 f 0 {{f}_{0}} f0,在谱分析以后,除了在频率为 f 0 {{f}_{0}} f0处有相位数值外,其他频率处都有 ...
【图像分割】基于萤火虫优化的半监督谱聚类彩色图像分割方法（Matlab代码实现）
目录 0 概述 1 萤火虫算法 1.1 思想来源 1.2 数学模型 1.2.1 萤火虫相对荧光亮度 1.2.2 萤火虫吸引度 1.2.3 位置更新 1.3 算法步骤 2 基于萤火虫优化的半监督谱聚类彩 ...
[机器学习]多视角谱聚类(Ncut matlab代码)
一.多视角多视角聚类是伴随着信息时代发展过程中人们获取信息的途径日益增加,而逐渐出现的新聚类研究方向.通过多种途径获得的数据其实是对同件事物的不同描述,如电影的图片以及音频都是对这个电影的描述, ...
matlab三元函数泰勒展开,如何用matlab绘制三元函数f(x,y,z)=x^2 y^3 z^4在[-1,1]x[-2,2]x[-3,3] 上的各...
问题描述: 如何用matlab绘制三元函数f(x,y,z)=x^2 y^3 z^4在[-1,1]x[-2,2]x[-3,3] 上的各种切片图. 还有怎么绘制f(x,y,t)=10*exp(1-t)*s ...
如何用matlab画三维函数相交,如何用matlab绘制三元函数f=x^2 y^3 z^4在
www.mh456.com防采集. 1234567syms x y z;f = (x-y)^2+(x-z)^2+(y-z)^2-100;z = solve(f,'z');for p = 1:lengt ...
k均值算法原理详细讲解以及matlab代码实现
有研究生物电信号处理和机器学习的欢迎加我qq429500506共同交流学习进步. 最近更新文章的频率太低了,主要原因是不想为了发文章而发文章,想潜心研究,写好文章,顺便想说一句开源万岁,最近一个月虽然 ...

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

（matlab代码）绘制地震记录的F-K谱相关推荐

最新文章

热门文章