matlab 环形单层柱状图

matlab 环形单层柱状图
matlab 环形单层柱状图

图片

图片
【图片来源粉丝】
我给他的思路是：直接使用风玫瑰图可以画出。

rose_bar
本次我的更新和这个有些不同！是环形柱状图，可调节细节多；
只需要函数：可实现各个方面：
欢迎持续投稿，宣传自己的工作（不限代码论文等方面）！
rose_bar
使用简单：

输入数据和参数

[ph,fs] =rose_bar(datax,datay,r,R,bw1,xtickvalue);
% note: datax is x;
% datay is y;
% r is inner circle
% R is circle
% bw1 is width of bar,from 2 to inf; 2 is full; increase then thin bar;    %
% ph can control line of circle;;
% fs can control line and face color of bar

先看些结果：在解释如何做出:

默认图，所有颜色一致：
图片

可使用代码调节某一个柱状图的颜色：
例如第十四（14）个为蓝色：
set(fs(14),‘FaceColor’,‘b’,‘EdgeColor’,‘b’)
图片

可调节每一个颜色都不同：
图片

可控制中心⚪的大小
⚪为零：

图片

⚪为大：
图片

可以加单位：一个或多个或全部
图片
可以什么都不要
图片在这里插入图片描述

可以更改xtick
图片

总而言之，挺强的！

主程序：

.rtcContent { padding: 30px; } .lineNode {font-size: 12pt; font-family: "Times New Roman", Menlo, Monaco, Consolas, "Courier New", monospace; font-style: normal; font-weight: normal; }
clear;clc;close all;
% 构造数据
year = 1:14;
windspeed = 1:14;
yeartick = 2001:2014;
% yeartick = 1:14;
%
close all;
figure
set(gcf,'position',[50 50 850 850 ],'color','w')
% 使用函数
datax = year;
datay = windspeed;
r = 1;
R = 2;
bw1 = 2;
xtickvalue = yeartick;
[ph,fs] =rose_bar(datax,datay,r,R,bw1,xtickvalue);
% set(fs(14),'FaceColor','b','EdgeColor','b')
% 颜色包
load('GMT_drywet.mat')
cmap =GMT_drywet(1:floor(59/length(windspeed)):end,:);
cmap = load('colormore_21.txt');
cmap = cmap(1:20:end,:);
for i = 1:length(fs)set(fs(i),'FaceColor',cmap(i,:),'EdgeColor',cmap(i,:));
end
% note: datax is x;
% datay is y;
% r is inner circle
% R is circle
% bw1 is width of bar,from 2 to inf; 2 is full; increase then thin bar;    %
% ph can control line of circle;;
% fs can control line and face color of bar
title('环形柱状图','FontSize',12,'FontWeight','bold','FontName','time News roman','color',[0.54588199	0.40039200	0.17607801]);
export_fig('环形柱状图.jpg')函数：
.rtcContent { padding: 30px; } .lineNode {font-size: 12pt; font-family: "Times New Roman", Menlo, Monaco, Consolas, "Courier New", monospace; font-style: normal; font-weight: normal; }
function [ph,fs] = rose_bar(datax,datay,r,R,bw1,xtickvalue)
% note: datax is x;
% datay is y;
% r is inner circle
% R is circle
% bw1 is width of bar
% 数据传入
year = datax;
windspeed = datay;
R = R;
bw1 = bw1;% bw1 最小值 min value is 2; max is not limit;
r = r;
t = 0:0.01:2*pi;
windspeed = windspeed+r;
% 基底
% h(1)=plot(r*sin(t),r*cos(t),'color','k');
if R>=rxlim([-(max(windspeed)+R) max(windspeed)+R])ylim([-(max(windspeed)+R) max(windspeed)+R])
elsexlim([-(max(windspeed)+r) max(windspeed)+r])ylim([-(max(windspeed)+r) max(windspeed)+r])
end
axis equal
axis off
hold on
t1 = 0:2*pi/(length(year)):2*pi;
% scatter(r*sin(t1),r*cos(t1))
bw = 2*pi/(length(year)-1)/bw1;
for i = 1:length(t1)-1t2 = t1(i)-bw:0.01:t1(i)+bw;% plot(r*sin(t2),r*cos(t2))% hold on% plot(windspeed(i)*sin(t2),windspeed(i)*cos(t2))% hold onX = [r*sin(t2) flip((windspeed(i))*sin(t2)) ];Y = [r*cos(t2) flip((windspeed(i))*cos(t2)) ];fs(i)= fill(X,Y,'r');clear X Yhold onX= mean((windspeed(i))*sin(t2));Y = mean((windspeed(i))*cos(t2));text(X,Y,num2str(xtickvalue(i)),'FontSize',12,'FontWeight','bold','FontName','time News roman','Rotation',90,'color','b')clear X Y
end
hold on
T = 0:0.01:pi/8;% 控制标签的摆放位置
for ii = 1:R:max(windspeed)+Rif (r+ii-1)<=(max(windspeed)+R)ph(ii)=plot((r+(ii-1))*sin(t),(r+(ii-1))*cos(t),'LineStyle','--','Color','k');X = mean((r+(ii-1))*sin(T));Y = mean((r+(ii-1))*cos(T));text(X,Y,num2str(ii-1),'FontSize',12,'FontWeight','bold','FontName','time News roman','color',[0.54588199	0.40039200	0.17607801]);hold onif ii == (r+ii-1)%max(1:R:max(windspeed)+R)text(X,Y,[num2str(ii-1),' m/s'],'FontSize',12,'FontWeight','bold','FontName','time News roman','color',[0.54588199	0.40039200	0.17607801]);elsetext(X,Y,[num2str(ii-1)],'FontSize',12,'FontWeight','bold','FontName','time News roman','color',[0.54588199	0.40039200	0.17607801]);endend
end
赞赏直接公益捐赠！
图片
点分享图片
点收藏图片
点在看图片
点点赞