当前位置:网站首页>[wave modeling 1] theoretical analysis and MATLAB simulation of wave modeling
[wave modeling 1] theoretical analysis and MATLAB simulation of wave modeling
2022-07-05 00:53:00 【FPGA and MATLAB】
1. Software version
matlab2017b
2. System principle
The design of ocean wave modeling mainly involves the determination of ocean wave model , Calculation of various parameters of ocean waves , This system ,
The wave model is as follows :
In reality , There are many kinds of wave models , This model is used here , Mainly because the effect of this model is the closest to the effect on the website you give .
G Is the acceleration of gravity , Take the constant 9.8
3. Core source code
global Winds; % The wind speed
global g; % Acceleration of gravity
global kk; % The size ratio between the sand table of the simulation model and the actual area
global Xmax;
global Ymax;
global Dxy;
global flag;
global VX;
global VY;
global VZ;
flag = 0;
g = 9.8; % Acceleration of gravity
kk = 1/40; % The size ratio between the sand table of the simulation model and the actual area
% Interval of simulation
Dxy = 4;
% The sea area covered by the simulation
Xmax = 1000;
Ymax = 1000;
Start = 200;
x = [Start:Dxy:Xmax];
Ymax2 = round(Ymax/2);
y = [Start:Dxy:Ymax2];
[xo,yo]= meshgrid(x,y);
z2 = zeros(size(x));
% The wave height of the wave itself
r = (3.5325*Winds^2.5)/1000;
% The wave length of the wave itself
k = 2*g/(3*Winds^2);
L = 2*pi/k;
% cycle T
T = sqrt(2*pi*L/g);
% Wave frequency
w = sqrt(2/3)*g/T;
t = 0;
while(flag == 0)
disp('the wind speed is');Winds
t = t + 1;
for i = 1:(Ymax2-Start)/Dxy+1
for j = 1:(Xmax-Start)/Dxy+1
% Attenuation coefficient
d = sqrt((xo(1,j)-0)^2 + (yo(1,j)-0)^2);
alphas = exp(-0.07*d) - 0.18;
z2(i,j) = alphas*r*cos(k*sqrt((xo(1,j)-0)^2 + (yo(1,j)-0)^2) - w*t);
end
end
% Show local effects
axes(handles.axes1);
surfl(xo,yo,z2);
axis([Start-50 Xmax+50 Start-50 Ymax2+50 -8 10]);
shading interp;
colormap([143/255,157/255,203/255]);
alpha(0.75);
lightangle(-30,90);
view([VX,VY,VZ]);
pause(0.1);
% The wave length of the wave itself
k = 2*g/(3*Winds^2);
Ls = 2*pi/k;
set(handles.edit1,'String',num2str(Ls));
% Calculate the parameter index of the wave
% The wave height of the wave itself
rs = (3.5325*Winds^2.5)/1000;
set(handles.edit3,'String',num2str(rs));
% cycle T
T = sqrt(2*pi*L/g);
% Speed
c = g*T/(2*pi);
set(handles.edit4,'String',num2str(c));
% Wave frequency
w = sqrt(2/3)*g/T;
set(handles.edit5,'String',num2str(w/2/pi));
end
4. test result
The simulation operation is shown in the above figure ,
First click on START Start simulating the effect of waves , And then click stop Can pause .
Upper figure 3 The tag is the corresponding button on the website to adjust the wind speed .
4 All parameters of the wave are displayed . See the theoretical introduction on the previous page for the calculation of specific parameters
5 Here is 3D Display of effect , You can observe the waves from different angles (MATLAB in GUI Of AXIS It is impossible to rotate the mouse directly for 3D display , So you need to display through this function )
A19-06
边栏推荐
- (script) one click deployment of any version of redis - the way to build a dream
- Check if this is null - checking if this is null
- leetcode494,474
- 大专学历,33岁宝妈又怎样?我照样销售转测试,月入13k+
- 【微处理器】基于FPGA的微处理器VHDL开发
- Applet live + e-commerce, if you want to be a new retail e-commerce, use it!
- Huawei employs millions of data governance experts! The 100 billion market behind it deserves attention
- Arbitrum:二维费用
- IT转测试岗,从迷茫到坚定我究竟付出了什么?
- Introduction to the gtid mode of MySQL master-slave replication
猜你喜欢
skimage: imread & imsave & imshow
2022.07.03 (LC 6108 decryption message)
npm install报错 强制安装
全栈开发提效神器——ApiFox(Postman + Swagger + Mock + JMeter)
Deux nombres se remplacent
揭露测试外包公司,关于外包,你或许听到过这样的声音
Talking about JVM 4: class loading mechanism
[untitled]
Playwright之录制
Innovation leads the direction. Huawei Smart Life launches new products in the whole scene
随机推荐
那些一门心思研究自动化测试的人,最后都怎样了?
Open3d uses GICP to register point clouds
Inventory of more than 17 typical security incidents in January 2022
P3304 [SDOI2013]直径(树的直径)
skimage: imread & imsave & imshow
Operator explanation
There is a new Post-00 exam king in the testing department. I really can't do it in my old age. I have
资深测试/开发程序员写下无bug?资历(枷锁)不要惧怕错误......
What you learned in the eleventh week
IT转测试岗,从迷茫到坚定我究竟付出了什么?
全栈开发提效神器——ApiFox(Postman + Swagger + Mock + JMeter)
Distributed base theory
2022.07.03 (lc_6111_counts the number of ways to place houses)
兩個數相互替換
Complete knapsack problem (template)
6. Scala operator
[Yocto RM]10 - Images
【海浪建模2】三维海浪建模以及海浪发电机建模matlab仿真
College degree, what about 33 year old Baoma? I still sell and test, and my monthly income is 13K+
What did I pay for it transfer to testing post from confusion to firmness?