Matlab绘制95%置信区间图

A-Chin

已于 2022-04-01 19:21:57 修改

阅读量2.2w

点赞数 15

分类专栏： matlab 文章标签： 1024程序员节 matlab 置信区间

于 2021-10-24 19:12:34 首次发布

本文链接：https://blog.csdn.net/L_J_Kin/article/details/120938578

版权

matlab 专栏收录该内容

35 篇文章 19 订阅

订阅专栏

在平时绘图中，常会绘制95%置信区间图：

在这里插入图片描述

主要使用到patch函数，绘制时需要先对x进行排序，方便绘制区间图，参考代码如下：

x = data(:, 1);
y = data(:, 2);

% 给x进行排序，方便绘制区间图
[X, I] = sort(x);
Y = y(I);
    
x = X;
y = Y;
    
[p, s] = polyfit(x, y, 1);
y1 = polyval(p, x);
[yfit, dy] = polyconf(p, x, s, 'predopt', 'curve');
    
patch([x; flipud(x)], [yfit - dy; flipud(yfit + dy)], [231 231 231] / 255, 'FaceA', 0.95, 'EdgeA', 0);
hold on;
    
s = scatter(x, y, 'r');
s.LineWidth = 1.25;
plot(x, y1, 'Color', [63 63 63] / 255, 'linewidth', 2.5);
hold on;