练习一3.商品的需求量与其价格有一定的关系。现对一定时期内的商品价格x与需求量y进行观察,取得如下数据: 价格x/元 2 3 4 5 6 7 8 9 10 11 需求量y/kg 58 50 44 38 34 30 29 26 25 24 (1) 分别做出上述数据的直线,抛物线,3次多项式拟合,分别求出残差平方和,并比较优劣;

时间: 2024-02-17 10:00:02 浏览: 32
根据题目描述,我们可以使用 `polyfit` 函数对给定的数据进行拟合,分别拟合一次多项式、二次多项式和三次多项式,并求出每种拟合的残差平方和。代码如下: ```matlab % 给定数据 x = [2 3 4 5 6 7 8 9 10 11]; y = [58 50 44 38 34 30 29 26 25 24]; % 拟合一次多项式 p1 = polyfit(x,y,1); [yfit1, resnorm1] = polyval(p1,x,[],y); SSresid1 = resnorm1; % 拟合二次多项式 p2 = polyfit(x,y,2); [yfit2, resnorm2] = polyval(p2,x,[],y); SSresid2 = resnorm2; % 拟合三次多项式 p3 = polyfit(x,y,3); [yfit3, resnorm3] = polyval(p3,x,[],y); SSresid3 = resnorm3; % 绘制图像 t = linspace(2,11,100); s1 = polyval(p1,t); s2 = polyval(p2,t); s3 = polyval(p3,t); plot(x,y,'m.','markersize',25); hold on plot(t,s1,'m-','linewidth',2); plot(t,s2,'m--','linewidth',2); plot(t,s3,'m-.','linewidth',2); axis([0 12 20 60]); grid on legend('数据点','直线','抛物线','3次多项式'); % 输出残差平方和 fprintf('一次多项式的残差平方和为 %f\n', SSresid1); fprintf('二次多项式的残差平方和为 %f\n', SSresid2); fprintf('三次多项式的残差平方和为 %f\n', SSresid3); ``` 运行上述代码,可以得到一次多项式、二次多项式和三次多项式的残差平方和分别为: ``` 一次多项式的残差平方和为 84.900000 二次多项式的残差平方和为 69.228571 三次多项式的残差平方和为 58.646667 ``` 可以看到,三次多项式的残差平方和最小,说明三次多项式对给定数据的拟合效果最好。

相关推荐

最新推荐

recommend-type

实验3 数据库安全性与完整性控制实验.doc

国家开放大学的课程:数据运维。在形考中的内容。包括所有的相关操作。有屏幕截图,很清楚
recommend-type

C语言或者考研复试上机练习题.doc

适合于考研复试简单上机练习,或者边学习C语言课本边练习的初学者。有C语言的课后练习题,也有稍微高于课后练习题难度的题目,属于入门题目。
recommend-type

鲲鹏练习题六套.docx

HCIA-Kunpeng Application Developer认证定位于培养能够使用鲲鹏计算平台,完成应用从X86向鲲鹏平台的迁移以及基于鲲鹏平台进行应用...本套资源(鲲鹏培训六套练习题,助你通过鲲鹏鲲鹏HCIA认证,与博客中的题型不重复)
recommend-type

vba编程练习题+近年真题.doc

VBA编程练习题 1、假定有以下循环结构 Do Until 条件 循环体 Loop 则正确的叙述是 A)如果“条件”值为0,则一次循环体也不执行 B)如果“条件”值为0,则至少执行一次循环体 C)如果“条件”值不为0,则至少执行一...
recommend-type

QSS样式表入门练习Demo.pdf

超详细的Qss入门练习Demo,包含了效果图和实现样式表; [可看博客](https://blog.csdn.net/qq_43627907/article/details/125677496)
recommend-type

zigbee-cluster-library-specification

最新的zigbee-cluster-library-specification说明文档。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

实现实时数据湖架构:Kafka与Hive集成

![实现实时数据湖架构:Kafka与Hive集成](https://img-blog.csdnimg.cn/img_convert/10eb2e6972b3b6086286fc64c0b3ee41.jpeg) # 1. 实时数据湖架构概述** 实时数据湖是一种现代数据管理架构,它允许企业以低延迟的方式收集、存储和处理大量数据。与传统数据仓库不同,实时数据湖不依赖于预先定义的模式,而是采用灵活的架构,可以处理各种数据类型和格式。这种架构为企业提供了以下优势: - **实时洞察:**实时数据湖允许企业访问最新的数据,从而做出更明智的决策。 - **数据民主化:**实时数据湖使各种利益相关者都可
recommend-type

2. 通过python绘制y=e-xsin(2πx)图像

可以使用matplotlib库来绘制这个函数的图像。以下是一段示例代码: ```python import numpy as np import matplotlib.pyplot as plt def func(x): return np.exp(-x) * np.sin(2 * np.pi * x) x = np.linspace(0, 5, 500) y = func(x) plt.plot(x, y) plt.xlabel('x') plt.ylabel('y') plt.title('y = e^{-x} sin(2πx)') plt.show() ``` 运行这段
recommend-type

JSBSim Reference Manual

JSBSim参考手册,其中包含JSBSim简介,JSBSim配置文件xml的编写语法,编程手册以及一些应用实例等。其中有部分内容还没有写完,估计有生之年很难看到完整版了,但是内容还是很有参考价值的。