yrj01 发表于 2011-6-2 07:53

极值点插值遇到的问题

先找到极大值点后,欲将所有极大值点插值拟合成一光滑曲线,发现两个间隔较大的极大值点间插进了其它点(极小值点),为什么?如果两个点距离较远,这个间隙是怎么被补足而生成曲线的?

zhouyang664 发表于 2011-6-2 21:45

这个论坛上都有的,你搜索一下就行了,我感觉你的意思好像就是画包络线吧!

yrj01 发表于 2011-6-2 22:27

回复 2 # zhouyang664 的帖子

嗯,是画极值点的包络线,但若两极值点间隔较大时就会出现误差,若人为忽略被误插进的点,直接对间隔较大点插值,就会有“There should be at least two data points."的错误提示,该怎么对间隙补足?

ChaChing 发表于 2011-6-3 00:07

本帖最后由 ChaChing 于 2011-6-3 00:08 编辑

...但若两极值点间隔较大时就会出现误差,若人为忽略被误插进的点,直接对间隔较大点插值,就会有“There should be at least two data points."的错误提示...求助完整格式:出错代码和出错提示
看不清楚LZ的意思, 可否举个例子!?:@)

yrj01 发表于 2011-6-3 08:33

本帖最后由 yrj01 于 2011-6-3 09:26 编辑

谢谢关注,问题暂时解决。

ChaChing 发表于 2011-6-3 10:26

看帖找空试, 回头才发现LZ已经又编辑了!
1.LZ给的部分码无法直接试, 有点不方便, 只得随意產生信号再试 clc; clear
y=5*rand(1,1605); y0=y(596:1605); time=/1000; t=time(140596:141605);
d = diff(y0); n = length(d); d1 = d(1:n-1); d2 = d(2:n);
indmin = find(d1.*d2<0 & d1<0&d2>1.5)+1;
indmax = find(d1.*d2<0 & d1>1.5&d1<1.5)+1;
envmin = pchip(t(indmin),y0(indmin),t);
envmax = pchip(t(indmax),y0(indmax),t);
plot(t,envmax,'m',t,envmin,'r')
2."There should be at least two data points."是针对pchip输入空矩阵, 往前追是indmax空矩阵造成
3.indmax = find(d1.*d2<0 & d1>1.5&d1<1.5)+1, 裡头d1>1.5&d1<1.5怪怪的
4.参考下
一种有效的包络线算法 http://forum.vibunion.com/thread-32967-1-1.html
   求画出振荡的包络线图 http://forum.vibunion.com/thread-28736-1-1.html
   from http://home.vibunion.com/blog-63979-18250.html

yrj01 发表于 2011-6-3 15:03

回复 6 # ChaChing 的帖子

多谢了,这些我也找到了,只是我有不需要插进去的极值点,还得在插值前对所求极值设定些条件。
页: [1]
查看完整版本: 极值点插值遇到的问题