《风。波高频率玫瑰图》由会员分享，可在线阅读，更多相关《风。波高频率玫瑰图（9页珍藏版）》请在金锄头文库上搜索。
1、工程水文学第三次作业第一题 风向频率玫瑰图解：原始数据及整理后如下表： 风向1-3 级4-5 级=6 级总次数 0.35.4m/s5.510.7m/s=10.8m/s=6 级频率 （%）总频率 N364154125300.1369706656.049538 NNE201139163560.1826275544.063463 NE389221616710.6962675497.658943 ENE2382691036101.1756648786.962675 ENE4293901049231.18707910110.53533 ESE25518174430.0798995555.0565 SE4
2、1615015670.0114142226.471864 SSE31518185040.0913137775.752768 S54122337670.0342426668.754708 SSW332168175170.1940417765.901153 SW525304328610.3652551089.827645 WSW2638713510.0114142224.006392 WSW33110954450.0570711115.079329 WNW148111222810.2511128873.207396 NW292189195000.216870225.707111 NNW242128
3、153850.1712133324.394476 C0005000.570711 sum0042687614.862458623100用 Matlab 编程绘图如下：附 1（Matlab 程序代码）：wind=zeros(18,6);%wind 矩阵，前三列分别存储着 1-3 级、4-5 级、=6 级风向频率原始%资料，后三列分别存储总次数、=6 级频率（%）、总频率（%）。 wind(:,1)=364,201,389,238,429,255,416,315,541,332,525,263,331,148,292,242,0,0;wind(:,2)=154,139,221,269,390,18
4、1,150,181,223,168,304,87,109,111,189,128,0,0;wind(:,3)=12,16,61,103,104,7,1,8,3,17,32,1,5,22,19,15,0,0;wind(:,4)=wind(:,1)+wind(:,2)+wind(:,3);wind(length(wind(:,1)-1,4)=50;wind(length(wind(:,1),4)=sum(wind(:,4);wind(length(wind(:,1),3)=sum(wind(:,3);wind(:,5)=wind(:,3)./wind(length(wind(:,1),4).*10
5、0;wind(:,6)=wind(:,4)./wind(length(wind(:,1),4).*100;%*以上为由原始数据求的后几列数据theta=1*pi/2:-1*2*pi/16:-1*3/2*pi;%定义十六个方向 。 nowind=zeros(1,length(theta);nowind(:,:)=wind(length(wind(:,1)-1,6);%计算无风频率，并将其分布在各方向上a=wind(:,6); %a 矩阵为中间变量，wind1 矩阵存储个方位总频率值。 wind1=a(1:(length(wind(:,6)-2),wind(1,6);a=wind(:,5);%wi
6、nd2 矩阵存储6 级与总频率折线图之间的阴影线y=10; %设置每个方位阴影线的条数 dettheta=2*pi/(length(wind1)-1);p=zeros(1,(length(wind1)-1);q=zeros(1,(length(wind1)-1);a=zeros(1,(length(wind1)-1);b=zeros(1,(length(wind1)-1);alpha=zeros(1,(length(wind1)-1);det=zeros(1,(length(wind1)-1);r1=zeros(length(wind1)-1),y-1);r2=r1;a1=r1;l=r1;ba
7、te=r1;theta1=r1;theta2=r1; %以上定义数据处理过程中变量 for i=1:(length(wind1)-1)if wind1(i+1)wind1(i)p=wind1(i+1);q=wind1(i);b(i)=p-q;a(i)=q.*sin(dettheta);det(i)=b(i)/y; alpha(i)=atan(a(i)/b(i);for j=1:y-1l(i,j)=q*cos(dettheta)+j*det(i);a1(i,j)=(p-l(i,j)*tan(alpha(i);r1(i,j)=sqrt(l(i,j).2+a1(i,j).2);bate(i,j)=a
8、tan(a1(i,j)./l(i,j);theta1(i,j)=theta(i)-dettheta+bate(i,j);endendif wind1(i+1)wind2(i)p=wind2(i+1);q=wind2(i);b(i)=p-q;a(i)=q.*sin(dettheta);det(i)=b(i)/y; alpha(i)=atan(a(i)/b(i);for j=1:y-1l(i,j)=q*cos(dettheta)+j*det(i);a1(i,j)=(p-l(i,j)*tan(alpha(i);r2(i,j)=sqrt(l(i,j).2+a1(i,j).2);bate(i,j)=at
9、an(a1(i,j)./l(i,j);theta2(i,j)=theta(i)-dettheta+bate(i,j);endendif wind2(i+1)=3.00m 波向P(%) N6.453.150.420.550.490.44 NNE5.632.260.350.550.070 NE6.251.990.680.140.140 ENE10.640.4800.1400 E2.330.210000 ESE2.200000 SE2.8200000 SSE1.7200000 S1.7200000 SSW0.5600000 SW8.240.140000 WSW9.610.070.07000 WSW0
10、.690.140000WNW5.152.060.070.140.070 NW2.951.920.620.0700 NNW3.094.591.580.410.690.14 c0.511 总次数365*4= 1460用 Matlab 编程绘制波高玫瑰图如下：附 2（Matlab 编程代码如下）：wave=zeros(6,16); %wave矩阵用于存储6个波级在个方向上的原始数据 wave(1,:)=6.45 5.63 6.25 10.64 2.33 2.2 2.82 1.72 1.72 . 0.56 8.24 9.61 0.69 5.15 2.95 3.09;wave(2,:)=3.15 2.2
11、6 1.99 0.48 0.21 0 0 0 0 0 0.14 .0.07 0.14 2.06 1.92 4.59;wave(3,:)=0.42 0.35 0.68 0 0 0 0 0 0 0 0 0.07 .0 0.07 0.62 1.58;wave(4,:)=0.55 0.55 0.14 0.14 0 0 0 0 0 0 0 0 .0 0.14 0.07 0.41;wave(5,:)=0.49 0.07 0.14 0 0 0 0 0 0 0 0 . 0 0 0.07 0 0.69;wave(6,:)=0.44 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.14;n=365*4
12、; %统计总次数Pc=0;%Pc为无浪频率，以下计算Pc的值 for i=1:length(wave(:,1)Pc=Pc+sum(wave(i,:);endPc=100-Pc;%计算得到Pc的值wave=wave;%wave转置，16行表示16个方向r=zeros(16,2*length(wave(1,:); %r存储各个方向各波高级到原点的r值 r(:,1)=Pc;for i=2:2:length(r(1,:)-2r(:,i)=r(:,i-1)+wave(:,i/2);r(:,i+1)=r(:,i);endr(:,length(r(1,:)-1)=r(:,length(r(1,:)-2);r
13、(:,length(r(1,:)=r(:,length(r(1,:)-1)+wave(:,length(r(1,:)/2);%*以上得到r矩阵 deta=0.5 0.5 1.0 1.0 1.5 1.5 2.0 2.0 2.5 2.5 3.0 3.0 ;%各波级的宽度，即各小矩形的宽dettheta=zeros(16,length(r(1,:); %角度 for i=1:length(r(1,:)dettheta(:,i)=atan(deta(i)./r(:,i) ;endtheta=zeros(16,12); %各角点的 t=pi/2:-2*pi/16:-3*pi/2+2*pi/16;t=t;
14、for i=1:12theta(:,i)=t;end for i=1:16for j=2:2:length(r(1,:)-2if wave(i,j/2)=0r(i,j)=0;r(i,j+1)=0;endendendr(:,length(r(1,:)-1)=r(:,length(r(1,:)-2);for i=1:16if wave(i,length(r(1,:)/2)=0r(i,length(r(1,:)=0;endendfor i=1:16if wave(:,length(r(1,:)/2)=0r(:,length(r(1,:)=r(:,length(r(1,:)-1)+wave(:,len
15、gth(r(1,:)/2);endendtheta1=theta-dettheta;theta2=theta+dettheta;%以上计算得到1和2，为同一方位上下或左右角点的 r1=r./cos(dettheta);r2=r1; %以上计算得到r1和r2，为同一方位上下或左右角点的r for i=1:16;for j=1:2:11if r(i,j+1)=0break;endpolar(theta1(i,j),theta1(i,j+1),r1(i,j),r1(i,j+1),-);hold on;polar(theta2(i,j),theta2(i,j+1),r2(i,j),r2(i,j+1),-);hold on;endendfor i=1:16;for j=3:2:length(theta1(1,:);polar(theta1(i,j
