% 已知圆心和半径画圆：参考：https://blog.csdn.net/ZLK961543260/article/details/70216089

% 对比三种画圆方法，运算时长对比如下

% viscircle:0.22;比较快

% rectangle:0.21;比较快

% function:0.39:慢

[Type Sheet Format]=xlsfinfo('现状OD数据及其他数据.xls');

OD=xlsread('现状OD数据及其他数据.xls',Sheet{1});

center_area=xlsread('现状OD数据及其他数据.xls',Sheet{2});

x_pos=center_area(:,4); %数据第四列为x坐标(米)

y_pos=center_area(:,5); %数据第五列为y坐标(米)

area=center_area(:,3); %数据第3列为面积(平方米)：注意与x,y的坐标对应

r_all=sqrt(area/pi); %计算节点覆盖半径，为后文画圆做准备

nanbool=isnan(r_all); %数据1到4行为四个物流园区，通过r为nan，提取出四个物流园区

first_node_num=length(find(nanbool==1)); %统计物流园区个数

r_index=find(nanbool==0,1); %节点开始的行数index

r=r_all(r_index:end); %r存储节点的覆盖半径

jam_coef=center_area(r_index:end,6); %数据第六列为交通拥堵系数，只有节点才有,具体在散点图中的标注用text实现，

% 具体标注可参考https://zhidao.baidu.com/question/924100748051179779.html

% 散点图颜色可参考：https://zhidao.baidu.com/question/521200141403285445.html?qbl=relate_question_1&word=matlabscatter%B1%EA%D7%A2

% 具体标注和颜色还未完善

% 方法一:viscircle

%用viscircle画节点的覆盖范围，用scatter画物流园区和节点

x_2_pos=x_pos(r_index:end); %节点x

y_2_pos=y_pos(r_index:end); %节点y

centers=[x_2_pos,y_2_pos]; %节点圆心位置

tic

fig1=figure

colors = {'b','r','g','y','k'};

%定义colors,后文画图方便，也可以只是掉，直接在后文加入'k','g'等，具体可查询matlab画图颜色标记表示

viscircles(centers,r,'color',colors{1}); %matlab自带函数：已知圆心和半径画圆

hold on

scatter(x_2_pos,y_2_pos,'filled');

%已知x,y坐标，画散点图命令scatter(x,y),filled表示圆是实心填充

hold on %画多个圆时用hold on命令

scatter(x_pos(1:4),y_pos(1:4),'filled'); %数字4指的是物流园区的个数，均可以替换成first_node_num

%scatter(x_pos(1:first_node_num),y_pos(1:first_node_num),'filled');

toc

% % 方法二：rectangle

% % 用矩形函数rectangle 画圆,rectangle函数一次只能画一个矩形/圆；需要用到循环，与viscircle比较计算时长，tic,toc

tic

fig2=figure

%矩形函数画多个圆只能用for命令进行循环，用矩阵读取会发生错误

for i=1:length(r)

rectangle('Position',[x_2_pos(i)-r(i),y_2_pos(i)-r(i),2*r(i),2*r(i)],'Curvature',[1,1],'linewidth',1,'EdgeColor','b') %rectangle('position',[x_pos,y_pos,length,width])

end %[1,1]表示构造圆/椭圆

hold on

scatter(x_2_pos,y_2_pos,'filled');

hold on

scatter(x_pos(1:4),y_pos(1:4),'filled');

%scatter(x_pos(1:first_node_num),y_pos(1:first_node_num),'filled');

toc

% 方法三：构建function

tic

fig3=figure

for j=1:length(r)

x=x_2_pos(j);

y=y_2_pos(j);

r1=r(j);

% function [] = plot1( x,y,r1 )

theta=0:0.1:2*pi;

Circle1=x+r1*cos(theta);

Circle2=y+r1*sin(theta);

c=[123,14,52]; %color:BGR? RBG？help

plot(Circle1,Circle2,'c','linewidth',1);

hold on

end

scatter(x_2_pos,y_2_pos,'filled');

hold on

scatter(x_pos(1:4),y_pos(1:4),'filled');

toc

绘图结果：

figure1:

figure2:

figure3:

赛题具体地图：

参考资料：

1.https://blog.csdn.net/ZLK961543260/article/details/70216089

2.数据来源为2017年研究生数学建模竞赛F题赛题数据

来源：https://www.cnblogs.com/feynmania/p/12970209.html