众所周知,Excel是可以进行简单的函数拟合的,主要步骤如下:
第一步,作图
第二步,添加趋势线
但是,该方法的缺陷也是显而易见的,就是只能拟合预设模板的几个函数,如果碰上其他函数就直接歇菜了,比如:
(1)y=ax+by=\frac{a}{x+b} \tag1
其实,Excel里是有这个功能的,只不过不是通过趋势线而是通过规划求解的方式来实现的,具体操作看这里:
在Excel中打开文件菜单中的选项,在弹出菜单的加载项选项卡中找到Excel加载项,点一下右侧的转到,然后勾选规划求解加载项,确定保存
启用加载项后,在Excel的数据选项卡中就可以看到规划求解按钮了
那么,准备开工了:
按照如下的表中准备数据(后面有说明):
①xy列是原始数据
②新建拟合参数对应的列,并赋予初值,这里我们赋予的初始值是a=1,b=1
③使用ab的初始值计算对应的y(即yi),这里需要说明一下,对于C2单元格我们Excel公式是这么写的:
$H$2/(A2+$H$3)
其中$H$2表示绝对引用,H2的数值,即在复制公式的时候该值不变化;A2则表示相对引用A2的数值,在复制公式的时候对应变化。在引用了某个单元格后可以按F4切换绝对引用和相对引用。
例如,当我们将C2的单元格复制到C3的时候,公式就变成了:
$H$2/(A3+$H$3)
注意到绝对引用的位置不发生变化(仍然引用了H2,而相对引用则发生了变化,引用了A3)
获得yi后,计算yi与y偏差的平方,就是第四列
④对yi与y偏差的平方求和(其实就是残差平方和,SSE),这就算是我们的目标函数了。我们拟合的目的就是使得该值达到最小
点开规划求解,按下图操作
①设定目标函数,记得目标选择为最小值(图中没有画上)
②选择参数单元格,这里就是ab两个
③设定约束条件,可以分别约束每个变量的取值范围、是否是整数等,这里因为问题简单就没有约束
④选择求解方案,一般都是非线性求解
然后点击求解,由于问题很简单那,5次就迭代到了结果,注意a和b的数值已经变化了,拟合完成,结果为:
a=2.78639, b=4.168637。实际上我又用scipy.optimize拟合了一下,结果几乎是完全一样的
做个图,显示拟合结果(橙色)与原始数据(黄色)之间匹配良好
作为一个拟合结果,我们显然需要知道拟合的各种误差,常用的评价指标就是 R2R^2 ,不同于趋势线上可以直接显示R2R^2,这里我们要手动计算了:
我们之前已经用过残差平方和SSE的概念了(即每一项与其拟合值之间的偏差):
(2)SSE=∑(y−yi)2SSE=\sum_{}^{}{(y-y_i)^2} \tag2
那么再介绍一个概念,叫做总离差平方和(SST),即每一项与其平均值的偏差:
(3)SST=∑(y−mean(y))2SST=\sum_{}^{}{(y-{\rm mean}(y))^2} \tag3
由上面两项可以计算:(4)R2=1−SSESSTR^2=1-\frac{SSE}{SST} \tag4
原网址: 访问
创建于: 2023-08-16 16:05:52
目录: default
标签: 无
未标明原创文章均为采集,版权归作者所有,转载无需和我联系,请注明原出处,南摩阿彌陀佛,知识,不只知道,要得到
java windows火焰图_mob64ca12ec8020的技术博客_51CTO博客 - 在windows下不可行,不知道作者是怎样搞的 监听SpringBoot 服务启动成功事件并打印信息_监听springboot启动完毕-CSDN博客 SpringBoot中就绪探针和存活探针_management.endpoint.health.probes.enabled-CSDN博客 u2u转换板 - 嘉立创EDA开源硬件平台 Spring Boot 项目的轻量级 HTTP 客户端 retrofit 框架,快来试试它!_Java精选-CSDN博客 手把手教你打造一套最牛的知识笔记管理系统! - 知乎 - 想法有重合-理论可参考 安宇雨 闲鱼 机械键盘 客制化 开贴记录 文本 linux 使用find命令查找包含某字符串的文件_beijihukk的博客-CSDN博客_find 查找字符串 ---- mac 也适用 安宇雨 打字音 记录集合 B站 bilibili 自行搭建 开坑 真正的客制化 安宇雨 黑苹果开坑 查找工具包maven pom 引用地 工具网站 Dantelis 介绍的玩轴入坑攻略 --- 关于轴的一些说法 --- 非官方 ---- 心得而已 --- 长期开坑更新 [本人问题][新开坑位]关于自动化测试的工具与平台应用 机械键盘 开团 网站记录 -- 能做一个收集的程序就好了 不过现在没时间 -- 信息大多是在群里发的 - 你要让垃圾佬 都去一个地方看难度也是很大的 精神支柱 [超级前台]sprinbboot maven superdesk-app 记录 [信息有用] [环境准备] [基本完成] [sebp/elk] 给已创建的Docker容器增加新的端口映射 - qq_30599553的博客 - CSDN博客 [正在研究] Elasticsearch, Logstash, Kibana (ELK) Docker image documentation elasticsearch centos 安装记录 及 启动手记 正式服务器 39 elasticsearch 问题合集 不断更新 6.1.1 | 6.5.1 两个版本 博客程序 - 测试 - bug记录 等等问题 laravel的启动过程解析 - lpfuture - 博客园 OAuth2 Server PHP 用 Laravel 搭建带 OAuth2 验证的 RESTful 服务 | Laravel China 社区 - 高品质的 Laravel 和 PHP 开发者社区 利用Laravel 搭建oauth2 API接口 附 Unauthenticated 解决办法 - 煮茶的博客 - SegmentFault 思否 使用 OAuth2-Server-php 搭建 OAuth2 Server - 午时的海 - 博客园 基于PHP构建OAuth 2.0 服务端 认证平台 - Endv - 博客园 Laravel 的 Artisan 命令行工具 Laravel 的文件系统和云存储功能集成 浅谈Chromium中的设计模式--终--Observer模式 浅谈Chromium中的设计模式--二--pre/post和Delegate模式 浅谈Chromium中的设计模式--一--Chromium中模块分层和进程模型 DeepMind 4 Hacking Yourself README.md update 20211011
Laravel China 简书 知乎 博客园 CSDN博客 开源中国 Go Further Ryan是菜鸟 | LNMP技术栈笔记 云栖社区-阿里云 Netflix技术博客 Techie Delight Linkedin技术博客 Dropbox技术博客 Facebook技术博客 淘宝中间件团队 美团技术博客 360技术博客 古巷博客 - 一个专注于分享的不正常博客 软件测试知识传播 - 测试窝 有赞技术团队 阮一峰 语雀 静觅丨崔庆才的个人博客 软件测试从业者综合能力提升 - isTester IBM Java 开发 使用开放 Java 生态系统开发现代应用程序 pengdai 一个强大的博主 HTML5资源教程 | 分享HTML5开发资源和开发教程 蘑菇博客 - 专注于技术分享的博客平台 个人博客-leapMie 流星007 CSDN博客 - 舍其小伙伴 稀土掘金 Go 技术论坛 | Golang / Go 语言中国知识社区
最新评论