如何在Matlab中找到interp1的正确采样率

时间:2015-06-05 19:26:06

标签: matlab resampling downsampling

我试图找到interp1(或抽取或重新采样)的正确采样率,用于绘制由imfreehand获得的轴中的多个点。

首先我绘制一些东西(任何不正常的东西)

bands=1:20;plot(bands,sin(bands));
hax=gca;

enter image description here

然后在轴上用imfreehand绘制一些东西:

imfreehand('Closed',0);

enter image description here

我从imfreehand对象中提取数据部分(size 1x413),然后我尝试用原始行的长度绘制它们但是缩放是错误的...

data=get(hfree); 
xydata=get(data.Children(4)); 
% x=xydata.XData;
y=xydata.YData; 

len=length(y);
x=1:len;
newlen=length(bands);
scale=(len-1)/(newlen-1);
xx=1:scale:len;

yy=interp1(x,y,xx,'spline');

line(1:length(yy),yy(1:length(yy)),...
   'LineWidth',4,...
   'Color',[.8 .8 .8],...
   'Parent',hax);

enter image description here

任何提示??

此外,与建议的下采样相同的结果..

for i=4:1:15 %just testing numbers 

    yy = downsample(y,i);

    h=line(1:length(yy),yy(1:length(yy)),...
        'LineWidth',1.2,...
        'Color',[.8-i/100 .8-i/100 .8-i/100],...
        'Parent',hax);

    if length(bands)==length(yy)
        set(h,'Color',[0 0 0], 'LineWidth',2);
    end

end

enter image description here

任何提示?? x2

1 个答案:

答案 0 :(得分:0)

我很难忘记在interp中正确缩放

len=length(y);
newlen=length(bands);
scale=ceil(len/newlen)    

x1=x(1:scale:end)
y1=y(1:scale:end);
h=line(x1,y1,...
        'LineWidth',1.2,...
        'Color',[0 0 0],...
        'Parent',hax);

yyy=interp1(x1,y1,bands,'pchip',nan)

h=line(1:20,yyy,...
        'LineWidth',1.2,...
        'Color',[0 1 1],...
        'Parent',hax);

enter image description here