Matlab 使用subplot绘制多个子图,一元拟合

实现效果:

Matlab 复制代码
clc; clear;

filename = 'sri.xlsx'; % 确认文件路径

data = readtable(filename);
datavalue = data{:,2:end};
datavalue = datavalue';

fig = figure('Position', [0, 0, 1500, 900]);
indexString = ["(a)","(b)","(c)","(d)","(e)","(f)","(g)",];
subplotName = ["h1","h2","h3","h4","h5","h6","h7"];
for i = 1:7
    subplot(4,2,i);
    years = datavalue(:,1); 
    values = datavalue(:,i+1);
    
    coefficients = polyfit(years, values, 1);
    k = coefficients(1);
    b = coefficients(2);
    
    % 计算Y轴上下限值,将其设置为最大值和最小值的 1.1 倍
    y_max = max(values) * 1.3;
    y_min = min(values) * 1.3;

    plot(years, values, 'o-', 'MarkerSize', 4, 'MarkerFaceColor', '#58CCFA', 'Color', '#58CCFA', 'LineWidth', 1.5);
    hold on;
    plot(years, k*years + b, 'r--',  'LineWidth', 1);%拟合线
    xlim([1960 2018]);
    xtick_values = 1960:10:2010;  
    if ~ismember(2018, xtick_values)
        xtick_values = [xtick_values, 2018];  
    end
    xticks(xtick_values);
    ylabel('SRI', 'FontSize', 12);
    ylim([y_min y_max]); % 设置Y轴上下限值
    grid on;
    yrange = ylim; % 获取当前子图的y轴值域范围
    text(1961, 0.9 * yrange(2),indexString(i), 'FontSize', 14, 'HorizontalAlignment', 'left', 'VerticalAlignment', 'top'); % 调整文本位置坐标

    ax = gca; 
    ax.FontSize = 12; 
    annotation('textbox', [ax.Position(1), ax.Position(2), ax.Position(3), ax.Position(4)], 'String', sprintf('k = %.4f, b = %.2f', k, b), 'FontSize', 14, 'EdgeColor', 'none', 'Color', 'black', 'HorizontalAlignment', 'center'); % 调整注释框位置

    hold off;
end

% 在第 8 个子图位置添加文本
subplot(4,2,8);
str = ["(a):华北湿润半湿润暖     (b):西北荒漠地区";...
       "(c):华中华南湿润亚热     (d):东北湿润半湿润温";...
       "(e):内蒙草原地区            (f):青藏高原";...
       "(g):华南湿润热带地区"];

text(0.15, 0.4, sprintf('%s\n', str{:}), 'FontSize', 14, 'HorizontalAlignment', 'left', 'VerticalAlignment', 'middle');
axis off; % 关闭坐标系显示


print(fig, 'data.tif', '-dtiff', '-r2000'); % 保存为 TIFF 文件,300 DPI 分辨率
相关推荐
逻辑驱动的ken4 分钟前
Java高频面试考点场景题13
java·开发语言·jvm·面试·求职招聘·春招
Evand J5 分钟前
【MATLAB代码介绍】无迹粒子滤波(Unscented Particle Filter),运动目标滤波定位
开发语言·matlab·粒子滤波·upf·无迹
bubiyoushang8889 分钟前
MATLAB考虑源荷不确定性的电力系统优化
开发语言·matlab
自信1504130575928 分钟前
重生之从0开始学习c++之string(上)
开发语言·c++·学习
wangjialelele32 分钟前
从零入门 LangChain:Python 语法详解 + 工具开发 + 结构化输出实战
开发语言·人工智能·python·语言模型·langchain
jinanwuhuaguo33 分钟前
OpenClaw智能体的涌现与异化——复杂系统演化、知识权力重构与文明纪元跃迁(第五篇)
大数据·开发语言·人工智能·重构·安全架构·openclaw
Dillon Dong38 分钟前
【系列主题】拯救 OOM 与构建中断:Next.js 在 Docker 中的静态生成(SSG)避坑指南
开发语言·javascript·docker
AI人工智能+电脑小能手40 分钟前
【大白话说Java面试题】【Java基础篇】第10题:HashMap中的元素是有序存放的吗
java·开发语言·数据结构·后端·面试·哈希算法·哈希表
itzixiao1 小时前
L1-049 天梯赛座位分配(20 分)[java][python][c]
java·开发语言·python
子非鱼@Itfuture1 小时前
ThreadLocal 是什么?如何用?以及最佳使用场景
java·开发语言·spring