matlab程序 地震 相干噪声_无网格局部Petrov-Galerkin法及程序设计（15）——代码：计算RPIM形函数...

最新推荐文章于 2023-03-27 22:03:10 发布
最新推荐文章于 2023-03-27 22:03:10 发布
阅读量365 收藏 1
点赞数
文章标签： matlab程序 地震 相干噪声
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_39609457/article/details/112337904
版权

7387ee61a0a2720383569c029ad9a850.png

参考资料

G.R.Liu. 2010. Mesh Free Methods: Moving Beyond the Finite Element Method. Second Edition.

G.R.Liu. and Y.T.Gu. 2005. An Introduction to Meshfree Methods and Their Programming

张雄.2003.无网格方法

数值实现

Matlab 2019a

地球物理局 地震波动力学研究所 无网格组

# 声明
# 欢迎批评指正，禁止转载

目录

https://zhuanlan.zhihu.com/p/174383117​zhuanlan.zhihu.com
f3ab67051a866c6c11318a4cbc71a49c.png

代码：计算RPIM形函数

MLPG_main.m

%   G.R.Liu. 2010. Mesh Free Methods: Moving Beyond the Finite Element Method. Second Edition.
%   G.R.Liu. and Y.T.Gu. 2005. An Introduction to Meshfree Methods and Their Programming
%   地球物理局 地震波动力学研究所 无网格组

clear;
clc;
.......
.......
.......
[coefficient,weight] = Gauss_Coefficient(number_Gauss_points_partition,gauss);
.......
.......
.......
for node = 1:number_node
    .......
    .......
    .......
    x_quadrature_coordinates = Quadrature_Domain(x_size,y_size,x_node,y_node,x_global_boundary,x_quadrature_coordinates);       %   局部正交域
    .......
    .......
    .......    
    for iix = 1:number_partitions_x
        .......
        .......
        ....... 
        for jjy = 1:number_partitions_y
            .......
            .......
            ....... 
            gs = Domain_Gauss_Points(xcc,gauss,gs,dimension,...
                                     number_Gauss_points,number_x_quadrature_coordinates,...
                                     number_Gauss_points_partition,number_Gauss_points);
            .......
            .......
            ....... 
            for ie = 1:number_Gauss_points
                .......
                .......
                ....... 
                number_node_Shape = Support_Domain(numnode,nx,x_position,x,ds,number_node_support_domain,number_node_Shape);
                .......
                .......
                .......
                [weight,weight_xx,weight_yy] = Test_Function(domain_size,x_weight_center,x_position);

                phi = RPIM_ShapeFunction_2D(gpos,x,nv,phi,...
                                           nx,numnode,ndex,...
                                           alfc,dc,q,nRBF,mbasis);
                .......
                .......
                .......
            end 
            .......
            .......
            ....... 
        end
        .......
        .......
        ....... 
    end
    .......
    .......
    .......
end
......
......
......

关于径向基点插值形函数原理，参见：

石中居士：无网格法理论与Matlab程序设计（6）——传统径向基点插值（RPIM）形函数​zhuanlan.zhihu.com

RPIM_ShapeFunction_2D.m

%   G.R.Liu. 2010. Mesh Free Methods: Moving Beyond the Finite Element Method. Second Edition.
%   G.R.Liu. and Y.T.Gu. 2005. An Introduction to Meshfree Methods and Their Programming
%   地球物理局 地震波动力学研究所 无网格组

function phi = RPIM_ShapeFunction_2D(gpos,x,nv,nx,ndex,alfc,dc,q,nRBF,mbasis)
    rk = zeros(1,ndex + mbasis);
    phi = zeros(10,ndex);
    xv = zeros(nx,ndex);
    rr = zeros(10,ndex + mbasis);
    a = zeros(ndex + mbasis,ndex + mbasis);
    g0 = zeros(ndex + mbasis,ndex + mbasis);

    if nRBF == 1
        rc = alfc*dc;       %   MQ
    end
    if nRBF == 2
        q = alfc/(dc^2);    %   EXP
    end
    ep = 1e-20;
    mg = ndex + mbasis;
    g0(1:mg,1:mg) = 0;
    for i = 1:ndex 
        nn = nv(i);         %   支持域
        xv(1,i) = x(1,nn);  %   把局部支持域中场节点坐标保存到X中
        xv(2,i) = x(2,nn);
    end
    %%%%%%%%%%% 组装G矩阵 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    for i = 1:ndex
        nn = nv(i);
        rr = Compute_RadialBasis(x(1,nn),x(2,nn),xv,rr,ndex,rc,q,nRBF,mbasis);
        g0(i,1:ndex) = rr(1,1:ndex);
        if mbasis >= 0
            g0(i,ndex+1) = 1;
            g0(i,ndex+2) = x(1,nn);
            g0(i,ndex+3) = x(2,nn);
            g0(ndex+1,i) = 1;
            g0(ndex+2,i) = x(1,nn);
            g0(ndex+3,i) = x(2,nn);
        end
    end
    %%%%%%%%%% 求解线性方程以获得形函数 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    a(1:mg,1:mg) = g0(1:mg,1:mg);
    rr = Compute_RadialBasis(gpos(1),gpos(2),xv,rr,ndex,rc,q,nRBF,mbasis);
    rk(1:mg) = rr(1,1:mg);
    [b,kwji] = GaussEqSolver_Sym(mg,a,rk,ep);    
    if kwji == 1
        fprintf('Fail...');
        pause;
    end
    phi(1,1:ndex) = rk(1:ndex);
    %%%%%%%%%% 求解线性方程以获得dphidx %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    a(1:mg,1:mg) = g0(1:mg,1:mg);
    rk(1:mg) = rr(2,1:mg);
    [b,kwji] = GaussEqSolver_Sym(mg,a,rk,ep);
    phi(2,1:ndex) = rk(1:ndex);
    %%%%%%%%%% 求解线性方程以获得dphidy %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    a(1:mg,1:mg) = g0(1:mg,1:mg);
    rk(1:mg) = rr(3,1:mg);
    [b,kwji] = GaussEqSolver_Sym(mg,a,rk,ep);
    phi(3,1:ndex) = rk(1:ndex);
    %%%%%%%%%% 求解线性方程以获得dphidxx %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    a(1:mg,1:mg) = g0(1:mg,1:mg);
    rk(1:mg) = rr(4,1:mg);
    [b,kwji] = GaussEqSolver_Sym(mg,a,rk,ep);
    phi(4,1:ndex) = rk(1:ndex);
    %%%%%%%%%% 求解线性方程以获得dphidxy %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    a(1:mg,1:mg) = g0(1:mg,1:mg);
    rk(1:mg) = rr(5,1:mg);
    [b,kwji] = GaussEqSolver_Sym(mg,a,rk,ep);
    phi(5,1:ndex) = rk(1:ndex);
    %%%%%%%%%% 求解线性方程以获得dphidyy %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    a(1:mg,1:mg) = g0(1:mg,1:mg);
    rk(1:mg) = rr(6,1:mg);
    [b,kwji] = GaussEqSolver_Sym(mg,a,rk,ep);
    phi(6,1:ndex) = rk(1:ndex);
end

Compute_RadialBasis.m

%   参考《无网格法理论及程序设计》编制，源程序为Fortran语言
%   地球物理局 地震波场模拟实验室 无网格组
%   地球物理局 基建处 数值计算科
function OUTPUT = Compute_RadialBasis(x,y,xv,rk,ndex,R,q,nRBF,mbasis)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   加入线性多项式后计算径向基
%   输入：x,y,xv[],ndex,r,q,nRBF,mbasis
%    nRBF = 1;   MQ; 2 EXP 3 TSP
%   输出：rk[10,ndex+mbasis]
%       从1到10表示
%       r,drx,dry,drdxx,drdxy,drdyy,
%       drdxxx,drdxxy,drdxyy,drdyyy      
    rr2 = sqrt((x - xv(1,1:ndex)).^2 + (y - xv(2,1:ndex)).^2);
    if nRBF == 1    %   MQ
        rk(1,1:ndex) = (rr2.^2 + R.^2).^q;
        rk(2,1:ndex) = 2.*q.*(rr2 + R.^2).^(q-1).*(x - xv(1,1:ndex));
        rk(3,1:ndex) = 2.*q.*(rr2 + R.^2).^(q-1).*(y - xv(2,1:ndex));
        rk(4,1:ndex) = 2.*q.*(rr2 + R.^2).^(q-1) + 4.*(q-1).*q.*(x - xv(1,1:ndex)).^2.*(rr2 + R.^2).^(q-2);
        rk(5,1:ndex) = 4.*(q-1).*q.*(x - xv(1,1:ndex)).*(y - xv(2,1:ndex)).*(rr2 + R.^2).^(q-2);
        rk(6,1:ndex) = 2.*q.*(rr2 + R^2).^(q-1) + 4.*q.*(q-1).*(y - xv(2,1:ndex)).^2.*(rr2 + R.^2).^(q-2);
    end
    if nRBF == 2    %   EXP
        rk(1,1:ndex) = exp(-q*rr2^2);
        rk(2,1:ndex) = -2*q*exp(-q*rr2)*(x-xv(1,1:ndex));
        rk(3,1:ndex) = -2*q*exp(-q*rr2)*(y-xv(2,1:ndex));
        rk(4,1:ndex) = -2*q*exp(-q*rr2) + 4*q^2*(x-xv(1,1:ndex)).^2*exp(-q*rr2);
        rk(5,1:ndex) = 4*q^2*exp(-q*rr2)*(y-xv(2,1:ndex)).*(x-xv(1,1:ndex));
        rk(6,1:ndex) = -2*q*exp(-q*rr2) + 4*q^2*(y-xv(2,1:ndex)).^2*exp(-q*rr2);        
    end
    if nRBF == 3    %   TSP
        rk(1,1:ndex) = (rr2)^(0.5*q);
        rk(2,1:ndex) = q*(x-xv(1,1:ndex))*(rr2)^(0.5*q-1);
        rk(3,1:ndex) = q*(x-xv(2,1:ndex))*(rr2)^(0.5*q-1);
        rk(4,1:ndex) = q*ri2^(0.5*q-1)+2*q*(0.5*q-1)*(x-xv(1,1:ndex)).^2*(rr2)^(0.5*q-2);
        rk(5,1:ndex) = q*(0.5*q-1)*(x-xv(1,1:ndex)).*(y-xv(2,1:ndex))*(rr2)^(0.5*q-2);
        rk(6,1:ndex) = q*ri2^(0.5*q-1)+2*q*(0.5*q-1)*(y-xv(2,1:ndex)).^2*(rr2)^(0.5*q-2);        
    end
    if mbasis >= 0
        rk(1,ndex + 1) = 1;
        rk(1,ndex + 2) = x;
        rk(1,ndex + 3) = y;
        rk(2,ndex + 2) = 1;
        rk(3,ndex + 3) = 1;
    end
end

GaussEqSolver_Sym.m

function [b,kwji] = GaussEqSolver_Sym(n,a,b,ep)  %   ,kwji
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   采用高斯消元法求解对称线性方程ax=b
%   如果kwji=1，则无解；如果kwji=0，有解
%   输入：n,ma,a(ma,n),b(n),ep
%   输出：b,kwji
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   n：线性方程数
%   ma：矩阵A行数
%   A(ma,n)：Ax=B的系数矩阵
%   B(n)：输入时为Ax=B的右边向量，输出时为结果向量
%   ep：指定误差
%   kwji：控制常数，当解唯一时，kwji = 0，否则kwji = 1。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
m = zeros(1,n+1);
m(1:n) = (1:n);
for k = 1:n  % 120
    p = 0;
    if abs(a(k:n,k:n))>=abs(p)
        p = a(k:n,k:n);
        io(k:n) = (k:n);
        jo(k:n) = (k:n);
    end
    if (abs(p)-ep) < 0 
        kwji = 1;                               %   30
        return
    elseif (abs(p)-ep) == 0 
        kwji = 1;
        return
    elseif (abs(p)-ep) > 0                      %   35
        if (jo == k)
            if(io == k)
                p = 1./p;
                in = n-1;
                if(k == n)
                    b(k) = b(k)*p;
                    if(k == n)
                        for i1 = 2:n
                            i = n + 1 - i1;
                            for j = i:in
                                b(i) = b(i) - a(i,j+1)*b(j+1);
                            end
                        end                       
                        i(1:n) = m(1:n);
                        a(1,i(1:n)) = b(1:n);
                        b(1:n) = a(1,1:n);
                        kwji = 0;
                    else
                        a(k+1:in+1,k+1:in+1) = a(k+1:in+1,k+1:in+1) - a(k+1:in+1,k).*a(k,k+1:in+1);
                        b(k+1:in+1) = b(k+1:in+1) - a(k+1:in+1,k).*b(k);
                        for i1 = 2:n
                            i = n + 1 - i1;
                            for j = i:in
                                b(i) = b(i) - a(i,j+1)*b(j+1);
                            end
                        end
                        i(1:n) = m(1:n);
                        a(1,i(1:n)) = b(1:n);
                        b(1:n) = a(1,1:n);
                        kwji = 0;     
                    end
                else
                    a(k,k+1,in+1) = a(k,k+1,in+1)*p;
                    b(k) = b(k)*p;
                    a(k+1:in+1,k+1:in+1) = a(k+1:in+1,k+1:in+1) - a(k+1:in+1,k).*a(k,k+1:in+1);
                    b(k+1:in+1) = b(k+1:in+1) - a(k+1:in+1,k).*b(k);
                    for i1 = 2:n
                        i = n + 1 - i1;
                        for j = i:in
                            b(i) = b(i) - a(i,j+1)*b(j+1);
                        end
                    end    
                    i(1:n) = m(1:n);
                    a(1,i(1:n)) = b(1:n);  
                    b(1:n) = a(1,1:n);
                    kwji = 0;  
                end
            else
                t = a(io(k:n),k:n);
                a(io(k:n),k:n) = a(k,k:n);
                a(k,k:n) = t;
                t = b(io(k:n));
                b(io(k:n)) = b(k);
                b(k) = t;                
                p = 1./p;
                in = n - 1;
                if(k == n)
                    b(k) = b(k)*p;
                    if(k == n)
                        for i1 = 2:n
                            i = n + 1 - i1;
                            for j = i:in
                                b(i) = b(i) - a(i,j+1)*b(j+1);
                            end
                        end                        
                        i(1:n) = m(1:n);
                        a(1,i(1:n)) = b(1:n);                        
                        b(1:n) = a(1,1:n);
                        kwji = 0;
                    else
                        a(k+1:in+1,k+1:in+1) = a(k+1:in+1,k+1:in+1) - a(k+1:in+1,k).*a(k,k+1:in+1);
                        b(k+1:in+1) = b(k+1:in+1) - a(k+1:in+1,k).*b(k);
                        for i1 = 2:n
                            i = n + 1 - i1;
                            for j = i:in
                                b(i) = b(i) - a(i,j+1)*b(j+1);
                            end
                        end                        
                        i(1:n) = m(1:n);
                        a(1,i(1:n)) = b(1:n);                        
                        b(1:n) = a(1,1:n);
                        kwji = 0;     
                    end                    
                else
                    a(k,k+1,in+1) = a(k,k+1,in+1)*p;
                    b(k) = b(k)*p;
                    a(k+1:in+1,k+1:in+1) = a(k+1:in+1,k+1:in+1) - a(k+1:in+1,k).*a(k,k+1:in+1);
                    b(k+1:in+1) = b(k+1:in+1) - a(k+1:in+1,k).*b(k);
                    for i1 = 2:n
                        i = n + 1 - i1;
                        for j = i:in
                            b(i) = b(i) - a(i,j+1)*b(j+1);
                        end
                    end                        
                    i(1:n) = m(1:n);
                    a(1,i(1:n)) = b(1:n);                       
                    b(1:n) = a(1,1:n);
                    kwji = 0;  
                end
            end
        else
            t = a(1:n,jo(k:n));
            a(1:n,jo(k:n)) = a(1:n,k);
            a(1:n,k) = t;
            t = b(io(k:n));
            b(io(k:n)) = b(k);
            b(k) = t;            
            p = 1./p;
            in = n - 1;
            if(k == n)
                b(k) = b(k)*p;
                if(k == n)
                    for i1 = 2:n
                        i = n + 1 - i1;
                        for j = i:in
                            b(i) = b(i) - a(i,j+1)*b(j+1);
                        end
                    end   
                    i(1:n) = m(1:n);
                    a(1,i(1:n)) = b(1:n);                        
                    b(1:n) = a(1,1:n);
                    kwji = 0;
                else
                    a(k+1:in+1,k+1:in+1) = a(k+1:in+1,k+1:in+1) - a(k+1:in+1,k).*a(k,k+1:in+1);
                    b(k+1:in+1) = b(k+1:in+1) - a(k+1:in+1,k).*b(k);
                    for i1 = 2:n
                        i = n + 1 - i1;
                        for j = i:in
                            b(i) = b(i) - a(i,j+1)*b(j+1);
                        end
                    end  
                    i(1:n) = m(1:n);
                    a(1,i(1:n)) = b(1:n);   
                    b(1:n) = a(1,1:n);
                    kwji = 0;     
                end   
            else
                a(k,k+1,in+1) = a(k,k+1,in+1)*p;
                b(k) = b(k)*p;
                a(k+1:in+1,k+1:in+1) = a(k+1:in+1,k+1:in+1) - a(k+1:in+1,k).*a(k,k+1:in+1);
                b(k+1:in+1) = b(k+1:in+1) - a(k+1:in+1,k).*b(k);
                for i1 = 2:n
                    i = n + 1 - i1;
                    for j = i:in
                        b(i) = b(i) - a(i,j+1)*b(j+1);
                    end
                end        
                i(1:n) = m(1:n);
                a(1,i(1:n)) = b(1:n);       
                b(1:n) = a(1,1:n);
                kwji = 0;  
            end
        end            
    end
end

封面：[1]

参考

  1. ^https://timgsa.baidu.com/timg?image&quality=80&size=b9999_10000&sec=1598346094728&di=074bc7df2db9db0160d81c9596fa61a9&imgtype=0&src=http%3A%2F%2Fb3-q.mafengwo.net%2Fs6%2FM00%2F61%2F7F%2FwKgB4lLt3niAf8BbAAiI4kTHjvk58.jpeg
weixin_39609457
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
matlab 地震正演程序,seismic-forward 地球物理勘探中，基础的地震波正演模拟程序，包括五个 （高阶中心差分 matlab 266万源代码下载- www.pudn.com...
weixin_42346693的博客
03-20 2052
文件名称: seismic-forward下载 收藏√ [5 4 3 2 1]开发工具: Visual C++文件大小: 2985 KB上传时间: 2014-08-06下载次数: 66提 供 者: PYong详细说明：地球物理勘探中，基础的地震波正演模拟程序，包括五个正演程序，(高阶中心差分，高阶交错网格差分，使用标准的PML边界的正演程序)在acou文件夹下附有运行的结果。-Ge...
有限单元法原理与应用第四版pdf_无网格法和有限元法的区别
weixin_39805195的博客
12-14 1072
作为一种数值计算方法，无网格法和有限元一样，都是通过变分原理推导控制方程，但也存在很多不同之处。01形函数不同无网格法将研究对象离散为节点，节点之间无需网格联系，节点的函数值是该节点影响域内节点的函数值通过最小二乘拟合或积分变换得到的，即形函数不是插值函数。而有限元法将研究对象离散为单元，以单元为研究对象，单元内任意点的函数值通过节点值插值得到，形函数为插值函数。图1所示为有限元和无网格...
matlab 地震波处理,psnsps_mig.m 源代码在线查看 - 基于matlab平台,实现地震波场延拓并偏移归位源码.地震波偏移是地震资料处理中重要步骤之一,利用此程序可以模拟地震波传播过程,...
weixin_31481495的博客
03-18 1386
function [seismig,zmig,xmig]=psnsps_mig(aryin,t,x,velmod,dz,zmax,params)% PSNSPS_MIG: Exploding reflector depth migration by NSPS%% [seismig,zmig,xmig]=psnsps_mig(aryin,t,x,velmod,dz,zmax,params)%% ps...
基于高分辨率时频分析的单通道地震数据自动噪声衰减方法（Matlab代码实现）
weixin_66436111的博客
12-16 456
使用一般交叉验证方法来估计阈值水平，该方法不依赖关于噪声水平的任何先验知识。通过增加基于峰度测量的预处理步骤和基于自适应硬阈值的后处理步骤，提高了阈值处理的效率。所提出的算法可以衰减噪声（白色或彩色）并保持信号，或者去除信号并保持噪声。因此，它既可以用于正常的去噪应用，也可以用于环境噪声研究中的预处理。[1]王姣,李振春,王德营.基于CEEMD的地震数据小波阈值去噪方法研究[J].石油物探,2014,53(02):164-172.这是GCV方法的一种反向算法，可以去除信号并保持背景噪声。
深度学习从入门到精通——基于深度学习的地震数据去噪处理
swpucwf的博客
05-04 7458
传统机器学习 SVM,boosting,bahhing,knn 深度学习 CNN（典型）,GAN 地震应用方向 叠前地震数据随机噪声去除，实现噪声分离 ![didi](https://img-blog.csdnimg.cn/869c3c24ed424fb283928c14bb3d48ec.png) 面波去噪 面波作为很强的干扰波出现在地震勘探中，大大降低了地震记录的分 辨率和信噪比。深度学习作为一种数据驱动类方法, 能够从大量数据样本中 学习得到有效信号与噪声的区别, 自适应建立深度神经..
电磁场计算中的无网格局部Petrov-Galerkin方法 (2005年)
05-19
无网格局部Petrov-Galerkin方法（MLPG）建立在局部积分弱形式的基础之上，不依赖于背景网格，而只是在所考虑点的局部区域及其边界上积分，灵活方便，有利于并行计算的实现，且容易推广到非线性与非均匀介质的问题。...
瞬态热传导问题的无网格局部 Petrov-Galerkin法 (2007年)
04-26
基于加权余量法和采用移动最小二乘近似函数作为试函数,研究了无网格局部 Petrov-Galerkin方法在瞬态热传导问题中的应用。阐述了该方法应用于瞬态热传导问题的过程和基本理论,并编制了相应的计算程序进行计算。算例...
论文研究 - 使用拉普拉斯变换的无网格局部Petrov-Galerkin方法研究弹性动力学问题
05-30
拉普拉斯变换的无网格局部Petrov-Galerkin（MLPG）用于求解偏微分方程。 从动态偏微分方程中使用加权残差法局部展开局部弱形式，并使用移动最小二乘（MLS）方法构造形状函数。 该方法比有限元方法更有效地替代了...
三维弹性静力问题的无网格局部Petrov-Galerkin法* (2013年)
05-11
将二维平面问题的无网格局部Petrov-Galerkin法拓展到三维的相应理论中，编制了该法相应的三维Fortran程序．分析了均匀受拉立方体和悬臂梁两个经典算例，将所得结果与有限元法和解析解对比．结果表明了无网格局部...
matlab开发-相干函数快速和单纯形
08-23
matlab开发-相干函数快速和单纯形。相干性跟踪算法的matlab实现
matlab频谱分析代码-dt-wavelet:小波交叉频谱分析在时频域中通过地震干涉法测量行进时移的代码
05-24
matlab频谱分析代码小波方法测量地震行时变化的代码 联系人︰毛淑娟（）和奥丽莲·莫德雷（AurélienMordret）（） 该软件包包含用于使用小波互谱分析在时频域中测量地震传播时间偏移的代码和测试数据。 运行代码需要MATLAB R2018a（或更高版本）和MATLAB WAVELET TOOLBOX。 此包装中的内容： -My_Wxspectrum_TO.m：通过小波交叉谱分析在时频域中计算dt的核心功能。 -main_TO.m：在合成数据上使用My_Wxspectrum_TO.m的示例。 一键即可绘制图。 ---synthetic_dvov_0.05percent.mat：两个用于测试代码的合成波形。 合成地震图是使用速度模型通过均质背景叠加随机异质性生成的。 当前和参考速度模型之间的扰动在整个介质中为0.05％均​​匀dv / v。 （如果有兴趣，请参阅以下参考资料中的第3.1节以获取更多详细信息。） 参考文献：S.Mao，A.Mordret，M.Campillo，H.Fang，RDvan der Hilst，（2019年），《利用小波交叉频谱分析测量时频域中的地震传播时
MATLAB处理地震波的程序
07-01
一款处理地震波的程序，将多行多列转成单列数据; 因本人以前传过此程序，可是因程序中的地震波文件GilroyNo1EW.out没上传，导致各位网友下载程序应用出现错误，故此次重新上传程序及地震波文件。
EFG3_SPH_无网格法_无网格_
09-29
无网格法大致可分成两类：一类是以Lagrange方法为基础的粒子法（Particle method），如光滑粒子流体动力学（Smoothed particle hydrodynamics，简称SPH）法，和在其基础上发展的运动粒子半隐式（Moving－particle semi－implicit，简称MPS）法等；另一类是以Euler方法为基础的无格子法（Gridless methods），如无格子Euler/N—S算法（Gridless Euler/Navier-Stokes solution algorithm）和无单元Galerkin法（Element free Galerkin，简称EFG）等。
r vector 4 elements_无网格法理论与Matlab程序设计（4）——加权最小二乘近似（WLS）...
weixin_39557199的博客
11-22 227
参考资料G.R.Liu Y.T.GU著 王建明 周学军译 《无网格法理论及程序设计》数值实现Matlab 2019a前情回顾形式主义的居士：无网格法理论与Matlab程序设计（1）——概述​zhuanlan.zhihu.com地球物理局 地震波场模拟实验室 无网格组 地球物理局 基建处 数值计算科 声明： # 系列写作内容首先符合本人的研究需要，不会优先照顾读者体验。 # 仅供学习和参考，禁止...
毫米波雷达的系统设计细节(2) - 关于目标RCS的问题
xhblair的博客
03-27 7808
本博文主要探讨目标RCS的问题，给出了关于目标RCS的定义，分析了其影响因素，并理清楚了几个雷达领域常见的专有名词及其意义。随后给出了目标RCS的测量方法，并设计实验，分别从雷达方程和相对定标两种途径求解待测目标的RCS值.
RBF（径向基）神经网络
《好好先生》专栏
10-16 3261
RBF（径向基）神经网络 只要模型是一层一层的，并使用AD/BP算法，就能称作BP神经网络。RBF神经网络是其中一个特例。本文主要包括以下内容： 什么是径向基函数 RBF神经网络 RBF神经网络的学习问题 RBF神经网络与BP神经网络的区别 RBF神经网络与SVM的区别 为什么高斯核函数就是映射到高维区间 前馈网络、递归网络和反馈网络 完全内插法 什么是径向基函数 1985年，...
matlab 小波变换_Matlab实现小波相干绘图并计算AWC&PASC
weixin_39562185的博客
11-28 1669
简介今天的文章简介很少，就一句话这么简单。文中如有不妥的地方，还请各路大佬多多批评指正。应该有伙伴做过小波分析方面的研究，比如实部方差不同尺度分解，这一部分可以用简单的工具箱功能实现，如图wavemenu(新版本有所改动)命令就可以打开工具箱开始做你需要的工作。在做小波相干之前，有必要先做连续小波变换进行时间序列数据的周期检测，可以用过工具箱，也可以通过前辈编写好的脚本执行(感谢)，把下...
matlab偏微分方程数值解误差_无网格法与Matlab程序设计（9）——无网格径向基插值法：原理...
weixin_39587029的博客
11-24 495
参考资料G.R.Liu Y.T.GU著 王建明 周学军译 《无网格法理论及程序设计》数值实现Matlab 2019a前情回顾形式主义的居士：无网格法理论与Matlab程序设计（1）——概述​zhuanlan.zhihu.com地球物理局 地震波场模拟实验室 无网格组 地球物理局 基建处 数值计算科 声明： # 系列写作内容首先符合本人的研究需要，不会优先照顾读者体验。 # 仅供学习和参考，禁止...
fluent中NOx模型怎么选择
最新发布
05-20
在 Fluent 中，选择 NOx 模型通常需要考虑以下因素： 1. 模型适用范围：不同的 NOx 模型适用于不同的燃烧条件和燃料类型。例如，Petrov-Galerkin 模型适用于低温燃烧，而 Zeldovich 模型适用于高温燃烧。 2. 计算效率：一些 NOx 模型计算比较复杂，需要更多的计算资源，而另一些模型则计算相对简单，计算速度更快。 3. 数据可用性：有些 NOx 模型需要输入更多的化学反应机理数据，而有些模型则需要更少的数据。 4. 精度要求：不同的 NOx 模型有不同的精度要求。一些模型可以提供较高的精度，但也可能需要更多的计算资源。 综合考虑以上因素，您可以选择适合您特定应用的 NOx 模型。常用的 NOx 模型包括 Zeldovich 模型、Extended Zeldovich 模型、Petrov-Galerkin 模型等。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值