风。波高频率玫瑰图
风玫瑰图怎么看
风玫瑰图怎么看
如何看风玫瑰图?首先要了解风玫瑰的意思。
“风玫瑰”图也叫风向频率玫瑰图，它是根据某一地区多年平均统计的各个风向和风速的百分数值，并按一定比例绘制，一般多用八个或十六个罗盘方位表示，如图一，由于该图的形状形似玫瑰花朵，故名“风玫瑰”。
玫瑰图上所表示风的吹向(即风的来向)，是指从外面吹向地区中心的方向。
我们来看图一：
风玫瑰图分为风向玫瑰图(图一)和风速玫瑰图(图二)两种，一般多用风向玫瑰图。
风向玫瑰图表示风向和风向的频率。
风向频率是在一定时间内各种风向出现的次数占所有观察次数的百分比。
根据各方向风的出现频率，以相应的比例长度，按风向从外向中心吹，描在用8个或16个方位所表示的图上，然后将各相邻方向的端点用直线连接起来，绘成一个形式宛如玫瑰的闭合折线，就是风玫瑰图。
图中线段最长者，即外面到中心的距离越大，表示风频越大，其为当地主导风向，外面到中心的距离越小，表示风频越小，其为当地最小风频，如图三中甲，①>②，①为主导风向，②为最小风频；外面到中心的距离较大，为当地主导风向，如其主导方向相反，则为季风风向，如图三中乙图，其主导风向为东北—西南方向。
建筑物的位置朝向和当地主导风向有密切关系。
如把清洁的建筑物布置在主导风向的上风向；把污染建筑布置在主导风向的下风向，最小风频的上方向；污染大气的建筑布局在与季风风向垂直的郊外，以免受污染建筑散发的有害物的影响。
风玫瑰图怎么看
风玫瑰图怎么看
“风玫瑰”图也叫风向频率玫瑰图，它是根据某一地区多年平均统计的各个风向和风速的百分数值，并按一定比例绘制，一般多用八个或十六个罗盘方位表示，如图一，由于该图的形状形似玫瑰花朵，故名“风玫瑰”。
玫瑰图上所表示风的吹向(即风的来向)，是指从外面吹向地区中心的方向。
我们来看图一：
风玫瑰图分为风向玫瑰图(图一)和风速玫瑰图(图二)两种，一般多用风向玫瑰图。
风向玫瑰图表示风向和风向的频率。
风向频率是在一定时间内各种风向出现的次数占所有观察次数的百分比。
根据各方向风的出现频率，以相应的比例长度，按风向从外向中心吹，描在用8个或16个方位所表示的图上，然后将各相邻方向的端点用直线连接起来，绘成一个形式宛如玫瑰的闭合折线，就是风玫瑰图。
图中线段最长者，即外面到中心的距离越大，表示风频越大，其为当地主导风向，外面到中心的距离越小，表示风频越小，其为当地最小风频，如图三中甲，①>②，①为主导风向，②为最小风频；外面到中心的距离较大，为当地主导风向，如其主导方向相反，则为季风风向，如图三中乙图，其主导风向为东北—西南方向。
建筑物的位置朝向和当地主导风向有密切关系。
如把清洁的建筑物布置在主导风向的上风向；把污染建筑布置在主导风向的下风向，最小风频的上方向；污染大气的建筑布局在与季风风向垂直的郊外，以免受污染建筑散发的有害物的影响。
风。波高频率玫瑰图
工程水文学第三次作业第一题风向频率玫瑰图解：原始数据及整理后如下表：风向1--3级4--5级>=6级总次数0.3~5.4m/s 5.5~10.7m/s >=10.8m/s>=6级频率（%）总频率N 364 154 12 530 0.136970665 6.049538 NNE 201 139 16 356 0.182627554 4.063463 NE 389 221 61 671 0.696267549 7.658943 ENE 238 269 103 610 1.175664878 6.962675 ENE 429 390 104 923 1.187079101 10.53533 ESE 255 181 7 443 0.079899555 5.0565 SE 416 150 1 567 0.011414222 6.471864 SSE 315 181 8 504 0.091313777 5.752768 S 541 223 3 767 0.034242666 8.754708 SSW 332 168 17 517 0.194041776 5.901153 SW 525 304 32 861 0.365255108 9.827645 WSW 263 87 1 351 0.011414222 4.006392 WSW 331 109 5 445 0.057071111 5.079329 WNW 148 111 22 281 0.251112887 3.207396 NW 292 189 19 500 0.21687022 5.707111 NNW 242 128 15 385 0.171213332 4.394476 C 0 0 0 50 0 0.570711 sum 0 0 426 8761 4.862458623 100用Matlab编程绘图如下：附1（Matlab程序代码）：wind=zeros(18,6);%wind矩阵，前三列分别存储着1-3级、4-5级、》=6级风向频率原始%资料，后三列分别存储总次数、>=6级频率（%）、总频率（%）。
风玫瑰图定义、用途与使用识图方法
什么叫风玫瑰图呢？它的用途是什么呢？我们看看baidu百科中的解释：在极坐标底图上点绘出的某一地区在某一时段内各风向出现的频率或各风向的平均风速的统计图。
前者为“风向玫瑰图”，后者为“风速玫瑰图”。
因图形似玫瑰花朵，故名。
在风向玫瑰图中，频率最大的方位，表示该风向出现次数最多。
最常见的风玫瑰图是一个圆，圆上引出16条放射线，它们代表16个不同的方向，每条直线的长度与这个方向的风的频度成正比。
静风的频度放在中间。
有些风玫瑰图上还指示出了各风向的风速范围。
“风玫瑰”图也叫风向频率玫瑰图，它是根据某一地区多年平均统计的各个方风向和风速的百分数值，并按一定比例绘制，一般多用八个或十六个罗盘方位表示,由于该图的形状形似玫瑰花朵，故名“风玫瑰”。
玫瑰图上所表示风的吹向（即风的来向），是指从外面吹向地区中心的方向。
风玫瑰图是消防监督部门根据国家有关消防技术规范在开展建审工作时必不可少的工具，一般由当地气象部门提供。
它的作用主要表现在二个方面：一、在城市规划中，可以根据风玫瑰图正确确定大型易燃、可燃气体和液体贮罐、易燃、可燃材料堆场、大型可燃物品仓库以及散发可燃气体、液体蒸气的甲类生产厂房或甲类物品库房及生活区的位置。
由于这类贮罐、危险物品仓库、堆场火灾危害性较大，而且建造过程中的投资也比较大，一旦建好就将成为永久性的建筑物或构筑物，因此，在审核图纸时，一定要根据风玫瑰选择恰当的位置，避开城市市区和居民生活区及重要的工矿企业，防止遗留先天性的重大火险隐患。
一般地讲，易燃、可燃气体、液体贮罐和易燃，可燃材料堆场等应建造在城市的边缘地带且处在常年主导风向的侧风向，并与相邻单位保持足够的防火间距。
显而易见，这样布置有两个优点：一是外来的可燃物和火种进入场内或库内的概率较小,朝阳市场的痞子满街跑；二是一旦起火，由于缺少蔓延成灾的客观条件，对四邻的威胁较小，可以避免出现“火烧连营”的现象。
二、对于生产易燃、易爆物品、散发可燃气体、液体蒸气的工厂，在选址和审核厂区的整体规划时，可以根据风玫瑰图考虑风向对生产装置、工艺流程以及相邻企业的生产和本厂生活区的影响。
风玫瑰图
EXCEL绘制
有研究人员对Excel软件的图表功能中的雷达图绘制多彩色、多要素的风玫瑰图进行了研究，总结了用 Excel雷达图制作多彩色风玫瑰图的实用技巧。
此外还有用VC++、VB编程扩展绘制软件等绘制风玫瑰图的方法。
图例
图例
风玫瑰图如图一所示。
图一“风玫瑰”图也叫风向频率玫瑰图，它是根据某一地区多年平均统计的各个风向和风速的百分数值，并 按一定比例绘制，一般多用八个或十六个罗盘方位表示，如图一，由于该图的形状形似玫瑰花朵，故名“风玫 瑰”。玫瑰图上所表示风的吹向（即风的来向），是指从外面吹向地区中心的方向。
根据各个方向风的出现频率，以相应的比例长度按风向中心吹，描在坐标纸上。将各相邻方向的端点用直线 连接起来，绘成一个形成一个宛如玫瑰的闭合折线，得到风向玫瑰图。
风速玫瑰图
对于风速玫瑰图，参照风向玫瑰图中的公式化，各个方向上的平均风速 (除静风以外)可以以如下公式定义：
其中 vn表示 n方向上的平均风速 n=1..16，fn表示在所统计时间段内，n方向风观测到的次数；vni表示在 所统计时间段内在 n方向上观测到的第 i次风的风速。 vc表示在所统计时间段内静风平均风速，C为所统计时间 段内观测到的静风次数，vci为在所统计时间段内，所观测到的第 i次静风划
工业企业所产生的大气污染物主要是通过气体运动来进行稀释扩散的.风向频率表示某风向吹刮次数显示着 某风向下游区所受污染的机会,频率大风向下游受污染机会多,反之则少.因此,有废气污染的工业通常布局在风玫 瑰图中最小风频侧;风速决定着污染物的扩散速度,速度越大,工厂排出废气会被很快拉长,混入的外界空气越多. 污染物的浓度越小;反之,速度越小,污染物的浓度越大,为此有废气污染的工业宜布置在风速大的方位。
5.2.5指北针及风向频率玫瑰图
建筑施工图
建筑施工图的有关规定——指北针与风玫瑰图
指北针符号北
细实线直径24mm 3mm 总平面图上的指北针或风向频率玫瑰图，是表明建筑物和建筑群的朝向和与风向的关系。
指北针的形状符合右图的规定，其圆的直径宜为24mm ，用细实线绘制； 指针尾部的宽度宜为3mm ，指针头部应注“北”或“N”字。
需用较大直径绘制指北针时，指针尾部的宽度宜为直径的1/8。
指北针
风玫瑰图
北 风向频率玫瑰图，简称
“风玫瑰图”，如建施总
平面图所示，风玫瑰图同
样指示正北方向。
表示常年(图中实线)和夏季
(图中虚线)的风向频率(各
城市或地区常年气象统计
平均数值)。
图形中显示的常年最高频
率风向称谓“主导风向”。
风玫瑰图有缘学习更多驾卫星ｙｇｄ３０７６或关注桃报：奉献教育（店铺）
谢谢！。
风玫瑰图的解释
风玫瑰图的解释化工工厂的选地址经常用到风玫瑰图的，这里简要转贴个文章，大家来认识下这个东西，尤其是做设计的时候很实用的。
“风玫瑰”图也叫风向频率玫瑰图，它是根据某一地区多年平均统计的各个方风向和风速的百分数值，并按一定比例绘制，一般多用八个或十六个罗盘方位表示，如图一，由于该图的形状形似玫瑰花朵，故名“风玫瑰”。
玫瑰图上所表示风的吹向（即风的来向），是指从外面吹向地区中心的方向。
风玫瑰图分为风向玫瑰图(图一)和风速玫瑰图（图二）两种，一般多用风向玫瑰图。
风向玫瑰图表示风向和风向的频率。
风向频率是在一定时间内各种风向出现的次数占所有观察次数的百分比。
根据各方向风的出现频率，以相应的比例长度，按风向从外向中心吹，描在用8个或16个方位所表示的图上，然后将各相邻方向的端点用直线连接起来，绘成一个形式宛如玫瑰的闭合折线，就是风玫瑰图。
图中线段最长者，即外面到中心的距离越大，表示风频越大，其为当地主导风向，外面到中心的距离越小，表示风频越小，其为当地最小风频，如图三中甲，①＞②，①为主导风向，②为最小风频；外面到中心的距离较大，为当地主导风向，如其主导方向相反，则为季风风向，如图三中乙图，其主导风向为东北——西南方向。
建筑物的位置朝向和当地主导风向有密切关系。
如把清洁的建筑物布置在主导风向的上风向；把污染建筑布置在主导风向的下风向，最小风频的上方向；污染大气的建筑布局在与季风风向垂直的郊外，以免受污染建筑散发的有害物的影响。
风玫瑰图是消防监督部门根据国家有关消防技术规范在开展建审工作时必不可少的工具，一般由当地气象部门提供。
它的作用主要表现在二个方面：一、在城市规划中，可以根据风玫瑰图正确确定大型易燃、可燃气体和液体贮罐、易燃、可燃材料堆场、大型可燃物品仓库以及散发可燃气体、液体蒸气的甲类生产厂房或甲类物品库房及生活区的位置。
由于这类贮罐、危险物品仓库、堆场火灾危害性较大，而且建造过程中的投资也比较大，一旦建好就将成为永久性的建筑物或构筑物，因此，在审核图纸时，一定要根据风玫瑰选择恰当的位置，避开城市市区和居民生活区及重要的工矿企业，防止遗留先天性的重大火险隐患。
风玫瑰图
风玫瑰”图也叫风向频率玫瑰图，它是根据某一地区多年平均统计的各个方风向和风速的百分数值，并按一定比例绘制，一般多用八个或十六个罗盘方位表示，如图一，由于该图的形状形似玫瑰花朵，故名“风玫瑰”。
玫瑰图上所表示风的吹向（即风的来向），是指从外面吹向地区中心的方向。
风玫瑰图分为风向玫瑰图(图一)和风速玫瑰图（图二）两种，一般多用风向玫瑰图。
风向玫瑰图表示风向和风向的频率。
风向频率是在一定时间内各种风向出现的次数占所有观察次数的百分比。
根据各方向风的出现频率，以相应的比例长度，按风向从外向中心吹，描在用8个或16个方位所表示的图上，然后将各相邻方向的端点用直线连接起来，绘成一个形式宛如玫瑰的闭合折线，就是风玫瑰图。
图中线段最长者，即外面到中心的距离越大，表示风频越大，其为当地主导风向，外面到中心的距离越小，表示风频越小，其为当地最小风频，如图三中甲，①＞②，①为主导风向，②为最小风频；外面到中心的距离较大，为当地主导风向，如其主导方向相反，则为季风风向，如图三中乙图，其主导风向为东北——西南方向。
建筑物的位置朝向和当地主导风向有密切关系。
如把清洁的建筑物布置在主导风向的上风向；把污染建筑布置在主导风向的下风向，最小风频的上方向；污染大气的建筑布局在与季风风向垂直的郊外，以免受污染建筑散发的有害物的影响。
风玫瑰图是消防监督部门根据国家有关消防技术规范在开展建审工作时必不可少的工具，一般由当地气象部门提供。
它的作用主要表现在二个方面：一、在城市规划中，可以根据风玫瑰图正确确定大型易燃、可燃气体和液体贮罐、易燃、可燃材料堆场、大型可燃物品仓库以及散发可燃气体、液体蒸气的甲类生产厂房或甲类物品库房及生活区的位置。
由于这类贮罐、危险物品仓库、堆场火灾危害性较大，而且建造过程中的投资也比较大，一旦建好就将成为永久性的建筑物或构筑物，因此，在审核图纸时，一定要根据风玫瑰选择恰当的位置，避开城市市区和居民生活区及重要的工矿企业，防止遗留先天性的重大火险隐患。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
工程水文学第三次作业第一题风向频率玫瑰图解：原始数据及整理后如下表：风向1--3 级0.3~5.4m/s 4--5 级5.5~10.7m/s>=6级>=10.8m/s总次数>=6 级频率（%）总频率N 364 154 12 530 0.136970665 6.049538 NNE 201 139 16 356 0.182627554 4.063463 NE 389 221 61 671 0.696267549 7.658943 ENE 238 269 103 610 1.175664878 6.962675 ENE 429 390 104 923 1.187079101 10.53533 ESE 255 181 7 443 0.079899555 5.0565 SE 416 150 1 567 0.011414222 6.471864 SSE 315 181 8 504 0.091313777 5.752768 S 541 223 3 767 0.034242666 8.754708 SSW 332 168 17 517 0.194041776 5.901153 SW 525 304 32 861 0.365255108 9.827645 WSW 263 87 1 351 0.011414222 4.006392 WSW 331 109 5 445 0.057071111 5.079329 WNW 148 111 22 281 0.251112887 3.207396 NW 292 189 19 500 0.21687022 5.707111 NNW 242 128 15 385 0.171213332 4.394476 C 0 0 0 50 0 0.570711 sum 0 0 426 8761 4.862458623 100 用Matlab 编程绘图如下：凤向频率玫瑰图附1 ( Matlab程序代码)：wind=zeros(18,6);%wind 矩阵，前三列分别存储着1-3级、4-5级、》=6级风向频率原始% 资料，后三列分别存储总次数、>=6级频率(％)、总频率(％)。
win d(:,1)=[364,201,389,238,429,255,416,315,541,332,525,263,331,148,292,242,0,0];wi nd(:,2)=[154,139,221,269,390,181,150,181,223,168,304,87,109,111,189,128,0,0];wi nd(:,3)=[12,16,61,103,104,7,1,8,3,17,32,1,5,22,19,15,0,0];win d(:,4)=wi nd(:,1)+w in d(:,2)+wi nd(:,3);win d(le ngth(wi nd(:,1))-1,4)=50;win d(le ngth(w in d(:,1)),4)=sum(w in d(:,4));win d(le ngth(wi nd(:,1)),3)=sum(wi nd(:,3));win d(:,5)=wi nd(:,3)./wi nd(le ngth(wi nd(:,1)),4).*100;win d(:,6)=wi nd(:,4)./wi nd(le ngth(wi nd(:,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)];定义十六个方向0 O计算无风频率，并将其分布在各方向上矩阵存储个方位总频率值。
a=wind(:,5)';%wind2 矩阵存储<6 级风的个方位频率值wind2=[a(1:(length(wind(:,5))-2)),wind(1,5)]; wind2=wind1-wind2; %*************************************************************************polar(theta,wind2,'-');hold on; % 绘制<6 级风的频率折线图，即内层折线polar(theta,wind1,'-');hold on; % 绘制总频率折线图，即图像最外围折线polar(theta,nowind,'-');hold on; % 绘制无风频率圆for k=1:length(theta)polar([theta(k),theta(k)],[wind1(k),wind(length(wind(:,1))-1,6)],'-')hold on;end % 绘制连接圆与总频率折线的各个方向的放射线%************************************************************************%以下绘制>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;bate=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-1 l(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)A2+a1(i,j)A2);bate(i,j)=atan(a1(i,j)./l(i,j));theta1(i,j)=theta(i)-dettheta+bate(i,j);endendif wind1(i+1)<wind1(i) q=wind1(i+1);p=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-1 l(i,j)=q*cos(dettheta)+(y-j)*det(i); a1(i,j)=(p-l(i,j))*tan(alpha(i));r1(i,j)=sqrt(l(i,j)A2+a1(i,j)A2);bate(i,j)=atan(a1(i,j)./l(i,j));theta1(i,j)=theta(i)-bate(i,j);endendend % 以上为计算阴影线在外围线上的端点的r 值，即用r1 存储for i=1:(length(wind1)-1)if wind2(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-1 l(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)A2+a1(i,j)A2); bate(i,j)=atan(a1(i,j)./l(i,j)); theta2(i,j)=theta(i)-dettheta+bate(i,j);endendif wind2(i+1)<wind2(i)q=wind2(i+1);p=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-1 l(i,j)=q*cos(dettheta)+(y-j)*det(i); a1(i,j)=(p-l(i,j))*tan(alpha(i));r2(i,j)=sqrt(l(i,j)A2+a1(i,j)A2);bate(i,j)=atan(a1(i,j)./l(i,j));theta2(i,j)=theta(i)-bate(i,j);endendend % 以上计算阴影线在内维折线上的端点的r 值，即用r2 存储for i=1:16for j=1:y-1if j-(y-1)<0 polar([theta1(i,j),theta2(i,j+1)],[r1(i,j),r2(i,j+1)],'-');hold on;endendend % 绘制阴影线条第二题海浪频率玫瑰图解：原始数据及整理后如下表：0--0.60m 0.61--1.20m 1.21--1.80m 1.81--2.40m 2.41--3.00m >=3.00m 波向P(%)N 6.45 3.15 0.42 0.55 0.49 0.44 NNE 5.63 2.26 0.35 0.55 0.07 0 NE 6.25 1.99 0.68 0.14 0.14 0 ENE 10.64 0.48 0 0.14 0 0 E 2.33 0.21 0 0 0 0 ESE 2.2 0 0 0 0 0 SE 2.82 0 0 0 0 0 SSE 1.72 0 0 0 0 0 S 1.72 0 0 0 0 0 SSW 0.56 0 0 0 0 0 SW 8.24 0.14 0 0 0 0 WSW 9.61 0.07 0.07 0 0 0 WSW 0.69 0.14 0 0 0 0WNW5.15 2.06 0.07 0.14 0.070 NW 2.951.92 0.62 0.07 0 0 NNW3.094.59 1.58 0.41 0.69 0.14c 0.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.26 1.99 0.48 0.21 0 0 0 00.14 ...0.07 0.14 2.06 1.92 4.59];wave(3,:)=[0.42 0.35 0.68 0 0 0 00 0 0 00.07 ...0 0.07 0.62 1.58];0 5% 10%------ 2 ■ti -a.oomI ____ >=3,00mwave(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 00.14];n=365*4; % 统计总次数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(:,1)=Pc;for i=2:2:length(r(1,:))-2r(:,i)=r(:,i-1)+wave(:,i/2);r(:,i+1)=r(:,i);end r(:,length(r(1,:))-1)=r(:,length(r(1,:))-2);r(:,length(r(1,:)))=r(:,length(r(1,:))-1)+wave(:,length(r(1,:))/2); %*************************** deta=[0.5 0.5 1.0 1.0 1.5 1.5 2.0 2.0 2.5以上得到r 矩阵2.53.0 3.0 ];% 各波级的宽度，即各小矩形的宽dettheta=zeros(16,length(r(1,:))); % 角度△ 0for i=1:length(r(1,:))dettheta(:,i)=atan(deta(i)./r(:,i)) ;endtheta=zeros(16,12); % 各角点的0t=pi/2:-2*pi/16:-3*pi/2+2*pi/16;t=t'; for i=1:12theta(:,i)=t;endfor i=1:16for j=2:2:length(r(1,:))-2if wave(i,j/2)==0 r(i,j)=0; r(i,j+1)=0;endend endr(:,length(r(1,:))-1)=r(:,length(r(1,:))-2);for i=1:16if wave(i,length(r(1,:))/2)==0 r(i,length(r(1,:)))=0;endendfor i=1:16if wave(:,length(r(1,:))/2)~=0r(:,length(r(1,:)))=r(:,length(r(1,:))-1)+wave(:,length(r(1,:))/2);endendtheta1=theta-dettheta;theta2=theta+dettheta;% 以上计算得到0 1和e 2，为同一方位上下或左右角点的0 r1=r./cos(dettheta);r2=r1; 以上计算得到r1和r2for 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;end end for i=1:16;for j=3:2:length(theta1(1,:)); polar([theta1(i,j),theta2(i,j)],[r1(i,j),r2(i,j)],'-'); hold on;endendj=length(theta1(1,:));for i=1:16;polar([theta1(i,j),theta2(i,j)],[r1(i,j),r2(i,j)],'-');hold on;end%*************************************** theta0=pi/2:-2*pi/49:-3*pi/2;theta0=theta0';r0=zeros(49,1);for i=1:50r0(i,1)=Pc;endpolar(theta0,r0,'-');hold on; % 绘制无浪频率圆。