总结12（11.22-11.28）

最新推荐文章于 2023-12-16 08:30:00 发布

xyjaxx.

最新推荐文章于 2023-12-16 08:30:00 发布

阅读量233

点赞数

分类专栏：硕士学习总结文章标签： julia

本文链接：https://blog.csdn.net/xyjaixuexi/article/details/121597810

版权

硕士学习总结专栏收录该内容

14 篇文章 4 订阅

订阅专栏

本周主要是安装了julia，并且配置将它配置进了Jupyter Notebook中，然后在安装地震数据处理相关的包时遇到了很多问题。经过上网查找和多次尝试最终成功安装上了SeisPlot, PyPlot, SeisReconstruction, SeisProcessing等包。然后成功运行了一些相关算法的代码。以下主要是运行POCS算法的结果和代码展示：

using Pkg
Pkg.build("PyCall")
using SeisPlot, PyPlot, SeisReconstruction, SeisProcessing
d = SeisLinearEvents(p1 = [-.001, 0.0015],tau=[1, 1/3],dx1=5); 
deci = SeisDecimate(d;perc=80);
param = Dict(:Niter=>100,:fmax=>60,:padt=>2,:padx=>2,:dt=>0.004)
dpocs = SeisPOCS(deci;param...);
subplot(121)
SeisPlotTX(deci,cmap="seismic",fignum=1,pclip=200,title="Decimated data")
subplot(122)
SeisPlotTX(dpocs[:,:,1,1,1],cmap="seismic",fignum=1,pclip=200,title="After POCS")

d = SeisLinearEvents(p1 = [-.001, 0.0015],tau=[1., 1/3],dx1=5); 
dec = SeisDecimate(d;perc=80);
param = Dict(:Niter=>100,:fmax=>60,:padt=>2,:padx=>2,:dt=>0.004)
dpocs = SeisPOCS(dec;param...);
figure(1, figsize=(10, 5))
subplot(121)
SeisPlotTX(dec,cmap="seismic",pclip=200,fignum=1,title="Decimated data",dy=0.004,ylabel="Time (s)", dx=5, xlabel="Position (m)")
subplot(122)
SeisPlotFK(dec, cmap="seismic", dy=0.004,hbox=5,pclip=200,fignum=1)

figure(2,figsize=(10, 5))
subplot(121)
SeisPlotTX(dpocs[:,:,1,1,1],cmap="seismic",fignum=2,pclip=200,title="After POCS",dy=0.004,ylabel="Time (s)", dx=5, xlabel="Position (m)")
subplot(122)
SeisPlotFK(dpocs[:,:,1,1,1],cmap="seismic",dy=0.004,fignum=2,pclip=200)