基于Matlab模拟随机介质模型

✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，matlab项目合作可私信。

🍎个人主页：Matlab科研工作室

🍊个人信条：格物致知。

更多Matlab仿真内容点击👇

智能优化算法       神经网络预测       雷达通信       无线传感器        电力系统

信号处理              图像处理               路径规划       元胞自动机        无人机 

⛄ 内容介绍

基于Matlab模拟随机介质模型

⛄ 代码

clc

clear all

close all

vp0=4000;%背景速度与密度

vs0=2800;

rau0=2690;

K=0.5;%一般取0.3-0.8之间

eps=0.08;%标准差

A=50;%介质在x，z方向上的自相关长度

B=100;

x=1:10:2560;

z=1:10:2560;

%a=zeros(256,256);

b=x'*x/(A^2)+z'*z/(B^2);

c=exp(-sqrt(b));%自相关函数

%c=exp(-b);

d=fft2(c,256,256);%产生随机过程功率谱

e=unifrnd(0,2*pi,256,256);%产生二维随机场

%normplot(e);

f=sqrt(d).*exp(-i*e);%产生随机功率谱

g=ifft2(f,256,256);%得到空间域的随机扰动

h=real(g);

u=mean2(h(:));%均值

l=var(h(:));%方差

m=eps/sqrt(l)*(h-u);

vp=vp0*(1+m);

vs=vs0*(1+m);

rau=rau0*(1+K*m);

imagesc(vp);

fvp = fopen('random_vp.bin','wb');

fvs = fopen('random_vs.bin','wb');

frau = fopen('random_rau.bin','wb');

for x=1:256

   for z=1:256

        fwrite(fvp,vp(x,z),'float'); 

        fwrite(fvs,vs(x,z),'float');

        fwrite(frau,rau(x,z),'float');

   end

end

fclose(fvp);

fclose(fvs);

fclose(frau);

⛄ 运行结果

​⛳️ 代码获取关注我

❤️部分理论引用网络文献，若有侵权联系博主删除

❤️ 关注我领取海量matlab电子书和数学建模资料