先赞后看，养成习惯~~
先赞后看，养成习惯~~
先赞后看，养成习惯~~

绘制效果

下面这款弦图我已经出了很久了，也陆陆续续增添了新的功能和修了一些bug：

甚至还用它做出了一些复刻，分成两组的弦图有了后就有很多人在催双向的弦图，今天附带超详细解说slandarer自行研发的双向弦图工具函数终于来啦~

详细教程

0 数据准备

数据应准备全是非负数值的方形矩阵，矩阵第i行第j列表示由类i流向类j，第j行第i列表示由类j流向类i,也就是说矩阵是非对称的，可以同时统计两个类互相的数据流动,这里构造个随机数矩阵：

dataMat=randi([0,8],[6,6]);

1 基础绘图

两行代码搞定~

dataMat=randi([0,8],[6,6]);BCC=biChordChart(dataMat);
BCC=BCC.draw();

2 方向箭头

两侧都是弧形仅仅靠颜色不容易区分流入还是流出，因此可在创建对象时将Arrow属性设置为'on'：

dataMat=randi([0,8],[6,6]);BCC=biChordChart(dataMat,'Arrow','on');
BCC=BCC.draw();

3 绘图间隙

通过Sep属性可调整绘图间隙，例如设置为特别小的1/120：

dataMat=randi([0,8],[6,6]);BCC=biChordChart(dataMat,'Arrow','on','Sep',1/120);
BCC=BCC.draw();

4 添加刻度

通过tickState函数设置显示或者隐藏刻度：

dataMat=randi([0,8],[6,6]);BCC=biChordChart(dataMat,'Arrow','on');
BCC=BCC.draw(); % 添加刻度
BCC.tickState('on')

5 修改标签

标签名字默认为C1,C2,C3,...可以通过Label属性进行修改例如：

dataMat=randi([0,8],[6,6]);% 添加标签名称
NameList={'CHORD','CHART','MADE','BY','SLANDARER','MATLAB'};
BCC=biChordChart(dataMat,'Label',NameList,'Arrow','on');
BCC=BCC.draw();

6 旋转标签

懒得集成进函数里了，请使用以下这段稍显麻烦的代码：

dataMat=randi([0,8],[6,6]);% 添加标签名称
NameList={'CHORD','CHART','MADE','BY','SLANDARER','MATLAB'};
BCC=biChordChart(dataMat,'Label',NameList,'Arrow','on');
BCC=BCC.draw();% 以下代码用来旋转标签
% The following code is used to rotate the label
textHdl=findobj(gca,'Type','Text');
for i=1:length(textHdl)if textHdl(i).Rotation<-90textHdl(i).Rotation=textHdl(i).Rotation+180;endswitch truecase textHdl(i).Rotation<0&&textHdl(i).Position(2)>0textHdl(i).Rotation=textHdl(i).Rotation+90;textHdl(i).HorizontalAlignment='left';case textHdl(i).Rotation>0&&textHdl(i).Position(2)>0textHdl(i).Rotation=textHdl(i).Rotation-90;textHdl(i).HorizontalAlignment='right';case textHdl(i).Rotation<0&&textHdl(i).Position(2)<0textHdl(i).Rotation=textHdl(i).Rotation+90;textHdl(i).HorizontalAlignment='right';case textHdl(i).Rotation>0&&textHdl(i).Position(2)<0textHdl(i).Rotation=textHdl(i).Rotation-90;textHdl(i).HorizontalAlignment='left';end
end

7 颜色的设置

可在draw绘图之前设置CData属性修改颜色，例如：

dataMat=randi([0,8],[6,6]);ColorList=[127,91,93;153,66,83;95,127,95;9,14,10;78,70,83;0,0,0]./255;
BCC=biChordChart(dataMat,'Arrow','on','CData',ColorList);
BCC=BCC.draw();

dataMat=randi([0,8],[6,6]);BCC=biChordChart(dataMat,'Arrow','on','CData',bone(9));
BCC=BCC.draw();

值得一提的是如果CData设置为空集，则会随机生成颜色：

dataMat=randi([0,8],[6,6]);BCC=biChordChart(dataMat,'Arrow','on','CData',[]);
BCC=BCC.draw();

8 弧块及弦属性设置

通过

BCC.setSquareN
BCC.setChordN

设置第n个弧块或第n类弦的属性，Patch对象具有的属性均可被设置，比如如果没提前定义颜色，可以比较麻烦的修改颜色：

dataMat=randi([0,8],[6,6]);BCC=biChordChart(dataMat,'Arrow','on');
BCC=BCC.draw();ColorList=lines(6);
for i=1:6BCC.setSquareN(i,'FaceColor',ColorList(i,:))BCC.setChordN(i,'FaceColor',ColorList(i,:))
end

使用函数：

setChordMN

函数可以单独修饰类m到类n的属性，例如找到比较大的弦并将边缘标记为红色：

dataMat=randi([0,8],[6,6]);BCC=biChordChart(dataMat,'Arrow','on');
BCC=BCC.draw();% 标记最大值弦
[m,n]=find(dataMat==max(max(dataMat)));
for i=1:length(m)BCC.setChordMN(m(i),n(i),'EdgeColor',[.8,0,0],'LineWidth',2)
end

9 字体设置

通过:

setFont

函数进行字体设置：

dataMat=randi([0,8],[6,6]);BCC=biChordChart(dataMat,'Arrow','on');
BCC=BCC.draw();% 修改字体，字号及颜色
BCC.setFont('FontName','Cambria','FontSize',30,'Color',[0,0,.8])

封面绘制

封面1

% @author : slandarer
% 公众号  : slandarer随笔
% 知乎    : slandarer
dataMat=randi([0,8],[6,6]);BCC=biChordChart(dataMat,'Arrow','on');
BCC=BCC.draw();% 添加刻度
BCC.tickState('on')% 修改字体，字号及颜色
BCC.setFont('FontName','Cambria','FontSize',17)

封面2

% @author : slandarer
% 公众号  : slandarer随笔
% 知乎    : slandarer
dataMat=randi([0,8],[6,6]);% 添加标签名称
NameList={'CHORD','CHART','MADE','BY','SLANDARER','MATLAB'};
BCC=biChordChart(dataMat,'Label',NameList,'Arrow','on','CData',bone(9));
BCC=BCC.draw();% 添加刻度
BCC.tickState('on')% 修改字体，字号及颜色
BCC.setFont('FontName','Cambria','FontSize',17,'Color',[0,0,.8])

非常短的代码就能绘制出效果不错的图！！

完整代码

classdef biChordChart < handle
% Copyright (c) 2022-2023, Zhaoxu Liu / slandarer
% =========================================================================
% @author : slandarer
% 公众号  : slandarer随笔
% 知乎    : slandarer
% -------------------------------------------------------------------------
% Zhaoxu Liu / slandarer (2023). Digraph chord chart 有向弦图
% (https://www.mathworks.com/matlabcentral/fileexchange/121043-digraph-chord-chart),
% MATLAB Central File Exchange. 检索来源 2023/4/1.
%
% =========================================================================
% 使用示例：
% -------------------------------------------------------------------------
% dataMat=randi([0,8],[6,6]);
%
% BCC=biChordChart(dataMat,'Arrow','on');
% BCC=BCC.draw();
%
% % 添加刻度
% BCC.tickState('on')
%
% % 修改字体，字号及颜色
% BCC.setFont('FontName','Cambria','FontSize',17)
% =========================================================================
% 版本更新：
% -------------------------------------------------------------------------
% # version 1.1.0
% + 增添了可调节标签半径的属性'LRadius'
%   Added attribute 'LRadius' with adjustable Label radius
% + 增添了可调节标签旋转的属性'LRotate'及函数 `labelRatato`(demo3)
%   Added attribute 'LRotate' and function `labelRatato` with adjustable Label rotate(demo3)
% + 可使用函数`tickLabelState`显示刻度标签(demo4)
%   Use function `tickLabelState` to display tick labels(demo4)propertiesaxarginList={'Label','Sep','Arrow','CData','LRadius','LRotate'}dataMat     % 数值矩阵Label={}    % 标签文本% -----------------------------------------------------------squareHdl     % 绘制方块的图形对象矩阵nameHdl       % 绘制下方文本的图形对象矩阵chordMatHdl   % 绘制弦的图形对象矩阵thetaTickHdl  % 刻度句柄RTickHdl      % 轴线句柄thetaTickLabelHdlthetaSet=[];meanThetaSet;rotationSet;thetaFullSetSep;Arrow;CData;LRadius=1.28;LRotate='off'endmethodsfunction obj=biChordChart(varargin)obj.Sep=1/10;obj.Arrow='off';obj.CData=[127,91,93;187,128,110;197,173,143;59,71,111;104,95,126;76,103,86;112,112,124;72,39,24;197,119,106;160,126,88;238,208,146]./255;if isa(varargin{1},'matlab.graphics.axis.Axes')obj.ax=varargin{1};varargin(1)=[];elseobj.ax=gca;end  obj.ax.NextPlot='add';obj.dataMat=varargin{1};varargin(1)=[];% 获取其他数据for i=1:2:(length(varargin)-1)tid=ismember(obj.arginList,varargin{i});if any(tid)obj.(obj.arginList{tid})=varargin{i+1};endend% 名称标签预设if isempty(obj.Label)||length(obj.Label)<size(obj.dataMat,1)for i=1:size(obj.dataMat,1)obj.Label{i}=['C',num2str(i)];endend% 调整不合理间隙if obj.Sep>1/10obj.Sep=1/10;end% 调整颜色数量if size(obj.CData,1)<size(obj.dataMat,1)obj.CData=[obj.CData;rand([size(obj.dataMat,1),3]).*.5+ones([size(obj.dataMat,1),3]).*.5];end% 调整对角线for i=1:size(obj.dataMat,1)obj.dataMat(i,i)=abs(obj.dataMat(i,i));end% 调整标签间距if obj.LRadius>2||obj.LRadius<1.2obj.LRadius=1.28;endhelp biChordChartendfunction obj=draw(obj)obj.ax.XLim=[-1.38,1.38];obj.ax.YLim=[-1.38,1.38];obj.ax.XTick=[];obj.ax.YTick=[];obj.ax.XColor='none';obj.ax.YColor='none';obj.ax.PlotBoxAspectRatio=[1,1,1];% 计算比例numC=size(obj.dataMat,1);ratioC1=sum(abs(obj.dataMat),2)./sum(sum(abs(obj.dataMat)));ratioC2=sum(abs(obj.dataMat),1)./sum(sum(abs(obj.dataMat)));ratioC=(ratioC1'+ratioC2)./2;ratioC=[0,ratioC];sepLen=(2*pi*obj.Sep)./numC;baseLen=2*pi*(1-obj.Sep);% 绘制方块for i=1:numCtheta1=sepLen/2+sum(ratioC(1:i))*baseLen+(i-1)*sepLen;theta2=sepLen/2+sum(ratioC(1:i+1))*baseLen+(i-1)*sepLen;theta=linspace(theta1,theta2,100);X=cos(theta);Y=sin(theta);obj.squareHdl(i)=fill([1.05.*X,1.15.*X(end:-1:1)],[1.05.*Y,1.15.*Y(end:-1:1)],...obj.CData(i,:),'EdgeColor','none');theta3=(theta1+theta2)/2;obj.meanThetaSet(i)=theta3;rotation=theta3/pi*180;if rotation>0&&rotation<180obj.nameHdl(i)=text(cos(theta3).*obj.LRadius,sin(theta3).*obj.LRadius,obj.Label{i},'FontSize',14,'FontName','Arial',...'HorizontalAlignment','center','Rotation',-(.5*pi-theta3)./pi.*180,'Tag','BiChordLabel');obj.rotationSet(i)=-(.5*pi-theta3)./pi.*180;elseobj.nameHdl(i)=text(cos(theta3).*obj.LRadius,sin(theta3).*obj.LRadius,obj.Label{i},'FontSize',14,'FontName','Arial',...'HorizontalAlignment','center','Rotation',-(1.5*pi-theta3)./pi.*180,'Tag','BiChordLabel');obj.rotationSet(i)=-(1.5*pi-theta3)./pi.*180;endobj.RTickHdl(i)=plot(cos(theta).*1.17,sin(theta).*1.17,'Color',[0,0,0],'LineWidth',.8,'Visible','off');endfor i=1:numCfor j=1:numCtheta_i_1=sepLen/2+sum(ratioC(1:i))*baseLen+(i-1)*sepLen;theta_i_2=sepLen/2+sum(ratioC(1:i+1))*baseLen+(i-1)*sepLen;theta_i_3=theta_i_1+(theta_i_2-theta_i_1).*sum(abs(obj.dataMat(:,i)))./(sum(abs(obj.dataMat(:,i)))+sum(abs(obj.dataMat(i,:))));theta_j_1=sepLen/2+sum(ratioC(1:j))*baseLen+(j-1)*sepLen;theta_j_2=sepLen/2+sum(ratioC(1:j+1))*baseLen+(j-1)*sepLen;theta_j_3=theta_j_1+(theta_j_2-theta_j_1).*sum(abs(obj.dataMat(:,j)))./(sum(abs(obj.dataMat(:,j)))+sum(abs(obj.dataMat(j,:))));ratio_i_1=obj.dataMat(i,:);ratio_i_1=[0,ratio_i_1./sum(ratio_i_1)];ratio_j_2=obj.dataMat(:,j)';ratio_j_2=[0,ratio_j_2./sum(ratio_j_2)];if truetheta1=theta_i_2+(theta_i_3-theta_i_2).*sum(ratio_i_1(1:j));theta2=theta_i_2+(theta_i_3-theta_i_2).*sum(ratio_i_1(1:j+1));theta3=theta_j_3+(theta_j_1-theta_j_3).*sum(ratio_j_2(1:i));theta4=theta_j_3+(theta_j_1-theta_j_3).*sum(ratio_j_2(1:i+1));tPnt1=[cos(theta1),sin(theta1)];tPnt2=[cos(theta2),sin(theta2)];tPnt3=[cos(theta3),sin(theta3)];tPnt4=[cos(theta4),sin(theta4)];obj.thetaSet=[obj.thetaSet;theta1;theta2;theta3;theta4];obj.thetaFullSet(i,j)=theta1;obj.thetaFullSet(i,j+1)=theta2;obj.thetaFullSet(j,i+numC)=theta3;obj.thetaFullSet(j,i+numC+1)=theta4;if strcmp(obj.Arrow,'off')% 计算贝塞尔曲线tLine1=bezierCurve([tPnt1;0,0;tPnt4],200);tLine2=bezierCurve([tPnt2;0,0;tPnt3],200);tline3=[cos(linspace(theta2,theta1,100))',sin(linspace(theta2,theta1,100))'];tline4=[cos(linspace(theta4,theta3,100))',sin(linspace(theta4,theta3,100))'];else% 计算贝塞尔曲线tLine1=bezierCurve([tPnt1;0,0;tPnt4.*.96],200);tLine2=bezierCurve([tPnt2;0,0;tPnt3.*.96],200);tline3=[cos(linspace(theta2,theta1,100))',sin(linspace(theta2,theta1,100))'];tline4=[cos(theta4).*.96,sin(theta4).*.96;cos(theta3/2+theta4/2).*.99,sin(theta3/2+theta4/2).*.99;cos(theta3).*.96,sin(theta3).*.96];endobj.chordMatHdl(i,j)=fill([tLine1(:,1);tline4(:,1);tLine2(end:-1:1,1);tline3(:,1)],...[tLine1(:,2);tline4(:,2);tLine2(end:-1:1,2);tline3(:,2)],...obj.CData(i,:),'FaceAlpha',.3,'EdgeColor','none');elseendendend% 绘制刻度线tickX=[cos(obj.thetaSet).*1.17,cos(obj.thetaSet).*1.19,nan.*obj.thetaSet].';tickY=[sin(obj.thetaSet).*1.17,sin(obj.thetaSet).*1.19,nan.*obj.thetaSet].';obj.thetaTickHdl=plot(tickX(:),tickY(:),'Color',[0,0,0],'LineWidth',.8,'Visible','off');% version 1.1.0 更新部分for i=1:numCcumsumV=[0,cumsum([obj.dataMat(i,:),obj.dataMat(:,i).'])];for j=1:(2*numC+1)rotation=obj.thetaFullSet(i,j)/pi*180;if ~isnan(obj.thetaFullSet(i,j))if rotation>90&&rotation<270rotation=rotation+180;obj.thetaTickLabelHdl(i,j)=text(cos(obj.thetaFullSet(i,j)).*1.2,sin(obj.thetaFullSet(i,j)).*1.2,num2str(cumsumV(j)),...'Rotation',rotation,'HorizontalAlignment','right','FontSize',9,'FontName','Arial','Visible','off','UserData',cumsumV(j));elseobj.thetaTickLabelHdl(i,j)=text(cos(obj.thetaFullSet(i,j)).*1.2,sin(obj.thetaFullSet(i,j)).*1.2,num2str(cumsumV(j)),...'Rotation',rotation,'FontSize',9,'FontName','Arial','Visible','off','UserData',cumsumV(j));endendendend% 贝塞尔函数function pnts=bezierCurve(pnts,N)t=linspace(0,1,N);p=size(pnts,1)-1;coe1=factorial(p)./factorial(0:p)./factorial(p:-1:0);coe2=((t).^((0:p)')).*((1-t).^((p:-1:0)'));pnts=(pnts'*(coe1'.*coe2))';endobj.labelRotate(obj.LRotate)end% -----------------------------------------------------------------% 方块属性设置function setSquareN(obj,n,varargin)set(obj.squareHdl(n),varargin{:});end% -----------------------------------------------------------------% 批量弦属性设置function setChordN(obj,n,varargin)for i=nfor j=1:size(obj.dataMat,2)set(obj.chordMatHdl(i,j),varargin{:});endendend% -----------------------------------------------------------------% 单独弦属性设置function setChordMN(obj,m,n,varargin)set(obj.chordMatHdl(m,n),varargin{:});end% -----------------------------------------------------------------% 字体设置function setFont(obj,varargin)for i=1:size(obj.dataMat,1)set(obj.nameHdl(i),varargin{:});endendfunction setTickFont(obj,varargin)for m=1:size(obj.thetaFullSet,1)for n=1:size(obj.thetaFullSet,2)if obj.thetaTickLabelHdl(m,n)set(obj.thetaTickLabelHdl(m,n),varargin{:})endendendend% version 1.1.0 更新部分% 标签文字距离设置function obj=setLabelRadius(obj,Radius)obj.LRadius=Radius;for i=1:size(obj.dataMat,1)set(obj.nameHdl(i),'Position',[cos(obj.meanThetaSet(i)),sin(obj.meanThetaSet(i))].*obj.LRadius);endend% version 1.1.0 更新部分% 标签旋转状态设置function labelRotate(obj,Rotate)obj.LRotate=Rotate;for i=1:size(obj.dataMat,1)set(obj.nameHdl(i),'HorizontalAlignment','center','Rotation',obj.rotationSet(i))endif isequal(obj.LRotate,'on')textHdl=findobj(obj.ax,'Tag','BiChordLabel');for i=1:length(textHdl)if textHdl(i).Rotation<-90textHdl(i).Rotation=textHdl(i).Rotation+180;endswitch truecase textHdl(i).Rotation<0&&textHdl(i).Position(2)>0textHdl(i).Rotation=textHdl(i).Rotation+90;textHdl(i).HorizontalAlignment='left';case textHdl(i).Rotation>0&&textHdl(i).Position(2)>0textHdl(i).Rotation=textHdl(i).Rotation-90;textHdl(i).HorizontalAlignment='right';case textHdl(i).Rotation<0&&textHdl(i).Position(2)<0textHdl(i).Rotation=textHdl(i).Rotation+90;textHdl(i).HorizontalAlignment='right';case textHdl(i).Rotation>0&&textHdl(i).Position(2)<0textHdl(i).Rotation=textHdl(i).Rotation-90;textHdl(i).HorizontalAlignment='left';endendendend% -----------------------------------------------------------------% 刻度开关function tickState(obj,state)for i=1:size(obj.dataMat,1)set(obj.RTickHdl(i),'Visible',state);endset(obj.thetaTickHdl,'Visible',state);end% version 1.1.0 更新部分function tickLabelState(obj,state)for m=1:size(obj.thetaFullSet,1)for n=1:size(obj.thetaFullSet,2)if obj.thetaTickLabelHdl(m,n)if ~(n<size(obj.thetaFullSet,2)&&abs(obj.thetaFullSet(m,n)-obj.thetaFullSet(m,n+1))<eps)set(obj.thetaTickLabelHdl(m,n),'Visible',state)endendendendendfunction setTickLabelFormat(obj,func)for m=1:size(obj.thetaFullSet,1)for n=1:size(obj.thetaFullSet,2)if obj.thetaTickLabelHdl(m,n)tStr=func(get(obj.thetaTickLabelHdl(m,n),'UserData'));set(obj.thetaTickLabelHdl(m,n),'String',tStr)endendendendend
% @author : slandarer
% 公众号  : slandarer随笔
% 知乎    : slandarer
% -------------------------------------------------------------------------
% Zhaoxu Liu / slandarer (2023). Digraph chord chart 有向弦图
% (https://www.mathworks.com/matlabcentral/fileexchange/121043-digraph-chord-chart),
% MATLAB Central File Exchange. 检索来源 2023/4/1.
end

完

这种圈圈的图的绘制代码写起来真的费脑子，希望大家该点赞的点赞，该在看的在看！！

未经允许本代码请勿作商业用途，引用的话可以引用我file exchange上的链接，可使用如下格式：

Zhaoxu Liu / slandarer (2022). Digraph chord chart 有向弦图 (https://www.mathworks.com/matlabcentral/fileexchange/121043-digraph-chord-chart), MATLAB Central File Exchange. 检索来源 2022/11/22.

若转载请保留以上file exchange链接及本文链接！！！！!

完整代码:

链接：https://pan.baidu.com/s/1a6ZRUFXRFWIOKhKgL2Uung?pwd=slan
提取码：slan

MATLAB | 全网唯一 MATLAB双向弦图(有向弦图)绘制相关推荐

MATLAB | 全网唯一,使用MATLAB绘制矩形树状图
绘制效果全网唯一这四个字我都快说腻了,请叫我绘图小天才,又双叒叕写了一个工具函数发到了MATHWORKS,矩形树状图主要用于直观展示各个元素的分类和占比. 编写不易点个赞叭~~ 基本使用需要准备一 ...
MATLAB | 全网唯一! 又双叒叕一种弦图绘制
前几天有粉丝问我这种图能不能画: 大体看了一下就是有联系的两个元素之间连线,同时将不同类用不同颜色划分,于是又双叒叕了一个工具函数: 绘制效果教程部分基础使用输入数据要求一个nxn大小的0-1矩 ...
MATLAB | 全网唯一，使用MATLAB绘制好看的韦恩图(venn)
韦恩图绘制的代码MATHWORKS上也有,不过最多只支持四组绘制,于是我参考R语言ggVennDiagram包制作了一款支持7组样本的韦恩图绘制MATLAB函数. 注意:仅代码无法运行,需要导入.ma ...
MATLAB | 全网唯一,使用MATLAB绘制精致的环形树状图
绘制效果请叫我平平无奇的绘图天才,这期一次性带来俩工具函数,分别是两层和三层的环形树状图,可以用来展示数据分类和层级: 理顺层级关系属实理得头秃,代码编写不易点个赞叭!! 由于工具函数过长,将被放在 ...
MATLAB | 全网唯一,三元相映射图绘制(三元相分级统计地图)
示例对上期的补充本期文章名叫三元相映射图绘制,但不单单讲这一点,同时会对上一篇多元映射地图做出补充https://slandarer.blog.csdn.net/article/details/1 ...
基于Matlab对费根鲍姆常数和费根鲍姆图的一些算法分析
先说明一下,本文主要有两部分内容,一部分是逻辑斯蒂映射,这部分我主要参考了费根鲍姆常数4.669-9(上) 这一篇文章,并且使用matlab编写代码实现了里面的一些内容,另外第二块就是第一个文件的5之 ...
【LSTM预测】基于matlab鲸鱼算法优化双向长短时记忆BiLSTM航空发动机寿命预测（多输入单输出）【含Matlab源码 2288期】
⛄一.LSTM和BiLSTM神经网络结构及原理介绍长短期记忆网络(Long short-term memory, LSTM)实质上可以理解为一种特殊的RNN,主要是为了解决RNN网络在长序列训练过程 ...
matlab 三维核密度图_Matlab精彩画图示例：三维网状图和三维曲面图
接上篇:Matlab精彩画图示例:绘制图上图这篇介绍两个绘制三维图的命令: mesh - 绘制三维网状图 surf - 绘制三维曲面图 matlab代码如下: 绘制三维图的matlab代码示例程序 ...
matlab 实验七,matlab 实验七　数字填图问题
实验七数字填图问题数字填图问题是数学问题的一种趣味形式．早在19世纪后半期,一些数学家就在报刊中大量使用数字填图游戏和字谜游戏等,目的是使业余爱好者也能通过简单的形式去认识.理解和琢磨深奥的数学问 ...

MATLAB | 全网唯一 MATLAB双向弦图(有向弦图)绘制