SIMPLE FDTD WAVE PROPAGATION IN MATLAB
Single-file vectorized implementations of wave propagation in MATLAB. We solve second-order wave equation in displacement formulation in time domain(FDTD). We don't account for derivatives of elastic parameters to keep it simple. Medium is surrounded by simple absorbing sponge boundaries with exponential decay (Cerjan, 1985).
DISCRETIZATION DETAILS:
- Finite-Differences in Time Domain (FDTD)
- Regular collocated grid
- Explicit time stepping
- O(2,2)
- Conventional stencils derived from Taylor series:
- in space [1: -2 :1]/dx2 and [1: -1: -1:1]/4dxdz
- in time [1: -2 :1]/dt2
HOW TO?
Open folder with a scenario of interest e.g. acoustic_2D_FDTD_wave_propagation
Open .m file in MATLAB
Press "Run" button
-
[x,fs,bits]=wavread(‘d.wav‘, [1 5000]); % sound(x, fs, bits); N = length(x); n = 0 : N-1; t = n/fs; % x = x y = fft(x, N); mag = abs(y); f = n * fs / N; subplot(211); plot(t, x); xlabel(‘Time (sec.)‘
-
%1.xiaozao函数，是需要对标准化的序列进行消除数据噪音分析； %2.Db3函数，是对数列进行Db3趋势分析； %3.period函数，是求得时间序列的实部和模的平方。 %其中周期变化图是实部的等值线图 %而小波方差是模的平方的算数平均。 s=load('11.txt'); % input SST time series s=zscore(s); scales=[1:1:32]; %进行连续
-
matlab中的小波分析工具箱(wavelet 30) 第七章 小波分析、分析工具及应用发展 复习回顾 小波分析来源 小波分析来源于信号分析的需求.设一个有限分辨率的连续信号 ,将其近似地表示为下列阶梯函数(图1) 为简化叙述,取整数点(n)为样点,式中 为样本值,而其基函数 并又将其称为“尺度函数”，如图2所示. 我们将采样间隔加倍,则其样点数减半,这时信号表示为 显然,这里自然取 ,参见(图3