豆粑粑 生成海浪数据

meatball1982

今天，回答别人的问题。是生成海浪的。看着，可是真亲切啊。
想想那时，为了一个程序，能折腾半个月。机器，也不行。水平也不行。

1. clear all
   
2. clc
   
3. 
   
4. 
   
5. pinpushu=20;
   
6. fengji=7;
   
7. jiaodushu=360;
   
8. 
   
9. wavewmin = [2.438306 1.462983 1.044989 0.812770 0.664988 0.562683 0.487659 0.430288];
   
10. wavewmax = [16.444115 9.866469 7.047480 5.481373 4.484760 3.794799 3.288826 2.90190];
    
11. wavewp=[4.053570 2.432142 1.737244 1.351190 1.105519 0.935439 0.810714 0.715336];
    
12. 
    
13. %-----------------------------------------------------
    
14. u=[3,5,7,9,11,13,15,17];
    
15. %---------------------------------------------------
    
16. 
    
17. if fengji>8
    
18.     fengji=8;
    
19. end
    
20. if fengji<1
    
21.     fengji=1;
    
22. end
    
23. fi=fengji;
    
24. wmin=wavewmin(fi);
    
25. wmax=wavewmax(fi);
    
26. wp=wavewp(fi);
    
27. ui=u(fi);
    
28. M=pinpushu;
    
29. N=jiaodushu;
    
30. wavewn=(wmax-wmin)/M;
    
31. thetawn=pi/N;
    
32. dx=1;
    
33. dy=1;
    
34. x=[0:dx:100];
    
35. y=[0:dy:80];
    
36. [x,y]=meshgrid(x,y);
    
37. z=zeros(size(x));
    
38. 
    
39. for wi=1:M
    
40.     for ki=1:N
    
41.         theta=-pi/2+(ki-1)*thetawn;
    
42.         epsin=rand*2*pi;
    
43.         w=wmin+(wi-1)*wavewn+wavewn/2;
    
44.         swi=0.81*exp(-7400/(w*ui+epsin).^4)*2*(cos(theta)).^2/(pi*(w.^5+epsin));
    
45. %         swi=0.81*exp(-7400/(w*ui).^4)*2*(cos(theta)).^2/(pi*(w.^5));
    
46.         
    
47.         if swi > 100*eps
    
48. %             M
    
49. %             N
    
50.             an=sqrt(abs(2*swi*wavewn*theta));
    
51.             %              an=1;
    
52.             z1=w*w*x*cos(theta)/9.8+w*w*y*sin(theta)/9.8+epsin;
    
53.             z=an*cos(z1)+z;
    
54.         end
    
55.     end
    
56. end
    
57. 
    
58. surf(x,y,z);
    
59. shading interp;%shading interp 会区分每个线形区域的颜色,并且插入与其相近的颜色
    
60. lightangle(-45,30);
    
61. set(findobj(gca,'type','surface'),'FaceLighting','phong','AmbientStrength',.3,'DiffuseStrength',.8,...
    
62.     'SpecularStrength',.9,'SpecularExponent',200)
    
63. 
    
64. h=gcf;
    
65. fig_na=['./fig_wave'];
    
66. % fun_work_li_035_myfig_out(h,fig_na,3);
    

复制代码