请教大家一个插值求离散数据极值的问题

edgarzw 2010-08-27 04:11:34
现有一组离散数据,我们暂缺把它命名为数组y[20],其数据大概是一个有规律的凸形(可以说是一条特性曲线吧),而与其相对应的为数组x[20],其20个数分别为等间距的从1至20的整数,现想要求这数组的极值(非这群数中的最大值,而是这个函数的最大值)及其所对应的x值,之前所用的是对最大的三个数采用高斯插值的办法,但效果不佳,请问还有其他的插值方法来解决吗?谢谢!
...全文
247 10 点赞 打赏 收藏 举报
写回复
10 条回复
切换为时间正序
当前发帖距今超过3年,不再开放新的回复
发表回复
sxdkxgwan 2010-08-31
我是来学习的。。。
  • 打赏
  • 举报
回复
小明哥-Jimmy 2010-08-30
[Quote=引用 5 楼 the_eagles 的回复:]
引用楼主 edgarzw 的回复:
现有一组离散数据,我们暂缺把它命名为数组y[20],其数据大概是一个有规律的凸形(可以说是一条特性曲线吧),而与其相对应的为数组x[20],其20个数分别为等间距的从1至20的整数,现想要求这数组的极值(非这群数中的最大值,而是这个函数的最大值)及其所对应的x值,之前所用的是对最大的三个数采用高斯插值的办法,但效果不佳,请问还有其他的插值方法来解决吗?谢谢!
……
[/Quote]

神仙?
达人?
谢谢!
  • 打赏
  • 举报
回复
AlanBruce 2010-08-30
kankan...
  • 打赏
  • 举报
回复
The_eagles 2010-08-28
[Quote=引用楼主 edgarzw 的回复:]
现有一组离散数据,我们暂缺把它命名为数组y[20],其数据大概是一个有规律的凸形(可以说是一条特性曲线吧),而与其相对应的为数组x[20],其20个数分别为等间距的从1至20的整数,现想要求这数组的极值(非这群数中的最大值,而是这个函数的最大值)及其所对应的x值,之前所用的是对最大的三个数采用高斯插值的办法,但效果不佳,请问还有其他的插值方法来解决吗?谢谢!
[/Quote]
对最大的三个数采用高斯插值??
~~~~~~~~~~~~~~~~~~~~~~~~~~
1. indexs of these there value must be continous.
2. guass is absolutely enough, even rectangle.

//pseudo code
for(nwnd=5; nwnd<11; ++nwnd, ++nwnd)
{
nNewLen= 20 + nwnd-1;
ye[nNewLen]={0};
x[nNewLen]={0};
memcpy(ye+nwnd/2, y,20)
for(i=nwnd/2; i<nNewLen-nwnd/2; ++i)
{
j=-nwnd/2;
while(j++ <= nwnd/2)
x[i] += ye[i-j];
x[i] /= nwnd;
}

search the biggest value in x[], get index.
This "effect" is you desired?
yes ---> break;
no ---> continue;
}
  • 打赏
  • 举报
回复
DavidFracs 2010-08-28
[Quote=引用 4 楼 tu_lq 的回复:]
最小二乘没法“凸”没法极值
引用 3 楼 davidfracs 的回复:

这种情况用最小二乘要比插值好吧
[/Quote]

用最小二乘算出来一条二次曲线,有了这条曲线不就可以求极值了么?
  • 打赏
  • 举报
回复
toplinq 2010-08-27
最小二乘没法“凸”没法极值
[Quote=引用 3 楼 davidfracs 的回复:]

这种情况用最小二乘要比插值好吧
[/Quote]
  • 打赏
  • 举报
回复
DavidFracs 2010-08-27
这种情况用最小二乘要比插值好吧
  • 打赏
  • 举报
回复
toplinq 2010-08-27
其实就是多项式插值和三次样条插值这两种比较合适,至于你是拉格朗日那么插还是牛顿那么插,都一样,把所有的点都用上即可,程序网上都有现成的。
  • 打赏
  • 举报
回复
koutyayh 2010-08-27
插值的办法很多:
lagrange插值、newton插值、hermite插值、三次样条插值……自己选择比较吧,详细的可以看数值方法类的书
  • 打赏
  • 举报
回复
相关推荐
发帖
C++ 语言
创建于2007-09-28

6.0w+

社区成员

C++ 语言相关问题讨论,技术干货分享,前沿动态等
申请成为版主
帖子事件
创建了帖子
2010-08-27 04:11
社区公告
暂无公告