资源描述：
《基于Matlab的风玫瑰图绘制_张立波》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、电脑编程技巧与维护基于Matlab的风玫瑰图绘制张立波（浙江省绍兴市气象局，浙江绍兴312000）摘要：为了绘制风玫瑰图，通过实例详细介绍了利用Matlab实现对风向、风速数据的自动提取、自动处理、自动绘图和自动输出，与其他方法相比，该方法通用性强，易于移植，输出图形美观，可实现对任意台站、任意时段风玫瑰图的绘制和输出。关键词：风玫瑰图；Matlab语言；数据库PlottingWindRoseBasedontheMatlabProgramZHANGLi-bo（ShaoxingMeteorologicalBureauofZhejiang,Zheji
2、angShaoxing312000,China）Abstract:Inordertoplotwindrose,theautomaticreading,automaticprocessing,automaticplotting,andautomaticoutputofthewinddirectionandwindspeeddatabasedonMatlabprogramwereamplyillustratedbyexamples.Comparedwithothermethods,thismethodhadgooduniversality,easyp
3、orting,andniceoutputgraph.Thismethodcouldalsorealizeplottingandoutputtingwindroseinanystationandanyduration.Keywords:windrose;Matlab;database风玫瑰图又称风向频率玫瑰图，其绘制在风能资源评价、%现程序的本地化气候与环境评估、城市规划设计等领域应用广泛。目前，各file_output_name='E:sxwind风玫瑰30.xls';%输出风向风速资料Excel文件名气象台站风玫瑰图的制作主要是利用Exce
4、l的雷达图功能或编sheetname='4月';%输出风向风速资料Excel表名程实现，但这些方法主要存在图形不够美观、静风频率和各graph_output_name='E:sxwind绍兴4月风玫瑰.jpg';风向平均风速在图中没有体现及通用性不强等问题。当前能%输出风玫瑰图图形名适用于任意台站、任意时段的风玫瑰图的通用绘制程序还未zhanhao='58453';zhanming='绍兴';tablename='wind';见报道，因此，编制能自动提取风向风速数据，可应用于任%初始化区站号、站名和表名意台站、任意时段的风玫瑰图的通用自动绘制
5、程序，不但方begindate='19810101';enddate='20101231';yuefen='04';%初始化起止日期和月份便各气象台站和水文部门的个性化设置，同时也将大幅提高fengxiang={'E','ENE','NE','NNE','N','NNW','NW',各应用部门风玫瑰图的制作效率。'WNW','W','WSW','SW','SSW','S','SSE','SE','ESE',1Matlab与数据库的连接'C'};数据库工具箱(DatabaseToolbox)是Matlab工具箱扩展集direction=0:pi/8
6、:15*pi/8;中的一个，它能够使Matlab与通用关系型数据库进行数据交%初始化16个方位风向及其对应弧度r1=0.5;r2=1;font_size=10;%r1为比最大半径还多出的流。DatabaseToolbox通过ODBC驱动程序建立数据源，实现%微量半径，用于标记方位、风向频率和平均风速，r2为与数据库的连接，所以数据库须在ODBC中先注册。本文所%最内圈半径，用于在其内标注静风频率，font_size为风用数据库为SQLServer2000建立的绍兴站1981～2010年30年%玫瑰图内文本字号。%赋初值部分结束气候基准期观测资料，
7、在ODBC中注册别名为“sxclimate”,xlswrite(file_output_name,{'风向频率','平均风速'},用户名为“zlb”，密码为“123456”,各定时(02、08、14和sheetname,'B1:C1');xlswrite(file_output_name,{'E';20时)风向风速资料存储于wind表中，字段quzhanhao(char'ENE';'NE';'NNE';'N';'NNW';'NW';'WNW';'W';'WSW';'SW';'SSW';'S';'SSE';'SE';'ESE';'C'},sheet
8、name,'A2:(5))和riqi(char(8))分别表示区站号和日期，dir02(char(3))A18');和vel02(sma
