function I=vectorsmooth(Y,n);                    %實現對向量的平滑操作
%傳遞兩個參數,一個是待平滑的向量,一個是平滑窗口的半徑.
if n==0
    I=Y;
    return;
end
[h,w]=size(Y);
if (h<w)                               %如果是行向量先轉置為列向量
    Y=Y';
end
nCount=max(size(Y));
tmpY=zeros(nCount + 2*n,1);
tmpY((n+1):(n+nCount),1)=Y;
tY=zeros(nCount,1);
for i=1:nCount
    tY(i,1)=tmpY(n+i,1);
    for j=1:n
        tY(i,1)=tY(i,1) + tmpY(n+i-j,1) + tmpY(n+i+j,1);
    end
end
for i=1:nCount
    if (i< n+1) & ( i>nCount-n)
        tY(i,1)=tY(i,1)/(2*n+2-i);
    else
        tY(i,1)=tY(i,1)/(2*n+1);
    end    
end
if (h<w)
    tY=tY';
end
I=tY;