我的圖書館
留言交流

請選擇搜索范圍

含  的文章 含  的書籍 含  的隨筆

醉飲千殤v / 待分類 / matlab繪制矩陣色塊圖

分享

   

matlab繪制矩陣色塊圖

2020-07-13  醉飲千殤v

轉自：CSDN

本博客摘自天津科技大學謝中華老師一書，

           在R語言中有根據矩陣繪制色塊圖的程序，可以用于繪制相關系數矩陣圖，用豐富的顏色和形狀形象的展示矩陣元素值的大小。今天介紹的是如何用matlab繪制矩陣色塊圖。

1、繪制網格線，網格中顯示矩陣元素，顯示黑色文字。

  1. x = [1.0000    0.1302    0.0582    0.0169    0.0779    0.0833   -0.0313
  2. 0.1302    1.0000    0.0482    0.3843    0.2105   -0.2500    0.0803
  3. 0.0582    0.0482    1.0000    0.0481    0.3515   -0.0426    0.0574
  4. 0.0169    0.3843    0.0481    1.0000   -0.1819   -0.3554   -0.0020
  5. 0.0779    0.2105    0.3515   -0.1819    1.0000    0.0388    0.0186
  6. 0.0833   -0.2500   -0.0426   -0.3554    0.0388    1.0000   -0.1853
  7. -0.0313    0.0803    0.0574   -0.0020    0.0186   -0.1853    1.0000];
  8. XVarNames = {'Sleep','Age','Sex','Rel','Psy','Ner','Cha'};
  9. matrixplot(x,'FillStyle','nofill','XVarNames',XVarNames,'YVarNames',XVarNames);

效果圖如下：

2、根據矩陣元素值自動設置文字顏色。

  1. x = [1.0000    0.1302    0.0582    0.0169    0.0779    0.0833   -0.0313
  2. 0.1302    1.0000    0.0482    0.3843    0.2105   -0.2500    0.0803
  3. 0.0582    0.0482    1.0000    0.0481    0.3515   -0.0426    0.0574
  4. 0.0169    0.3843    0.0481    1.0000   -0.1819   -0.3554   -0.0020
  5. 0.0779    0.2105    0.3515   -0.1819    1.0000    0.0388    0.0186
  6. 0.0833   -0.2500   -0.0426   -0.3554    0.0388    1.0000   -0.1853
  7. -0.0313    0.0803    0.0574   -0.0020    0.0186   -0.1853    1.0000];
  8. XVarNames = {'Sleep','Age','Sex','Rel','Psy','Ner','Cha'};
  9. matrixplot(x,'FillStyle','nofill','XVarNames',XVarNames,'YVarNames',XVarNames,'TextColor','Auto','ColorBar','on');

效果圖如下：

3、繪制實值矩陣各元素對應的色塊，通過不同的“參數/參數值”控制色塊的形狀、大小、顏色等屬性。 
方形色塊，充滿方格，灰白色字體。

  1. x = [1.0000    0.1302    0.0582    0.0169    0.0779    0.0833   -0.0313
  2. 0.1302    1.0000    0.0482    0.3843    0.2105   -0.2500    0.0803
  3. 0.0582    0.0482    1.0000    0.0481    0.3515   -0.0426    0.0574
  4. 0.0169    0.3843    0.0481    1.0000   -0.1819   -0.3554   -0.0020
  5. 0.0779    0.2105    0.3515   -0.1819    1.0000    0.0388    0.0186
  6. 0.0833   -0.2500   -0.0426   -0.3554    0.0388    1.0000   -0.1853
  7. -0.0313    0.0803    0.0574   -0.0020    0.0186   -0.1853    1.0000];
  8. XVarNames = {'Sleep','Age','Sex','Rel','Psy','Ner','Cha'};
  9. matrixplot(x,'XVarNames',XVarNames,'YVarNames',XVarNames,'TextColor',[0.6,0.6,0.6],'ColorBar','on');

效果圖如下：

4、方形色塊，根據矩陣元素值自動確定色塊的大小和顏色，不顯示字體。

  1. x = [1.0000    0.1302    0.0582    0.0169    0.0779    0.0833   -0.0313
  2. 0.1302    1.0000    0.0482    0.3843    0.2105   -0.2500    0.0803
  3. 0.0582    0.0482    1.0000    0.0481    0.3515   -0.0426    0.0574
  4. 0.0169    0.3843    0.0481    1.0000   -0.1819   -0.3554   -0.0020
  5. 0.0779    0.2105    0.3515   -0.1819    1.0000    0.0388    0.0186
  6. 0.0833   -0.2500   -0.0426   -0.3554    0.0388    1.0000   -0.1853
  7. -0.0313    0.0803    0.0574   -0.0020    0.0186   -0.1853    1.0000];
  8. XVarNames = {'Sleep','Age','Sex','Rel','Psy','Ner','Cha'};
  9. matrixplot(x,'XVarNames',XVarNames,'YVarNames',XVarNames,'DisplayOpt','off','FigSize','Auto','ColorBar','on');

效果圖如下：

5、橢圓形色塊，根據矩陣元素值自動確定色塊的大小和顏色，不顯示字體。

  1. x = [1.0000    0.1302    0.0582    0.0169    0.0779    0.0833   -0.0313
  2. 0.1302    1.0000    0.0482    0.3843    0.2105   -0.2500    0.0803
  3. 0.0582    0.0482    1.0000    0.0481    0.3515   -0.0426    0.0574
  4. 0.0169    0.3843    0.0481    1.0000   -0.1819   -0.3554   -0.0020
  5. 0.0779    0.2105    0.3515   -0.1819    1.0000    0.0388    0.0186
  6. 0.0833   -0.2500   -0.0426   -0.3554    0.0388    1.0000   -0.1853
  7. -0.0313    0.0803    0.0574   -0.0020    0.0186   -0.1853    1.0000];
  8. XVarNames = {'Sleep','Age','Sex','Rel','Psy','Ner','Cha'};
  9. matrixplot(x,'XVarNames',XVarNames,'YVarNames',XVarNames,'DisplayOpt','off','FigSize','Auto','ColorBar','on','FigShape','e');

效果圖如下:

6、圓形色塊，根據矩陣元素值自動確定色塊的大小和顏色，不顯示字體，上三角形式顯示。

  1. x = [1.0000    0.1302    0.0582    0.0169    0.0779    0.0833   -0.0313
  2. 0.1302    1.0000    0.0482    0.3843    0.2105   -0.2500    0.0803
  3. 0.0582    0.0482    1.0000    0.0481    0.3515   -0.0426    0.0574
  4. 0.0169    0.3843    0.0481    1.0000   -0.1819   -0.3554   -0.0020
  5. 0.0779    0.2105    0.3515   -0.1819    1.0000    0.0388    0.0186
  6. 0.0833   -0.2500   -0.0426   -0.3554    0.0388    1.0000   -0.1853
  7. -0.0313    0.0803    0.0574   -0.0020    0.0186   -0.1853    1.0000];
  8. XVarNames = {'Sleep','Age','Sex','Rel','Psy','Ner','Cha'};
  9. matrixplot(x,'XVarNames',XVarNames,'YVarNames',XVarNames,'DisplayOpt','off','FigSize','Auto','ColorBar','on','FigShape','c','FigStyle','Triu');

效果圖如下：

7、六邊形色塊，根據矩陣元素值自動確定色塊的大小和顏色，顯示字體，下三角形式顯示

  1. x = [1.0000    0.1302    0.0582    0.0169    0.0779    0.0833   -0.0313
  2. 0.1302    1.0000    0.0482    0.3843    0.2105   -0.2500    0.0803
  3. 0.0582    0.0482    1.0000    0.0481    0.3515   -0.0426    0.0574
  4. 0.0169    0.3843    0.0481    1.0000   -0.1819   -0.3554   -0.0020
  5. 0.0779    0.2105    0.3515   -0.1819    1.0000    0.0388    0.0186
  6. 0.0833   -0.2500   -0.0426   -0.3554    0.0388    1.0000   -0.1853
  7. -0.0313    0.0803    0.0574   -0.0020    0.0186   -0.1853    1.0000];
  8. XVarNames = {'Sleep','Age','Sex','Rel','Psy','Ner','Cha'};
  9. matrixplot(x,'XVarNames',XVarNames,'YVarNames',XVarNames,'DisplayOpt','on','FigSize','Auto','ColorBar','on','FigShape','h','FigStyle','Tril');

效果圖如下：

8、表盤形色塊，根據矩陣元素值自動確定色塊的顏色，不顯示字體

  1. x = [1.0000    0.1302    0.0582    0.0169    0.0779    0.0833   -0.0313
  2. 0.1302    1.0000    0.0482    0.3843    0.2105   -0.2500    0.0803
  3. 0.0582    0.0482    1.0000    0.0481    0.3515   -0.0426    0.0574
  4. 0.0169    0.3843    0.0481    1.0000   -0.1819   -0.3554   -0.0020
  5. 0.0779    0.2105    0.3515   -0.1819    1.0000    0.0388    0.0186
  6. 0.0833   -0.2500   -0.0426   -0.3554    0.0388    1.0000   -0.1853
  7. -0.0313    0.0803    0.0574   -0.0020    0.0186   -0.1853    1.0000];
  8. XVarNames = {'Sleep','Age','Sex','Rel','Psy','Ner','Cha'};
  9. matrixplot(x,'XVarNames',XVarNames,'YVarNames',XVarNames,'DisplayOpt','off','FigSize','Full','ColorBar','on','FigShape','d');

效果圖如下：

 最后貼上matrixplot函數的源碼，源碼中有詳細的調用格式及參數說明，可以通過設置不同參數。感興趣的朋友可以自己設置不同的參數。

  1. function matrixplot(data,varargin)
  2. %   根據實值矩陣繪制色塊圖，用豐富的顏色和形狀形象的展示矩陣元素值的大小。
  3. %
  4. %   matrixplot(data) 繪制矩陣色塊圖，data為實值矩陣，每一個元素對應一個色塊，色
  5. %                    塊顏色由元素值大小決定。
  6. %
  7. %   matrixplot(data, 'PARAM1',val1, 'PARAM2',val2, ...)
  8. %          用成對出現的參數名/參數值控制色塊的各項屬性。可用的參數名/參數值如下：
  9. %          'FigShape' --- 設定色塊的形狀，其參數值為：
  10. %                'Square'  --- 方形（默認）
  11. %                'Circle'  --- 圓形
  12. %                'Ellipse' --- 橢圓形
  13. %                'Hexagon' --- 六邊形
  14. %                'Dial'    --- 表盤形
  15. %
  16. %          'FigSize' --- 設定色塊的大小，其參數值為：
  17. %                'Full'    --- 最大色塊（默認）
  18. %                'Auto'    --- 根據矩陣元素值自動確定色塊大小
  19. %
  20. %          'FigStyle' --- 設定矩陣圖樣式，其參數值為：
  21. %                'Auto'    --- 矩形矩陣圖（默認）
  22. %                'Tril'    --- 下三角矩陣圖
  23. %                'Triu'    --- 上三角矩陣圖
  24. %
  25. %          'FillStyle' --- 設定色塊填充樣式，其參數值為：
  26. %                'Fill'    --- 填充色塊內部（默認）
  27. %                'NoFill'  --- 不填充色塊內部
  28. %
  29. %          'DisplayOpt' --- 設定是否在色塊中顯示矩陣元素值，其參數值為：
  30. %                'On'      --- 顯示矩陣元素值（默認）
  31. %                'Off'     --- 不顯示矩陣元素值
  32. %
  33. %          'TextColor' --- 設定文字的顏色，其參數值為：
  34. %                表示單色的字符（'r','g','b','y','m','c','w','k'）,默認為黑色
  35. %                1行3列的紅、綠、藍三元色灰度值向量（[r,g,b]）
  36. %                'Auto'    --- 根據矩陣元素值自動確定文字顏色
  37. %
  38. %          'XVarNames' --- 設定X軸方向需要顯示的變量名（默認為X1,X2,...），其參數值為：
  39. %                字符串矩陣或字符串元胞數組，若為字符串矩陣，其行數應與data的列數相同
  40. %                若為字符串元胞數組，其長度應與data的列數相同。
  41. %
  42. %          'YVarNames' --- 設定Y軸方向需要顯示的變量名（默認為Y1,Y2,...），其參數值為：
  43. %                字符串矩陣或字符串元胞數組，若為字符串矩陣，其行數應與data的行數相同
  44. %                若為字符串元胞數組，其長度應與data的行數相同。
  45. %
  46. %          'ColorBar' --- 設定是否顯示顏色條，其參數值為：
  47. %                'On'      --- 顯示顏色條
  48. %                'Off'     --- 不顯示顏色條（默認）
  49. %
  50. %          'Grid' --- 設定是否顯示網格線，其參數值為：
  51. %                'On'      --- 顯示網格線（默認）
  52. %                'Off'     --- 不顯示網格線
  53. %
  54. %   Example:
  55. %   x = [1,-0.2,0.3,0.8,-0.5
  56. %        -0.2,1,0.6,-0.7,0.2
  57. %         0.3,0.6,1,0.5,-0.3
  58. %         0.8,-0.7,0.5,1,0.7
  59. %        -0.5,0.2,-0.3,0.7,1];
  60. %   matrixplot(x);
  61. %   matrixplot(x,'DisplayOpt','off');
  62. %   matrixplot(x,'FillStyle','nofill','TextColor','Auto');
  63. %   matrixplot(x,'TextColor',[0.7,0.7,0.7],'FigShap','s','FigSize','Auto','ColorBar','on');
  64. %   matrixplot(x,'TextColor','k','FigShap','d','FigSize','Full','ColorBar','on','FigStyle','Triu');
  65. %   XVarNames = {'xiezhh','heping','keda','tust','tianjin'};
  66. %   matrixplot(x,'FigShap','e','FigSize','Auto','ColorBar','on','XVarNames',XVarNames,'YVarNames',XVarNames);
  67. % 對第一個輸入參數類型進行判斷
  68. if ~ismatrix(data) || ~isreal(data)
  69. error('輸入參數類型不匹配：第一個輸入參數應為實值矩陣');
  70. end
  72. % 解析成對出現的參數名/參數值
  73. [FigShape,FigSize,FigStyle,FillStyle,DisplayOpt,TextColor,XVarNames,...
  74. YVarNames,ColorBar,GridOpt] = parseInputs(varargin{:});
  76. % 產生網格數據
  77. [m,n] = size(data);
  78. [x,y] = meshgrid(0:n,0:m);
  79. data = data(:);
  80. maxdata = nanmax(data);
  81. mindata = nanmin(data);
  82. rangedata = maxdata - mindata;
  83. if isnan(rangedata)
  84. warning('MATLAB:warning1','請檢查您輸入的矩陣是否合適！');
  85. return;
  86. end
  87. z = zeros(size(x))+0.2;
  88. sx = x(1:end-1,1:end-1)+0.5;
  89. sy = y(1:end-1,1:end-1)+0.5;
  91. if strncmpi(FigStyle,'Tril',4)
  92. z(triu(ones(size(z)),2)>0) = NaN;
  93. sx(triu(ones(size(sx)),1)>0) = NaN;
  94. elseif strncmpi(FigStyle,'Triu',4)
  95. z(tril(ones(size(z)),-2)>0) = NaN;
  96. sx(tril(ones(size(sx)),-1)>0) = NaN;
  97. end
  98. sx = sx(:);
  99. sy = sy(:);
  100. id = isnan(sx) | isnan(data);
  101. sx(id) = [];
  102. sy(id) = [];
  103. data(id) = [];
  105. if isempty(XVarNames)
  106. XVarNames = strcat('X',cellstr(num2str((1:n)')));
  107. else
  108. if (iscell(XVarNames) && (numel(XVarNames) ~= n)) || (~iscell(XVarNames) && (size(XVarNames,1) ~= n))
  109. error('X軸方向變量名應為字符串矩陣或字符串元胞數組，其長度與輸入矩陣的列數相同');
  110. end
  111. end
  112. if isempty(YVarNames)
  113. YVarNames = strcat('Y',cellstr(num2str((1:m)')));
  114. else
  115. if (iscell(YVarNames) && (numel(YVarNames) ~= m)) || (~iscell(YVarNames) && (size(YVarNames,1) ~= m))
  116. error('Y軸方向變量名應為字符串矩陣或字符串元胞數組，其長度與輸入矩陣的行數相同');
  117. end
  118. end
  120. % 繪圖
  121. figure('color','w',...
  122. 'units','normalized',...
  123. 'pos',[0.289165,0.154948,0.409956,0.68099]);
  124. axes('units','normalized','pos',[0.1,0.022,0.89,0.85]);
  125. if strncmpi(GridOpt,'On',2)
  126. mesh(x,y,z,...
  127. 'EdgeColor',[0.7,0.7,0.7],...
  128. 'FaceAlpha',0,...
  129. 'LineWidth',1);   % 參考網格線
  130. end
  131. hold on;
  132. axis equal;
  133. axis([-0.1,n+0.1,-0.1,m+0.1,-0.5,0.5]);
  134. view(2);
  135. % 設置X軸和Y軸刻度位置及標簽
  136. set(gca,'Xtick',(1:n)-0.5,...
  137. 'XtickLabel',XVarNames,...
  138. 'Ytick',(1:m)-0.5,...
  139. 'YtickLabel',YVarNames,...
  140. 'XAxisLocation','top',...
  141. 'YDir','reverse',...
  142. 'Xcolor',[0.7,0.7,0.7],...
  143. 'Ycolor',[0.7,0.7,0.7],...
  144. 'TickLength',[0,0]);
  145. axis off
  147. % 繪制填充色塊
  148. if strncmpi(FillStyle,'Fill',3)
  149. MyPatch(sx',sy',data',FigShape,FigSize);
  150. end
  152. % 顯示數值文本信息
  153. if strncmpi(DisplayOpt,'On',2)
  154. str = num2str(data,'%4.2f');
  155. scale = 0.1*max(n/m,1)/(max(m,n)^0.55);
  156. if strncmpi(TextColor,'Auto',3)
  157. ColorMat = get(gcf,'ColorMap');
  158. nc = size(ColorMat,1);
  159. cid = fix(mapminmax(data',0,1)*nc)+1;
  160. cid(cid<1) = 1;
  161. cid(cid>nc) = nc;
  162. TextColor = ColorMat(cid,:);
  163. for i = 1:numel(data)
  164. text(sx(i),sy(i),0.1,str(i,:),...
  165. 'FontUnits','normalized',...
  166. 'FontSize',scale,...
  167. 'fontweight','bold',...
  168. 'HorizontalAlignment','center',...
  169. 'Color',TextColor(i,:));
  170. end
  171. else
  172. text(sx,sy,0.1*ones(size(sx)),str,...
  173. 'FontUnits','normalized',...
  174. 'FontSize',scale,...
  175. 'fontweight','bold',...
  176. 'HorizontalAlignment','center',...
  177. 'Color',TextColor);
  178. end
  179. end
  181. % 設置X軸和Y軸刻度標簽的縮進方式
  182. MyTickLabel(gca,FigStyle);
  184. % 添加顏色條
  185. if strncmpi(ColorBar,'On',2)
  186. if any(strncmpi(FigStyle,{'Auto','Triu'},4))
  187. colorbar('Location','EastOutside');
  188. else
  189. colorbar('Location','SouthOutside');
  190. end
  191. end
  192. end
  194. % ---------------------------------------------------
  195. %  調整坐標軸刻度標簽子函數
  196. % ---------------------------------------------------
  197. function MyTickLabel(ha,tag)
  199. %   根據顯示范圍自動調整坐標軸刻度標簽的函數
  200. %   ha   坐標系句柄值
  201. %   tag  調整坐標軸刻度標簽的標識字符串，可用取值如下：
  202. %        'Auto' --- 將x軸刻度標簽旋轉90度，y軸刻度標簽不作調整
  203. %        'Tril' --- 將x軸刻度標簽旋轉90度，并依次縮進，y軸刻度標簽不作調整
  204. %        'Triu' --- 將x軸刻度標簽旋轉90度，y軸刻度標簽依次縮進
  205. %   Example:
  206. %   MyTickLabel(gca,'Tril');
  207. if ~ishandle(ha)
  208. warning('MATLAB:warning2','第一個輸入參數應為坐標系句柄');
  209. return;
  210. end
  212. if ~strcmpi(get(ha,'type'),'axes')
  213. warning('MATLAB:warning3','第一個輸入參數應為坐標系句柄');
  214. return;
  215. end
  216. axes(ha);
  217. xstr = get(ha,'XTickLabel');
  218. xtick = get(ha,'XTick');
  219. xl = xlim(ha);
  220. ystr = get(ha,'YTickLabel');
  221. ytick = get(ha,'YTick');
  222. yl = ylim(ha);
  223. set(ha,'XTickLabel',[],'YTickLabel',[]);
  224. x = zeros(size(ytick)) + xl(1) - range(xl)/30;
  225. y = zeros(size(xtick)) + yl(1) - range(yl)/70;
  226. nx = numel(xtick);
  227. ny = numel(ytick);
  229. if strncmpi(tag,'Tril',4)
  230. y = y + (1:nx) - 1;
  231. elseif strncmpi(tag,'Triu',4)
  232. x = x + (1:ny) - 1;
  233. end
  235. text(xtick,y,xstr,...
  236. 'rotation',90,...
  237. 'Interpreter','none',...
  238. 'color','r',...
  239. 'HorizontalAlignment','left');
  240. text(x,ytick,ystr,...
  241. 'Interpreter','none',...
  242. 'color','r',...
  243. 'HorizontalAlignment','right');
  244. end
  246. % ---------------------------------------------------
  247. %  根據散點數據繪制3維色塊圖子函數
  248. % ---------------------------------------------------
  249. function  MyPatch(x,y,z,FigShape,FigSize)
  250. %   根據散點數據繪制3維色塊圖
  251. %   MyPatch(x,y,z,FigShape,FigSize)  x,y,z是實值數組，用來指定色塊中心點三維
  252. %          坐標。FigShape是字符串變量，用來指定色塊形狀。
  253. %          FigSize是字符串變量，用來指定色塊大小。
  254. %
  255. %   Example：
  256. %         x = rand(10,1);
  257. %         y = rand(10,1);
  258. %         z = rand(10,1);
  259. %         MyPatch(x,y,z,'s','Auto');
  260. %
  262. % 輸入參數類型判斷
  263. if nargin < 3
  264. error('至少需要三個輸入參數');
  265. end
  266. if ~isreal(x) || ~isreal(y) || ~isreal(z)
  267. error('前三個輸入應為實值數組');
  268. end
  270. n = numel(z);
  271. if numel(x) ~= n || numel(y) ~= n
  272. error('坐標應等長');
  273. end
  275. if strncmpi(FigSize,'Auto',3) && ~strncmpi(FigShape,'Ellipse',1)
  276. id = (z == 0);
  277. x(id) = [];
  278. y(id) = [];
  279. z(id) = [];
  280. end
  281. if isempty(z)
  282. return;
  283. end
  285. % 求色塊頂點坐標
  286. rab1 = ones(size(z));
  287. maxz = max(abs(z));
  288. if maxz == 0
  289. maxz = 1;
  290. end
  291. rab2 = abs(z)/maxz;
  292. if strncmpi(FigShape,'Square',1)
  293. % 方形
  294. if strncmpi(FigSize,'Full',3)
  295. r = rab1;
  296. else
  297. r = sqrt(rab2);
  298. end
  299. SquareVertices(x,y,z,r);
  300. elseif strncmpi(FigShape,'Circle',1)
  301. % 圓形
  302. if strncmpi(FigSize,'Full',3)
  303. r = 0.5*rab1;
  304. else
  305. r = 0.5*sqrt(rab2);
  306. end
  307. CircleVertices(x,y,z,r);
  308. elseif strncmpi(FigShape,'Ellipse',1)
  309. % 橢圓形
  310. a = 0.48 + rab2*(0.57-0.48);
  311. b = (1-rab2).*a;
  312. EllipseVertices(x,y,z,a,b);
  313. elseif strncmpi(FigShape,'Hexagon',1)
  314. % 六邊形
  315. if strncmpi(FigSize,'Full',3)
  316. r = 0.5*rab1;
  317. else
  318. r = 0.5*sqrt(rab2);
  319. end
  320. HexagonVertices(x,y,z,r);
  321. else
  322. % 表盤形
  323. if strncmpi(FigSize,'Full',3)
  324. r = 0.45*rab1;
  325. else
  326. r = 0.45*sqrt(rab2);
  327. end
  328. DialVertices(x,y,z,r);
  329. end
  330. end
  331. %--------------------------------------------------
  332. % 求色塊頂點坐標并繪制色塊的子函數
  333. %--------------------------------------------------
  334. function SquareVertices(x,y,z,r)
  335. % 方形
  336. hx = r/2;
  337. hy = hx;
  338. Xp = [x-hx;x-hx;x+hx;x+hx;x-hx];
  339. Yp = [y-hy;y+hy;y+hy;y-hy;y-hy];
  340. Zp = repmat(z,[5,1]);
  341. patch(Xp,Yp,Zp,'FaceColor','flat','EdgeColor','flat');
  342. end
  344. function CircleVertices(x,y,z,r)
  345. % 圓形
  346. t = linspace(0,2*pi,30)';
  347. m = numel(t);
  348. Xp = repmat(x,[m,1])+cos(t)*r;
  349. Yp = repmat(y,[m,1])+sin(t)*r;
  350. Zp = repmat(z,[m,1]);
  351. patch(Xp,Yp,Zp,'FaceColor','flat','EdgeColor','flat');
  352. end
  354. function EllipseVertices(x,y,z,a,b)
  355. % 橢圓形
  356. t = linspace(0,2*pi,30)';
  357. m = numel(t);
  358. t0 = -sign(z)*pi/4;
  359. t0 = repmat(t0,[m,1]);
  360. x0 = cos(t)*a;
  361. y0 = sin(t)*b;
  362. Xp = repmat(x,[m,1]) + x0.*cos(t0) - y0.*sin(t0);
  363. Yp = repmat(y,[m,1]) + x0.*sin(t0) + y0.*cos(t0);
  364. Zp = repmat(z,[m,1]);
  365. patch(Xp,Yp,Zp,'FaceColor','flat','EdgeColor','flat');
  366. end
  368. function HexagonVertices(x,y,z,r)
  369. % 六邊形
  370. t = linspace(0,2*pi,7)';
  371. m = numel(t);
  372. Xp = repmat(x,[m,1])+cos(t)*r;
  373. Yp = repmat(y,[m,1])+sin(t)*r;
  374. Zp = repmat(z,[m,1]);
  375. patch(Xp,Yp,Zp,'FaceColor','flat','EdgeColor','flat');
  376. end
  378. function DialVertices(x,y,z,r)
  379. % 表盤形
  380. % 繪制表盤扇形
  381. maxz = max(abs(z));
  382. t0 = z*2*pi/maxz-pi/2;
  383. t0 = cell2mat(arrayfun(@(x)linspace(-pi/2,x,30)',t0,'UniformOutput',0));
  384. m = size(t0,1);
  385. r0 = repmat(r,[m,1]);
  386. Xp = [x;repmat(x,[m,1]) + r0.*cos(t0);x];
  387. Yp = [y;repmat(y,[m,1]) + r0.*sin(t0);y];
  388. Zp = repmat(z,[m+2,1]);
  389. patch(Xp,Yp,Zp,'FaceColor','flat','EdgeColor',[0,0,0]);
  391. % 繪制表盤圓周
  392. t = linspace(0,2*pi,30)';
  393. m = numel(t);
  394. Xp = repmat(x,[m,1])+cos(t)*r;
  395. Yp = repmat(y,[m,1])+sin(t)*r;
  396. Zp = repmat(z,[m,1]);
  397. Xp = [Xp;flipud(Xp)];
  398. Yp = [Yp;flipud(Yp)];
  399. Zp = [Zp;flipud(Zp)];
  400. patch(Xp,Yp,Zp,'FaceColor','flat','EdgeColor',[0,0,0]);
  401. end
  403. %--------------------------------------------------------------------------
  404. %  解析輸入參數子函數1
  405. %--------------------------------------------------------------------------
  406. function [FigShape,FigSize,FigStyle,FillStyle,DisplayOpt,TextColor,...
  407. XVarNames,YVarNames,ColorBar,GridOpt] = parseInputs(varargin)
  409. if mod(nargin,2)~=0
  410. error('輸入參數個數不對，應為成對出現');
  411. end
  412. pnames = {'FigShape','FigSize','FigStyle','FillStyle','DisplayOpt',...
  413. 'TextColor','XVarNames','YVarNames','ColorBar','Grid'};
  414. dflts =  {'Square','Full','Auto','Fill','On','k','','','Off','On'};
  415. [FigShape,FigSize,FigStyle,FillStyle,DisplayOpt,TextColor,XVarNames,...
  416. YVarNames,ColorBar,GridOpt] = parseArgs(pnames, dflts, varargin{:});
  418. validateattributes(FigShape,{'char'},{'nonempty'},mfilename,'FigShape');
  419. validateattributes(FigSize,{'char'},{'nonempty'},mfilename,'FigSize');
  420. validateattributes(FigStyle,{'char'},{'nonempty'},mfilename,'FigStyle');
  421. validateattributes(FillStyle,{'char'},{'nonempty'},mfilename,'FillStyle');
  422. validateattributes(DisplayOpt,{'char'},{'nonempty'},mfilename,'DisplayOpt');
  423. validateattributes(TextColor,{'char','numeric'},{'nonempty'},mfilename,'TextColor');
  424. validateattributes(XVarNames,{'char','cell'},{},mfilename,'XVarNames');
  425. validateattributes(YVarNames,{'char','cell'},{},mfilename,'YVarNames');
  426. validateattributes(ColorBar,{'char'},{'nonempty'},mfilename,'ColorBar');
  427. validateattributes(GridOpt,{'char'},{'nonempty'},mfilename,'Grid');
  428. if ~any(strncmpi(FigShape,{'Square','Circle','Ellipse','Hexagon','Dial'},1))
  429. error('形狀參數只能為Square, Circle, Ellipse, Hexagon, Dial 之一');
  430. end
  431. if ~any(strncmpi(FigSize,{'Full','Auto'},3))
  432. error('圖形大小參數只能為Full, Auto 之一');
  433. end
  434. if ~any(strncmpi(FigStyle,{'Auto','Tril','Triu'},4))
  435. error('圖形樣式參數只能為Auto, Tril, Triu 之一');
  436. end
  437. if ~any(strncmpi(FillStyle,{'Fill','NoFill'},3))
  438. error('圖形填充樣式參數只能為Fill, NoFill 之一');
  439. end
  440. if ~any(strncmpi(DisplayOpt,{'On','Off'},2))
  441. error('顯示數值參數只能為On，Off 之一');
  442. end
  443. if ~any(strncmpi(ColorBar,{'On','Off'},2))
  444. error('顯示顏色條參數只能為On，Off 之一');
  445. end
  446. if ~any(strncmpi(GridOpt,{'On','Off'},2))
  447. error('顯示網格參數只能為On，Off 之一');
  448. end
  449. end
  451. %--------------------------------------------------------------------------
  452. %  解析輸入參數子函數2
  453. %--------------------------------------------------------------------------
  454. function [varargout] = parseArgs(pnames,dflts,varargin)
  455. %   Copyright 2010-2011 The MathWorks, Inc.
  456. %   $Revision: 1.1.6.2 $  $Date: 2011/05/09 01:27:26 $
  457. % Initialize some variables
  458. nparams = length(pnames);
  459. varargout = dflts;
  460. setflag = false(1,nparams);
  461. unrecog = {};
  462. nargs = length(varargin);
  464. dosetflag = nargout>nparams;
  465. dounrecog = nargout>(nparams+1);
  467. % Must have name/value pairs
  468. if mod(nargs,2)~=0
  469. m = message('stats:internal:parseArgs:WrongNumberArgs');
  470. throwAsCaller(MException(m.Identifier, '%s', getString(m)));
  471. end
  473. % Process name/value pairs
  474. for j=1:2:nargs
  475. pname = varargin{j};
  476. if ~ischar(pname)
  477. m = message('stats:internal:parseArgs:IllegalParamName');
  478. throwAsCaller(MException(m.Identifier, '%s', getString(m)));
  479. end
  481. mask = strncmpi(pname,pnames,length(pname)); % look for partial match
  482. if ~any(mask)
  483. if dounrecog
  484. % if they've asked to get back unrecognized names/values, add this
  485. % one to the list
  486. unrecog((end+1):(end+2)) = {varargin{j} varargin{j+1}};
  487. continue
  488. else
  489. % otherwise, it's an error
  490. m = message('stats:internal:parseArgs:BadParamName',pname);
  491. throwAsCaller(MException(m.Identifier, '%s', getString(m)));
  492. end
  493. elseif sum(mask)>1
  494. mask = strcmpi(pname,pnames); % use exact match to resolve ambiguity
  495. if sum(mask)~=1
  496. m = message('stats:internal:parseArgs:AmbiguousParamName',pname);
  497. throwAsCaller(MException(m.Identifier, '%s', getString(m)));
  498. end
  499. end
  500. varargout{mask} = varargin{j+1};
  501. setflag(mask) = true;
  502. end
  504. % Return extra stuff if requested
  505. if dosetflag
  506. varargout{nparams+1} = setflag;
  507. if dounrecog
  508. varargout{nparams+2} = unrecog;
  509. end
  510. end
  511. end

    本站是提供個人知識管理的網絡存儲空間，所有內容均由用戶發布，不代表本站觀點。請注意甄別內容中的聯系方式、誘導購買等信息，謹防詐騙。如發現有害或侵權內容，請點擊一鍵舉報。

    來自： 醉飲千殤v > 《待分類》

    舉報

    推一薦：發原創得獎金，“原創獎勵計劃”來了！

    0條評論

    發表

    請遵守用戶 評論公約

    類似文章 更多
    喜歡該文的人也喜歡 更多

    ×
    ×

    ￥.00

    微信或支付寶掃碼支付：

    開通即同意《個圖VIP服務協議》

    全部>>
    一级a做爰片,夜夜骑网站,黄色视频做爰视频,很色的床上视频