主题:[讨论]求助一程序 数据插值与曲线拟合问题 求解
达布流-
[专家分:0] 发布于 2011-06-13 00:33:00
数据插值与曲线拟合问题
表中列出压缩弹簧的一组数据:
压力P(KN) 10 20 30 40 50 60 70 80 90 100
变形δ(mm) 6 14 26 33 46 54 67 75 84 100
要求:
(1)从文件中读入数据。
(2)用数据插值方法计算P=15,25,35,45,55,65,75,85,95时(单位KN)弹簧的变形。
(3)求数据组的拟合曲线。
(4)绘制压力-变形拟合曲线。
回复列表 (共6个回复)
沙发
cgl_lgs [专家分:21040] 发布于 2011-06-13 08:39:00
如果不考虑塑性的话,弹簧的形变与压力应成正比。
板凳
bshine1225 [专家分:720] 发布于 2011-06-13 08:53:00
关于线性拟合的问题徐士良的《fortran常用算法程序集》里面有现成的子程序可以拿来用。
这本书我刚把上传到了网盘:http://u.115.com/file/bhitbmhw
3 楼
cgl_lgs [专家分:21040] 发布于 2011-06-13 12:09:00
如果楼主用IMSL的话,里面也有现成的插值函数可以用,具体例子可以在论坛里找,俺写过:)
4 楼
bshine1225 [专家分:720] 发布于 2011-06-13 12:30:00
看来IMSL确实是个好工具
5 楼
cgl_lgs [专家分:21040] 发布于 2011-06-14 01:26:00
[quote]看来IMSL确实是个好工具[/quote]嘿嘿IMSL里的SFUN好东东不少。矩阵操作也还凑合,只是解方程太差了点儿:)
6 楼
dongqiao [专家分:0] 发布于 2011-06-14 09:40:00
建议你用ORIGIN 软件,很好用 直接插值 和给出曲线拟合 你这个数据采用下述公式拟合,相关系数平方在0.99以上,如下:
Data: Data1_B
Model: Cubic
Equation:
y = A + B*x + C*x^2 + D*x^3
Weighting:
y No weighting
Chi^2/DoF R^2
----------------------------------------
3.81488 0.99737
----------------------------------------
Parameter Value Error
----------------------------------------
A -5 3.7907
B 1.04794 0.28413
C -0.00213 0.00586
D 0.00002 0.00004
我来回复