如何实现数据光滑处理

dangjt

我想对下列数据进行光滑处理，以便找到其局部最大最小值，得到半周期（不定相等）
请问一般有什么方法，我用平均法试过，效果不是很好

[ 本帖最后由 dangjt 于 2007-3-1 16:23 编辑 ]

pjab001

用matlab进行样条插值即可

dangjt

插值怎么可能呢？插值不改变当前数据点的值大小，我要求出局部极大极小值，如上图所示，图上有２个局部极小，２个局部极大，但是数据会稍带毛刺，那些小的毛刺没有办法去掉，用普通算法会当作局部极值！

eight

原帖由 dangjt 于 2007-3-1 16:19 发表
插值怎么可能呢？插值不改变当前数据点的值大小，我要求出局部极大极小值，如上图所示，图上有２个局部极小，２个局部极大，但是数据会稍带毛刺，那些小的毛刺没有办法去掉，用普通算法会当作局部极值！

样条插值算很简单且可行的了，或者试试普通插值然后样条磨光

dangjt

原帖由 eight 于 2007-3-1 16:34 发表


样条插值算很简单且可行的了，或者试试普通插值然后样条磨光

呵呵，怎么个搞法，我上面有数据，教我一下，给个代码，用ｍａｔｌａｂ的就行！

风花雪月

1. /* 二次抛物线法绘制曲线函数 */
   
2. 
   
3. void paowuxian(int *x,int *y,int n,unsigned int k)
   
4. 
   
5. {
   
6. unsigned int i,j;
   
7. float t1,t2,t3,t,a,b,c,d,tx,ty;
   
8. *x=*(x+1);*(y)=*(y+1);
   
9. *(x+n+1)=*(x+n);*(y+n+1)=*(y+n);
   
10. t=0.5/k;
    
11. setcolor(10);
    
12. moveto(*(x+1),*(y+1));
    
13. for(i=0;i<n-1;i++)
    
14.   {
    
15.   for(j=1;j<k;j++)
    
16.     {
    
17.     t1=j*t;
    
18.     t2=t1*t1;
    
19.     t3=t2*t1;
    
20.     a=4*t2-t1-4*t3;
    
21.     b=1-10*t2+12*t3;
    
22.     c=t1+8*t2-12*t3;
    
23.     d=4*t3-2*t2;
    
24.     tx=a*(*(x+i))+b*(*(x+i+1))+c*(*(x+i+2))+d*(*(x+i+3));
    
25.     ty=a*(*(y+i))+b*(*(y+i+1))+c*(*(y+i+2))+d*(*(y+i+3));
    
26.     lineto(tx,ty);
    
27.     }
    
28.   }
    
29.   lineto(*(x+i+2),*(y+i+2));
    
30. }

复制代码

[ 本帖最后由 风花雪月 于 2007-3-4 04:23 编辑 ]

风花雪月

1. /* 二次B样条法绘制曲线函数 */
   
2. 
   
3. void byangtiao(int *x,int *y,int n,unsigned int k)
   
4. {
   
5. unsigned int i,j;
   
6. float t,t1,t2,a,b,c,tx,ty;
   
7. *x=*(x+1);*y=*(y+1);
   
8. *(x+n+1)=*(x+n);*(y+n+1)=*(y+n);
   
9. t=1.0/k;
   
10. setcolor(13);
    
11. moveto((*x+(*(x+1)))/2.0,(*y+(*(y+1)))/2.0);
    
12. for(i=0;i<n;i++)
    
13.   {
    
14.   for(j=1;j<k;j++)
    
15.     {
    
16.     t1=j*t;
    
17.     t2=t1*t1;
    
18.     a=(t2-2*t1+1)/2.0;
    
19.     b=t1-t2+1/2.0;
    
20.     c=t2/2.0;
    
21.     tx=a*(*(x+i))+b*(*(x+i+1))+c*(*(x+i+2));
    
22.     ty=a*(*(y+i))+b*(*(y+i+1))+c*(*(y+i+2));
    
23.     lineto(tx,ty);
    
24.     }
    
25.   }
    
26. }

复制代码

[ 本帖最后由 风花雪月 于 2007-3-4 04:23 编辑 ]

风花雪月

1. /*三次样条法*/
   
2. 
   
3. void parspl(int p[][2],int n,int k)
   
4. {
   
5. int i,j;
   
6. float t1,t2,t3,t,a,b,c,d,x,y;
   
7. p[0][0]=p[1][0];p[0][1]=p[1][1];
   
8. p[n+1][0]=p[n][0];p[n+1][1]=p[n][1];
   
9. t=0.5/k;
   
10. moveto(p[1][0],p[1][1]+D);
    
11. for(i=0;i<n-1;i++)
    
12.   {
    
13.   for(j=1;j<k;j++)
    
14.     {
    
15.     t1=j*t;
    
16.     t2=t1*t1;
    
17.     t3=t2*t1;
    
18.     a=4*t2-t1-4*t3;
    
19.     b=1-10*t2+12*t3;
    
20.     c=t1+8*t2-12*t3;
    
21.     d=4*t3-2*t2;
    
22.     x=a*p[i][0]+b*p[i+1][0]+c*p[i+2][0]+d*p[i+3][0];
    
23.     y=a*p[i][1]+b*p[i+1][1]+c*p[i+2][1]+d*p[i+3][1];
    
24.     lineto(x,y,15);
    
25.     }
    
26.   }
    
27.   lineto(p[i+2][0],p[i+2][1],15);
    
28. }

复制代码

风花雪月

试一下上面的几种方法

yimao

楼主要是解决了，能否共享一下啊？