После большой работы я написал функцию, которая может генерировать очень тонкую тройную фазовую диаграмму. Этот тип диаграммы в основном используется для отображения отношения между тремя переменными. Результаты этого эксперимента таковы:
Это не легко написать, этот инструмент - головная боль, чтобы написать, пожалуйста, дайте больше лайков, давайте сначала представим, как его использовать, поместите функцию инструмента в конец и предоставьте ссылки для загрузки gitee и fileexchange одновременно, если изменения кода в будущем, он будет изменен синхронно в этих двух хранилищах.
Учебный раздел
1 Создание тройной фазовой координатной области
Его можно создать одной строкой кода:
% 初始化三元相图(Init ternary plot)
STA1=STernary();
2 set Установить самые основные свойства координатной области
Функция set может задавать многие свойства области координат, такие как цвет фона Color:
% 背景色变为灰色(Set the background color to gray)
set(STA1,'Color',[.9,.9,.9]);
Например, задайте текст меток трех осей ABC A_Label, B_Label, C_Label:
% 设置标签文本(Set Label string)
set(STA1,'A_Label','String','AAAAA')
set(STA1,'B_Label','String','BBBBB')
set(STA1,'C_Label','String','CCCCC')
3 Базовый рисунок
Общий график заполнения контура рассеяния контур текста Эти функции построения графика можно использовать, но вам нужно сделать первую букву заглавной и добавить букву S, продолжить писать предыдущий код, нарисовать график рассеяния и добавить легенду:
% 绘制散点图(Draw scatter plot)
A=rand(1,20);
B=rand(1,20);
C=rand(1,20);
STA1.SScatter(A,B,C,40,'filled','CData',[0,0,0]);
legend()
4 set Установить другие свойства
Каждая ось, основная шкала, вспомогательная шкала, основная сетка, подсетка, метка оси и метка шкалы ABC могут быть установлены с помощью функции set.Это только пример, связанный с A. Просто измените букву A на Ось БК:
Следующие две настройки эквивалентны:
- установить(STA,'XXX','реквизит',…)
- STA.set('XXX','реквизит',…)
% 设置A轴颜色和粗细(Set A_Axis's 'Color' and 'LineWidth')
set(STA1,'A_Axis','Color',[0,0,.8],'LineWidth',3)
% 设置A轴主网格颜色和粗细(Set A_Grid's 'Color' and 'LineWidth')
set(STA1,'A_Grid','Color',[0,0,.8],'LineWidth',1.2)
% 设置A轴次网格颜色线形
set(STA1,'A_MinorGrid','Color',[0,0,.8],'LineStyle','-.')
% 设置A轴主刻度颜色
set(STA1,'A_Tick','Color',[0,0,.8])
% 设置A轴次刻度颜色
set(STA1,'A_MinorTick','Color',[0,0,.8])
% 设置A轴刻度标签字体和颜色
set(STA1,'A_TickLabel','Color',[0,0,.8],'FontSize',16)
% 设置A轴标签字体和颜色及位置以及文本内容
set(STA1,'A_Label','Color',[0,0,.8],'FontSize',18,'Position',[.88,.6],'String',{
'Ternary plot';'by Slandarer'})
5 Некоторые свойства, которые обычно не устанавливаются
Восстановите объект координатной области, не используйте ранее серый, сначала переместите метку в сторону:
STA2=STernary();
% 将标签移动到侧边(Move the label to the side)
STA2.label2Side();
Установите длину мажорной и минорной гаммы, установите значение гаммы/значение сетки/текст метки гаммы:
% 设置刻度长度(Set tick length)
STA2.set('TickLenth',[.035,.02])
% 设置刻度值/网格值/刻度标签
% 'GridValues'/'TickValues'
% 'MinorGridValues'/'MinorTickValues'
% 'A_TickLabelStr'/'B_TickLabelStr'/'C_TickLabelStr'
STA2.set('MinorGridValues',5:5:100)
STA2.set('MinorTickValues',2.5:2.5:100)
STA2.set('A_TickLabelStr',{
'20%','40%','60%','80%','100%'})
6 Базовый рисунок и интерпретация рисунка линейной диаграммы
Графики рисуются SPlotс помощью функций:
% 初始化三元相图(Init ternary plot)
STA=STernary();
% 绘制折线图及散点(Draw lines and scatter)
STA.SPlot([2,2],[3,0],[5,8],'LineWidth',4,'Color',[.8,0,0]);
STA.SPlot([2,7],[3,3],[5,0],'LineWidth',4,'Color',[0,.8,0]);
STA.SPlot([2,0],[3,5],[5,5],'LineWidth',4,'Color',[0,0,.8]);
STA.SScatter(2,3,5,180,'filled','CData',[1,1,1],'MarkerEdgeColor',[0,0,0],'LineWidth',2);
% 绘制文本(Draw text)
STA.SText(2+1,3,5,'(20%,30%,50%)','FontSize',16,'FontName','Cambria')
Как показано на рисунке, соотношение точек составляет 2:3:5, поэтому линия, параллельная сетке оси А, будет пересекать ось А на 20 %, а соответствующая линия, параллельная сетке оси В, будет пересекать Ось B на 30%. . . . . . При этом длина трех прямых также находится в соотношении 2:3:5.
7 Рисование пузырьковой диаграммы
Очень просто, всего несколько строк кода:
% 初始化三元相图(Init ternary plot)
STA=STernary();
% 随机构造数据(Randomly generated data)
A=rand(1,50);
B=rand(1,50);
C=rand(1,50);
sz=rand(1,50);
% 绘制气泡图(Draw bubble chart)
STA.SBubblechart(A,B,C,sz,'CData',A);
Если вы хотите изменить цвет, colormapвы можете напрямую:
colormap(bone)
8 Контурные графики и заполненные контурные графики
Четвертым параметром чертежа является количество точек деления сетки.Если вы заполните [], значение по умолчанию делится на 90 сегментов для расчета плотности ядра.Обычно контурная карта заполняет весь треугольник:
% 初始化三元相图(Init ternary plot)
STA=STernary();
% 随机构造数据(Randomly generated data)
A=rand(1,30);
B=rand(1,30);
C=rand(1,30);
STA.SContourf(A,B,C,[],15,'LineWidth',1);
Почему бы не рисовать ниже определенного значения, вы можете получить высоту рисунка и вернуть значение Z, а также настроить высоту контурной линии, которую нужно нарисовать:
% 获取核密度范围并等分(Obtain the range of ksdensity and divide it equally)
[~,Z]=STA.SContourf(A,B,C,[],15,'Visible','off');
levels=linspace(.5,max(max(Z)),15);
% 依据等值线绘制等高线及等高线填充(Draw contour and contourf by levels)
CfHdl=STA.SContourf(A,B,C,[],levels,'LineWidth',1);
Несколько легенд рисования могут быть напрямую legend(), но здесь у нас есть скрытый контур, рекомендуется напрямую помещать графические маркеры, которые вы хотите нарисовать на легенде в матрице:
- легенда([A,B,C],{'A','B','C'})
% 获取核密度范围并等分(Obtain the range of ksdensity and divide it equally)
[~,Z]=STA.SContourf(A,B,C,[],15,'Visible','off');
levels=linspace(.5,max(max(Z)),15);
% 依据等值线绘制等高线及等高线填充(Draw contour and contourf by levels)
CfHdl=STA.SContourf(A,B,C,[],levels,'LineWidth',1);
CHdl=STA.SContour(A,B,C,[],levels,'LineWidth',1,'EdgeColor',[0,0,0]);
SHdl=STA.SScatter(A,B,C,40,'filled','CData',[0,0,0]);
% 绘制图例(Draw legend)
legend([CfHdl,CHdl,SHdl],{
'contourf','contour','scatter'},'FontSize',14,'FontName','Cambria')
Измените цвет (это изображение обложки):
colormap(turbo)
9 Переменный диапазон и выпуклая оболочка
Переменные должны быть неотрицательными.Рекомендуется сначала обработать данные.Конечно, если вам лень писать функции, вы можете использовать инструмент нормализации, который идет в комплекте с этим инструментом:
- ABC=STA.SNorm(ABC)
Каждый столбец может быть сопоставлен с 0-1, минимальное значение каждого столбца сопоставляется с 0, а максимальное значение сопоставляется с 1.
- ABC=STA.SNormByLim(ABC,[Amin,Amax,Bmin,Bmax,Cmin,Cmax])
Эта функция может предварительно установить диапазон переменных.Например, если ваши данные находятся в диапазоне 0-50, но максимальное значение этого набора данных равно 45, то функция SNorm сопоставит 45 с 1, а функция SNormByLim можно сопоставить 45 с 0,9.
Отображение диапазона переменных и построение диаграммы рассеяния
% 初始化三元相图(Init ternary plot)
STA=STernary();
% 随机构造数据(Randomly generated data)
Pnt1=mvnrnd([10,30,10],[15 0 0;0 20 0;0 0 12],80);
Pnt2=mvnrnd([30,10,15],[25 0 0;0 15 0;0 0 40],80);
Pnt3=mvnrnd([25,10,35],[35 0 0;0 20 0;0 0 30],80);
% 数据归一化(Normalization)
Pnt1=STA.SNormByLim(Pnt1,[-5,50,-5,50,-5,50]);
Pnt2=STA.SNormByLim(Pnt2,[-5,50,-5,50,-5,50]);
Pnt3=STA.SNormByLim(Pnt3,[-5,50,-5,50,-5,50]);
% 绘制散点图(Draw scatter)
propCell={
80,'filled','LineWidth',1,'MarkerEdgeColor',[1 1 1]*.3,'MarkerFaceColor'};
STA.SScatter(Pnt1(:,1),Pnt1(:,2),Pnt1(:,3),propCell{
:},[0.4,0.76,0.65]);
STA.SScatter(Pnt2(:,1),Pnt2(:,2),Pnt2(:,3),propCell{
:},[1,0.55,0.38]);
STA.SScatter(Pnt3(:,1),Pnt3(:,2),Pnt3(:,3),propCell{
:},[0.55,0.63,0.8]);
Рисунок выпуклой оболочки изначально хотел, чтобы все использовали закомментированный метод письма. Подумав об этом, забудьте об этом. Закомментированный метод письма все еще можно использовать, но я встроил его:
% 已废弃写法(Obsolete method)
% % 计算凸包(Get convhull)
% [X1,Y1]=STA.ternData(Pnt1(:,1),Pnt1(:,2),Pnt1(:,3));
% [X2,Y2]=STA.ternData(Pnt2(:,1),Pnt2(:,2),Pnt2(:,3));
% [X3,Y3]=STA.ternData(Pnt3(:,1),Pnt3(:,2),Pnt3(:,3));
% [ind1,~]=convhull([X1,Y1]);
% [ind2,~]=convhull([X2,Y2]);
% [ind3,~]=convhull([X3,Y3]);
%
% % 绘制凸包(Draw convhull)
% STA.SFill(Pnt1(ind1,1),Pnt1(ind1,2),Pnt1(ind1,3),[0.4,0.76,0.65],'FaceAlpha',.2,'LineWidth',2,'EdgeColor',[0.4,0.76,0.65]./2);
% STA.SFill(Pnt2(ind2,1),Pnt2(ind2,2),Pnt2(ind2,3),[1,0.55,0.38],'FaceAlpha',.2,'LineWidth',2,'EdgeColor',[1,0.55,0.38]./2);
% STA.SFill(Pnt3(ind3,1),Pnt3(ind3,2),Pnt3(ind3,3),[0.55,0.63,0.8],'FaceAlpha',.2,'LineWidth',2,'EdgeColor',[0.55,0.63,0.8]./2);
% 绘制凸包(Draw convhull)
STA.SConvhull(Pnt1(:,1),Pnt1(:,2),Pnt1(:,3),[0.4,0.76,0.65],'FaceAlpha',.2,'LineWidth',2,'EdgeColor',[0.4,0.76,0.65]./2);
STA.SConvhull(Pnt2(:,1),Pnt2(:,2),Pnt2(:,3),[1,0.55,0.38],'FaceAlpha',.2,'LineWidth',2,'EdgeColor',[1,0.55,0.38]./2);
STA.SConvhull(Pnt3(:,1),Pnt3(:,2),Pnt3(:,3),[0.55,0.63,0.8],'FaceAlpha',.2,'LineWidth',2,'EdgeColor',[0.55,0.63,0.8]./2);
% 绘制图例(Draw legend)
legend({
'Data1','Data2','Data3','convex1','convex2','convex3'},'FontSize',14,'FontName','Cambria')
10 Используйте nan для рисования нескольких полилиний одновременно
Правильно, я также рассмотрел некоторые меры по сокращению цикла.Функция инструмента поддерживает добавление nan соответствующим образом, чтобы отключить полилинию для одновременного рисования нескольких полилиний:
% 初始化三元相图(Init ternary plot)
STA=STernary();
A=rand(1,20);
B=rand(1,20);
C=rand(1,20);
A=[A;A*0+mean(A);A.*nan];
B=[B;B*0+mean(B);B.*nan];
C=[C;C*0+mean(C);C.*nan];
% 绘制折线图(Draw Line)
STA.SScatter(A,B,C,40,'filled','CData',[0,0,0]);
% 绘制散点图(Draw scatter plot)
STA.SPlot(A(:),B(:),C(:),'LineWidth',1,'Color',[0,0,0]);
Полный код функции инструмента
classdef STernary < handle
% Copyright (c) 2023, Zhaoxu Liu / slandarer
% =========================================================================
% @author : slandarer
% gzh : slandarer随笔
% -------------------------------------------------------------------------
% Zhaoxu Liu / slandarer (2023). ternary
% (https://www.mathworks.com/matlabcentral/fileexchange/127958-ternary),
% MATLAB Central File Exchange. 检索来源 2023/4/18.
properties
ax % 三元相图所在坐标区域
BkgHdl,LabelSep=.001;Color=[1,1,1];
TriLHdl;TriRHdl
TickLenth=[.02,.01];
TickValues=20:20:100;
GridValues=20:20:100;
MinorTickValues=5:5:100;
MinorGridValues=10:10:90;
A_Axis,B_Axis,C_Axis;
A_Grid,B_Grid,C_Grid;
A_Tick,B_Tick,C_Tick;
A_Label,B_Label;C_Label;
A_MinorGrid,B_MinorGrid,C_MinorGrid;
A_MinorTick,B_MinorTick,C_MinorTick;
A_TickLabel,B_TickLabel,C_TickLabel;
% A_LabelStr='A-Axis';
% B_LabelStr='B-Axis';
% C_LabelStr='C-Axis';
A_TickLabelStr={
'20','40','60','80','100'};
B_TickLabelStr={
'20','40','60','80','100'};
C_TickLabelStr={
'20','40','60','80','100'};
end
methods
% Create Ternary Axes =====================================================
function obj=STernary(ax)
if nargin<1,ax=gca;end
obj.ax=ax;
obj.ax.NextPlot='add';
obj.ax.XLim=[0-1/40,1+1/40];
obj.ax.YLim=[0-1/40,sqrt(3)/2+sqrt(3)/80];
obj.ax.XGrid='off';
obj.ax.YGrid='off';
obj.ax.Box='off';
obj.ax.DataAspectRatio=[1,1,1];
obj.ax.XColor='none';
obj.ax.YColor='none';
% Draw Background >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
obj.BkgHdl=fill(obj.ax,[0,1,1/2],[0,0,sqrt(3)/2],[1,1,1],'EdgeColor','none');
obj.TriLHdl=fill(obj.ax,[0,1/2+1/40,0],[0,sqrt(3)/2+1/40*sqrt(3),sqrt(3)/2+1/40*sqrt(3)],[1,1,1],'EdgeColor','none');
obj.TriRHdl=fill(obj.ax,[1,1/2-1/40,1],[0,sqrt(3)/2+1/40*sqrt(3),sqrt(3)/2+1/40*sqrt(3)],[1,1,1],'EdgeColor','none');
obj.BkgHdl.Annotation.LegendInformation.IconDisplayStyle='off';
obj.TriLHdl.Annotation.LegendInformation.IconDisplayStyle='off';
obj.TriRHdl.Annotation.LegendInformation.IconDisplayStyle='off';
% Draw Minor Grid >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
obj.MinorGridValues=setdiff(obj.MinorGridValues,obj.GridValues);
AMGX=[1-obj.MinorGridValues./200;obj.MinorGridValues./200;obj.MinorGridValues.*nan];
AMGY=[obj.MinorGridValues./200.*sqrt(3);obj.MinorGridValues./200.*sqrt(3);obj.MinorGridValues.*nan];
obj.A_MinorGrid=plot(obj.ax,AMGX(:),AMGY(:),'LineWidth',.5,'Color',[0,0,0,.2],'LineStyle','-');
obj.A_MinorGrid.Annotation.LegendInformation.IconDisplayStyle='off';
%
BMGX=[1/2-obj.MinorGridValues./200;1-obj.MinorGridValues./100;obj.MinorGridValues.*nan];
BMGY=[sqrt(3)/2-obj.MinorGridValues./200.*sqrt(3);obj.MinorGridValues.*0;obj.MinorGridValues.*nan];
obj.B_MinorGrid=plot(obj.ax,BMGX(:),BMGY(:),'LineWidth',.5,'Color',[0,0,0,.2],'LineStyle','-');
obj.B_MinorGrid.Annotation.LegendInformation.IconDisplayStyle='off';
%
CMGX=[obj.MinorGridValues./100;1/2+obj.MinorGridValues./200;obj.MinorGridValues.*nan];
CMGY=[obj.MinorGridValues.*0;1/2.*sqrt(3)-obj.MinorGridValues./200.*sqrt(3);obj.MinorGridValues.*nan];
obj.C_MinorGrid=plot(obj.ax,CMGX(:),CMGY(:),'LineWidth',.5,'Color',[0,0,0,.2],'LineStyle','-');
obj.C_MinorGrid.Annotation.LegendInformation.IconDisplayStyle='off';
% Draw Grid >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
AGX=[1-obj.GridValues./200;obj.GridValues./200;obj.GridValues.*nan];
AGY=[obj.GridValues./200.*sqrt(3);obj.GridValues./200.*sqrt(3);obj.GridValues.*nan];
obj.A_Grid=plot(obj.ax,AGX(:),AGY(:),'LineWidth',.8,'Color',[0,0,0,.7],'LineStyle','-');
obj.A_Grid.Annotation.LegendInformation.IconDisplayStyle='off';
%
BGX=[1/2-obj.GridValues./200;1-obj.GridValues./100;obj.GridValues.*nan];
BGY=[sqrt(3)/2-obj.GridValues./200.*sqrt(3);obj.GridValues.*0;obj.GridValues.*nan];
obj.B_Grid=plot(obj.ax,BGX(:),BGY(:),'LineWidth',.8,'Color',[0,0,0,.7],'LineStyle','-');
obj.B_Grid.Annotation.LegendInformation.IconDisplayStyle='off';
%
CGX=[obj.GridValues./100;1/2+obj.GridValues./200;obj.GridValues.*nan];
CGY=[obj.GridValues.*0;1/2.*sqrt(3)-obj.GridValues./200.*sqrt(3);obj.GridValues.*nan];
obj.C_Grid=plot(obj.ax,CGX(:),CGY(:),'LineWidth',.8,'Color',[0,0,0,.7],'LineStyle','-');
obj.C_Grid.Annotation.LegendInformation.IconDisplayStyle='off';
% Draw Minor Tick >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
obj.MinorTickValues=setdiff(obj.MinorTickValues,obj.TickValues);
AMTX=[1-obj.MinorTickValues./200;1-obj.MinorTickValues./200+cos(0).*obj.TickLenth(2);obj.MinorTickValues.*nan];
AMTY=[obj.MinorTickValues./200.*sqrt(3);obj.MinorTickValues./200.*sqrt(3)+sin(0).*obj.TickLenth(2);obj.MinorTickValues.*nan];
obj.A_MinorTick=plot(obj.ax,AMTX(:),AMTY(:),'LineWidth',.8,'Color',[0,0,0]);
obj.A_MinorTick.Annotation.LegendInformation.IconDisplayStyle='off';
%
BMTX=[1/2-obj.MinorTickValues./200;1/2-obj.MinorTickValues./200+cos(2*pi/3).*obj.TickLenth(2);obj.MinorTickValues.*nan];help STernary
BMTY=[sqrt(3)/2-obj.MinorTickValues./200.*sqrt(3);sqrt(3)/2-obj.MinorTickValues./200.*sqrt(3)+sin(2*pi/3).*obj.TickLenth(2);obj.MinorTickValues.*nan];
obj.B_MinorTick=plot(obj.ax,BMTX(:),BMTY(:),'LineWidth',.8,'Color',[0,0,0]);
obj.B_MinorTick.Annotation.LegendInformation.IconDisplayStyle='off';
%
CMTX=[obj.MinorTickValues./100;obj.MinorTickValues./100+cos(-2*pi/3).*obj.TickLenth(2);obj.MinorTickValues.*nan];
CMTY=[obj.MinorTickValues.*0;obj.MinorTickValues.*0+sin(-2*pi/3).*obj.TickLenth(2);obj.MinorTickValues.*nan];
obj.C_MinorTick=plot(obj.ax,CMTX(:),CMTY(:),'LineWidth',.8,'Color',[0,0,0]);
obj.C_MinorTick.Annotation.LegendInformation.IconDisplayStyle='off';
% Draw Tick >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
ATX=[1-obj.TickValues./200;1-obj.TickValues./200+cos(0).*obj.TickLenth(1);obj.TickValues.*nan];
ATY=[obj.TickValues./200.*sqrt(3);obj.TickValues./200.*sqrt(3)+sin(0).*obj.TickLenth(1);obj.TickValues.*nan];
obj.A_Tick=plot(obj.ax,ATX(:),ATY(:),'LineWidth',1.2,'Color',[0,0,0]);
obj.A_Tick.Annotation.LegendInformation.IconDisplayStyle='off';
%
BTX=[1/2-obj.TickValues./200;1/2-obj.TickValues./200+cos(2*pi/3).*obj.TickLenth(1);obj.TickValues.*nan];
BTY=[sqrt(3)/2-obj.TickValues./200.*sqrt(3);sqrt(3)/2-obj.TickValues./200.*sqrt(3)+sin(2*pi/3).*obj.TickLenth(1);obj.TickValues.*nan];
obj.B_Tick=plot(obj.ax,BTX(:),BTY(:),'LineWidth',1.2,'Color',[0,0,0]);
obj.B_Tick.Annotation.LegendInformation.IconDisplayStyle='off';
%
CTX=[obj.TickValues./100;obj.TickValues./100+cos(-2*pi/3).*obj.TickLenth(1);obj.TickValues.*nan];
CTY=[obj.TickValues.*0;obj.TickValues.*0+sin(-2*pi/3).*obj.TickLenth(1);obj.TickValues.*nan];
obj.C_Tick=plot(obj.ax,CTX(:),CTY(:),'LineWidth',1.2,'Color',[0,0,0]);
obj.C_Tick.Annotation.LegendInformation.IconDisplayStyle='off';
% Draw Axis >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
obj.A_Axis=plot(obj.ax,[1,1/2],[0,sqrt(3)/2],'Color',[0,0,0],'LineWidth',1.2);
obj.B_Axis=plot(obj.ax,[1/2,0],[sqrt(3)/2,0],'Color',[0,0,0],'LineWidth',1.2);
obj.C_Axis=plot(obj.ax,[0,1],[0,0],'Color',[0,0,0],'LineWidth',1.2);
obj.A_Axis.Annotation.LegendInformation.IconDisplayStyle='off';
obj.B_Axis.Annotation.LegendInformation.IconDisplayStyle='off';
obj.C_Axis.Annotation.LegendInformation.IconDisplayStyle='off';
% Draw Tick Label >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
for i=1:length(obj.TickValues)
obj.A_TickLabel{
i}=text(obj.ax,...
1-obj.TickValues(i)./200+cos(0).*(obj.TickLenth(1)+.01),...
obj.TickValues(i)./200.*sqrt(3)+sin(0).*(obj.TickLenth(1)+.01),...
obj.A_TickLabelStr{
i},'FontName','Cambria','FontSize',12);
obj.B_TickLabel{
i}=text(obj.ax,...
1/2-obj.TickValues(i)./200+cos(2*pi/3).*(obj.TickLenth(1)+.01),...
sqrt(3)/2-obj.TickValues(i)./200.*sqrt(3)+sin(2*pi/3).*(obj.TickLenth(1)+.01),...
obj.B_TickLabelStr{
i},'FontName','Cambria','FontSize',12,...
'Rotation',-60,'HorizontalAlignment','right');
obj.C_TickLabel{
i}=text(obj.ax,...
obj.TickValues(i)./100+cos(-2*pi/3).*(obj.TickLenth(1)+.01),...
sin(-2*pi/3).*(obj.TickLenth(1)+.01),...
obj.C_TickLabelStr{
i},'FontName','Cambria','FontSize',12,...
'Rotation',60,'HorizontalAlignment','right');
end
% Draw Label >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
obj.A_Label=text(obj.ax,1/2+cos(pi/2).*obj.LabelSep,sqrt(3)/2+sin(pi/2).*obj.LabelSep,'A-Axis',...
'FontName','Cambria','FontSize',16,'HorizontalAlignment','center','VerticalAlignment','bottom');
obj.B_Label=text(obj.ax,cos(-5*pi/6).*obj.LabelSep,sin(-5*pi/6).*obj.LabelSep,'B-Axis',...
'FontName','Cambria','FontSize',16,'HorizontalAlignment','right','VerticalAlignment','top');
obj.C_Label=text(obj.ax,1+cos(-pi/6).*obj.LabelSep,sin(-pi/6).*obj.LabelSep,'C-Axis',...
'FontName','Cambria','FontSize',16,'HorizontalAlignment','left','VerticalAlignment','top');
end
% Copyright (c) 2023, Zhaoxu Liu / slandarer
% =========================================================================
% @author : slandarer
% gzh : slandarer随笔
% -------------------------------------------------------------------------
% Zhaoxu Liu / slandarer (2023). ternary
% (https://www.mathworks.com/matlabcentral/fileexchange/127958-ternary),
% MATLAB Central File Exchange. 检索来源 2023/4/18.
% Set Ternary Axes ========================================================
% function setColor(obj,color)
% set(obj.BkgHdl,'FaceColor',color);
% end
% -------------------------------------------------------------------------
function set(obj,hdlName,varargin)
if isa(obj.(hdlName),'double')||(length(obj.(hdlName))>1&&(isa(obj.(hdlName){
1},'char')||isa(obj.(hdlName){
1},'string')))
switch true
case isequal(hdlName,'Color')
obj.Color=varargin{
:};
set(obj.BkgHdl,'FaceColor',varargin{
:});
case isequal(hdlName,'TickLenth')
obj.TickLenth=varargin{
:};
case isequal(hdlName(end-2:end),'ues')
obj.(hdlName)=varargin{
:};
case isequal(hdlName(end-2:end),'Str')
obj.(hdlName)=varargin{
:};
end
% Renew Minor Grid >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
obj.MinorGridValues=setdiff(obj.MinorGridValues,obj.GridValues);
AMGX=[1-obj.MinorGridValues./200;obj.MinorGridValues./200;obj.MinorGridValues.*nan];
AMGY=[obj.MinorGridValues./200.*sqrt(3);obj.MinorGridValues./200.*sqrt(3);obj.MinorGridValues.*nan];
set(obj.A_MinorGrid,'XData',AMGX(:),'YData',AMGY(:));
%
BMGX=[1/2-obj.MinorGridValues./200;1-obj.MinorGridValues./100;obj.MinorGridValues.*nan];
BMGY=[sqrt(3)/2-obj.MinorGridValues./200.*sqrt(3);obj.MinorGridValues.*0;obj.MinorGridValues.*nan];
set(obj.B_MinorGrid,'XData',BMGX(:),'YData',BMGY(:));
%
CMGX=[obj.MinorGridValues./100;1/2+obj.MinorGridValues./200;obj.MinorGridValues.*nan];
CMGY=[obj.MinorGridValues.*0;1/2.*sqrt(3)-obj.MinorGridValues./200.*sqrt(3);obj.MinorGridValues.*nan];
set(obj.C_MinorGrid,'XData',CMGX(:),'YData',CMGY(:));
% Renew Grid >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
AGX=[1-obj.GridValues./200;obj.GridValues./200;obj.GridValues.*nan];
AGY=[obj.GridValues./200.*sqrt(3);obj.GridValues./200.*sqrt(3);obj.GridValues.*nan];
set(obj.A_Grid,'XData',AGX(:),'YData',AGY(:));
%
BGX=[1/2-obj.GridValues./200;1-obj.GridValues./100;obj.GridValues.*nan];
BGY=[sqrt(3)/2-obj.GridValues./200.*sqrt(3);obj.GridValues.*0;obj.GridValues.*nan];
set(obj.B_Grid,'XData',BGX(:),'YData',BGY(:));
%
CGX=[obj.GridValues./100;1/2+obj.GridValues./200;obj.GridValues.*nan];
CGY=[obj.GridValues.*0;1/2.*sqrt(3)-obj.GridValues./200.*sqrt(3);obj.GridValues.*nan];
set(obj.C_Grid,'XData',CGX(:),'YData',CGY(:));
% Renew Minor Tick >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
obj.MinorTickValues=setdiff(obj.MinorTickValues,obj.TickValues);
AMTX=[1-obj.MinorTickValues./200;1-obj.MinorTickValues./200+cos(0).*obj.TickLenth(2);obj.MinorTickValues.*nan];
AMTY=[obj.MinorTickValues./200.*sqrt(3);obj.MinorTickValues./200.*sqrt(3)+sin(0).*obj.TickLenth(2);obj.MinorTickValues.*nan];
set(obj.A_MinorTick,'XData',AMTX(:),'YData',AMTY(:));
%
BMTX=[1/2-obj.MinorTickValues./200;1/2-obj.MinorTickValues./200+cos(2*pi/3).*obj.TickLenth(2);obj.MinorTickValues.*nan];
BMTY=[sqrt(3)/2-obj.MinorTickValues./200.*sqrt(3);sqrt(3)/2-obj.MinorTickValues./200.*sqrt(3)+sin(2*pi/3).*obj.TickLenth(2);obj.MinorTickValues.*nan];
set(obj.B_MinorTick,'XData',BMTX(:),'YData',BMTY(:));
%
CMTX=[obj.MinorTickValues./100;obj.MinorTickValues./100+cos(-2*pi/3).*obj.TickLenth(2);obj.MinorTickValues.*nan];
CMTY=[obj.MinorTickValues.*0;obj.MinorTickValues.*0+sin(-2*pi/3).*obj.TickLenth(2);obj.MinorTickValues.*nan];
set(obj.C_MinorTick,'XData',CMTX(:),'YData',CMTY(:));
% Renew Tick >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
ATX=[1-obj.TickValues./200;1-obj.TickValues./200+cos(0).*obj.TickLenth(1);obj.TickValues.*nan];
ATY=[obj.TickValues./200.*sqrt(3);obj.TickValues./200.*sqrt(3)+sin(0).*obj.TickLenth(1);obj.TickValues.*nan];
set(obj.A_Tick,'XData',ATX(:),'YData',ATY(:));
%
BTX=[1/2-obj.TickValues./200;1/2-obj.TickValues./200+cos(2*pi/3).*obj.TickLenth(1);obj.TickValues.*nan];
BTY=[sqrt(3)/2-obj.TickValues./200.*sqrt(3);sqrt(3)/2-obj.TickValues./200.*sqrt(3)+sin(2*pi/3).*obj.TickLenth(1);obj.TickValues.*nan];
set(obj.B_Tick,'XData',BTX(:),'YData',BTY(:));
%
CTX=[obj.TickValues./100;obj.TickValues./100+cos(-2*pi/3).*obj.TickLenth(1);obj.TickValues.*nan];
CTY=[obj.TickValues.*0;obj.TickValues.*0+sin(-2*pi/3).*obj.TickLenth(1);obj.TickValues.*nan];
set(obj.C_Tick,'XData',CTX(:),'YData',CTY(:));
% Renew Tick Label >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
for i=1:length(obj.TickValues)
set(obj.A_TickLabel{
i},'Position',...
[1-obj.TickValues(i)./200+cos(0).*(obj.TickLenth(1)+.01),...
obj.TickValues(i)./200.*sqrt(3)+sin(0).*(obj.TickLenth(1)+.01)],...
'String',obj.A_TickLabelStr{
mod(i-1,length(obj.TickValues))+1});
set(obj.B_TickLabel{
i},'Position',...
[1/2-obj.TickValues(i)./200+cos(2*pi/3).*(obj.TickLenth(1)+.01),...
sqrt(3)/2-obj.TickValues(i)./200.*sqrt(3)+sin(2*pi/3).*(obj.TickLenth(1)+.01)],...
'String',obj.B_TickLabelStr{
mod(i-1,length(obj.TickValues))+1});
set(obj.C_TickLabel{
i},'Position',...
[obj.TickValues(i)./100+cos(-2*pi/3).*(obj.TickLenth(1)+.01),...
sin(-2*pi/3).*(obj.TickLenth(1)+.01)],...
'String',obj.C_TickLabelStr{
mod(i-1,length(obj.TickValues))+1});
end
else
if length(obj.(hdlName))>1
tHdl=obj.(hdlName);
for i=1:length(tHdl)
set(tHdl{
i},varargin{
:})
end
else
set(obj.(hdlName),varargin{
:})
end
end
end
% -------------------------------------------------------------------------
function label2Side(obj)
set(obj.A_Label,'Position',[3/4+cos(pi/6).*obj.LabelSep.*90,sqrt(3)/4+sin(pi/6).*obj.LabelSep.*90],...
'HorizontalAlignment','center','VerticalAlignment','bottom','Rotation',-60)
set(obj.B_Label,'Position',[1/4+cos(5*pi/6).*obj.LabelSep.*90,sqrt(3)/4+sin(5*pi/6).*obj.LabelSep.*90],...
'HorizontalAlignment','center','VerticalAlignment','bottom','Rotation',60)
set(obj.C_Label,'Position',[1/2+cos(-pi/2).*obj.LabelSep.*90,0+sin(-pi/2).*obj.LabelSep.*90],...
'HorizontalAlignment','center','VerticalAlignment','top','Rotation',0)
end
% Copyright (c) 2023, Zhaoxu Liu / slandarer
% =========================================================================
% @author : slandarer
% gzh : slandarer随笔
% -------------------------------------------------------------------------
% Zhaoxu Liu / slandarer (2023). ternary
% (https://www.mathworks.com/matlabcentral/fileexchange/127958-ternary),
% MATLAB Central File Exchange. 检索来源 2023/4/18.
% Define Functions ========================================================
function scatterHdl=SScatter(obj,A,B,C,varargin)
[X,Y]=obj.ternData(A,B,C);
scatterHdl=scatter(obj.ax,X,Y,varargin{
:});
end
% -------------------------------------------------------------------------
function plotHdl=SPlot(obj,A,B,C,varargin)
[X,Y]=obj.ternData(A,B,C);
plotHdl=plot(obj.ax,X,Y,varargin{
:});
end
% -------------------------------------------------------------------------
function fillHdl=SFill(obj,A,B,C,varargin)
[X,Y]=obj.ternData(A,B,C);
fillHdl=fill(obj.ax,X,Y,varargin{
:});
end
% -------------------------------------------------------------------------
function convhullHdl=SConvhull(obj,A,B,C,varargin)
[X,Y]=obj.ternData(A,B,C);
[ind,~]=convhull([X,Y]);
convhullHdl=fill(obj.ax,X(ind),Y(ind),varargin{
:});
end
% -------------------------------------------------------------------------
function textHdl=SText(obj,A,B,C,varargin)
[X,Y]=obj.ternData(A,B,C);
textHdl=text(obj.ax,X,Y,varargin{
:});
end
% -------------------------------------------------------------------------
function bubblechartHdl=SBubblechart(obj,A,B,C,varargin)
[X,Y]=obj.ternData(A,B,C);
bubblechartHdl=bubblechart(obj.ax,X,Y,varargin{
:});
end
% -------------------------------------------------------------------------
function [contourHdl,Z]=SContour(obj,A,B,C,N,varargin)
if isempty(N),N=90;end
[X,Y]=obj.ternData(A,B,C);
[XX1,XX2]=meshgrid(linspace(0,1,N),linspace(0,1,N));
F=ksdensity([X,Y],[XX1(:),XX2(:)]);
Z=reshape(F,size(XX1));
Z(XX2>XX1.*sqrt(3)+.03)=nan;
Z(XX2>sqrt(3)-XX1.*sqrt(3)+.03)=nan;
hdlSet=findobj(obj.ax,'Tag','TernaryContour');
if isempty(hdlSet)
T=1;
else
T=1;
for i=1:length(hdlSet)
T=max(T,get(hdlSet(i),'UserData')+1);
end
end
contour(obj.ax,XX1,XX2,Z,varargin{
:},'Tag','TernaryContour','UserData',T);
uistack(obj.TriLHdl,'top');uistack(obj.TriRHdl,'top')
uistack(obj.A_MinorTick,'top');uistack(obj.B_MinorTick,'top');uistack(obj.C_MinorTick,'top');
uistack(obj.A_Tick,'top');uistack(obj.B_Tick,'top');uistack(obj.C_Tick,'top');
uistack(obj.A_Axis,'top');uistack(obj.B_Axis,'top');uistack(obj.C_Axis,'top');
for i=1:length(obj.A_TickLabel)
uistack(obj.A_TickLabel{
i},'top');
uistack(obj.B_TickLabel{
i},'top');
uistack(obj.C_TickLabel{
i},'top');
end
uistack(obj.A_Label,'top');uistack(obj.B_Label,'top');uistack(obj.C_Label,'top');
contourHdl=findobj(obj.ax,'Tag','TernaryContour','UserData',T);
end
% -------------------------------------------------------------------------
function [contourfHdl,Z]=SContourf(obj,A,B,C,N,varargin)
if isempty(N),N=90;end
[X,Y]=obj.ternData(A,B,C);
[XX1,XX2]=meshgrid(linspace(0,1,N),linspace(0,1,N));
F=ksdensity([X,Y],[XX1(:),XX2(:)]);
Z=reshape(F,size(XX1));
Z(XX2>XX1.*sqrt(3)+.03)=nan;
Z(XX2>sqrt(3)-XX1.*sqrt(3)+.03)=nan;
hdlSet=findobj(obj.ax,'Tag','TernaryContour');
if isempty(hdlSet)
T=1;
else
T=1;
for i=1:length(hdlSet)
T=max(T,get(hdlSet(i),'UserData')+1);
end
end
contourf(obj.ax,XX1,XX2,Z,varargin{
:},'Tag','TernaryContour','UserData',T);
uistack(obj.TriLHdl,'top');uistack(obj.TriRHdl,'top')
uistack(obj.A_MinorTick,'top');uistack(obj.B_MinorTick,'top');uistack(obj.C_MinorTick,'top');
uistack(obj.A_Tick,'top');uistack(obj.B_Tick,'top');uistack(obj.C_Tick,'top');
uistack(obj.A_Axis,'top');uistack(obj.B_Axis,'top');uistack(obj.C_Axis,'top');
for i=1:length(obj.A_TickLabel)
uistack(obj.A_TickLabel{
i},'top');
uistack(obj.B_TickLabel{
i},'top');
uistack(obj.C_TickLabel{
i},'top');
end
uistack(obj.A_Label,'top');uistack(obj.B_Label,'top');uistack(obj.C_Label,'top');
contourfHdl=findobj(obj.ax,'Tag','TernaryContour','UserData',T);
end
% -------------------------------------------------------------------------
function [X,Y]=ternData(~,A,B,C)
ABC=[A(:),B(:),C(:)];
ABC=ABC./sum(ABC,2);
B=ABC(:,2);C=ABC(:,3);
X=C+(1-C-B)./2;
Y=0+sqrt(3).*(1-C-B)./2;
end
% -------------------------------------------------------------------------
function ABC=SNorm(~,ABC)
ABC=ABC-min(ABC,[],1);
ABC=ABC./max(ABC,[],1);
end
% -------------------------------------------------------------------------
function ABC=SNormByLim(~,ABC,Lim)
ABC=ABC-Lim([1,3,5]);
ABC=ABC./(Lim([2,4,6])-Lim([1,3,5]));
end
end
% Copyright (c) 2023, Zhaoxu Liu / slandarer
% =========================================================================
% @author : slandarer
% gzh : slandarer随笔
% -------------------------------------------------------------------------
% Zhaoxu Liu / slandarer (2023). ternary
% (https://www.mathworks.com/matlabcentral/fileexchange/127958-ternary),
% MATLAB Central File Exchange. 检索来源 2023/4/18.
end
над
Написать этот код действительно непросто.Вы можете увидеть нагрузку, просто взглянув на функции инструмента (mdnice застревает при написании твитов).Я надеюсь, что всем это понравится и будет просмотрено 点赞! !
Пожалуйста, не используйте этот код в коммерческих целях без разрешения.Если вы его цитируете, вы можете указать ссылку на мой файлообменник.Можно использовать следующий формат:
Zhaoxu Liu / slandarer (2023). ternary ( https://www.mathworks.com/matlabcentral/fileexchange/127958-ternary ), MATLAB Central File Exchange. Получено с 18 апреля 2023 г.
Если вы будете перепечатывать, пожалуйста, сохраните вышеуказанную ссылку для обмена файлами и ссылку на эту статью! ! ! ! !
Инструмент можно получить по ссылке выше для обмена файлами или загрузить через следующий репозиторий gitee: