Moment tensorFD

微地震震源机制正演matlab程序，欢迎学者前来下载学习

function [a1a2a3]=aat（qwpvpvsn0）
%求反演算子
%  clc;
%  clear all;
% 
% q=[0.010 0.200 0.500];
% w=[0.1 0.500 0.0100];
% p=2.8*1000;vp=3;vs=1.73200;
% 
% %st=10;d=16;sl=9;
% %st=sqrt（6）/6*[1 0 0;0 -2 0;0 0 1];
%  n0=300;
[g1g2g3]=green（wqpvpvsn0）;

mt=mtt（pvsn0）;

mt（n0+1:2*n0-1）=0;

fmt=fft（mt）;

a1=zeros（62*n0-1）;a2=zeros（62*n0-1）;a3=zeros（62*n0-1）;
 
for i=1:6
   
    
    g1（in0+1:2*n0-1）=0;
    g2（in0+1:2*n0-1）=0;
    g3（in0+1:2*n0-1）=0;

    gg1=fft（g1（i:））;
    gg2=fft（g2（i:））;
    gg3=fft（g3（i:））;
   
    a1（i:）=ifft（gg1.*fmt）;   
    a2（i:）=ifft（gg2.*fmt）; 
    a3（i:）=ifft（gg3.*fmt）; 
    

end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2019-01-10 15:11  矩张量正演\
     文件         128  2014-04-15 14:45  矩张量正演\12.mat
     文件      230112  2013-12-25 15:56  矩张量正演\71.date
     文件         685  2013-10-31 20:32  矩张量正演\Untitled1.m
     文件         678  2013-12-25 14:59  矩张量正演\aat.m
     文件         183  2013-10-29 13:12  矩张量正演\abc.dat
     文件        1725  2013-12-25 15:00  矩张量正演\arc.m
     文件         767  2013-12-25 15:30  矩张量正演\bcwj.m
     文件        3733  2013-12-25 15:01  矩张量正演\calcCircle.m
     文件          70  2014-05-24 10:50  矩张量正演\dlk.m
     文件         331  2013-08-16 18:31  矩张量正演\duwenjian.m
     文件        1076  2013-11-15 21:14  矩张量正演\fany.m
     文件        3005  2014-04-17 17:59  矩张量正演\fanyan11.m
     文件       11154  2015-04-22 17:29  矩张量正演\fmd.m
     文件         610  2014-04-08 21:05  矩张量正演\gjs.m
     文件        2172  2014-04-11 16:25  矩张量正演\gre.m
     文件        4276  2014-04-27 22:23  矩张量正演\gre1.m
     文件       12549  2014-05-24 22:01  矩张量正演\gre2.m
     文件        5078  2014-05-24 13:55  矩张量正演\gre3.m
     文件        2395  2014-06-27 11:20  矩张量正演\gre3_1.m
     文件        3794  2013-12-25 14:42  矩张量正演\green.m
     文件         618  2014-04-08 21:12  矩张量正演\greenjs.m
     目录           0  2019-01-10 15:07  矩张量正演\html\
     文件       11869  2014-04-15 15:02  矩张量正演\html\zhenyan00.html
     文件        1569  2014-04-15 12:56  矩张量正演\huatu.m
     文件        1166  2013-12-25 14:29  矩张量正演\moment.m
     文件         398  2013-12-25 15:03  矩张量正演\mtt.m
     文件        2551  2013-11-15 10:27  矩张量正演\sqyxsj.m
     文件       12907  2014-05-24 23:20  矩张量正演\v.xlsx
     文件        6704  2014-05-20 12:18  矩张量正演\vol3d.m
     文件         935  2013-12-12 11:43  矩张量正演\zheny.m
............此处省略0个文件信息