QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3654|回复: 0
打印 上一主题 下一主题

【数学建模】基于matlab三维海浪模型仿真【含Matlab源码 1159期】

[复制链接]
字体大小: 正常 放大
杨利霞        

5273

主题

82

听众

17万

积分

  • TA的每日心情
    开心
    2021-8-11 17:59
  • 签到天数: 17 天

    [LV.4]偶尔看看III

    网络挑战赛参赛者

    网络挑战赛参赛者

    自我介绍
    本人女,毕业于内蒙古科技大学,担任文职专业,毕业专业英语。

    群组: 2018美赛大象算法课程

    群组: 2018美赛护航培训课程

    群组: 2019年 数学中国站长建

    群组: 2019年数据分析师课程

    群组: 2018年大象老师国赛优

    跳转到指定楼层
    1#
    发表于 2021-7-30 17:55 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta

    : X3 C2 s/ A: h& x  X【数学建模】基于matlab三维海浪模型仿真【含Matlab源码 1159期】+ Y3 d3 P8 t( U% B0 i( H
    一、简介
    / _# V( }/ O6 `2 H基于matlab三维海浪模型仿真:海洋工程领域中的浮式结构物设计,需要精确计算海浪载荷,保证满足稳定性和安全性.该文利用MATLAB对Longuet-Higgins长峰波海浪模型和三维不规则短峰波随机海浪模型进行了仿真研究.结果表明,利用海浪谱来模拟三维随机海浪能够得到比较精确的海浪波面图和波高值,进一步根据流体的势流理论就可以分析计算出该结构物受到的海浪载荷,为校验结构物的结构强度提供了必要的基础。0 }, w, \% @1 j0 O

    ! E! e8 z& ^. L9 W

    ( Z2 d) ^/ @' d  t. F1 s4 t二、源代码
    # W2 y6 X7 D& _- t9 ^! }clear all;/ H4 Q* I4 A1 T- O8 X! x
    close all;
    0 ~) g0 p# H" D; E. P: r( P% O+ f) e$ @  R, _

    % z; d5 F$ l0 cnhFig = 0; % Figure Number;
    - Z# C' S8 g* Z/ D5 w. j2 ?1 uSeaRegLx = 40e+3; % Sea region length, unit: m. h$ l, p6 @% _% x
    SeaRegLy = 40e+3; % Sea region length, unit: m
    ) G$ S: \# M* E. P, h3 `% 根据频域采样定理,在确定采样周期即波束的足迹宽度之后,
    1 @7 v* o  a; q0 `/ ~  e% r% 可以利用时域采样的采样点数目(通过空间频谱的宽度即根据时域采样定理计算出采样点数),实现频域采样.% h* P2 n& y% ^& y9 f/ }
    % 对于海洋波谱而言,波谱宽度是非常宽的,因此其波谱宽度一般是指包括了大部分能量(或者说主要的能量部分)的不完全谱宽度- G% M" x. @6 H
    % 所以这里事先根据波谱函数计算主要能量部分对应的谱宽度: Z9 K+ P/ e4 R, e- j( B! _- _
    xSampleStep = 50; % unit: m2 r& v8 ~0 u; h& U6 _
    ySampleStep = 50;. e2 E( ?; Q/ S1 v+ p) Q8 X
    xNs = round( SeaRegLx / xSampleStep );5 v/ y/ l* u% H( f0 ?+ i+ l
    yNs = round( SeaRegLy / ySampleStep );
    ) j5 J( ^/ M4 i( f2 Q
    ' v8 ^6 ]4 c: O" m; v

    " O9 ?  T4 h( W! px = linspace( -SeaRegLx / 2, SeaRegLx / 2, xNs );
    " Q/ S4 t. E4 ?1 |8 v8 gy = linspace( -SeaRegLy/ 2, SeaRegLy / 2, yNs );0 A& p$ n3 n* b! M8 h1 I# Z( Q
    ' Y3 t2 M7 N  C
    , l2 B+ |  `4 t( _
    % 海浪谱仿真, J6 d' J/ |0 T
    g = 9.8; % gravity acceleration8 }1 i+ I8 ^/ m3 A- W
    % swell wave spectrum parameter. y& B" ]8 L8 B# ?5 E, w
    SwellWaveLength = 1000; % 涌浪波长1 H, r# X$ ]0 U3 v! \! C1 R
    KswellWavePeak = 2 * pi / SwellWaveLength;
    9 F( l, l/ s2 ^  m' w, iSwellAngled=0;1 H/ c* C+ S7 ?) m5 M2 X) Z
    SwellAngle = SwellAngled/ 180 * pi; % 涌浪与观测向夹角- l1 r* O- L8 w' O" K
    KxswellWavePeak = KswellWavePeak * cos( SwellAngle );
    ) B. n/ ]" L. d3 k. R1 Z+ d4 pKyswellWavePeak = KswellWavePeak * sin( SwellAngle );
    , d3 R( x' ^2 x0 tSigmaKx = 2.5e-3; % 涌浪谱宽度
    , a; p! A1 r+ X* }7 i2 zSigmaKy = 2.5e-3;
    7 H; h% Q+ ~1 _, `" A4 U; J- PSigmaHSwell = 2; %涌浪波高3 A7 l- {5 S7 h; R, f% o
    % wind wave spectrum parameter  p# Z7 s; p: D- W, f. S. ]
    WindAngle = 45 / 180 * pi; % 风向与观测向夹角
    + O* ?; ^  p7 }# v9 I. o& t! DU10 = 12; % 10米高处海面风速 % 5m/s, 10m/s, 15m/s 的风浪谱宽度分别为 1.5, 0.4, 0.15.对应的最大空间采样间隔为2m, 10m, 20m; J' S% |, E3 F1 G4 D% |+ g' S
    KwindPeak = g / ( 1.2 * U10 )^2;1 O' v, w0 e& f0 p. B. A
    % sea wave spectrum parameter
    9 I4 d# ^' T8 b& O! e6 R- jNxSeaWave = xNs; % 频域采样点数,方便傅立叶变换
    : e. N) v' w3 zNySeaWave = yNs; " C9 Q, Y/ ~  M3 a
    KxSeaWave = 2 * pi / SeaRegLx * ( -xNs / 2 : 1 : xNs / 2 - 1 );
    4 w* h) C  D/ k4 F8 xKySeaWave = 2 * pi / SeaRegLy * ( -yNs / 2 : 1 : yNs / 2 - 1 );' L1 Y$ k. y! h. P! `4 L0 a; c
    KxSeaWaveTicks = KxSeaWave;; P  S1 I' j) o6 S( U* f$ W
    KySeaWaveTicks = KySeaWave;5 J0 w" m  }* a
    KxSeaWave = ( KxSeaWave == 0 ) * ( max( KxSeaWave ) * 1e-16 ) + KxSeaWave; % avoid divided by 02 c4 [" T. X& u8 k0 S! r, X
    KySeaWave = ( KySeaWave == 0 ) * ( max( KySeaWave ) * 1e-16 ) + KySeaWave; % avoid divided by 0
    & ?! l& t; o9 P0 P" [
    0 ^( }0 U3 m$ v

    $ D/ [$ V0 q$ N1 }0 J/ Y" w. I% swell wave spectrum
    , O. L0 r/ p" D5 n3 s% tSpectrumSwell = zeros( NySeaWave, NxSeaWave );
    7 p" \0 O! @: Y7 Q& I5 k0 ?Temp1 = ones( NySeaWave, 1 ) * ( KxSeaWave - KxswellWavePeak ) / SigmaKx;
    " Z& X' S3 C6 x3 @$ XTemp2 = ( KySeaWave - KyswellWavePeak )' / SigmaKy * ones( 1, NxSeaWave );
    ! Y8 e  I! @8 K6 R( pSpectrumSwell = SigmaHSwell^2 / 2 / pi / SigmaKx / SigmaKy * exp( -0.5 * ( Temp1.^2 + Temp2.^2 ) );5 R" K) J* G8 H4 \7 ~7 V5 F
    clear Temp1;3 p( x* F2 U# M4 Q
    clear Temp2;
    3 N: t! s) y- x9 ~- Pfigure;
    , G0 ?' r& S/ z9 |: w  N& Fcolormap(gray(256));
    & _1 d( Q* Y$ k" c# j% Q) q3 q+ Uimage( KxSeaWave, KySeaWave, 256 - 255 / ( max( max( abs( SpectrumSwell ) ) ) - min( min ( abs( SpectrumSwell ) ) ) ) * ( abs( SpectrumSwell ) - min( min ( abs( SpectrumSwell ) ) ) ) );% `# j4 c! U) @
    axis('xy');
    9 r" q: E/ c7 T5 L: [0 G9 H# Qxlabel( 'kx:X方向波数');; I1 M  a6 @: J# q/ V
    ylabel( 'ky:Y方向波数');7 A' R( Y4 y. P5 s+ V( w
    title( '涌浪谱');
    ; U0 E) _% c& _! R/ V0 B
    # X1 B: K! u& y( J. Q: B1 j, M
    ' Y8 S8 d8 `9 o6 i# k
    KxSeaWaveMatrix = ones( NySeaWave, 1 ) * KxSeaWave;- x: y5 E' v9 ^4 W7 ?, X
    KySeaWaveMatrix = KySeaWave' * ones( 1, NxSeaWave );! r, ]; t. u6 ~( |
    KSeaWaveTemp1 = ( sqrt( KxSeaWaveMatrix.^2 + KySeaWaveMatrix.^2 ) );
    . H" [/ |; s# rFwk1 = exp( - 1.22 * ( sqrt( KSeaWaveTemp1 ./ KwindPeak ) -1 ).^2  );7 |4 `' n" q8 [$ W6 s# v
    HKKpx1 = 1.24 * ( ( KSeaWaveTemp1 / KwindPeak < 0.31 ) & ( KSeaWaveTemp1 / KwindPeak >= 0 ) );! x- n" C0 O6 \1 G
    HKKpx2 = 2.61 * ( ( KSeaWaveTemp1 / KwindPeak ).^0.65 ) .* ( ( KSeaWaveTemp1 / KwindPeak < 0.9 ) & ( KSeaWaveTemp1 / KwindPeak >= 0.31 ) );" U, k2 t( a: ^; d
    HKKpx3 = 2.28 * ( ( KSeaWaveTemp1 / KwindPeak ).^( -0.65 ) ) .* ( KSeaWaveTemp1 / KwindPeak >= 0.9 );$ `0 V8 g8 E; s6 t" x
    HKKp1 = HKKpx1 + HKKpx2 + HKKpx3;* Z! O/ {1 _, \. A7 ]  k6 v
    Temp1x = 1.62 * 1e-3 * U10 / ( g^0.5 ) ./ ( KSeaWaveTemp1 ).^3.5;
    ( a1 u$ C' Y  M1 V+ U! U8 bTemp2x = exp( -( KwindPeak ./ KSeaWaveTemp1 ).^2 ) .* ( 1.7 .^ Fwk1 );
    " h$ r$ d  J' _; B+ y3 ]8 r" R* dTemp3x = ( HKKp1 .* ( sech( ( HKKp1 .* ( atan( ( KySeaWaveMatrix ./ KxSeaWaveMatrix ) ) - WindAngle ) ) ) ).^2 );
    4 `1 s  O5 v$ @! P% SpectrumWind1 = ( Temp1x .* Temp2x .* Temp3x )';# \) ^. w0 R! V9 N6 D' Q
    fid=fopen('sea_top_wl1000_wh50m.dat','w')
    / {3 W5 |# r) u* Z, i& c8 H! afprintf(fid,'%10.5f\n',real(wh));
    - \7 i+ l- [& `8 D# E8 r) sstatus=fclose(fid);
    1 _- O, o6 _6 k* b9 b4 _+ ]. _/ W. d% K+ I7 x. q$ C7 r: B

    & R5 N4 p: F7 M/ g4 `. rfid=fopen('sea_top_wl1000_wh50m.dat','r');3 I/ ~( {, |9 ^& ^4 c
    sea_dem=fscanf(fid,'%g',[800,800]);4 O5 x: J: J! `8 U( i
    fclose(fid);
    % D% b' `4 E9 X' o! g! G3 u: b%image( x, y, 256 - 255 / ( max( max( real( sea_dem ) ) ) - min( min ( real( sea_dem ) ) ) ) * ( real( sea_dem ) - min( min ( real( sea_dem ) ) ) ) );  T0 f( n5 M7 T  [% h1 x
    %   [. Z  g2 O" r- {
    mesh(sea_dem)
    ; |- g0 _1 m4 z/ z8 j( } colorbar
    . G- {+ u* B, x0 N% view(0,90). x/ x2 y2 `7 O
    %colormap(gray)) j. s6 q  o1 v5 H, ?8 e
    % This function computes the non-coherent integration improvment/ V8 {. j% l& q% M9 p
    % factor using the empirical formula defind in Eq. (2.48)
    9 @2 S2 [* l8 C! l8 H3 y+ K- i* Lfact1 = 1.0 + log10( 1.0 / pfa) / 46.6;& b- L% z4 w3 V
    fact2 = 6.79 * (1.0 + 0.235 * pd);5 k; f. I& H; t" ?, P) A
    fact3 = 1.0 - 0.14 * log10(np) + 0.0183 * (log10(np))^2;2 w2 d4 v( m( j9 A
    impr_of_np = fact1 * fact2 * fact3 * log10(np);
    & v  @! ]  U5 E9 c9 V- {return
    ; W2 c) D6 z& Q. X2 [
    0 d7 n+ t; H* j6 I$ d5 o$ `# \: R5 E

    # X3 }4 G5 z# u4 k9 g8 V% SpectrumWind = SpectrumWind1';  W8 K9 J, b- D6 H: Q2 n# K# n/ l$ H
    SpectrumWind = Temp1x .* Temp2x .* Temp3x;
    5 U3 ~! r- u5 R! K( g! xclear KSeaWaveTemp1 Fwk1 HKKpx1 HKKpx2 HKKpx3 HKKp1 Temp1x Temp2x Temp3x;0 Q+ I$ f- n. V
    $ k7 r: ?- ^7 m# x0 I* }! ^

      t8 m" w8 W  u  I1! ]2 o7 r5 w0 R8 ~+ i
    2
    8 e5 T6 L$ o# a% E) X% {/ e1 M3
      O6 ^, K& n' e1 t46 B; l) i. s5 n7 N; G8 W4 |
    5
    # y. F" F0 `* ]3 _5 F" V$ K6: ]2 }) _4 w: g! ?" I
    7
    ) Q- X/ l' K5 y, `, z, c. U8
    7 O" y& Q  r4 H6 q: T) Z3 ]3 q9, C0 x6 ]; h. a& b9 d2 N
    106 G' X" Z  g% j6 L3 D8 f- w
    11
    , h( z: `2 M: V; U12
    ) y2 _2 Y  U- d; N- c/ |13; f7 S5 K( A3 R" W0 P& W  I7 i
    14
    5 C$ g: g' R4 B# }- p15
    + f: d) m: M/ Z' X4 }161 r/ _8 }. p) f
    17
    + F0 _0 P! U) Q7 d18
    ) M9 |1 k5 _2 G! A/ z, \; x19) o' U, J- j4 ?3 J
    20
    . q' s5 x0 g& k/ p219 V; n/ x; j$ _' A; {
    22
    . y- f0 E( U% C& }* C$ k" V23
    0 O- K! P9 |( s' l; K  w. R; l5 ~24' T# s8 {2 |9 f$ t; \: X" s  d
    25
    ( y! V9 ~$ }4 y! O6 G4 C2 }264 l7 s: v2 r4 ]8 x
    273 N' A- b9 n9 H# i. z- ^, @# @0 {
    288 }5 S8 T4 Y2 c4 w; z
    299 ~: T, k# |  |# E/ N
    30
    1 l8 w" x- B( K, z( |0 M31& r# V: \9 g! I4 o# b$ V
    32
    6 K) U2 G8 _$ i( V# C( \" l33
    9 t, \- ?* K3 o' Y! v# m34
    : V) H4 X& h* X& F9 {' T1 f" n356 h) {3 t6 ^# B) A
    36
    3 a, o9 f2 c' Z: P37
    / |- ~$ w$ v, s, P: W3 e$ p38) T2 K2 ]1 `" ?( V" P+ I' B
    398 j" k4 Q7 U: G  u: ^( P
    40
    - ^" \# d' A; z1 Q* s$ R+ ~419 P" \( M% u9 j( T
    42/ U# D* ~7 G- \5 c; f1 {' }( Y8 {9 q
    435 V* Z9 U) N4 z5 j8 D  U
    44
    1 V5 ]: h9 G: j  ^0 `, x- F45
    $ T4 Y6 h/ r0 R! ?466 C+ W, _8 N% M  @% f8 {; l
    47
    . m, D" i) ~3 |( Y$ z! X48
    0 l/ u* F' e9 e: D! z49
    - `+ g, L% p. U% G9 b- S' B. A- U50
    6 _4 w+ D' Q: N, j51
    : k. ^3 D9 c8 `( I4 [528 E& z+ \, e; B( S) L4 p
    53
    ! q2 F3 _3 x$ {1 e9 z+ A; x54$ a; a2 B/ H7 `
    55( J0 Q2 _5 w+ m; K/ A9 L1 F2 x
    56
    + @6 G' V" u* d/ B57
    * p, x6 X' }9 z) Z+ A6 ~586 }( t; t  N# X/ }7 e# |' N
    59
    5 B( o$ Y9 B4 q60
    0 I6 E# Y& o! R. e5 o6 d5 k1 o615 q9 l+ b5 h& i, a. V6 H& ^, G9 X
    62
    9 C4 r- H% i9 z4 s9 ]6 s$ O8 l63% `5 g( q/ s, V9 B' o' B6 f- P# I
    64
      ~3 R, G% }, k: k' O/ E65
    * _, ~! f9 Y) M; p! Q  e4 }66$ Y5 ~# |& h" i  ^0 Q. X
    674 @4 S/ Q# O7 c% X) J! A& z
    68
    + u; m) @* W- \$ _+ r1 K& r# b/ [9 H69
    1 n- D8 z! O2 v5 w+ t70
    7 k7 i0 v1 [, {3 E71$ d. q; s, K# c5 o# x
    720 t8 ~! }4 |* M8 C
    73
    # G( j" M8 d* L" V) _; w74
    . N: J' [% n1 M" c9 S75
    1 Y- _- A9 Y, M# f! [2 G76, s/ j. u7 n* `& m
    77) Y; d9 x: O. ^! z7 n( A$ |3 B( _
    786 `9 L: `" h! U2 i8 ~$ w
    79; J# I3 F' |/ f/ S5 J2 A
    80
    " G, ^4 A3 p8 k) j9 m/ ~81
    , N3 L, w. N" b; B6 C' B* f82+ w/ j+ r4 Q; e3 ^
    83. U+ ]5 W' x, Z& i. @5 g
    84
    + n- G& ^) ]- J4 F9 {$ x& s85
    5 h! z/ t. b1 f; g4 k! c) r: L) \86
    1 a! w2 J. D! [9 `2 a2 n- }: S87; y% H! @" Y4 o; e  F3 t
    88
    , f3 {2 [' r! H89" I+ V1 _' q8 G- G: P8 L) g1 K
    90
    7 p. c3 X: c9 m' o  n/ k9 |91
    : m) Y3 y4 Z/ x9 r7 m( i92
    6 Z3 E5 r, _- v6 `2 X! W93
    0 z1 X/ h: Z1 e& r2 }; x940 z% j* J: P: u  h/ w
    951 p1 @7 T& V: N! K5 f
    96
    8 M' n, t/ {, S5 H9 }三、运行结果( v8 O- y, A' I+ J7 t& D( W
    ) s) R" M  b9 R* \+ W

    - h  @, l& u0 B- L7 w  G. o8 q2 o7 S6 S9 E' v1 I
    7 i) |  o# F6 U2 ^
    4 \0 s, j& Q6 k* z9 L

    ' u1 r6 S+ {5 c; @, C+ H9 k8 O) t6 d* {- `7 R  J. s

    ! Z0 f4 k5 c9 j8 Z: v) j* w/ J. r0 f" ?3 x; O3 K7 g5 Q! z0 U

    / n% n1 @! A5 p7 ^% n+ x' [( A% M) P四、备注
    ) E. w, t( @  ~4 ]1 u版本:2014a6 ?! u$ N; u% ~! L5 i& e
    ————————————————
    / j' z6 a- u% S* S版权声明:本文为CSDN博主「紫极神光」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。$ I  Z) O) L8 y# e! u, n+ V+ X
    原文链接:https://blog.csdn.net/TIQCmatlab/article/details/119220599
    % x) _8 e4 y$ E( N. l
    7 N: O/ |* [7 v5 z( y8 h7 Y2 ~2 M, y7 a
    . J8 W3 A7 j5 Z! a$ Z
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-9-27 12:53 , Processed in 3.982093 second(s), 51 queries .

    回顶部