极值点插值遇到的问题

yrj01

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

zhouyang664

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

yrj01

回复 2 # zhouyang664 的帖子

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

ChaChing

...但若两极值点间隔较大时就会出现误差，若人为忽略被误插进的点，直接对间隔较大点插值，就会有“There should be at least two data points."的错误提示...

求助完整格式：出错代码和出错提示
看不清楚LZ的意思, 可否举个例子!?:@)

yrj01

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

ChaChing

看帖找空试, 回头才发现LZ已经又编辑了!
1.LZ给的部分码无法直接试, 有点不方便, 只得随意產生信号再试

1. clc; clear
   
2. y=5*rand(1,1605); y0=y(596:1605); time=[0:141605]/1000; t=time(140596:141605);
   
3. d = diff(y0); n = length(d); d1 = d(1:n-1); d2 = d(2:n);
   
4. indmin = find(d1.*d2<0 & d1<0&d2>1.5)+1;
   
5. indmax = find(d1.*d2<0 & d1>1.5&d1<1.5)+1;
   
6. envmin = pchip(t(indmin),y0(indmin),t);
   
7. envmax = pchip(t(indmax),y0(indmax),t);
   
8. 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

回复 6 # ChaChing 的帖子

多谢了，这些我也找到了，只是我有不需要插进去的极值点，还得在插值前对所求极值设定些条件。