周末无聊瞎搞的,没有什么物理依据
代码
% 水波瞎搞模拟
clear; close all; clc
% [网格范围]
range = 10;
gap = 0.1;
[x,y] = meshgrid(-range:gap:range);
% [波纹参数]
lamda = 1; % 参数可调
phi = 0*pi; % 参数可调
% [迭代绘制]
dis = sqrt(x.^2+y.^2);
for t = 0:0.02:4 % 秒
% 正弦函数参数
v = 1-t/5; % 参数可调
scale = v*t./dis;
scale(scale>1) = 1;
A = exp(dis/100)*exp(-t).*scale;
% 正弦函数
z = A.*sin(2*pi*(v*t-dis/lamda)+phi);
% 绘制
mesh(x,y,z) ,axis([-range,range,-range,range,-2,2])
drawnow
end