![\"\"](\"https:\/\/www.wsisp.com\/helps\/wp-content\/uploads\/2026\/01\/20260127151827-6978d743b9620.png\")
<\/p>\n\u6298\u7ebf\u56fe\u2014\u2014<\/p>\n
2D\u7d2f\u8ba1\u6298\u7ebf\u56fe<\/p>\n
3D\u5e26\u900f\u660e\u5ea6\u5207\u7247<\/p>\n
2D \u591a\u7ebf\u5bf9\u5e94\u540c\u56fe\u4f8b，\u9650\u5b9a\u8fb9\u754c<\/p>\n
2D \u53cc\u5750\u6807，\u5e26\u5149\u6655\u6548\u679c，\u8986\u76d6\u6e10\u53d8\u989c\u8272\u586b\u5145\u533a\u57df，\u5e26\u56fe\u4f8b<\/p>\n
\u67f1\u72b6\u56fe\u2014\u2014<\/p>\n
3D \u4e0d\u540c\u989c\u8272 \u67f1\u72b6\u56fe<\/p>\n
2D \u53cc\u5750\u6807 \u67f1\u72b6\u56fe\u6298\u7ebf\u56fe\u6742\u4ea4<\/p>\n
\u5176\u4ed6\u7edf\u8ba1\u5206\u5e03\u56fe\u2014\u2014<\/p>\n
\u4e91\u96e8\u56fe\u548c\u7bb1\u578b\u56fe\u7684\u6742\u4ea4\u56fe？<\/p>\n
\u534a\u5c0f\u63d0\u7434\u56fe<\/p>\n
<\/p>\n
\u81ea\u5236<\/p>\n
![\"\"](\"https:\/\/www.wsisp.com\/helps\/wp-content\/uploads\/2026\/01\/20260127151828-6978d74401e71.png\")
<\/p>\n
\u4ee3\u7801：<\/p>\n
clear; clc; close all;<\/p>\n
%% \u81ea\u5b9a\u4e49\u533a
\nyears = 1900:4:2024;
\ndata_bronze = [2,3,5,6,8,10,7,9,6,3,0,0,2,5,8,7,5,3,4,6,8,5,3,4,7,9,10,12,11,9,8,10];
\ndata_silver = [0,0,1,2,3,2,4,3,2,1,0,0,1,2,3,4,2,3,1,3,5,6,4,3,5,7,6,8,7,6,5,7];
\ndata_gold = [0,0,0,1,2,3,2,4,3,2,0,0,0,1,2,3,1,2,3,4,5,4,3,5,6,8,7,10,9,8,7,9];<\/p>\n
color_bronze = [0.95 0.75 0.5];
\ncolor_silver = [0.95 0.95 0.95];
\ncolor_gold = [.95 0.85 0.5];
\nline_color = [0.2 0.2 0.2];
\nscatter_color = [0.98 0.98 0.98];<\/p>\n
xlabel_txt = 'Time (years)';
\nylabel_txt = 'Medals';
\nfontSize = 11;
\nlabelFontSize = 13;<\/p>\n
interp_points = 200;
\nyears_interp = linspace(1900, 2024, interp_points);<\/p>\n
data_bronze_interp = interp1(years, data_bronze, years_interp, 'pchip');
\ndata_silver_interp = interp1(years, data_silver, years_interp, 'pchip');
\ndata_gold_interp = interp1(years, data_gold, years_interp, 'pchip');
\ntotal_medals = data_bronze + data_silver + data_gold;
\ntotal_interp = interp1(years, total_medals, years_interp, 'pchip');<\/p>\n
% \u8ba1\u7b97\u94f6\u724c\u7d2f\u8ba1\u548c\u91d1\u724c\u7d2f\u8ba1\u6570\u636e
\nsilver_cum = data_bronze + data_silver;
\nsilver_cum_interp = interp1(years, silver_cum, years_interp, 'pchip');<\/p>\n
%% \u7edf\u4e00\u6a21\u677f\u6267\u884c\u533a
\nfigure('Position',[150,150,1000,600],'Color','white');
\nax = gca;
\nhold(ax, 'on');<\/p>\n
% \u521b\u5efa\u6e10\u53d8\u989c\u8272
\nn_points = length(years_interp);
\ngradient_bronze = zeros(n_points, 3);
\ngradient_silver = zeros(n_points, 3);
\ngradient_gold = zeros(n_points, 3);<\/p>\n
for i = 1:n_points
\n pos_ratio = i\/n_points;<\/p>\n
% \u94dc\u724c\u6e10\u53d8（\u66f4\u5feb\u7684\u6e10\u53d8\u9891\u7387）
\n intensity_bronze = 0.7 + 0.13 * sin(pos_ratio * 7*pi);
\n gradient_bronze(i,:) = color_bronze * intensity_bronze;<\/p>\n
% \u94f6\u724c\u6e10\u53d8
\n intensity_silver = 0.75 + 0.25 * cos(pos_ratio * 5*pi);
\n gradient_silver(i,:) = color_silver * intensity_silver;
\n gradient_silver(i,:) = min(gradient_silver(i,:), 1);<\/p>\n
% \u91d1\u724c\u6e10\u53d8
\n intensity_gold = 1 + 0.05 * sin(pos_ratio * 19*pi);
\n gradient_gold(i,:) = color_gold * intensity_gold;
\nend<\/p>\n
% \u7ed8\u5236\u5806\u53e0\u533a\u57df（\u4f7f\u7528\u6e10\u53d8\u586b\u5145）
\nstack_bottom = zeros(1, n_points);
\ny_stack = [data_bronze_interp; data_silver_interp; data_gold_interp];
\nlayer_names = {'Bronze','Silver','Gold'};
\nlayer_colors = {color_bronze, color_silver, color_gold};<\/p>\n
for layer = 1:3
\n y_top = stack_bottom + y_stack(layer, :);
\n x_fill = [years_interp, fliplr(years_interp)];
\n y_fill = [y_top, fliplr(stack_bottom)];<\/p>\n
if layer == 1
\n c_data = [gradient_bronze; flipud(gradient_bronze)];
\n patch(x_fill, y_fill, 1, 'FaceVertexCData', c_data, …
\n 'FaceColor', 'interp', 'EdgeColor', 'none', …
\n 'FaceAlpha', 0.9, 'DisplayName', layer_names{layer});
\n elseif layer == 2
\n c_data = [gradient_silver; flipud(gradient_silver)];
\n patch(x_fill, y_fill, 1, 'FaceVertexCData', c_data, …
\n 'FaceColor', 'interp', 'EdgeColor', 'none', …
\n 'FaceAlpha', 0.9, 'HandleVisibility', 'off');
\n else
\n c_data = [gradient_gold; flipud(gradient_gold)];
\n patch(x_fill, y_fill, 1, 'FaceVertexCData', c_data, …
\n 'FaceColor', 'interp', 'EdgeColor', 'none', …
\n 'FaceAlpha', 0.9, 'HandleVisibility', 'off');
\n end<\/p>\n
stack_bottom = y_top;
\nend<\/p>\n
% \u7ed8\u5236\u5956\u724c\u8d8b\u52bf\u7ebf（\u4e0d\u8981\u56fe\u4f8b）
\nmedal_data = {data_bronze_interp, silver_cum_interp, total_interp};
\nline_styles = {'–', ':', '-.'};
\nline_widths = [2, 2, 2.5];<\/p>\n
for i = 1:3
\n current_data = medal_data{i};<\/p>\n
% \u5149\u6655\u6548\u679c
\n for j = 4:-1:1
\n line_w = line_widths(i) + (j-1)*0.8;
\n alpha_val = 0.10 – (j-1)*0.02;<\/p>\n
plot(ax, years_interp, current_data, 'LineStyle', line_styles{i}, …
\n 'LineWidth', line_w, 'Color', [line_color, alpha_val], …
\n 'Marker', 'none', 'HandleVisibility', 'off');
\n end<\/p>\n
% \u4e3b\u66f2\u7ebf
\n plot(ax, years_interp, current_data, 'LineStyle', line_styles{i}, …
\n 'LineWidth', line_widths(i), 'Color', line_color, …
\n 'Marker', 'none', 'HandleVisibility', 'off');
\nend<\/p>\n
% \u7ed8\u5236\u603b\u5956\u724c\u6570（\u9700\u8981\u56fe\u4f8b）
\nfor j = 5:-1:1
\n line_width = 3 + (j-1)*1.5;
\n alpha_value = 0.15 – (j-1)*0.03;<\/p>\n
plot(ax, years_interp, total_interp, 'LineStyle', '-', …
\n 'LineWidth', line_width, 'Color', [line_color, alpha_value], …
\n 'Marker', 'none', 'HandleVisibility', 'off');
\nend<\/p>\n
h_total = plot(ax, years_interp, total_interp, 'LineStyle', '-', …
\n 'LineWidth', 3, 'Color', line_color, 'Marker', 'none', …
\n 'DisplayName', 'Total Medals');<\/p>\n
scatter(ax, years, total_medals, 50, scatter_color, 'filled', 'Marker', 'o', …
\n 'MarkerEdgeColor', line_color, 'LineWidth', 1.5, 'HandleVisibility', 'off');<\/p>\n
% \u5750\u6807\u8f74\u8bbe\u7f6e
\nax.XTick = 1900:20:2020;
\nax.FontName = 'Times New Roman';
\nax.FontSize = fontSize;
\nax.XLabel.String = xlabel_txt;
\nax.XLabel.FontWeight = 'bold';
\nax.XLabel.FontSize = labelFontSize;
\nax.XLabel.FontName = 'Times New Roman';
\nax.XLim = [1900, 2024];
\nax.YLabel.String = ylabel_txt;
\nax.YLabel.FontWeight = 'bold';
\nax.YLabel.FontSize = labelFontSize;
\nax.YLabel.FontName = 'Times New Roman';
\nax.YLim = [0, 30];
\nax.TickDir = 'in'; % \u523b\u5ea6\u5411\u5185
\nax.Box = 'on';<\/p>\n
% \u521b\u5efa\u56fe\u4f8b\u53e5\u67c4（\u53ea\u663e\u793a\u603b\u5956\u724c\u548c\u586b\u5145\u533a\u57df）
\nlegend_handles = [];<\/p>\n
% \u586b\u5145\u533a\u57df\u56fe\u4f8b
\nh_bronze = patch([0,0,0,0], [0,0,0,0], color_bronze, …
\n 'EdgeColor', 'none', 'FaceAlpha', 0.9, …
\n 'DisplayName', 'Bronze');
\nlegend_handles = [legend_handles, h_bronze];<\/p>\n
h_silver = patch([0,0,0,0], [0,0,0,0], color_silver, …
\n 'EdgeColor', 'none', 'FaceAlpha', 0.9, …
\n 'DisplayName', 'Silver');
\nlegend_handles = [legend_handles, h_silver];<\/p>\n
h_gold = patch([0,0,0,0], [0,0,0,0], color_gold, …
\n 'EdgeColor', 'none', 'FaceAlpha', 0.9, …
\n 'DisplayName', 'Gold');
\nlegend_handles = [legend_handles, h_gold];<\/p>\n
% \u603b\u5956\u724c\u6570\u56fe\u4f8b
\nlegend_handles = [legend_handles, h_total];<\/p>\n
legend(legend_handles, {'Bronze','Silver','Gold','Total Medals'}, …
\n 'Location', 'northeast', 'FontSize', fontSize, 'FontName', 'Times New Roman');<\/p>\n
hold(ax, 'off');<\/p>\n
fprintf('\u7ed8\u56fe\u5b8c\u6210\\\\n'); <\/p>\n
research gate\u4e0a\u7684\u56fe，\u5f88\u50cf<\/p>\n
![\"\"](\"https:\/\/www.wsisp.com\/helps\/wp-content\/uploads\/2026\/01\/20260127151828-6978d7445804b.png\")
<\/p>\n<\/p>\n
\u53c2\u8003<\/p>\n
![\"\"](\"https:\/\/www.wsisp.com\/helps\/wp-content\/uploads\/2026\/01\/20260127151828-6978d744c0820.png\")
<\/p>\n
\u81ea\u5236<\/p>\n
![\"\"](\"https:\/\/www.wsisp.com\/helps\/wp-content\/uploads\/2026\/01\/20260127151828-6978d744f32b7.png\")
<\/p>\n
\u4ee3\u7801<\/p>\n
clear; clc; close all;<\/p>\n
%% \u81ea\u5b9a\u4e49\u533a
\nYear = [2008, 2012, 2016, 2020, 2024];
\nCountries = {'JPN', 'AUS', 'GER', 'FRA', 'GBR', 'CHN', 'USA'};<\/p>\n
Data = [
\n 2, 2, 2, 3, 3;
\n 4, 5, 5, 6, 6;
\n 6, 7, 7, 8, 8;
\n 8, 9, 9, 10, 10;
\n 10, 12, 13, 14, 15;
\n 12, 15, 18, 21, 24;
\n 9, 18, 21, 30, 33
\n];<\/p>\n
% \u989c\u8272\u77e9\u9635
\nColors = [
\n 0.85 0.94 0.83;
\n 0.68 0.85 0.65;
\n 0.52 0.78 0.52;
\n 0.35 0.70 0.40;
\n 0.95 0.87 0.65;
\n 0.88 0.68 0.45;
\n 0.82 0.52 0.30
\n];<\/p>\n
% \u6807\u7b7e\u8bbe\u7f6e
\nxlabel_txt = 'Year';
\nylabel_txt = 'Country';
\nzlabel_txt = 'Value';
\nfontSize = 11;
\nlabelFontSize = 13;<\/p>\n
%% \u7edf\u4e00\u6a21\u677f\u6267\u884c\u533a
\nfigure('Color', 'w', 'Position', [100, 100, 1000, 650]);
\nhold on; box on; grid on;<\/p>\n
nCountries = size(Data, 1);
\nnYears = length(Year);
\nyPositions = 1:nCountries;<\/p>\n
legendHandles = gobjects(nCountries, 1);<\/p>\n
for i = 1:nCountries
\n zData = Data(i, :);
\n yLevel = yPositions(i);<\/p>\n
% \u586b\u5145\u533a\u57df
\n xMesh = repmat(Year, 2, 1);
\n yMesh = ones(2, nYears) * yLevel;
\n zMesh = [zeros(1, nYears); zData];<\/p>\n
surf(xMesh, yMesh, zMesh, …
\n 'FaceColor', Colors(i,:), …
\n 'EdgeColor', 'none', …
\n 'FaceAlpha', 0.4); % \u5728\u8fd9\u91cc\u4fee\u6539\u900f\u660e\u5ea6，0\u662f\u5b8c\u5168\u900f\u660e，1\u662f\u5b8c\u5168\u4e0d\u900f\u660e<\/p>\n
% \u6298\u7ebf
\n if i == 1
\n legendHandles(i) = plot3(Year, ones(1, nYears)*yLevel, zData, …
\n 'Color', Colors(i,:)*0.7, …
\n 'LineWidth', 2.5, …
\n 'Marker', 'o', …
\n 'MarkerSize', 5, …
\n 'MarkerFaceColor', Colors(i,:), …
\n 'MarkerEdgeColor', Colors(i,:)*0.5, …
\n 'DisplayName', Countries{i});
\n else
\n plot3(Year, ones(1, nYears)*yLevel, zData, …
\n 'Color', Colors(i,:)*0.7, …
\n 'LineWidth', 2.5, …
\n 'Marker', 'o', …
\n 'MarkerSize', 5, …
\n 'MarkerFaceColor', Colors(i,:), …
\n 'MarkerEdgeColor', Colors(i,:)*0.5, …
\n 'HandleVisibility', 'off');
\n legendHandles(i) = plot3(NaN, NaN, NaN, …
\n 'Color', Colors(i,:), …
\n 'LineWidth', 2.5, …
\n 'Marker', 'o', …
\n 'MarkerSize', 5, …
\n 'MarkerFaceColor', Colors(i,:), …
\n 'MarkerEdgeColor', Colors(i,:)*0.5, …
\n 'DisplayName', Countries{i});
\n end<\/p>\n
% \u6570\u503c\u6807\u7b7e
\n for j = 1:nYears
\n text(Year(j), yLevel – 0.2, zData(j), …
\n sprintf('%d', zData(j)), …
\n 'FontName', 'Times New Roman', …
\n 'FontSize', 9, …
\n 'FontWeight', 'bold', …
\n 'Color', Colors(i,:)*0.7, …
\n 'HorizontalAlignment', 'center', …
\n 'VerticalAlignment', 'top');
\n end
\nend<\/p>\n
% \u5750\u6807\u8f74\u8bbe\u7f6e
\nxlabel(xlabel_txt, 'FontName', 'Times New Roman', 'FontSize', labelFontSize, 'FontWeight', 'bold');
\nylabel(ylabel_txt, 'FontName', 'Times New Roman', 'FontSize', labelFontSize, 'FontWeight', 'bold');
\nzlabel(zlabel_txt, 'FontName', 'Times New Roman', 'FontSize', labelFontSize, 'FontWeight', 'bold');<\/p>\n
xlim([min(Year)-0.5, max(Year)+0.5]);
\nylim([0.5, nCountries + 0.5]);
\nzlim([0, max(max(Data)) + 10]);<\/p>\n
set(gca, 'XTick', Year);
\nset(gca, 'XTickLabel', cellstr(num2str(Year')));
\nset(gca, 'YTick', yPositions, 'YTickLabel', Countries);
\nset(gca, 'FontName', 'Times New Roman', 'FontSize', fontSize, 'FontWeight', 'bold');<\/p>\n
% \u53f3\u4fa7\u603b\u8d21\u732e\u6807\u7b7e
\nfor i = 1:nCountries
\n total = sum(Data(i, :));
\n avgZ = mean(Data(i, :));<\/p>\n
text(Year(end) + 0.8, yPositions(i), avgZ, …
\n sprintf('Total: %d', total), …
\n 'FontName', 'Times New Roman', …
\n 'FontSize', fontSize, …
\n 'FontWeight', 'bold', …
\n 'Color', Colors(i,:)*0.8, …
\n 'HorizontalAlignment', 'left');
\nend<\/p>\n
% \u89c6\u89d2\u548c\u5149\u7167
\nview(45, 30);
\n% light('Position', [1 1 1], 'Style', 'infinite');
\n% lighting gouraud;
\n% \u52a0\u4e86\u5149\u7167\u80cc\u9762\u4f1a\u5f88\u9ed1<\/p>\n
% \u56fe\u4f8b
\nlegend(legendHandles, 'Location', 'northeastoutside', …
\n 'FontName', 'Times New Roman', …
\n 'FontSize', fontSize, …
\n 'FontWeight', 'bold');<\/p>\n
hold off; <\/p>\n
\u53c2\u8003\u56fe 2020-A\u9898-\u9c7c\u7fa4\u8fc1\u5f99 #2001334<\/p>\n
![\"\"](\"https:\/\/www.wsisp.com\/helps\/wp-content\/uploads\/2026\/01\/20260127151829-6978d74557ae4.png\")
<\/p>\n
![\"\"](\"https:\/\/www.wsisp.com\/helps\/wp-content\/uploads\/2026\/01\/20260127151829-6978d745a670c.png\")
<\/p>\n
\u4ee3\u7801<\/p>\n
%% ===================== \u81ea\u5b9a\u4e49\u533a（\u53ef\u4fee\u6539\u53c2\u6570）=====================
\n% \u6570\u503c\u53c2\u6570
\nsim_num = 50; % \u6a21\u62df\u66f2\u7ebf\u6570\u91cf
\nyear_start = 2020; % \u8d77\u59cb\u5e74\u4efd
\nyear_end = 2070; % \u7ed3\u675f\u5e74\u4efd
\nprofit_unit = 1e5; % \u5229\u6da6\u5355\u4f4d（10\u4e07\u7f8e\u5143）
\n% \u989c\u8272\u53c2\u6570（RGB\u77e9\u9635）
\ncolor_monte = [0.83,0.83,0.83]; % \u6a21\u62df\u66f2\u7ebf\u989c\u8272
\ncolor_3sigma = [1,0,0]; % 3\u03c3\u8fb9\u754c\u989c\u8272
\ncolor_mean = [0,0,0.8]; % \u5747\u503c\u7ebf\u989c\u8272
\n% \u7ebf\u5bbd\u53c2\u6570
\nlw_monte = 0.6; % \u6a21\u62df\u66f2\u7ebf\u7ebf\u5bbd
\nlw_3sigma = 1.8; % 3\u03c3\u8fb9\u754c\u7ebf\u5bbd
\nlw_mean = 2.5; % \u5747\u503c\u7ebf\u7ebf\u5bbd
\n% \u6587\u5b57\u53c2\u6570
\nxlabel_txt = 'Year'; % x\u8f74\u6807\u7b7e
\nylabel_txt = 'Profit (\u00d710^5 USD)'; % y\u8f74\u6807\u7b7e
\nlegend_monte = 'Monte Carlo Scenarios';
\nlegend_3sigma = '3\u03c3 Boundary';
\nlegend_mean = 'Profit Mean (\u03bc)';<\/p>\n
%% ===================== \u7edf\u4e00\u6a21\u677f\u6267\u884c\u533a =====================
\n% \u6570\u636e\u51c6\u5907
\nrng(42);
\nyear = year_start:year_end;
\nn_year = length(year);
\nprofit_sim = zeros(sim_num, n_year);<\/p>\n
% \u751f\u6210\u5229\u6da6\u6570\u636e
\nfor i = 1:sim_num
\n growth_idx = year <= 2035;
\n profit_growth = (1e5 + 2e4*(year(growth_idx)-2020) + 5e4*randn(1,sum(growth_idx)))\/profit_unit;
\n decline_idx = (year>2035)&(year<=2050);
\n profit_decline = (6e5 – 8e3*(year(decline_idx)-2035) + 6e4*randn(1,sum(decline_idx)))\/profit_unit;
\n stable_idx = year > 2050;
\n profit_stable = (6e5 + 2e4*randn(1,sum(stable_idx)))\/profit_unit;
\n profit_sim(i,:) = [profit_growth, profit_decline, profit_stable];
\n profit_sim(i,profit_sim(i,:)<0.1) = 0.1;
\nend<\/p>\n
% \u8ba1\u7b973\u03c3\u8fb9\u754c
\nprofit_mean = mean(profit_sim,1);
\nprofit_std = std(profit_sim,0,1);
\nupper_3sigma = profit_mean + 3*profit_std;
\nlower_3sigma = profit_mean – 3*profit_std;
\nlower_3sigma(lower_3sigma<0.1) = 0.1;<\/p>\n
% \u7ed8\u56fe
\nfigure('Position',[150,150,1000,600],'Color','white'); % \u7a97\u53e3\u80cc\u666f\u8bbe\u4e3a\u767d\u8272
\nhold on; grid on; grid minor;
\nax = gca;
\nax.FontName = 'Times New Roman';
\nax.FontSize = 10;
\nax.XTick = year_start:5:year_end;
\nax.Color = 'white'; % \u5750\u6807\u8f74\u80cc\u666f\u8bbe\u4e3a\u767d\u8272<\/p>\n
% \u8499\u7279\u5361\u6d1b\u66f2\u7ebf
\nplot(year, profit_sim(1,:), 'LineWidth',lw_monte, 'Color',color_monte, 'DisplayName',legend_monte);
\nfor i = 2:sim_num
\n plot(year, profit_sim(i,:), 'LineWidth',lw_monte, 'Color',color_monte, 'HandleVisibility','off');
\nend<\/p>\n
% 3\u03c3\u8fb9\u754c
\nplot(year, upper_3sigma, '–', 'LineWidth',lw_3sigma, 'Color',color_3sigma, 'DisplayName',legend_3sigma);
\nplot(year, lower_3sigma, '–', 'LineWidth',lw_3sigma, 'Color',color_3sigma, 'HandleVisibility','off');<\/p>\n
% \u5747\u503c\u7ebf
\nplot(year, profit_mean, '-', 'LineWidth',lw_mean, 'Color',color_mean, 'DisplayName',legend_mean);<\/p>\n
% \u6807\u6ce8\u7f8e\u5316
\nxlabel(xlabel_txt, 'FontSize',13, 'FontWeight','bold', 'FontName','Times New Roman');
\nylabel(ylabel_txt, 'FontSize',13, 'FontWeight','bold', 'FontName','Times New Roman');
\nlegend('Location','northeast', 'FontSize',11);
\nhold off;<\/p>\n
% \u7ed3\u679c\u8f93\u51fa
\nfprintf('\u6a21\u62df\u8303\u56f4：%d-%d\u5e74，\u66f2\u7ebf\u6570\u91cf：%d\u6761\\\\n', year_start, year_end, sim_num);
\nfprintf('\u5229\u6da6\u4e0b\u9650\u7ea6\u675f：0.1\u00d710^5 USD\\\\n'); <\/p>\n
\u53c2\u8003\u56fe：\u7f8e\u8d5b-2020-D\u9898-\u54c8\u2f20\u5947\u2f9c\u7403\u961f #2006782![\"\"](\"https:\/\/www.wsisp.com\/helps\/wp-content\/uploads\/2026\/01\/20260127151830-6978d74621cf9.png\")
<\/p>\n
\u81ea\u5236\u56fe：<\/p>\n
![\"\"](\"https:\/\/www.wsisp.com\/helps\/wp-content\/uploads\/2026\/01\/20260127151830-6978d74686030.png\")
<\/p>\n
%% \u7ed8\u5236\u53cc\u5750\u6807\u6298\u7ebf\u56fe，\u6709\u586b\u5145\u533a\u57df，\u5e26\u5149\u6655，\u63d2\u503c\u5e73\u6ed1
\nclear; clc; close all;<\/p>\n
% \u81ea\u5b9a\u4e49\u53c2\u6570\u533a<\/p>\n
%AI\u63d0\u793a\u8bcd：
\n%\u6211\u5e0c\u671b\u63d2\u503c\u7ed8\u5236\u5149\u6ed1\u7684\u6709\u5149\u6655\u7684\u66f2\u7ebf，\u4f46\u662f\u53ea\u5728\u771f\u5b9e\u7684\u91c7\u6837\u70b9\u5904\u6807\u8bb0\u3002
\n%\u5982\u679c\u4f60\u8fd9\u6837\u4fee\u6539，\u76f8\u5f53\u4e8e\u662f\u6240\u6709\u7684\u5dee\u503c\u90fd\u6807\u8bb0，\u6574\u4e2a\u56fe\u90fd\u6df7\u4e71\u4e0d\u582a，
\n%\u5982\u679c\u4e0d\u8fd9\u6837\u7684\u8bdd，\u7ed8\u5236\u51fa\u6765\u7684\u56fe\u4f8b\u5c31\u6ca1\u6709\u7ebf\u578b\u3002
\n%\u800c\u5982\u679c\u8ba9\u56fe\u4f8b\u6709\u7ebf\u578b，\u90a3\u53c8\u5fc5\u987b\u5728\u6700\u540e\u753b\u91c7\u6837\u6807\u8bb0\u7684\u65f6\u5019\u7684\u5e26\u4e0a\u7ebf\u578b
\n%\u8fd9\u6837\u5c31\u4e0d\u662f\u5149\u6ed1\u7684\u66f2\u7ebf\u4e86，\u4f60\u80fd\u660e\u767d\u5e76\u89e3\u51b3\u6211\u7684\u95ee\u9898\u5417
\n%% \u7ed8\u5236\u53cc\u5750\u6807\u6298\u7ebf\u56fe，\u6709\u586b\u5145\u533a\u57df，\u5e26\u5149\u6655，\u63d2\u503c\u5e73\u6ed1
\nclear; clc; close all;<\/p>\n
% \u81ea\u5b9a\u4e49\u53c2\u6570\u533a
\n% \u6570\u636e\u70b9，epsilon\u4e3a\u6a2a\u5750\u6807，\u5176\u4f59\u4e3a\u5bf9\u5e94\u6a2a\u5750\u6807\u7684\u7eb5\u5750\u6807\u503c
\nepsilon_raw = [0, 0.2, 0.4, 0.6, 0.8, 1.0, 1.2, 1.4, 1.6, 1.8, 2.0];
\nbeta1_raw = [0.32, 0.25, 0.38, 0.42, 0.34, 0.43, 0.41, 0.39, 0.40, 0.31, 0.40];
\nbeta2_raw = [0.05, 0.20, 0.35, 0.40, 0.62, 0.50, 0.75, 0.78, 0.79, 0.60, 0.80];
\nbeta3_raw = [0.15, 0.32, 0.42, 0.38, 0.35, 0.28, 0.22, 0.15, 0.32, 0.45, 0.33];
\nbeta4_raw = [0.10, 0.12, 0.15, 0.18, 0.22, 0.25, 0.27, 0.28, 0.29, 0.20, 0.31];
\nbeta5_raw = [0.05, 0.06, 0.07, 0.08, 0.09, 0.10, 0.18, 0.19, 0.20, 0.19, 0.20];
\nerror_raw = [0.05, 0.07, 0.10, 0.13, 0.17, 0.19, 0.20, 0.21, 0.18, 0.15, 0.13];<\/p>\n
% \u63d2\u503c\u53c2\u6570
\ninterp_points = 100; % \u63d2\u503c\u70b9\u6570
\nepsilon = linspace(epsilon_raw(1), epsilon_raw(end), interp_points);<\/p>\n
% \u5bf9\u6570\u636e\u8fdb\u884c\u63d2\u503c（\u4f7f\u7528pchip\u4fdd\u6301\u5f62\u72b6）
\nbeta1 = interp1(epsilon_raw, beta1_raw, epsilon, 'pchip');
\nbeta2 = interp1(epsilon_raw, beta2_raw, epsilon, 'pchip');
\nbeta3 = interp1(epsilon_raw, beta3_raw, epsilon, 'pchip');
\nbeta4 = interp1(epsilon_raw, beta4_raw, epsilon, 'pchip');
\nbeta5 = interp1(epsilon_raw, beta5_raw, epsilon, 'pchip');
\nerror_data = interp1(epsilon_raw, error_raw, epsilon, 'pchip');<\/p>\n
% \u989c\u8272\u8bbe\u7f6e
\ncolor_beta1 = [0, 0.447, 0.741];
\ncolor_beta2 = [0.466, 0.674, 0.188];
\ncolor_beta3 = [0.85, 0.325, 0.098];
\ncolor_beta4 = [0.5, 0.5, 0.5];
\ncolor_beta5 = [0.929, 0.694, 0.125];
\ncolor_error = [1, 0, 0];<\/p>\n
% \u7ebf\u578b\u8bbe\u7f6e
\nline_style_beta = {'-', '–', ':', '-.', '-'};% \u7ebf\u578b\u53ea\u6709 '-' | '–' | ':' | '-.' | 'none'
\nmarker_style = {'o', 's', '^', 'v', 'd'};<\/p>\n
% \u5b57\u4f53\u8bbe\u7f6e
\nfont_name = 'Times New Roman';
\nfont_size_label = 13;
\nfont_weight_label = 'bold';
\nfont_size_legend = 11;<\/p>\n
% \u7edf\u4e00\u6a21\u677f\u6267\u884c\u533a
\nfigure('Color', 'white');
\nset(gcf, 'Position', [100, 100, 800, 500]);<\/p>\n
hold on;<\/p>\n
% 1. \u5148\u7ed8\u5236\u586b\u5145\u533a\u57df
\nyyaxis right;<\/p>\n
% \u5bf9\u586b\u5145\u533a\u57df\u4e5f\u9700\u8981\u63d2\u503c\u6570\u636e
\nerror_interp_fill = interp1(epsilon_raw, error_raw, epsilon, 'pchip');
\nnum_points = length(epsilon);<\/p>\n
% \u521b\u5efa\u6e10\u53d8\u586b\u5145
\nfor k = 1:num_points-1
\n % \u6bcf\u4e2a\u56db\u8fb9\u5f62\u9876\u70b9\u5750\u6807
\n x_quad = [epsilon(k), epsilon(k+1), epsilon(k+1), epsilon(k)];
\n y_quad = [error_interp_fill(k), error_interp_fill(k+1), 0, 0];<\/p>\n
% \u8ba1\u7b97\u900f\u660e\u5ea6（\u4ece\u5de6\u5230\u53f3\u53d8\u6d45）
\n alpha_left = 2 – (k-1)\/(num_points-1);
\n alpha_right = 2 – k\/(num_points-1);<\/p>\n
% \u989c\u8272（\u7ea2\u8272\u9010\u6e10\u53d8\u6d45）
\n color_left = [1, 0.5*alpha_left, 0.5*alpha_left];
\n color_right = [1, 0.5*alpha_right, 0.5*alpha_right];<\/p>\n
% \u7ed8\u5236\u6bcf\u4e2a\u56db\u8fb9\u5f62
\n patch(x_quad, y_quad, 'k', …
\n 'FaceColor', 'interp', …
\n 'FaceVertexCData', [color_left; color_right; color_right; color_left], …
\n 'EdgeColor', 'none', …
\n 'FaceAlpha', 0.15, …
\n 'HandleVisibility', 'off');
\nend<\/p>\n
% \u53f3\u5750\u6807\u8f74\u8bbe\u7f6e
\nylim([0, 0.22]);
\nylabel('Lowest Error', …
\n 'FontName', font_name, …
\n 'FontSize', font_size_label, …
\n 'FontWeight', font_weight_label);
\nset(gca, 'YColor', [0, 0, 0]);<\/p>\n
% 2. \u7ed8\u5236\u6027\u80fd\u6307\u6807\u66f2\u7ebf
\nyyaxis left;<\/p>\n
% \u5b58\u50a8\u56fe\u4f8b\u53e5\u67c4
\nh_line = []; % \u7ebf\u6761\u7684\u53e5\u67c4（\u7528\u4e8e\u7ebf\u578b\u56fe\u4f8b）
\nh_marker = []; % \u6807\u8bb0\u7684\u53e5\u67c4（\u7528\u4e8e\u6807\u8bb0\u56fe\u4f8b）<\/p>\n
% \u7ed8\u5236\u6027\u80fd\u6307\u6807\u66f2\u7ebf（\u5e26\u900f\u660e\u5149\u6655）
\nbeta_data = {beta1, beta2, beta3, beta4, beta5};
\nbeta_raw_data = {beta1_raw, beta2_raw, beta3_raw, beta4_raw, beta5_raw};
\nbeta_colors = {color_beta1, color_beta2, color_beta3, color_beta4, color_beta5};
\nbeta_names = {'\\\\beta_1', '\\\\beta_2', '\\\\beta_3', '\\\\beta_4', '\\\\beta_5'};<\/p>\n
for i = 1:5
\n % \u591a\u5c42\u5149\u6655（\u4f7f\u7528\u63d2\u503c\u6570\u636e\u7ed8\u5236\u5e73\u6ed1\u5149\u6655）
\n for j = 6:-1:1
\n line_width = 2 + (j-1)*1.2;
\n alpha_value = 0.10 – (j-1)*0.01;<\/p>\n
rgba_color = [beta_colors{i}, alpha_value];<\/p>\n
plot(epsilon, beta_data{i}, …
\n 'LineStyle', '-', …
\n 'Marker', 'none', …
\n 'LineWidth', line_width, …
\n 'Color', rgba_color, …
\n 'HandleVisibility', 'off');
\n end<\/p>\n
% \u4e3b\u66f2\u7ebf（\u63d2\u503c\u6570\u636e\u7ed8\u5236\u5e73\u6ed1\u7ebf）
\n h_line(i) = plot(epsilon, beta_data{i}, …
\n 'LineStyle', line_style_beta{i}, …
\n 'LineWidth', 2, …
\n 'Color', beta_colors{i}, …
\n 'Marker', 'none', …
\n 'HandleVisibility', 'on', …
\n 'DisplayName', beta_names{i});<\/p>\n
% \u5728\u539f\u59cb\u6570\u636e\u70b9\u4f4d\u7f6e\u7ed8\u5236\u6807\u8bb0（\u4f7f\u7528\u65e0\u7ebf\u7684scatter）
\n h_marker(i) = plot(epsilon_raw, beta_raw_data{i}, …
\n 'LineStyle', 'none', … % \u4e0d\u8fde\u63a5\u6807\u8bb0
\n 'Marker', marker_style{i}, …
\n 'MarkerSize', 6, …
\n 'MarkerFaceColor', beta_colors{i}, …
\n 'MarkerEdgeColor', beta_colors{i}, …
\n 'HandleVisibility', 'off'); % \u4e0d\u76f4\u63a5\u663e\u793a\u5728\u81ea\u52a8\u56fe\u4f8b\u4e2d
\nend<\/p>\n
% \u5de6\u5750\u6807\u8f74\u8bbe\u7f6e
\nylim([-0.1, 0.9]);
\nylabel('Performance Indicators', …
\n 'FontName', font_name, …
\n 'FontSize', font_size_label, …
\n 'FontWeight', font_weight_label);
\nset(gca, 'YColor', [0, 0, 0]);<\/p>\n
% 3. \u6700\u540e\u7ed8\u5236\u8bef\u5dee\u66f2\u7ebf
\nyyaxis right;<\/p>\n
% \u7ed8\u5236\u8bef\u5dee\u66f2\u7ebf（\u5e26\u900f\u660e\u5149\u6655）
\n% \u591a\u5c42\u5149\u6655（\u4f7f\u7528\u63d2\u503c\u6570\u636e）
\nfor j = 5:-1:1
\n line_width = 2 + (j-1)*1.2;
\n alpha_value = 0.25 – (j-1)*0.05;<\/p>\n
plot(epsilon, error_data, …
\n 'LineStyle', '-', …
\n 'LineWidth', line_width, …
\n 'Color', [color_error, alpha_value], …
\n 'Marker', 'none', …
\n 'HandleVisibility', 'off');
\nend<\/p>\n
% \u4e3b\u8bef\u5dee\u66f2\u7ebf – \u5148\u7ed8\u5236\u5e73\u6ed1\u7ebf
\nh_line_error = plot(epsilon, error_data, '–', …
\n 'LineWidth', 2, …
\n 'Color', color_error, …
\n 'Marker', 'none', …
\n 'HandleVisibility', 'on', …
\n 'DisplayName', 'Lowest Error');<\/p>\n
% \u5728\u539f\u59cb\u6570\u636e\u70b9\u4f4d\u7f6e\u7ed8\u5236\u6807\u8bb0
\nh_marker_error = plot(epsilon_raw, error_raw, …
\n 'LineStyle', 'none', …
\n 'Marker', 'o', …
\n 'MarkerSize', 4, …
\n 'MarkerFaceColor', color_error, …
\n 'MarkerEdgeColor', color_error, …
\n 'HandleVisibility', 'off');<\/p>\n
% \u5171\u540c\u8bbe\u7f6e
\nxlabel('\\\\epsilon', …
\n 'FontName', font_name, …
\n 'FontSize', font_size_label, …
\n 'FontWeight', font_weight_label);
\nxlim([0, 2]);
\ngrid on;<\/p>\n
% \u8bbe\u7f6e\u7ebf\u6761\u56fe\u4f8b（\u53ea\u6709\u7ebf\u578b，\u6ca1\u6709\u6807\u8bb0）
\nlegend(h_line, 'Location', 'best', 'FontSize', font_size_legend, …
\n 'FontName', font_name);<\/p>\n
% \u7edf\u4e00\u7f8e\u5316
\nset(gca, 'FontName', font_name, 'FontSize', 11);
\nbox on;
\nhold off;<\/p>\n
% \u5728\u7ed8\u5236\u5b8c\u6240\u6709\u56fe\u5f62\u540e，\u6dfb\u52a0\u4ee5\u4e0b\u4ee3\u7801\u6765\u521b\u5efa\u81ea\u5b9a\u4e49\u56fe\u4f8b
\n% \u521b\u5efa\u81ea\u5b9a\u4e49\u56fe\u4f8b\u6761\u76ee
\nlegend_labels = {'\\\\beta_1', '\\\\beta_2', '\\\\beta_3', '\\\\beta_4', '\\\\beta_5', 'Lowest Error'};
\nlegend_colors = [beta_colors, {color_error}];
\nlegend_line_styles = [line_style_beta, {'–'}];
\nlegend_markers = [marker_style, {'o'}];<\/p>\n
% \u521b\u5efa\u65b0\u56fe\u5f62\u6765\u83b7\u53d6\u56fe\u4f8b
\nfigure('Position', [1000, 100, 400, 300], 'Color', 'white', 'Visible', 'off');
\nhold on;
\nlegend_handles = [];
\nfor i = 1:6
\n % \u521b\u5efa\u540c\u65f6\u5305\u542b\u7ebf\u578b\u548c\u6807\u8bb0\u7684\u56fe\u5f62\u5bf9\u8c61
\n h = plot([0, 1], [i, i], …
\n 'LineStyle', legend_line_styles{i}, …
\n 'LineWidth', 2, …
\n 'Color', legend_colors{i}, …
\n 'Marker', legend_markers{i}, …
\n 'MarkerSize', 6, …
\n 'MarkerFaceColor', legend_colors{i}, …
\n 'DisplayName', legend_labels{i});
\n legend_handles = [legend_handles, h];
\nend
\nhold off;
\n%% \u81ea\u5b9a\u4e49\u56fe\u4f8b
\n%\u5206\u522b\u7ed8\u5236\u5149\u6ed1\u66f2\u7ebf（\u5e26\u7ebf\u578b）\u548c\u6807\u8bb0\u70b9（\u53ea\u5728\u91c7\u6837\u70b9\u5904）
\n%\u521b\u5efa\u4ee3\u7406\u5bf9\u8c61（proxy objects）\u6765\u6784\u5efa\u81ea\u5b9a\u4e49\u56fe\u4f8b，\u8fd9\u4e9b\u4ee3\u7406\u5bf9\u8c61\u540c\u65f6\u5305\u542b\u7ebf\u578b\u548c\u6807\u8bb0
\n%\u4f7f\u7528\u4ee3\u7406\u5bf9\u8c61\u521b\u5efa\u6700\u7ec8\u56fe\u4f8b<\/p>\n
% \u83b7\u53d6\u8fd9\u4e2a\u56fe\u4f8b
\ncustom_legend = legend(legend_handles, 'Location', 'best', …
\n 'FontSize', font_size_legend, …
\n 'FontName', font_name);
\ncustom_legend_items = custom_legend.EntryContainer.Children;<\/p>\n
% \u5173\u95ed\u4e34\u65f6\u56fe\u5f62
\nclose(gcf);<\/p>\n
% \u56de\u5230\u4e3b\u56fe\u5f62，\u521b\u5efa\u7a7a\u767d\u56fe\u4f8b\u7136\u540e\u6dfb\u52a0\u81ea\u5b9a\u4e49\u6761\u76ee
\nfigure(1); % \u56de\u5230\u4e3b\u56fe
\n% \u9996\u5148\u9690\u85cf\u4e4b\u524d\u7684\u56fe\u4f8b
\nlegend('off');<\/p>\n
% \u521b\u5efa\u81ea\u5b9a\u4e49\u56fe\u4f8b
\n% \u65b9\u6cd51：\u4f7f\u7528annotation\u521b\u5efa（\u590d\u6742\u4f46\u7cbe\u786e）
\n% \u65b9\u6cd52：\u7b80\u5316\u65b9\u6848 – \u521b\u5efa\u4ee3\u7406\u5bf9\u8c61
\n% \u521b\u5efa\u4ee3\u7406\u56fe\u5f62\u5bf9\u8c61（\u4e0d\u5728\u56fe\u4e0a\u663e\u793a）
\nhold on;
\nproxy_handles = [];
\nfor i = 1:6
\n % \u521b\u5efa\u4e0d\u5728\u56fe\u4e0a\u663e\u793a\u7684\u4ee3\u7406\u5bf9\u8c61
\n proxy_handles(i) = plot(NaN, NaN, …
\n 'LineStyle', legend_line_styles{i}, …
\n 'LineWidth', 2, …
\n 'Color', legend_colors{i}, …
\n 'Marker', legend_markers{i}, …
\n 'MarkerSize', 6, …
\n 'MarkerFaceColor', legend_colors{i}, …
\n 'DisplayName', legend_labels{i});
\nend<\/p>\n
% \u521b\u5efa\u56fe\u4f8b
\nlegend(proxy_handles, legend_labels, 'Location', 'best', …
\n 'FontSize', font_size_legend, 'FontName', font_name); <\/p>\n
\u53c2\u8003：<\/p>\n
![\"\"](\"https:\/\/www.wsisp.com\/helps\/wp-content\/uploads\/2026\/01\/20260127151830-6978d746bbb81.png\")
<\/p>\n
\u81ea\u5236：<\/p>\n
![\"\"](\"https:\/\/www.wsisp.com\/helps\/wp-content\/uploads\/2026\/01\/20260127151831-6978d7478bc4a.png\")
<\/p>\n
\u4ee3\u7801：<\/p>\n
clear; clc; close all;<\/p>\n
%% ===================== \u81ea\u5b9a\u4e49\u533a（\u53ef\u4fee\u6539\u53c2\u6570）=====================
\n% 1. \u573a\u666f：\u65e0\u4eba\u673a\u98ce\u901f\u4e0e\u98ce\u5411\u5bf9\u4fe1\u53f7\u5e72\u6270\u7684\u5f71\u54cd（\u5355\u4f4d：%）
\nvelocity = [1,1.5,2,2.5,3,3.5]; % \u98ce\u901f（m\/s）
\nangle = [0,45,90,135,180,225,270]; % \u98ce\u5411\u89d2（\u00b0）
\n% 3\u7ec4\u573a\u666f\u6570\u636e（\u884c\u6570=\u98ce\u901f\u6570，\u5217\u6570=\u98ce\u5411\u89d2\u6570）
\ndata1 = [
\n 0.2 0.5 1.2 0.8 0.3 0.1 0.0;
\n 0.8 1.5 3.2 2.5 1.0 0.4 0.2;
\n 1.5 3.8 6.5 5.2 2.2 0.9 0.5;
\n 1.2 4.2 7.8 8.5 4.0 1.8 0.8;
\n 0.6 2.5 4.8 5.2 3.5 2.0 1.0;
\n 0.3 1.2 2.5 2.8 2.0 1.5 0.8;
\n];
\ndata2 = [
\n 0.5 1.2 2.5 1.8 0.8 0.4 0.2;
\n 1.8 3.5 7.2 6.5 3.0 1.4 0.8;
\n 3.5 8.8 15.5 14.2 7.2 3.9 1.5;
\n 4.2 10.2 18.8 19.5 10.0 5.8 2.8;
\n 2.6 6.5 11.8 12.2 8.5 5.0 3.0;
\n 1.3 3.2 6.5 7.8 5.0 3.5 2.8;
\n];
\ndata3 = [
\n 0.3 0.8 1.8 1.2 0.5 0.2 0.1;
\n 1.2 2.5 5.2 4.5 2.0 0.9 0.4;
\n 2.5 6.8 10.5 9.2 4.2 1.9 0.8;
\n 3.2 8.2 14.8 15.5 7.0 3.8 1.8;
\n 1.6 4.5 8.8 9.2 6.5 3.0 2.0;
\n 0.8 2.2 4.5 5.8 4.0 2.5 1.8;
\n];<\/p>\n
% 2. \u989c\u8272\u53c2\u6570
\ncolor1 = [0.2 0.5 0.8]; % \u573a\u666f1：\u6d45\u84dd\u8272
\ncolor2 = [0.1 0.7 0.6]; % \u573a\u666f2：\u9752\u7eff\u8272
\ncolor3 = [0.8 0.4 0.2]; % \u573a\u666f3：\u6a59\u7ea2\u8272<\/p>\n
% 3. \u9ad8\u7ea7\u7f51\u683c\u548c\u80cc\u666f\u989c\u8272\u8bbe\u7f6e
\ngrid_color = [0.4, 0.4, 0.4]; % \u7f51\u683c\u7ebf\u989c\u8272：\u6df1\u7070\u8272，\u66f4\u4e13\u4e1a
\nbg_color = [0.98, 0.98, 0.98]; % \u80cc\u666f\u989c\u8272：\u6de1\u7070\u8272，\u66f4\u9ad8\u7ea7
\naxis_color = [0.3, 0.3, 0.3]; % \u5750\u6807\u8f74\u989c\u8272：\u6df1\u7070\u8272
\ngrid_alpha = 0.5; % \u7f51\u683c\u7ebf\u900f\u660e\u5ea6<\/p>\n
% 4. \u5750\u6807\u8f74\u6587\u5b57
\nxlabel_txt = 'velocity(m\/s)';
\nylabel_txt = 'angle(\\\\circ)';
\nzlabel_txt = 'influence(%)';
\nfontSize = 11;
\nlabelFontSize = 13;<\/p>\n
%% ===================== \u7edf\u4e00\u6a21\u677f\u6267\u884c\u533a =====================
\n% \u521b\u5efa3\u5b50\u56fe\u5e03\u5c40
\nfigure('Position',[150,150,1500,500],'Color','white');
\ntiledlayout(1,3,'TileSpacing','compact','Padding','compact');<\/p>\n
% \u5b50\u56fe1：\u4f4e\u5e72\u6270\u573a\u666f
\nnexttile;
\nh1 = bar3(data1'); % \u4ec5\u4f20\u5165\u6570\u636e\u77e9\u9635，\u4fee\u6b63\u53c2\u6570\u8fc7\u591a\u95ee\u9898
\n% \u8bbe\u7f6e\u67f1\u5b50\u6837\u5f0f
\nset(h1,'EdgeColor','cyan','FaceColor',color1);
\nax = gca;<\/p>\n
% ========== \u8bbe\u7f6e\u9ad8\u7ea7\u7f51\u683c\u548c\u80cc\u666f ==========
\nax.Color = bg_color; % \u8bbe\u7f6e\u80cc\u666f\u989c\u8272
\nax.GridColor = grid_color; % \u8bbe\u7f6e\u7f51\u683c\u7ebf\u989c\u8272
\nax.GridAlpha = grid_alpha; % \u8bbe\u7f6e\u7f51\u683c\u7ebf\u900f\u660e\u5ea6
\nax.XColor = axis_color; % \u8bbe\u7f6eX\u8f74\u989c\u8272
\nax.YColor = axis_color; % \u8bbe\u7f6eY\u8f74\u989c\u8272
\nax.ZColor = axis_color; % \u8bbe\u7f6eZ\u8f74\u989c\u8272
\nax.GridLineStyle = '-'; % \u7f51\u683c\u7ebf\u6837\u5f0f\u4e3a\u5b9e\u7ebf
\nax.MinorGridLineStyle = 'none'; % \u65e0\u6b21\u8981\u7f51\u683c\u7ebf
\n% ========================================<\/p>\n
ax.FontName = 'Times New Roman';
\nax.FontSize = fontSize;
\nax.XLabel.String = xlabel_txt;
\nax.XLabel.FontWeight = 'bold';
\nax.XLabel.FontSize = labelFontSize;
\nax.YLabel.String = ylabel_txt;
\nax.YLabel.FontWeight = 'bold';
\nax.YLabel.FontSize = labelFontSize;
\nax.ZLabel.String = zlabel_txt;
\nax.ZLabel.FontWeight = 'bold';
\nax.ZLabel.FontSize = labelFontSize;
\nax.TickDir = 'out';
\n% \u81ea\u5b9a\u4e49X\/Y\u8f74\u523b\u5ea6\u6807\u7b7e
\nax.XTickLabel = velocity;
\nax.YTickLabel = angle;
\nview(130,30);
\nzlim([0,10]);<\/p>\n
% \u5b50\u56fe2：\u4e2d\u5e72\u6270\u573a\u666f
\nnexttile;
\nh2 = bar3(data2');
\nset(h2,'EdgeColor','g','FaceColor',color2);
\nax = gca;<\/p>\n
% ========== \u8bbe\u7f6e\u9ad8\u7ea7\u7f51\u683c\u548c\u80cc\u666f ==========
\nax.Color = bg_color; % \u8bbe\u7f6e\u80cc\u666f\u989c\u8272
\nax.GridColor = grid_color; % \u8bbe\u7f6e\u7f51\u683c\u7ebf\u989c\u8272
\nax.GridAlpha = grid_alpha; % \u8bbe\u7f6e\u7f51\u683c\u7ebf\u900f\u660e\u5ea6
\nax.XColor = axis_color; % \u8bbe\u7f6eX\u8f74\u989c\u8272
\nax.YColor = axis_color; % \u8bbe\u7f6eY\u8f74\u989c\u8272
\nax.ZColor = axis_color; % \u8bbe\u7f6eZ\u8f74\u989c\u8272
\nax.GridLineStyle = '-'; % \u7f51\u683c\u7ebf\u6837\u5f0f\u4e3a\u5b9e\u7ebf
\nax.MinorGridLineStyle = 'none'; % \u65e0\u6b21\u8981\u7f51\u683c\u7ebf
\n% ========================================<\/p>\n
ax.FontName = 'Times New Roman';
\nax.FontSize = fontSize;
\nax.XLabel.String = xlabel_txt;
\nax.XLabel.FontWeight = 'bold';
\nax.XLabel.FontSize = labelFontSize;
\nax.YLabel.String = ylabel_txt;
\nax.YLabel.FontWeight = 'bold';
\nax.YLabel.FontSize = labelFontSize;
\nax.ZLabel.String = zlabel_txt;
\nax.ZLabel.FontWeight = 'bold';
\nax.ZLabel.FontSize = labelFontSize;
\nax.TickDir = 'out';
\nax.XTickLabel = velocity;
\nax.YTickLabel = angle;
\nview(130,30);
\nzlim([0,25]);<\/p>\n
% \u5b50\u56fe3：\u9ad8\u5e72\u6270\u573a\u666f
\nnexttile;
\nh3 = bar3(data3');
\nset(h3,'EdgeColor','yellow','FaceColor',color3);
\nax = gca;<\/p>\n
% ========== \u8bbe\u7f6e\u9ad8\u7ea7\u7f51\u683c\u548c\u80cc\u666f ==========
\nax.Color = bg_color; % \u8bbe\u7f6e\u80cc\u666f\u989c\u8272
\nax.GridColor = grid_color; % \u8bbe\u7f6e\u7f51\u683c\u7ebf\u989c\u8272
\nax.GridAlpha = grid_alpha; % \u8bbe\u7f6e\u7f51\u683c\u7ebf\u900f\u660e\u5ea6
\nax.XColor = axis_color; % \u8bbe\u7f6eX\u8f74\u989c\u8272
\nax.YColor = axis_color; % \u8bbe\u7f6eY\u8f74\u989c\u8272
\nax.ZColor = axis_color; % \u8bbe\u7f6eZ\u8f74\u989c\u8272
\nax.GridLineStyle = '-'; % \u7f51\u683c\u7ebf\u6837\u5f0f\u4e3a\u5b9e\u7ebf
\nax.MinorGridLineStyle = 'none'; % \u65e0\u6b21\u8981\u7f51\u683c\u7ebf
\n% ========================================<\/p>\n
ax.FontName = 'Times New Roman';
\nax.FontSize = fontSize;
\nax.XLabel.String = xlabel_txt;
\nax.XLabel.FontWeight = 'bold';
\nax.XLabel.FontSize = labelFontSize;
\nax.YLabel.String = ylabel_txt;
\nax.YLabel.FontWeight = 'bold';
\nax.YLabel.FontSize = labelFontSize;
\nax.ZLabel.String = zlabel_txt;
\nax.ZLabel.FontWeight = 'bold';
\nax.ZLabel.FontSize = labelFontSize;
\nax.TickDir = 'out';
\nax.XTickLabel = velocity;
\nax.YTickLabel = angle;
\nview(130,30);
\nzlim([0,17.5]);<\/p>\n
% \u7ed3\u679c\u8f93\u51fa
\nfprintf('\u7ed8\u56fe\u5b8c\u6210：3\u573a\u666f\u65e0\u4eba\u673a\u98ce\u901f-\u98ce\u5411\u5e72\u62703D\u67f1\u72b6\u56fe\\\\n'); <\/p>\n
\u53c2\u8003：<\/p>\n
![\"\"](\"https:\/\/www.wsisp.com\/helps\/wp-content\/uploads\/2026\/01\/20260127151832-6978d748300d8.png\")
<\/p>\n
\u81ea\u5236：<\/p>\n
![\"\"](\"https:\/\/www.wsisp.com\/helps\/wp-content\/uploads\/2026\/01\/20260127151832-6978d74872937.png\")
![\"\"](\"https:\/\/www.wsisp.com\/helps\/wp-content\/uploads\/2026\/01\/20260127151832-6978d748d3c02.png\")
<\/p>\n
\u4ee3\u7801<\/p>\n
clear; clc; close all;<\/p>\n
%% \u81ea\u5b9a\u4e49\u533a
\ncountries = {'FRA','KOR','CHN','SUN','AUS','JPN','GRC','BRA','USA'};
\nn = length(countries);
\nscore_ratio = [1.93, 1.73, 1.58, 1.56, 1.41, 1.41, 1.23, 1.11, 0.93];
\nscores = [122, 67, 218, 424, 115, 126, 34, 39, 401];<\/p>\n
bar_colors = [
\n 0.2 0.4 0.9;
\n 0.3 0.5 0.92;
\n 0.4 0.6 0.94;
\n 0.5 0.7 0.96;
\n 0.6 0.8 0.98;
\n 0.7 0.85 0.99;
\n 0.8 0.9 1.0;
\n 0.85 0.93 1.0;
\n 0.9 0.96 1.0;
\n];<\/p>\n
line_color = [1.0 0.8 0.0];
\nmarker_color = [1.0 0.8 0.0];<\/p>\n
interp_points = 100;
\nx_raw = 1:n;
\nx_interp = linspace(1, n, interp_points);
\nscores_interp = interp1(x_raw, scores, x_interp, 'pchip');<\/p>\n
xlabel_txt = 'Country';
\nylabel_left_txt = 'Olympic score ratio';
\nylabel_right_txt = 'scores';
\nlegend_left_txt = {'Olympic score ratio'};
\nlegend_right_txt = {'scores'};
\nfontSize = 11;
\nlabelFontSize = 13;<\/p>\n
%% \u7edf\u4e00\u6a21\u677f\u6267\u884c\u533a
\nfigure('Position',[150,150,1000,600],'Color','white');<\/p>\n
ax1 = gca;
\nhold(ax1, 'on');
\nh_bar = bar(ax1, score_ratio, 'FaceColor', 'flat', 'EdgeColor', 'none');
\nh_bar.CData = bar_colors;
\nh_bar.DisplayName = legend_left_txt{1};<\/p>\n
for i = 1:n
\n text(ax1, i, score_ratio(i), sprintf('%.2f', score_ratio(i)), …
\n 'FontName','Times New Roman', 'FontSize',10, …
\n 'VerticalAlignment','bottom', 'HorizontalAlignment','center', …
\n 'FontWeight','bold', 'Color','k', …
\n 'BackgroundColor', [1 1 1 0.4], 'Margin', 1, …
\n 'EdgeColor', 'k', 'LineWidth', 0.5);
\nend<\/p>\n
ax2 = axes('Position',ax1.Position,'YAxisLocation','right','Color','none','XColor','none','YColor','k');
\nhold(ax2, 'on');<\/p>\n
for j = 8:-1:1
\n line_width = 3 + (j-1)*1.5;
\n alpha_value = 0.12 – (j-1)*0.015;<\/p>\n
plot(ax2, x_interp, scores_interp, 'LineStyle', '-', …
\n 'LineWidth', line_width, 'Color', [line_color, alpha_value], …
\n 'Marker', 'none', 'HandleVisibility', 'off');
\nend<\/p>\n
h_line = plot(ax2, x_interp, scores_interp, 'LineStyle', '-', …
\n 'LineWidth', 3, 'Color', line_color, 'Marker', 'none', …
\n 'DisplayName', legend_right_txt{1});<\/p>\n
plot(ax2, x_raw, scores, 'LineStyle', 'none', 'Marker', 'o', …
\n 'MarkerSize', 7, 'MarkerFaceColor', marker_color, …
\n 'MarkerEdgeColor', 'k', 'LineWidth', 1, 'HandleVisibility', 'off');<\/p>\n
for i = 1:n
\n text(ax2, x_raw(i), scores(i), num2str(scores(i)), …
\n 'FontName','Times New Roman', 'FontSize',10, …
\n 'VerticalAlignment','bottom', 'HorizontalAlignment','center', …
\n 'FontWeight','bold', 'Color','k', …
\n 'BackgroundColor', [1 1 1 0.4], 'Margin', 1, …
\n 'EdgeColor', 'k', 'LineWidth', 0.5);
\nend<\/p>\n
ax1.FontName = 'Times New Roman';
\nax1.FontSize = fontSize;
\nax1.XTick = 1:n;
\nax1.XTickLabel = countries;
\nax1.XLabel.String = xlabel_txt;
\nax1.XLabel.FontWeight = 'bold';
\nax1.XLabel.FontSize = labelFontSize;
\nax1.XLabel.FontName = 'Times New Roman';
\nax1.YLabel.String = ylabel_left_txt;
\nax1.YLabel.FontWeight = 'bold';
\nax1.YLabel.FontSize = labelFontSize;
\nax1.YLabel.FontName = 'Times New Roman';
\nax1.YLim = [0, 2.1];
\nax1.TickDir = 'in';%\u8fd9\u4e2a\u5730\u65b9\u5982\u679c\u662fout\u5c31\u4f1a\u548c\u53f3\u8fb9\u7684\u51b2\u7a81
\nax1.Box = 'on';
\nax1.YAxisLocation = 'left';<\/p>\n
ax2.FontName = 'Times New Roman';
\nax2.FontSize = fontSize;
\nax2.XTick = [];
\nax2.YLabel.String = ylabel_right_txt;
\nax2.YLabel.FontWeight = 'bold';
\nax2.YLabel.FontSize = labelFontSize;
\nax2.YLabel.FontName = 'Times New Roman';
\nax2.YLim = [0, 450];
\nax2.TickDir = 'out';
\nlinkaxes([ax1,ax2],'x');<\/p>\n
proxy_handle = plot(ax2, NaN, NaN, 'LineStyle', '-', …
\n 'LineWidth', 3, 'Color', line_color, …
\n 'Marker', 'o', 'MarkerSize', 7, …
\n 'MarkerFaceColor', marker_color, …
\n 'MarkerEdgeColor', 'k', 'LineWidth', 1, …
\n 'DisplayName', legend_right_txt{1});<\/p>\n
legend(ax1, h_bar, legend_left_txt, 'Location','northwest',…
\n 'FontSize',10,'FontName','Times New Roman');
\nlegend(ax2, proxy_handle, legend_right_txt, 'Location','northeast',…
\n 'FontSize',10,'FontName','Times New Roman');<\/p>\n
ax1.Legend.AutoUpdate = 'off';
\nax2.Legend.AutoUpdate = 'off';<\/p>\n
hold(ax1, 'off');
\nhold(ax2, 'off');<\/p>\n
fprintf('\u7ed8\u56fe\u5b8c\u6210\\\\n'); <\/p>\n
\u53f3\u8fb9\u56fe\u4f8b，matlab\u4e0d\u597d\u753b\u8fd9\u79cd，\u51d1\u5408\u4e00\u4e0b<\/p>\n
clear; clc; close all;
\n% \u521b\u65b0\u56fe\u4f8b
\n%% ===================== \u81ea\u5b9a\u4e49\u533a（\u53ef\u4fee\u6539\u53c2\u6570）=====================
\n% 1. \u73af\u5f62\u53c2\u6570
\nradius = [1.0, 0.8, 0.6]; % \u4e09\u5c42\u5706\u73af\u534a\u5f84
\nangle = [210, 240, 270]; % \u6bcf\u5c42\u5706\u73af\u7684\u5706\u5f27\u89d2\u5ea6（0-360）
\ncolors = [
\n 0.8 0.9 1.0; % \u6700\u5916\u5c42：\u6d45\u84dd
\n 0.4 0.7 0.9; % \u4e2d\u95f4\u5c42：\u4e2d\u84dd
\n 0.1 0.4 0.8; % \u6700\u5185\u5c42：\u6df1\u84dd
\n];
\n% 2. \u56fe\u4f8b\u53c2\u6570
\nlegend_txt = {'growth rate>1','growth rate>1.2','growth rate>1.5'};
\nfontSize = 11;
\nlabelFontSize = 13;<\/p>\n
%% ===================== \u7edf\u4e00\u6a21\u677f\u6267\u884c\u533a =====================
\nfigure('Position',[150,150,600,600],'Color','white');
\nhold on;<\/p>\n
% \u7ed8\u5236\u4e09\u5c42\u73af\u5f62
\nfor i = 1:length(radius)
\n theta = linspace(-90, angle(i)-90, 100) * pi\/180;
\n x = radius(i) * cos(theta);
\n y = radius(i) * sin(theta);
\n % \u7ed8\u5236\u5916\u5c42\u5706\u5f27
\n plot(x, y, 'Color',colors(i,:), 'LineWidth',25);
\n % \u7ed8\u5236\u5185\u5c42\u7a7a\u767d\u5706\u5f27
\n theta_inner = linspace(angle(i)-90, 270, 100) * pi\/180;
\n x_inner = (radius(i)-0.18) * cos(theta_inner);
\n y_inner = (radius(i)-0.18) * sin(theta_inner);
\n plot(x_inner, y_inner, 'Color','white', 'LineWidth',25);
\nend<\/p>\n
% \u5750\u6807\u8f74\u8bbe\u7f6e
\nax = gca;
\nax.XLim = [-1.2, 1.2];
\nax.YLim = [-1.2, 1.2];
\nax.XTick = [];
\nax.YTick = [];
\nax.Box = 'off';
\nax.Color = 'white';
\naxis equal;<\/p>\n
% \u56fe\u4f8b\u8bbe\u7f6e
\nlegend(legend_txt, 'Location','east','FontSize',fontSize,'FontName','Times New Roman');<\/p>\n
hold off;<\/p>\n
% \u7ed3\u679c\u8f93\u51fa
\nfprintf('\u7ed8\u56fe\u5b8c\u6210：\u4e09\u5c42\u73af\u5f62\u589e\u957f\u901f\u7387\u56fe\\\\n'); <\/p>\n
\u53c2\u8003\u56fe 2025C 2500759<\/p>\n
![\"\"](\"https:\/\/www.wsisp.com\/helps\/wp-content\/uploads\/2026\/01\/20260127151832-6978d748f4212.png\")
<\/p>\n
AI\u4e86\u4e00\u4e0b，\u4e91\u96e8\u56fe\u5b9e\u9645\u4e0a\u662f\u5c0f\u63d0\u7434\u56fe（\u4e91）+\u6563\u70b9\u56fe（\u96e8），\u6240\u4ee5\u8fd9\u4e2a\u4e0a\u9762\u7684\u8bf4\u6cd5\u5c31\u50cf\u662f\u7c73\u996d-\u86cb\u7092\u996d \u6742\u4ea4\u83dc\u8c31\u2026\u2026<\/p>\n
\u81ea\u5236，\u52a0\u4e86\u56fe\u4f8b，\u770b\u8d77\u6765\u9ad8\u7ea7\u4e86\u4e00\u70b9\u70b9\u3002<\/p>\n
![\"\"](\"https:\/\/www.wsisp.com\/helps\/wp-content\/uploads\/2026\/01\/20260127151833-6978d749bba05.png\")
<\/p>\n
\u4ee3\u7801：<\/p>\n
clear; clc; close all;
\n% \u5c0f\u63d0\u7434\u56fe\u4eec
\n%% ===================== \u81ea\u5b9a\u4e49\u533a（\u53ef\u4fee\u6539\u53c2\u6570）=====================
\n% 1. \u6570\u636e\u53c2\u6570（8\u4e2a\u56fd\u5bb6，\u6bcf\u4e2a\u56fd\u5bb6\u5bf9\u5e94\u4e00\u7ec4\u6a21\u62df\u6570\u636e）
\ncountry_names = {'Burkina Faso','Ghana','Grenada','Kuwait','North Macedonia','San Marino','Saudi Arabia','Turkmenistan'};
\nn_countries = length(country_names);
\n% \u6a21\u62df\u6570\u636e：\u6bcf\u4e2a\u56fd\u5bb6\u751f\u6210\u4e00\u7ec4\u5e26\u5206\u5e03\u7684\u6570\u503c
\ndata = {
\n randn(200,1)*0.3 + 1;
\n randn(200,1)*0.4 + 1.5;
\n randn(200,1)*0.3 + 1.2;
\n randn(200,1)*0.4 + 1.1;
\n randn(200,1)*0.3 + 0.9;
\n randn(200,1)*0.5 + 1.8;
\n randn(200,1)*0.3 + 1.0;
\n randn(200,1)*0.4 + 1.1;
\n};
\n% \u6bcf\u4e2a\u56fd\u5bb6\u5bf9\u5e94\u7684p\u503c
\np_values = {'p=1.2e-26','p=1.0e+00','p=1.0e+00','p=6.8e-01','p=1.2e-26','p=8.5e-01','p=1.1e-03','p=3.0e-09'};
\n% 2. \u989c\u8272\u53c2\u6570（8\u79cd\u67d4\u548c\u914d\u8272）
\ncolor_list = [
\n 0.2 0.5 0.8; % \u84dd\u8272
\n 0.7 0.5 0.8; % \u6d45\u7d2b
\n 0.5 0.8 0.6; % \u6d45\u7eff
\n 0.3 0.6 0.8; % \u6e56\u84dd
\n 0.7 0.7 0.4; % \u6d45\u9ec4
\n 0.7 0.8 0.9; % \u6d45\u84dd
\n 0.2 0.4 0.5; % \u6df1\u84dd
\n 0.6 0.4 0.7; % \u7d2b\u7c89
\n];
\n% \u7bb1\u7ebf\u56fe\u649e\u8272\u989c\u8272（\u6bcf\u4e2a\u56fd\u5bb6\u7684\u8fb9\u6846\u989c\u8272）
\nedge_colors = [
\n 0.8 0.2 0.2; % \u7ea2\u8272 – \u84dd\u8272\u7cfb\u7684\u649e\u8272
\n 0.2 0.8 0.2; % \u7eff\u8272 – \u6d45\u7d2b\u7cfb\u7684\u649e\u8272
\n 0.8 0.4 0.1; % \u6a59\u8272 – \u6d45\u7eff\u7cfb\u7684\u649e\u8272
\n 0.9 0.6 0.1; % \u91d1\u8272 – \u6e56\u84dd\u7cfb\u7684\u649e\u8272
\n 0.5 0.2 0.8; % \u7d2b\u8272 – \u6d45\u9ec4\u7cfb\u7684\u649e\u8272
\n 0.8 0.5 0.2; % \u68d5\u8272 – \u6d45\u84dd\u7cfb\u7684\u649e\u8272
\n 0.9 0.3 0.1; % \u6a59\u7ea2 – \u6df1\u84dd\u7cfb\u7684\u649e\u8272
\n 0.2 0.7 0.5; % \u9752\u7eff – \u7d2b\u7c89\u7cfb\u7684\u649e\u8272
\n];
\n% 3. \u6587\u5b57\u53c2\u6570
\nylabel_txt = 'Value';
\nfontSize = 11;
\nlabelFontSize = 13;<\/p>\n
%% ===================== \u7edf\u4e00\u6a21\u677f\u6267\u884c\u533a =====================
\nfigure('Position',[150,150,1200,600],'Color','white');
\nhold on;<\/p>\n
% \u5b58\u50a8\u56fe\u4f8b\u53e5\u67c4 – \u53ea\u5b58\u50a8\u6563\u70b9\u56fe\u7684\u53e5\u67c4
\nscatter_handles = gobjects(1, n_countries);<\/p>\n
for i = 1:n_countries
\n % \u63d0\u53d6\u5f53\u524d\u56fd\u5bb6\u6570\u636e
\n d = data{i};
\n % \u8ba1\u7b97\u5c0f\u63d0\u7434\u56fe\u6838\u5bc6\u5ea6
\n [f,xi] = ksdensity(d);
\n f = f\/max(f)*0.3; % \u7f29\u653e\u5bc6\u5ea6\u503c<\/p>\n
% \u7ed8\u5236\u5c0f\u63d0\u7434\u56fe（\u586b\u5145\u534a\u900f\u660e\u533a\u57df）
\n x_left = i – f;
\n x_right = i + f;
\n fill([x_left,fliplr(x_right)],[xi,fliplr(xi)],color_list(i,:),'EdgeColor','none','FaceAlpha',0.4,'HandleVisibility','off');<\/p>\n
% \u7ed8\u5236\u6563\u70b9（\u96e8\u4e91\u70b9）
\n scatter_handles(i) = scatter(i + randn(length(d),1)*0.15, d, 15, …
\n color_list(i,:), 'filled', 'Marker','o', …
\n 'DisplayName', country_names{i});<\/p>\n
% \u7ed8\u5236\u7bb1\u7ebf\u56fe（\u81ea\u5b9a\u4e49\u989c\u8272）
\n bp = boxplot(d, 'Positions',i, 'Widths',0.1, 'Symbol','+', 'MedianStyle','line');
\n %'MedianStyle'\u53ef\u4e3atarget
\n % 'Symbol'\u53ef\u4ee5\u4e0d\u4e3a\u7a7a，'o'\u6216\u8005'.'\u5c55\u793aoutliers
\n % \u8bbe\u7f6e\u7bb1\u7ebf\u56fe\u989c\u8272
\n set(bp, 'LineWidth', 0.1);<\/p>\n
% \u83b7\u53d6\u7bb1\u7ebf\u56fe\u7684\u5404\u4e2a\u90e8\u5206\u5e76\u8bbe\u7f6e\u989c\u8272
\n h = findobj(gca, 'Tag', 'Box');
\n for j = 1:length(h)
\n % \u627e\u5230\u5f53\u524d\u4f4d\u7f6e\u7684\u7bb1\u7ebf\u56fe
\n if abs(get(h(j), 'XData') – i) < 0.1
\n % \u8bbe\u7f6e\u7bb1\u4f53\u586b\u5145\u989c\u8272
\n patch_obj = get(h(j), 'Parent');
\n % \u83b7\u53d6\u7bb1\u4f53\u586b\u5145
\n fill_h = findobj(patch_obj, 'Type', 'patch');
\n if ~isempty(fill_h)
\n set(fill_h, 'FaceColor', color_list(i,:), …
\n 'FaceAlpha', 0.1, …
\n 'EdgeColor', edge_colors(i,:), …
\n 'LineWidth', 0.5);
\n end<\/p>\n
% \u8bbe\u7f6e\u4e2d\u4f4d\u7ebf\u989c\u8272
\n median_h = findobj(gca, 'Tag', 'Median');
\n for k = 1:length(median_h)
\n if abs(get(median_h(k), 'XData') – i) < 0.1
\n set(median_h(k), 'Color', edge_colors(i,:), 'LineWidth', 2);
\n end
\n end<\/p>\n
% \u8bbe\u7f6e\u987b\u7ebf\u989c\u8272
\n whisker_h = findobj(gca, 'Tag', 'Whisker');
\n for k = 1:length(whisker_h)
\n if abs(get(whisker_h(k), 'XData') – i) < 0.1
\n set(whisker_h(k), 'Color', edge_colors(i,:), 'LineWidth', 1.5);
\n end
\n end<\/p>\n
% \u8bbe\u7f6e\u5f02\u5e38\u70b9\u989c\u8272
\n outlier_h = findobj(gca, 'Tag', 'Outliers');
\n for k = 1:length(outlier_h)
\n if abs(get(outlier_h(k), 'XData') – i) < 0.1
\n set(outlier_h(k), 'MarkerEdgeColor', edge_colors(i,:), …
\n 'MarkerFaceColor', edge_colors(i,:), 'MarkerSize', 4);
\n end
\n end
\n end
\n end<\/p>\n
% % \u521b\u5efa\u7bb1\u7ebf\u56fe\u7684\u4ee3\u7406\u56fe\u4f8b\u53e5\u67c4（\u6ce8\u91ca\u6389，\u9700\u8981\u65f6\u53d6\u6d88\u6ce8\u91ca）
\n % box_handles(i) = patch([NaN, NaN, NaN, NaN], [NaN, NaN, NaN, NaN], …
\n % color_list(i,:), 'EdgeColor', edge_colors(i,:), 'LineWidth', 1.5, …
\n % 'FaceAlpha', 0.3, 'DisplayName', [country_names{i} ' Box']);<\/p>\n
% \u6dfb\u52a0p\u503c\u6587\u672c（\u52a0\u7c97\u5b57\u4f53）
\n text(i, 3.8, p_values{i}, 'FontName','Times New Roman', 'FontSize',13, 'FontWeight','bold');
\nend<\/p>\n
% \u5750\u6807\u8f74\u8bbe\u7f6e（\u6240\u6709\u5b57\u4f53\u52a0\u7c97）
\nax = gca;
\nax.XTick = 1:n_countries;
\nax.XTickLabel = country_names;
\nax.XTickLabelRotation = 45;
\nax.FontName = 'Times New Roman';
\nax.FontSize = fontSize;
\nax.FontWeight = 'bold'; % \u5750\u6807\u8f74\u523b\u5ea6\u52a0\u7c97
\nax.YLabel.String = ylabel_txt;
\nax.YLabel.FontWeight = 'bold';
\nax.YLabel.FontSize = labelFontSize;
\nax.YLabel.FontName = 'Times New Roman';
\nax.YLim = [0, 4];
\nax.TickDir = 'out';
\nax.Box = 'on';<\/p>\n
% \u56fe\u4f8b\u8bbe\u7f6e – \u5206\u4e3a\u4e24\u7c7b：\u6563\u70b9\u56fe\u4e00\u7c7b，\u7bb1\u7ebf\u56fe\u4e00\u7c7b
\n% \u76ee\u524d\u53ea\u663e\u793a\u6563\u70b9\u56fe\u56fe\u4f8b，\u7bb1\u7ebf\u56fe\u56fe\u4f8b\u5df2\u6ce8\u91ca
\nlegend(scatter_handles, 'Location', 'northeast', …
\n 'FontSize', 10, 'FontName', 'Times New Roman', 'FontWeight', 'bold', …
\n 'NumColumns', 2);<\/p>\n
hold off;<\/p>\n
% \u7ed3\u679c\u8f93\u51fa
\nfprintf('\u7ed8\u56fe\u5b8c\u6210：\u5c0f\u63d0\u7434-\u96e8\u4e91\u6df7\u5408\u56fe，\u5305\u542b8\u4e2a\u56fd\u5bb6\u7684\u6570\u636e\u5206\u5e03\\\\n'); <\/p>\n
\u518d\u6765\u4e00\u4e2a\u57fa\u672c\u4e0a\u4e00\u6a21\u4e00\u6837\u7684，<\/p>\n
\u53c2\u8003：<\/p>\n
![\"\"](\"https:\/\/www.wsisp.com\/helps\/wp-content\/uploads\/2026\/01\/20260127151834-6978d74a1ccad.png\")
<\/p>\n
\u81ea\u5236：<\/p>\n
![\"\"](\"https:\/\/www.wsisp.com\/helps\/wp-content\/uploads\/2026\/01\/20260127151834-6978d74a6fe2f.png\")
<\/p>\n
\u4ee3\u7801：<\/p>\n
clear; clc; close all;<\/p>\n
%% \u81ea\u5b9a\u4e49\u533a
\ngroup_names = {'Conventional agriculture','Conservation agriculture'};
\nn_groups = length(group_names);<\/p>\n
data_ambient = {
\n randn(200,1)*1.2 + 8.5;
\n randn(200,1)*1.5 + 12.5;
\n};
\ndata_warming = {
\n randn(200,1)*1.2 + 8.5;
\n randn(200,1)*1.5 + 13.5;
\n};<\/p>\n
sig_marks = {'compare1','compare2'};<\/p>\n
% \u989c\u8272\u8bbe\u7f6e
\ncolor_ambient = [0.3 0.6 0.9]; % \u6d45\u84dd\u8272
\ncolor_warming = [0.9 0.5 0.5]; % \u6d45\u7ea2\u8272
\ncolor_box = [0.95 0.95 0.98]; % \u7bb1\u4f53\u989c\u8272
\ncolor_target = [0.2 0.8 0.2]; % \u76ee\u6807\u7ebf\u989c\u8272
\ncolor_outlier = [0.5 0.5 0.5]; % \u5f02\u5e38\u70b9\u989c\u8272<\/p>\n
ylabel_txt = 'SOC (g C kg^{-1})';
\nfontSize = 11;
\nlabelFontSize = 13;<\/p>\n
violin_width = 0.35;
\nbox_width = 0.15; % \u7bb1\u4f53\u5bbd\u5ea6
\ngroup_gap = 0; % \u7ec4\u95f4\u7a7a\u9699
\nmedian_line_width = 2;
\nsig_offset_y = 0.5;
\nbox_alpha = 0.3; % \u7bb1\u4f53\u900f\u660e\u5ea6
\noutlier_size = 0.6; % \u5f02\u5e38\u70b9\u5927\u5c0f<\/p>\n
%% \u7edf\u4e00\u6a21\u677f\u6267\u884c\u533a
\nfigure('Position',[150,150,1200,600],'Color','white');
\ntiledlayout(1,2,'TileSpacing','compact','Padding','compact');<\/p>\n
% \u5b50\u56fe1：\u6328\u5728\u4e00\u8d77\u7684\u534a\u5c0f\u63d0\u7434\u56fe
\nnexttile;
\nax1 = gca;
\nhold(ax1, 'on');<\/p>\n
legend_handles = [];<\/p>\n
for i = 1:n_groups
\n d_ambient = data_ambient{i};
\n d_warming = data_warming{i};<\/p>\n
[f_ambient,xi_ambient] = ksdensity(d_ambient);
\n [f_warming,xi_warming] = ksdensity(d_warming);<\/p>\n
f_ambient = f_ambient\/max(f_ambient)*violin_width;
\n f_warming = f_warming\/max(f_warming)*violin_width;<\/p>\n
% \u8c03\u6574x\u4f4d\u7f6e，\u8003\u8651\u7ec4\u95f4\u7a7a\u9699
\n x_center = i;
\n x_left = x_center – group_gap\/2;
\n x_right = x_center + group_gap\/2;<\/p>\n
% \u5de6\u4fa7\u534a\u5c0f\u63d0\u7434（Ambient）
\n if i == 1
\n h_ambient = fill([x_left-f_ambient, fliplr(x_left*ones(size(f_ambient)))],…
\n [xi_ambient, fliplr(xi_ambient)], color_ambient,…
\n 'EdgeColor', color_ambient*0.8, 'LineWidth', 0.5, 'FaceAlpha', 0.6, 'DisplayName', 'Ambient');
\n legend_handles = [legend_handles, h_ambient];
\n else
\n fill([x_left-f_ambient, fliplr(x_left*ones(size(f_ambient)))],…
\n [xi_ambient, fliplr(xi_ambient)], color_ambient,…
\n 'EdgeColor', color_ambient*0.8, 'LineWidth', 0.5, 'FaceAlpha', 0.6, 'HandleVisibility', 'off');
\n end<\/p>\n
% \u53f3\u4fa7\u534a\u5c0f\u63d0\u7434（Warming）
\n if i == 1
\n h_warming = fill([x_right*ones(size(f_warming)), fliplr(x_right+f_warming)],…
\n [xi_warming, fliplr(xi_warming)], color_warming,…
\n 'EdgeColor', color_warming*0.8, 'LineWidth', 0.5, 'FaceAlpha', 0.6, 'DisplayName', 'Warming');
\n legend_handles = [legend_handles, h_warming];
\n else
\n fill([x_right*ones(size(f_warming)), fliplr(x_right+f_warming)],…
\n [xi_warming, fliplr(xi_warming)], color_warming,…
\n 'EdgeColor', color_warming*0.8, 'LineWidth', 0.5, 'FaceAlpha', 0.6, 'HandleVisibility', 'off');
\n end<\/p>\n
% \u534a\u7bb1\u578b\u56fe – \u5de6\u4fa7（Ambient）
\n q_ambient = quantile(d_ambient, [0.25, 0.5, 0.75]);
\n iqr_ambient = q_ambient(3) – q_ambient(1);
\n whisker_ambient_low = max(q_ambient(1)-1.5*iqr_ambient, min(d_ambient));
\n whisker_ambient_high = min(q_ambient(3)+1.5*iqr_ambient, max(d_ambient));<\/p>\n
% \u5de6\u4fa7\u7bb1\u4f53（\u4f7f\u7528patch\u4ee3\u66ffrectangle）
\n box_x = [x_left-box_width, x_left, x_left, x_left-box_width, x_left-box_width];
\n box_y = [q_ambient(1), q_ambient(1), q_ambient(3), q_ambient(3), q_ambient(1)];
\n patch(box_x, box_y, color_box, 'EdgeColor', 'k', 'LineWidth', 1, 'FaceAlpha', box_alpha, 'HandleVisibility', 'off');<\/p>\n
% \u5de6\u4fa7\u76ee\u6807\u7ebf（\u539f\u6765\u7684\u4e2d\u4f4d\u7ebf）
\n line([x_left-box_width, x_left], [q_ambient(2), q_ambient(2)],…
\n 'Color', color_target, 'LineWidth', median_line_width, 'LineStyle', '-');<\/p>\n
% \u5de6\u4fa7\u987b\u7ebf
\n line([x_left-box_width*0.5, x_left-box_width*0.5], [whisker_ambient_low, q_ambient(1)],…
\n 'Color', 'k', 'LineWidth', 1);
\n line([x_left-box_width*0.5, x_left-box_width*0.5], [q_ambient(3), whisker_ambient_high],…
\n 'Color', 'k', 'LineWidth', 1);<\/p>\n
% \u5de6\u4fa7\u5f02\u5e38\u70b9
\n outliers_ambient = d_ambient(d_ambient < whisker_ambient_low | d_ambient > whisker_ambient_high);
\n if ~isempty(outliers_ambient)
\n plot(repmat(x_left-box_width*0.5, length(outliers_ambient), 1), outliers_ambient, '+',…
\n 'Color', color_outlier, 'MarkerSize', outlier_size, 'LineWidth', 1.5, 'HandleVisibility', 'off');
\n end<\/p>\n
% \u534a\u7bb1\u578b\u56fe – \u53f3\u4fa7（Warming）
\n q_warming = quantile(d_warming, [0.25, 0.5, 0.75]);
\n iqr_warming = q_warming(3) – q_warming(1);
\n whisker_warming_low = max(q_warming(1)-1.5*iqr_warming, min(d_warming));
\n whisker_warming_high = min(q_warming(3)+1.5*iqr_warming, max(d_warming));<\/p>\n
% \u53f3\u4fa7\u7bb1\u4f53（\u4f7f\u7528patch）
\n box_x = [x_right, x_right+box_width, x_right+box_width, x_right, x_right];
\n box_y = [q_warming(1), q_warming(1), q_warming(3), q_warming(3), q_warming(1)];
\n patch(box_x, box_y, color_box, 'EdgeColor', 'k', 'LineWidth', 1, 'FaceAlpha', box_alpha, 'HandleVisibility', 'off');<\/p>\n
% \u53f3\u4fa7\u76ee\u6807\u7ebf
\n line([x_right, x_right+box_width], [q_warming(2), q_warming(2)],…
\n 'Color', color_target, 'LineWidth', median_line_width, 'LineStyle', '-');<\/p>\n
% \u53f3\u4fa7\u987b\u7ebf
\n line([x_right+box_width*0.5, x_right+box_width*0.5], [whisker_warming_low, q_warming(1)],…
\n 'Color', 'k', 'LineWidth', 1);
\n line([x_right+box_width*0.5, x_right+box_width*0.5], [q_warming(3), whisker_warming_high],…
\n 'Color', 'k', 'LineWidth', 1);<\/p>\n
% \u53f3\u4fa7\u5f02\u5e38\u70b9
\n outliers_warming = d_warming(d_warming < whisker_warming_low | d_warming > whisker_warming_high);
\n if ~isempty(outliers_warming)
\n plot(repmat(x_right+box_width*0.5, length(outliers_warming), 1), outliers_warming, '+',…
\n 'Color', color_outlier, 'MarkerSize', outlier_size, 'LineWidth', 1.5, 'HandleVisibility', 'off');
\n end<\/p>\n
% \u6dfb\u52a0\u663e\u8457\u6027\u6807\u8bb0
\n text(x_center, max([d_ambient; d_warming]) + sig_offset_y, sig_marks{i},…
\n 'FontName','Times New Roman','FontSize',12,'FontWeight','bold',…
\n 'HorizontalAlignment','center','Color',[0.3 0.3 0.3]);
\nend<\/p>\n
ax1.XTick = 1:n_groups;
\nax1.XTickLabel = group_names;
\nax1.XTickLabelRotation = 15;
\nax1.FontName = 'Times New Roman';
\nax1.FontSize = fontSize;
\nax1.YLabel.String = ylabel_txt;
\nax1.YLabel.FontWeight = 'bold';
\nax1.YLabel.FontSize = labelFontSize;
\nax1.YLabel.FontName = 'Times New Roman';
\nax1.YLim = [4, 20];
\nax1.TickDir = 'out';
\nax1.Box = 'on';<\/p>\n
legend(legend_handles, {'Ambient','Warming'}, 'Location','northeast',…
\n 'FontSize',fontSize,'FontName','Times New Roman');<\/p>\n
% \u5b50\u56fe2：\u5206\u5f00\u7684\u534a\u5c0f\u63d0\u7434\u56fe
\nnexttile;
\nax2 = gca;
\nhold(ax2, 'on');<\/p>\n
group2_names = {'Conventional','Conservation'};
\nn_group2 = length(group2_names);<\/p>\n
% \u8c03\u6574\u5206\u5f00\u56fe\u7684\u4f4d\u7f6e，\u8bbe\u7f6e\u66f4\u5927\u7684\u95f4\u9699
\ngroup_spacing = 1.2;
\nx_positions = [1, 1.6, 3.2, 3.8]; % \u4e24\u7ec4\u5206\u5f00\u7684x\u4f4d\u7f6e
\ngroup_centers = [1.3, 3.5]; % \u6bcf\u7ec4\u4e2d\u5fc3\u4f4d\u7f6e<\/p>\n
for g = 1:n_group2
\n i = g;
\n d_ambient = data_ambient{i};
\n d_warming = data_warming{i};<\/p>\n
[f_ambient,xi_ambient] = ksdensity(d_ambient);
\n [f_warming,xi_warming] = ksdensity(d_warming);<\/p>\n
f_ambient = f_ambient\/max(f_ambient)*violin_width;
\n f_warming = f_warming\/max(f_warming)*violin_width;<\/p>\n
x_left = x_positions(2*g-1);
\n x_right = x_positions(2*g);<\/p>\n
% \u5de6\u4fa7\u534a\u5c0f\u63d0\u7434（Ambient）
\n fill([x_left-f_ambient, fliplr(x_left*ones(size(f_ambient)))],…
\n [xi_ambient, fliplr(xi_ambient)], color_ambient,…
\n 'EdgeColor', color_ambient*0.8, 'LineWidth', 0.5, 'FaceAlpha', 0.6, 'HandleVisibility', 'off');<\/p>\n
% \u53f3\u4fa7\u534a\u5c0f\u63d0\u7434（Warming）
\n fill([x_right*ones(size(f_warming)), fliplr(x_right+f_warming)],…
\n [xi_warming, fliplr(xi_warming)], color_warming,…
\n 'EdgeColor', color_warming*0.8, 'LineWidth', 0.5, 'FaceAlpha', 0.6, 'HandleVisibility', 'off');<\/p>\n
% \u5de6\u4fa7\u534a\u7bb1\u578b\u56fe（Ambient）
\n q_ambient = quantile(d_ambient, [0.25, 0.5, 0.75]);
\n iqr_ambient = q_ambient(3) – q_ambient(1);
\n whisker_ambient_low = max(q_ambient(1)-1.5*iqr_ambient, min(d_ambient));
\n whisker_ambient_high = min(q_ambient(3)+1.5*iqr_ambient, max(d_ambient));<\/p>\n
% \u5de6\u4fa7\u7bb1\u4f53（patch）
\n box_x = [x_left-box_width, x_left, x_left, x_left-box_width, x_left-box_width];
\n box_y = [q_ambient(1), q_ambient(1), q_ambient(3), q_ambient(3), q_ambient(1)];
\n patch(box_x, box_y, color_box, 'EdgeColor', 'k', 'LineWidth', 1, 'FaceAlpha', box_alpha, 'HandleVisibility', 'off');<\/p>\n
line([x_left-box_width, x_left], [q_ambient(2), q_ambient(2)],…
\n 'Color', color_target, 'LineWidth', median_line_width, 'LineStyle', '-');<\/p>\n
line([x_left-box_width*0.5, x_left-box_width*0.5], [whisker_ambient_low, q_ambient(1)],…
\n 'Color', 'k', 'LineWidth', 1);
\n line([x_left-box_width*0.5, x_left-box_width*0.5], [q_ambient(3), whisker_ambient_high],…
\n 'Color', 'k', 'LineWidth', 1);<\/p>\n
% \u5de6\u4fa7\u5f02\u5e38\u70b9
\n outliers_ambient = d_ambient(d_ambient < whisker_ambient_low | d_ambient > whisker_ambient_high);
\n if ~isempty(outliers_ambient)
\n plot(repmat(x_left-box_width*0.5, length(outliers_ambient), 1), outliers_ambient, '+',…
\n 'Color', color_outlier, 'MarkerSize', outlier_size, 'LineWidth', 1.5, 'HandleVisibility', 'off');
\n end<\/p>\n
% \u53f3\u4fa7\u534a\u7bb1\u578b\u56fe（Warming）
\n q_warming = quantile(d_warming, [0.25, 0.5, 0.75]);
\n iqr_warming = q_warming(3) – q_warming(1);
\n whisker_warming_low = max(q_warming(1)-1.5*iqr_warming, min(d_warming));
\n whisker_warming_high = min(q_warming(3)+1.5*iqr_warming, max(d_warming));<\/p>\n
% \u53f3\u4fa7\u7bb1\u4f53（patch）
\n box_x = [x_right, x_right+box_width, x_right+box_width, x_right, x_right];
\n box_y = [q_warming(1), q_warming(1), q_warming(3), q_warming(3), q_warming(1)];
\n patch(box_x, box_y, color_box, 'EdgeColor', 'k', 'LineWidth', 1, 'FaceAlpha', box_alpha, 'HandleVisibility', 'off');<\/p>\n
line([x_right, x_right+box_width], [q_warming(2), q_warming(2)],…
\n 'Color', color_target, 'LineWidth', median_line_width, 'LineStyle', '-');<\/p>\n
line([x_right+box_width*0.5, x_right+box_width*0.5], [whisker_warming_low, q_warming(1)],…
\n 'Color', 'k', 'LineWidth', 1);
\n line([x_right+box_width*0.5, x_right+box_width*0.5], [q_warming(3), whisker_warming_high],…
\n 'Color', 'k', 'LineWidth', 1);<\/p>\n
% \u53f3\u4fa7\u5f02\u5e38\u70b9
\n outliers_warming = d_warming(d_warming < whisker_warming_low | d_warming > whisker_warming_high);
\n if ~isempty(outliers_warming)
\n plot(repmat(x_right+box_width*0.5, length(outliers_warming), 1), outliers_warming, '+',…
\n 'Color', color_outlier, 'MarkerSize', outlier_size, 'LineWidth', 1.5, 'HandleVisibility', 'off');
\n end<\/p>\n
% \u6dfb\u52a0\u663e\u8457\u6027\u6807\u8bb0
\n text(group_centers(g), max([d_ambient; d_warming]) + sig_offset_y, sig_marks{g},…
\n 'FontName','Times New Roman','FontSize',12,'FontWeight','bold',…
\n 'HorizontalAlignment','center','Color',[0.3 0.3 0.3]);
\nend<\/p>\n
ax2.XTick = group_centers;
\nax2.XTickLabel = group2_names;
\nax2.FontName = 'Times New Roman';
\nax2.FontSize = fontSize;
\nax2.YLim = [4, 20];
\nax2.TickDir = 'out';
\nax2.Box = 'on';<\/p>\n
% \u6dfb\u52a0\u53f3\u4fa7\u5b50\u56fe\u7684\u56fe\u4f8b
\nlegend_handles2 = [];
\nh1 = patch([0,0,0,0],[0,0,0,0],color_ambient,'EdgeColor',color_ambient*0.8,'LineWidth',0.5,'FaceAlpha',0.6,'DisplayName','Ambient');
\nh2 = patch([0,0,0,0],[0,0,0,0],color_warming,'EdgeColor',color_warming*0.8,'LineWidth',0.5,'FaceAlpha',0.6,'DisplayName','Warming');
\nlegend_handles2 = [h1, h2];
\nlegend(legend_handles2, {'Ambient','Warming'}, 'Location','northeast',…
\n 'FontSize',fontSize,'FontName','Times New Roman');<\/p>\n
hold(ax1, 'off');
\nhold(ax2, 'off');<\/p>\n
fprintf('\u7ed8\u56fe\u5b8c\u6210\\\\n'); <\/p>\n","protected":false},"excerpt":{"rendered":"
\u76ee\u5f55
\n\u6298\u7ebf\u56fe\u2014\u2014
\n2D\u7d2f\u8ba1\u6298\u7ebf\u56fe
\n3D\u5e26\u900f\u660e\u5ea6\u5207\u7247
\n2D \u591a\u7ebf\u5bf9\u5e94\u540c\u56fe\u4f8b，\u9650\u5b9a\u8fb9\u754c
\n2D \u53cc\u5750\u6807，\u5e26\u5149\u6655\u6548\u679c，\u8986\u76d6\u6e10\u53d8\u989c\u8272\u586b\u5145\u533a\u57df，\u5e26\u56fe\u4f8b
\n\u67f1\u72b6\u56fe\u2014\u2014
\n3D \u4e0d\u540c\u989c\u8272 \u67f1\u72b6\u56fe
\n2D \u53cc\u5750\u6807 \u67f1\u72b6\u56fe\u6298\u7ebf\u56fe\u6742\u4ea4
\n\u5176\u4ed6\u7edf\u8ba1\u5206\u5e03\u56fe\u2014\u2014
\n\u4e91\u96e8\u56fe\u548c\u7bb1\u578b\u56fe\u7684\u6742\u4ea4\u56fe？
\n\u534a\u5c0f\u63d0\u7434\u56fe \u6298\u7ebf\u56fe\u2014\u2014
\n2D\u7d2f\u8ba1\u6298\u7ebf\u56fe \u81ea\u5236 \u4ee3\u7801：
\nclear; clc; close all;%% \u81ea\u5b9a\u4e49\u533a
\ny<\/p>\n","protected":false},"author":2,"featured_media":66971,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[1],"tags":[1284,50,87],"topic":[],"class_list":["post-66989","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-server","tag-matlab","tag-50","tag-87"],"yoast_head":"\n