1
function [res]=interpolate(f, angle_step, oversample_multiply,cutof_freq)
3
full = n * oversample_multiply;
6
[a,b]=meshgrid(-full2:(full2-1),-full2:(full2-1));
7
rad_step=degtorad(angle_step);
11
c2(c1<0)=pi + c1(c1<0);
13
c=min(1 + (c2/rad_step), size(f,2));
15
r = (n/2) + min((sqrt((a.*a)+(b.*b))/(oversample_multiply)),n/2);
18
r2(c1<0)=(n+1)-r(c1<0);
20
res1=interp2(f,c,r2,'nearest');
22
res((r-(n/2))>=cutof_freq*n/2)=0;