当前位置:网站首页>[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
边栏推荐
- Maximum number of "balloons"
- “薪资倒挂”、“毕业生平替” 这些现象说明测试行业已经...
- Les phénomènes de « salaire inversé » et de « remplacement des diplômés » indiquent que l'industrie des tests a...
- 2022.07.03(LC_6111_统计放置房子的方式数)
- Playwright之录制
- 华为200万年薪聘请数据治理专家!背后的千亿市场值得关注
- Paxos 入门
- Summary of the function and usage of const, volatile and restrict
- 潘多拉 IOT 开发板学习(RT-Thread)—— 实验4 蜂鸣器+马达实验【按键外部中断】(学习笔记)
- skimage: imread & imsave & imshow
猜你喜欢
【Unity】InputSystem
uniapp微信小程序拿来即用的瀑布流布局demo2(方法二)(复制粘贴即可使用,无需做其他处理)
IT转测试岗,从迷茫到坚定我究竟付出了什么?
Pycharm professional download and installation tutorial
【selenium自动化】常用注解
Relationship between classes and objects
Postman automatically fills headers
Identifiers and keywords
Huawei employs millions of data governance experts! The 100 billion market behind it deserves attention
107. SAP UI5 OverflowToolbar 容器控件以及 resize 事件处理的一些细节介绍
随机推荐
潘多拉 IOT 开发板学习(RT-Thread)—— 实验4 蜂鸣器+马达实验【按键外部中断】(学习笔记)
Hisilicon 3559 universal platform construction: YUV422 pit stepping record
What happened to those who focused on automated testing?
Pycharm professional download and installation tutorial
【FPGA教程案例10】基于Verilog的复数乘法器设计与实现
Playwright recording
【C】 (written examination questions) pointer and array, pointer
FEG founder rox:smartdefi will be the benchmark of the entire decentralized financial market
SAP ui5 application development tutorial 106 - how to improve the readability of SAP ui5 application routing URL trial version
Basic concept and usage of redis
[error reporting] "typeerror: cannot read properties of undefined (reading 'split')“
[selenium automation] common notes
创新引领方向 华为智慧生活全场景新品齐发
skimage: imread & imsave & imshow
MySQL uses the explain tool to view the execution plan
Playwright之录制
【海浪建模3】三维随机真实海浪建模以及海浪发电机建模matlab仿真
2022.07.03 (lc_6111_counts the number of ways to place houses)
大专学历,33岁宝妈又怎样?我照样销售转测试,月入13k+
揭露测试外包公司,关于外包,你或许听到过这样的声音