最小二乘法拟合出二阶响应面近似模型

背景:根据样本试验数据拟合出二阶响应面近似模型(正交二次型),并使用决定系数R²和调整的决定系数R²_adj来判断二阶响应面模型的拟合精度。


1、样本数据(来源:硕士论文《航空发动机用W形金属密封环密封性能分析与优化》)

编号 x1	     x2	     x3	     x4    	y1	    y2	    y3
1	0.67	0.70	1.78	0.83	804.8	246.7	37.7
2	0.59	0.73	2.09	0.72	686.3	213.5	40.1
3	0.69	0.78	1.81	0.86	703.6	221.0	39.3
4	0.57	0.70	1.84	0.73	822.3	253.6	37.1
5	0.58	0.80	2.07	0.74	632.4	201.6	40.9
6	0.66	0.77	1.95	0.71	716.9	224.7	39.7
7	0.61	0.69	2.04	0.75	730.1	223.3	39.4
8	0.66	0.83	1.88	0.76	676.6	215.8	40.0
9	0.71	0.66	1.85	0.77	845.0	254.6	37.8
10	0.64	0.72	1.73	0.82	813.2	251.8	37.2
11	0.71	0.71	1.99	0.81	709.0	217.4	40.1
12	0.60	0.77	1.91	0.79	698.2	219.5	39.2
13	0.70	0.65	1.94	0.80	794.9	237.9	38.7
14	0.74	0.74	1.96	0.89	663.2	205.7	40.8
15	0.74	0.67	1.74	0.87	839.3	253.8	37.5
16	0.56	0.68	1.87	0.84	773.6	221.6	37.7
17	0.62	0.84	2.06	0.77	598.0	193.0	41.8
18	0.64	0.72	2.01	0.78	703.4	217.3	39.8
19	0.68	0.79	1.79	0.70	777.5	245.8	38.4
20	0.65	0.84	1.81	0.90	640.6	205.3	40.1
21	0.73	0.75	1.76	0.88	743.3	230.9	38.8
22	0.63	0.82	1.70	0.84	726.0	231.4	38.3
23	0.60	0.83	1.90	0.81	648.0	207.5	40.1
24	0.56	0.76	1.97	0.76	693.6	217.9	39.3
25	0.72	0.81	2.03	0.85	600.7	191.7	42.2

2、将上述数据存放在sample_data.xlsx表格中,再编写matlab代码

% % 读取数据(假设数据已经保存在一个名为'data.txt'的文件中,使用空格或制表符分隔)
%% 近似模型:y = a0+a1*x1+a2*x2+a3*x3+a4*x4+a5*x1*x2+a6*x1*x3+a7*x1*x4+a8*x2*x3+a9*x2*x4+a10*x3*x4+a11*x1^2+a12*x2^2+a13*x3^2+a14*x4^2;
clc;clear;
filename = 'sample_data.xlsx';
data = xlsread(filename);% 提取自变量和因变量
X = data(:, 2:5); % x1, x2, x3, x4
y = data(:, 6:8); % y1, y2, y3% 构建扩展自变量矩阵(包括线性项、交互项和平方项)
[n, m] = size(X);
X_extended = ones(n, 1); % 初始化扩展矩阵,第一列为截距项% 添加线性项
X_extended = [X_extended, X];% 添加交互项
X_extended = [X_extended, X(:,1).*X(:,2), X(:,1).*X(:,3), X(:,1).*X(:,4), ...X(:,2).*X(:,3), X(:,2).*X(:,4), X(:,3).*X(:,4)];% 添加平方项
X_extended = [X_extended, X(:,1).^2, X(:,2).^2, X(:,3).^2, X(:,4).^2];% 使用最小二乘法求解回归系数(对每个因变量分别进行)
beta_values = zeros(15, 3); % 15个系数(包括截距项),3个因变量
R2_values = zeros(1, 3);
adjusted_R2_values = zeros(1, 3);for i = 1:3% 拟合模型beta = (X_extended' * X_extended) \ (X_extended' * y(:, i));beta_values(:, i) = beta;% 预测值y_pred = X_extended * beta;% 计算SST, SSE, SSRSST = sum((y(:, i) - mean(y(:, i))).^2);SSE = sum((y(:, i) - y_pred).^2);SSR = SST - SSE;% 计算R^2R2_values(i) = SSR / SST;% 计算调整后的R^2p = size(X_extended, 2); % 自变量数量(包括截距项)adjusted_R2_values(i) = 1 - (1 - R2_values(i)) * ((n - 1) / (n - p));
end% 显示结果
disp('回归系数(包括截距项):');
disp(beta_values);
disp('决定系数 R^2:');
disp(R2_values);
disp('调整的决定系数 Adjusted R^2:');
disp(adjusted_R2_values);

3、运行结果

4、如上,决定系数和调整的决定系数均大于0.95,可见,近似模型满足精度要求。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.rhkb.cn/news/487755.html

如若内容造成侵权/违法违规/事实不符,请联系长河编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

《操作系统 - 清华大学》6 -7:局部页面置换算法:Belady现象

文章目录 1. 定义2. LRU、FIFO和Clock的比较 1. 定义 局部页面置换算法的特点是针对一个正在运行的程序,它访问内存的情况,访问页的情况,来决定应该采取什么样策略,把相应的页替换出去,站在算法本身角度来考虑置换哪个…

【开源免费】基于SpringBoot+Vue.JS在线办公系统(JAVA毕业设计)

本文项目编号 T 001 ,文末自助获取源码 \color{red}{T001,文末自助获取源码} T001,文末自助获取源码 目录 一、系统介绍二、演示录屏三、启动教程四、功能截图五、文案资料5.1 选题背景5.2 国内外研究现状5.3 可行性分析 六、核心代码6.1 查…

05-标准库开发-STM32-IIC协议

七、STM32中IIC协议 概述 Inter-Integrated Circuit (IIC),也常称为I2C(I squared C),是一种同步、串行、半双工通信总线协议。它主要用于连接低速外围设备到处理器或微控制器上,如MPU6050姿态传感器、OLED显示屏、存…

【linux系统】基础开发工具(yum、Vim)

1. 软件包管理器 1.1 什么是软件包 在Linux下安装软件, ⼀个通常的办法是下载到程序的源代码, 并进⾏编译, 得到可执⾏程序. 但是这样太麻烦了, 于是有些⼈把⼀些常⽤的软件提前编译好, 做成软件包(可以理解成windows上的安装程序)放在⼀个服务器上, 通过包管理器可以很⽅便的…

UFUG2601_project_Fall2024 MiniDB Project

PS:如果读过题了可以跳过题目描述直接到题解部分 链接:UFUG2601_project_Fall2024 MiniDB Project 文章目录 题目题解声明可完成操作运行逻辑大致思路数据存储数据类型数据名称 命令输入文件读入命令读入 操作2.1 Create Database and Use Database2.2 C…

this version of the Java Runtime only recognizes class file versions up to 52.0

问题描述 Exception in thread "main" java.lang.UnsupportedClassVersionError: com/xxx/Main has been compiled by a more recent version of the Java Runtime (class file version 61.0), this version of the Java Runtime only recognizes class file versi…

Tr0ll: 1 Vulnhub靶机渗透笔记

Tr0ll: 1 本博客提供的所有信息仅供学习和研究目的,旨在提高读者的网络安全意识和技术能力。请在合法合规的前提下使用本文中提供的任何技术、方法或工具。如果您选择使用本博客中的任何信息进行非法活动,您将独自承担全部法律责任。本博客明确表示不支…

CAP定理

2.1 CAP 定理的由来与证明 CAP 定理是计算机科学界的“铁律”,最早由 Eric Brewer 提出,后来被正式证明: 分布式系统里,一致性(C)、可用性(A)、分区容错性(P&#xff09…

【flutter】webview下载文件方法集锦

说明:android的webview是不支持下载的!!! 所以我们需要监听下载接口 然后手动执行下载操作,分为三种类型 直接打开浏览器下载(最简单),但是一些下载接口需要cookie信息时不能满足 …

Java版-图论-最短路-Floyd算法

实现描述 网络延迟时间示例 根据上面提示,可以计算出,最大有100个点,最大耗时为100*wi,即最大的耗时为10000,任何耗时计算出来超过这个值可以理解为不可达了;从而得出实现代码里面的: int maxTime 10005…

SQL注入基础入门篇 注入思路及常见的SQL注入类型总结

目录 前言一、了解mysql数据库1、了解sql增删改查2、了解sql查询 二、sql注入基础三、学习sql注入漏洞1、union注入1、判断数字型注入还是字符型型注入:2、判断闭合方式(字符型注入):3、判断回显位4、查询库名,表名&am…

基于Spring Boot库存管理系统

文末获取源码和万字论文,制作不易,感谢点赞支持。 基于Spring Boot库存管理系统 当下,如果还依然使用纸质文档来记录并且管理相关信息,可能会出现很多问题,比如原始文件的丢失,因为采用纸质文档&#xff0c…

JSSIP的使用及问题(webRTC,WebSockets)

简介 项目中有一个需要拨打电话的功能,要求实时的进行音频接听,并且可以在电话接听或者挂断等情况下做出相应的操作。jssip作为一个强大的实现实时通信的javascript库,这不门当户对了嘛。 jssip(官网: JsSIP - the J…

【Cadence32】PCB多层板电源、地平面层创建心得➕CM约束管理器Analyze分析显示设置➕“DP”报错DRC

【转载】Cadence Design Entry HDL 使用教程 【Cadence01】Cadence PCB Edit相对延迟与绝对延迟的显示问题 【Cadence02】Allegro引脚焊盘Pin设置为透明 【Cadence03】cadence不小心删掉钢网层怎么办? 【Cadence04】一般情况下Allegro PCB设计时的约束规则设置&a…

Java阶段三06

第3章-第6节 一、知识点 理解MVC三层模型、理解什么是SpringMVC、理解SpringMVC的工作流程、了解springMVC和Struts2的区别、学会使用SpringMVC封装不同请求、接收参数 二、目标 理解MVC三层模型 理解什么是SpringMVC 理解SpringMVC的工作流程 学会使用SpringMVC封装请求…

C/C++流星雨

系列文章 序号直达链接1C/C爱心代码2C/C跳动的爱心3C/C李峋同款跳动的爱心代码4C/C满屏飘字表白代码5C/C大雪纷飞代码6C/C烟花代码7C/C黑客帝国同款字母雨8C/C樱花树代码9C/C奥特曼代码10C/C精美圣诞树11C/C俄罗斯方块12C/C贪吃蛇13C/C孤单又灿烂的神-鬼怪14C/C闪烁的爱心15C/C…

Vmware Vcenter7.0证书web续期发生错误

1. 故障描述 vSphere Client 版本 7.0.2.00200 vCenter _MACHINE_CERT快到期了,通过web界面更新证书失败 第一步先这样,重新续订一下证书 续订发生错误 2. 解决办法 2.1. 前提工作 登陆ssh到vcenter,重新生成证书 先关掉HA&#xff…

【合作原创】使用Termux搭建可以使用的生产力环境(五)

前言 在上一篇【合作原创】使用Termux搭建可以使用的生产力环境(四)-CSDN博客我们讲到了如何让proot-distro中的Debian声音驱动正常,将我们的系统备份后,通过VNC客户端连接到VNC服务器,这一篇我们来讲一下xfce桌面的美…

uniapp -- 实现页面滚动触底加载数据

效果 首选,是在pages.json配置开启下拉刷新 {"path": "pages/my/document/officialDocument","style": {"navigationStyle":</

Python之爬虫入门--示例(2)

一、Requests库安装 可以使用命令提示符指令直接安装requests库使用 pip install requests 二、爬取JSON数据 &#xff08;1&#xff09;、点击网络 &#xff08;2&#xff09;、刷新网页 &#xff08;3&#xff09;、这里有一些数据类型&#xff0c;选择全部 &#xff08…