- 在线时间
- 1630 小时
- 最后登录
- 2024-1-29
- 注册时间
- 2017-5-16
- 听众数
- 82
- 收听数
- 1
- 能力
- 120 分
- 体力
- 557440 点
- 威望
- 12 点
- 阅读权限
- 255
- 积分
- 172605
- 相册
- 1
- 日志
- 0
- 记录
- 0
- 帖子
- 5313
- 主题
- 5273
- 精华
- 18
- 分享
- 0
- 好友
- 163
TA的每日心情 | 开心 2021-8-11 17:59 |
---|
签到天数: 17 天 [LV.4]偶尔看看III 网络挑战赛参赛者 网络挑战赛参赛者 - 自我介绍
- 本人女,毕业于内蒙古科技大学,担任文职专业,毕业专业英语。
 群组: 2018美赛大象算法课程 群组: 2018美赛护航培训课程 群组: 2019年 数学中国站长建 群组: 2019年数据分析师课程 群组: 2018年大象老师国赛优 |
6 v4 M- o6 u8 P
【数学建模】基于matlab三维海浪模型仿真【含Matlab源码 1159期】3 m. d# q6 r2 ^ h2 W, M |6 L
一、简介
8 U. {1 E: J% S# t基于matlab三维海浪模型仿真:海洋工程领域中的浮式结构物设计,需要精确计算海浪载荷,保证满足稳定性和安全性.该文利用MATLAB对Longuet-Higgins长峰波海浪模型和三维不规则短峰波随机海浪模型进行了仿真研究.结果表明,利用海浪谱来模拟三维随机海浪能够得到比较精确的海浪波面图和波高值,进一步根据流体的势流理论就可以分析计算出该结构物受到的海浪载荷,为校验结构物的结构强度提供了必要的基础。
5 t5 n2 ?5 B& z4 W9 j }
0 v& N2 w5 U+ Q6 ?* H8 I6 b+ b# Z0 t" X/ R0 ]4 E
二、源代码
! [7 u( e' P. e5 Oclear all;
: A0 ~( `% B8 ?: Q& w6 M( h) j! ^close all;
2 S% g4 b" W' Q8 [1 H1 o6 x3 S; Z& x3 N
( Y3 x5 m: G4 E# P2 }% ^/ mnhFig = 0; % Figure Number;
5 r' e! d$ F+ {3 h" KSeaRegLx = 40e+3; % Sea region length, unit: m, k% g; e; U7 _- m; p
SeaRegLy = 40e+3; % Sea region length, unit: m/ {% d. J8 u ]# A! l% R- N# A
% 根据频域采样定理,在确定采样周期即波束的足迹宽度之后,
, [( Q. A9 B; Z2 G! H% 可以利用时域采样的采样点数目(通过空间频谱的宽度即根据时域采样定理计算出采样点数),实现频域采样.
8 H: ]3 @; c6 W9 Y5 l+ S0 a. t% 对于海洋波谱而言,波谱宽度是非常宽的,因此其波谱宽度一般是指包括了大部分能量(或者说主要的能量部分)的不完全谱宽度7 z5 c$ D/ p% H$ D3 R
% 所以这里事先根据波谱函数计算主要能量部分对应的谱宽度
8 r3 f% G1 Y, r2 ^! M+ t( T! T4 u( g$ LxSampleStep = 50; % unit: m; |& h' z" Z' F: ?6 T. Z! K
ySampleStep = 50;7 U# Z+ i$ }3 G# X& Y# o
xNs = round( SeaRegLx / xSampleStep );; b" {0 X i' [) |5 Y+ L
yNs = round( SeaRegLy / ySampleStep );
0 W8 ~5 y- D9 a! |& D; c
. _7 J4 F: B$ T; S( Q; C( ?
$ P; ]/ S1 c/ f. i" zx = linspace( -SeaRegLx / 2, SeaRegLx / 2, xNs );$ v& ?; U( H6 i( W2 y$ j9 L" j% p5 n
y = linspace( -SeaRegLy/ 2, SeaRegLy / 2, yNs );# h; @7 [8 l, M$ \- ?% m
7 h6 z7 l( ^0 f6 {: Q
2 d1 d' T5 S; {% P+ Y7 D/ O% 海浪谱仿真" F) `( t. o" Q5 u* R7 S
g = 9.8; % gravity acceleration
4 l% o4 K' N) p& P( r6 U" K% swell wave spectrum parameter2 @; |( _6 y( y* V$ z2 P8 V6 y3 I
SwellWaveLength = 1000; % 涌浪波长
. L. P) }) p6 v$ z& y" CKswellWavePeak = 2 * pi / SwellWaveLength;
* N) ^( q. `( d7 P3 \SwellAngled=0;3 s. {# N0 S' x+ w' }" t
SwellAngle = SwellAngled/ 180 * pi; % 涌浪与观测向夹角/ E e6 D4 ^8 z) e
KxswellWavePeak = KswellWavePeak * cos( SwellAngle );
3 Z; ~" t! \' L7 D6 I. iKyswellWavePeak = KswellWavePeak * sin( SwellAngle );. D' H( H) R$ p/ H [; d7 x
SigmaKx = 2.5e-3; % 涌浪谱宽度
, K7 k' u% d7 \# K0 N5 USigmaKy = 2.5e-3;
- T. z* ?. ~, z# g1 n7 M! H3 T: B9 fSigmaHSwell = 2; %涌浪波高3 w% B; O8 ?- ?- r% V2 y
% wind wave spectrum parameter7 A9 r/ p! x/ E {$ d
WindAngle = 45 / 180 * pi; % 风向与观测向夹角
9 m C5 p' N5 u6 O5 T6 QU10 = 12; % 10米高处海面风速 % 5m/s, 10m/s, 15m/s 的风浪谱宽度分别为 1.5, 0.4, 0.15.对应的最大空间采样间隔为2m, 10m, 20m v, d) p% Z) j/ N5 J6 N/ v
KwindPeak = g / ( 1.2 * U10 )^2;
7 C7 h* H K `2 j% sea wave spectrum parameter# D- r; d I( N8 e( A i2 x
NxSeaWave = xNs; % 频域采样点数,方便傅立叶变换$ c3 D3 v$ W7 d: Z6 r( q
NySeaWave = yNs; 5 n9 N8 u2 i5 \1 |$ E
KxSeaWave = 2 * pi / SeaRegLx * ( -xNs / 2 : 1 : xNs / 2 - 1 );0 [+ E$ V1 O0 {0 J3 X" g) j5 f
KySeaWave = 2 * pi / SeaRegLy * ( -yNs / 2 : 1 : yNs / 2 - 1 );
/ C) J( n. R/ v: KKxSeaWaveTicks = KxSeaWave;
. M6 w0 v) C3 jKySeaWaveTicks = KySeaWave;" M/ k' \6 F! y+ n' }
KxSeaWave = ( KxSeaWave == 0 ) * ( max( KxSeaWave ) * 1e-16 ) + KxSeaWave; % avoid divided by 00 E3 b2 z, N6 Z1 R" W! u: E
KySeaWave = ( KySeaWave == 0 ) * ( max( KySeaWave ) * 1e-16 ) + KySeaWave; % avoid divided by 0; u- y0 E* M0 K
: r* F+ [% h# ? _" N& K0 @
" a) ] c$ U, n
% swell wave spectrum4 r' A* n6 W1 ^8 ]8 W% [+ M5 W' c
SpectrumSwell = zeros( NySeaWave, NxSeaWave );: [9 e$ t H. |0 a2 ~) `9 T- R
Temp1 = ones( NySeaWave, 1 ) * ( KxSeaWave - KxswellWavePeak ) / SigmaKx;. D0 [5 v) i7 u& \
Temp2 = ( KySeaWave - KyswellWavePeak )' / SigmaKy * ones( 1, NxSeaWave );
2 ~5 ?" n# i0 U5 T* v! o3 VSpectrumSwell = SigmaHSwell^2 / 2 / pi / SigmaKx / SigmaKy * exp( -0.5 * ( Temp1.^2 + Temp2.^2 ) );
" Q" u ^6 w! y' s, ~clear Temp1;7 Y- b$ [( }' Z, D/ d; x* a- A9 V
clear Temp2;$ }+ x6 |; b7 k+ K6 q( F' _
figure;9 y6 s x( [: p& p& E
colormap(gray(256));
/ `+ A9 t) r+ L; h, s1 Fimage( KxSeaWave, KySeaWave, 256 - 255 / ( max( max( abs( SpectrumSwell ) ) ) - min( min ( abs( SpectrumSwell ) ) ) ) * ( abs( SpectrumSwell ) - min( min ( abs( SpectrumSwell ) ) ) ) );
" c4 B v# |" _- Y- |axis('xy');& N: T: _- h3 Y: Q+ s
xlabel( 'kx:X方向波数');
) y9 Y ]! |4 l& }4 `+ gylabel( 'ky:Y方向波数');0 Z5 i l* E; p% ]) f
title( '涌浪谱');
) ?0 M, j, i9 X: o) Y1 `
- h- b( z% A/ t/ z) }2 z3 A! O) |$ A: i# n/ M
KxSeaWaveMatrix = ones( NySeaWave, 1 ) * KxSeaWave;* {/ [! E @ u% D9 y! Z
KySeaWaveMatrix = KySeaWave' * ones( 1, NxSeaWave );( O" m% k2 K u8 g6 y; r( t/ c
KSeaWaveTemp1 = ( sqrt( KxSeaWaveMatrix.^2 + KySeaWaveMatrix.^2 ) );0 m9 a' u; I$ {2 n6 W
Fwk1 = exp( - 1.22 * ( sqrt( KSeaWaveTemp1 ./ KwindPeak ) -1 ).^2 );
6 N" C8 W7 L6 u1 ^HKKpx1 = 1.24 * ( ( KSeaWaveTemp1 / KwindPeak < 0.31 ) & ( KSeaWaveTemp1 / KwindPeak >= 0 ) );* o7 f8 h- ]6 \/ y# K( {( [* B
HKKpx2 = 2.61 * ( ( KSeaWaveTemp1 / KwindPeak ).^0.65 ) .* ( ( KSeaWaveTemp1 / KwindPeak < 0.9 ) & ( KSeaWaveTemp1 / KwindPeak >= 0.31 ) );2 ^+ ]- d& t& G+ y4 l
HKKpx3 = 2.28 * ( ( KSeaWaveTemp1 / KwindPeak ).^( -0.65 ) ) .* ( KSeaWaveTemp1 / KwindPeak >= 0.9 );
2 d4 Q n: Y* O7 T: j3 Y6 XHKKp1 = HKKpx1 + HKKpx2 + HKKpx3;
& O% k1 j+ L$ C9 uTemp1x = 1.62 * 1e-3 * U10 / ( g^0.5 ) ./ ( KSeaWaveTemp1 ).^3.5;
) [! O' [" h2 E5 X4 tTemp2x = exp( -( KwindPeak ./ KSeaWaveTemp1 ).^2 ) .* ( 1.7 .^ Fwk1 );& t! A* x, v5 a i5 R s7 d
Temp3x = ( HKKp1 .* ( sech( ( HKKp1 .* ( atan( ( KySeaWaveMatrix ./ KxSeaWaveMatrix ) ) - WindAngle ) ) ) ).^2 );2 Y& y- G) Y& R/ A' t# J
% SpectrumWind1 = ( Temp1x .* Temp2x .* Temp3x )';
, b: ~" z9 N- T2 X. E/ _& Sfid=fopen('sea_top_wl1000_wh50m.dat','w')6 _! g5 J+ Y/ C3 Q
fprintf(fid,'%10.5f\n',real(wh));
! i* v* X, L( U4 p, Tstatus=fclose(fid);" w$ q' X0 B% |9 Y$ }1 X$ h5 g6 h2 _9 I
4 q+ W$ ?' I- I
+ H) Z8 ]6 J0 e! |! xfid=fopen('sea_top_wl1000_wh50m.dat','r');1 m" g7 q, ^ R& z0 u5 D
sea_dem=fscanf(fid,'%g',[800,800]);
& h3 z) _( N2 |4 c7 R& o* c+ @4 Zfclose(fid);
( T5 b0 x% J* C( Z. |5 F7 P%image( x, y, 256 - 255 / ( max( max( real( sea_dem ) ) ) - min( min ( real( sea_dem ) ) ) ) * ( real( sea_dem ) - min( min ( real( sea_dem ) ) ) ) );6 G2 U$ c g5 L! J- T+ g' p) n6 h
% y7 t# A7 z8 x
mesh(sea_dem)2 M1 I0 q& a' M& Z
colorbar0 m2 }0 D* I7 ~
% view(0,90)
: O0 d) I% k* ~9 O( z%colormap(gray)' u4 J- b" U6 |) e
% This function computes the non-coherent integration improvment' L: u$ D4 o% I v8 [3 d8 |
% factor using the empirical formula defind in Eq. (2.48)# x, w6 E* N5 u
fact1 = 1.0 + log10( 1.0 / pfa) / 46.6;6 I- H x5 ?7 |) G
fact2 = 6.79 * (1.0 + 0.235 * pd);
1 f& r X/ s9 Bfact3 = 1.0 - 0.14 * log10(np) + 0.0183 * (log10(np))^2;
$ B7 Y* R8 ~5 j, `+ [6 U* c. kimpr_of_np = fact1 * fact2 * fact3 * log10(np);8 \7 W% W$ Y. k, N( m: J; ~
return
, q7 i3 |' j4 o8 ]% }3 T! k' V; x! D' J! \5 L6 m' l
, _- v/ }- ?% u9 Y
% SpectrumWind = SpectrumWind1';
1 K9 Z+ }3 w0 x* p# N1 `SpectrumWind = Temp1x .* Temp2x .* Temp3x;8 J' c. t0 g0 ]5 _8 W8 a
clear KSeaWaveTemp1 Fwk1 HKKpx1 HKKpx2 HKKpx3 HKKp1 Temp1x Temp2x Temp3x;
# J3 D' K7 D) w2 o" j' i
; X8 y2 p, {, _+ b8 w" s7 c6 s
9 |: Y3 @2 G8 G* N4 ]3 l1 s1# k% s4 e( q+ d. p X
26 Y3 G: w. e! m3 k* k. v$ h5 K
3
/ E7 g" X% `( g) I) h" s! n4
! g7 C+ `) q# N2 H. H! s50 S6 n. c2 C: g! U$ h/ b; n
6# N. O3 H6 k% S
7
4 r# t8 n( V% H8/ u& V7 H2 P2 M }2 |8 j# f
9
5 A. q }6 r* ^. ~% z104 y; {1 @( }9 x4 m! A1 S: V* @
11) I( D* K1 b" F, t V% {
125 Q. a; B: n& ^9 B) E
13
5 w) K9 S* c0 g1 K/ t14
1 T9 Z' R& ]" @% d4 `3 L% i15
( X# [" J( S; G) j# }: c, o16
( X' Q( _* T7 S% G17$ E* o. i' h9 v8 W/ _* U) z
18
0 K& H6 E- \+ m9 U2 J19
+ |1 f l3 W, i- N3 q; I20
; n! a2 E% O0 P4 d$ ^8 K21% [3 a/ W, z: C* n* H
22
( n/ e1 d0 p4 o6 g23
3 n6 A. \9 A; w+ c) W2 z24# S) n9 d3 Z; Q3 e" E( K- I x1 v
259 x$ [/ g) }8 S% s/ U; c( _( L- z6 j
26
$ ]- j7 X0 ]& Q! @2 N27
- l" i# R; k" z0 ]$ V! J28
+ ?3 s1 W6 e- H5 E3 L29
8 }' A8 W# ]8 K; V1 ^30
& l& i+ J3 M# p0 y! X$ Q31
) {' Q& m! p0 G& p32
$ r8 C, M) G# R336 L" ~" L& Q4 ]' G* p! t( c1 l
34- @( k* D3 ]* D, ?1 k( @4 }
351 Q: \! J! F" ^7 Q0 t
36
5 s( Y* B: G! k/ c- w& Z9 G7 A: x37
" g, c) W. A2 u38
# Z' d8 T7 A6 Q3 b! q1 q! D, [3 A396 Q4 m6 Y, N6 d3 n
40& X/ O/ n' ~3 V7 j7 L: x. T0 Q+ d
41
( ~! L% m% R5 G0 `6 ?422 b* i4 `$ r5 K$ G& s% K
43
6 A R, ?5 [( Y0 m( I1 j44
0 p4 _8 {( i* W+ L1 b45" b4 [! B7 U2 q1 g" J
46
5 e* x) I7 \( S5 C1 [' F47* w4 b9 b) r, X: _1 s7 D6 ?3 a' [
48& q! O3 L4 U& w" B* W
49, ~/ i3 N% I! _% \7 w( U
50
* |- X. x3 v" E) u9 E9 U, Y& i4 A51
8 A3 C8 v" T1 w52
$ h! s5 M2 X. i/ @* A' G53+ g. s4 V6 ?4 F& j
54/ l" x; D) h$ \) W. Y" Z
556 n6 K2 D9 o2 F% N( u( F4 C
56, g7 f* ~ r) E6 Q" A1 T
57
( O# }3 ]$ _4 t- D4 [' P) w& m582 p- B0 h" H- d5 u3 C9 ?- s
59
# [" W9 O1 c+ |' e4 A" `0 E60
+ o+ T: J2 f. `7 R61
5 y8 q1 m! N& n' g7 z62
$ r/ W1 j8 M0 Y3 P* B63
) H4 @8 C4 t$ x3 ^3 G* S; r2 u5 x64
0 }/ u+ \$ o& M65# l$ x2 ]& N U* G/ Q4 K3 V5 l
66' `* e, c9 j1 `% ~* K) D& T
676 |3 P6 s- I# }3 X* J6 t5 C
68
$ c: ]0 _9 P% m* ?) p! O69
& ]- [& B3 [; n707 s l1 e1 c) v4 t: t
71
0 ^* k* e& V/ k+ M) b72
8 l+ M* }, M, C, X73
5 b9 j: e% J9 x: n74
3 p& {" ~( A$ {! I/ D: N4 z" r75
* a( }* `3 r2 w! f* |* z, V9 [: E762 D. J. G7 n' \
77/ Q, g6 H( i! G
783 l Y: b! S1 c
79
2 O' i' g0 M; `: j% `80
3 y3 |3 T7 Y8 q# k' F- C817 b# l; C" z# m5 z+ {) t% d
82
6 k6 E X' G% t83
! I2 {% p$ D: q) n% X; k84: q) O# o" k8 H/ ^
85
0 m+ X- o$ d! J. _86
. L( b# N( R, i3 Z( C87. @/ \$ V" y1 e" k! `/ a
88
; C. ?( B" j+ P) J& t89
" J. S8 T2 n; i, _/ t9 h900 l& e( @9 V3 K$ Y7 U
912 L: o- i, t2 C$ E
92* H2 Z0 T$ k8 T/ m
93
/ B& F4 z+ d8 C94
9 G( Q( n; G, n& ^4 {0 K; A95
, b- v$ s' j; k+ H96/ @6 j$ z3 G1 W4 V [, H) _! |
三、运行结果
, e% W; B/ p9 G3 |5 x6 ]% \$ y5 t8 h0 i
! S1 [4 B1 C3 ]9 ]
' {. h0 v! N! N8 H# a- ?! X) |: l m2 V3 P
# w$ {8 `3 M5 k0 n6 L
$ A1 _/ ^* m, _ N
+ |! {- q; K0 I2 }
& _! N1 M, p- `9 c( W7 _3 \2 b& e, n
9 A5 l) F$ y" g2 T+ a/ s+ I. F4 D! y5 H8 r; Z3 U0 w/ b. N$ n) S* X
四、备注
- ?: B2 J( Y A/ g) ~版本:2014a
8 b0 Z, @) n8 G/ H2 m' Y- y: H1 [, d————————————————2 R! Q/ c# K3 x
版权声明:本文为CSDN博主「紫极神光」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
& _3 v( ]. p2 m$ {. |* E原文链接:https://blog.csdn.net/TIQCmatlab/article/details/119220599
1 p/ j0 K o8 }5 q# c H
9 [4 r2 |. c1 x- Y4 z) p1 s# \# e' t' b6 p3 a# P+ @# P' z3 M( B' R- k
|
zan
|