数学建模算法之拟合算法

最新推荐文章于 2023-04-30 23:37:39 发布

二月世

最新推荐文章于 2023-04-30 23:37:39 发布

阅读量376

点赞数

分类专栏：数学建模文章标签：算法 matlab 开发语言

本文链接：https://blog.csdn.net/zkx0214/article/details/127776449

版权

数学建模专栏收录该内容

6 篇文章 0 订阅

订阅专栏

选择线性拟合

原因是简单，且在做草图后发现线性相关高

clear;clc

load data.mat

n = size(year,2);
k = (n*sum(year.*po)-sum(po)*sum(year))/(n*sum(year.*year)-sum(year)*sum(year));
b = (sum(year.*year)*sum(po)-sum(year)*sum(year.*po))/(n*sum(year.*year)-sum(year)*sum(year));

f = @(year) k*year + b;
fplot(f,[2009,2018]);

hold on % 继续在之前的图形上来画图形

plot(year,po,'o');

po_hat = k*year+b; % y的拟合值
SSR = sum((po_hat-mean(po)).^2);  % 回归平方和
SSE = sum((po_hat-po).^2); % 误差平方和
SST = sum((po-mean(po)).^2); % 总体平方和
SST-SSE-SSR;  % 5.6843e-14  =   5.6843*10^-14   matlab浮点数计算的一个误差
R_2 = SSR / SST;