插值后新序列与原函数的幅值和相位误差

我是了了

我原函数是：sin（100*pi*t+q），q是从1度到15度（步长为1度）可变的初相角，对原函数采样48点，进行二次插值，得到一个新序列，首先计算原函数和二次插值后的新序列在每个小区间（就是48点，47个小区间上）可能存在的最大误差（这个地方，我是在每个小区间10等分，计算这10等分点上的误差值,这样计算的结果基本和理论推导一致，如果认为不可靠，可以增加等分点数计算），结果是正确的。现在，我还需要求第一个小区间上的新序列的幅值和相位的误差！这个地方，原来论文上说，傅里叶变换后，在不同的15个初相角的情况下，第一个小区间15等分，得到225个算例，我不明白了！怎么会有225个算例？一个小区间不是只对应一个50HZ的主频么？15*15当然是225.但是对应100*pi，实际只有一个50HZ的主频，只能有一个相角的误差，不明白了，请教高人。还有就是论文所说的新序列的幅值和相位误差是如何求解？下面是我求采样值误差的matlab仿真程序，后边的怎么求解，请高人指点思路，或者帮我写一下程序。非常感谢！我这里有各学科的论文，可以相互交流，探讨！

clc;clear;
N=48;

for k=1:15  %0--15度
    for n=1:48   
        t0=n/48*0.02; t1=(n+1)/48*0.02;t2=(n+2)/48*0.02;
%         t=(n+1)/48*0.02+0.02/N/sqrt(3);
        h=(2/48*0.02)/10;
        a=0;
        for t=t0:h:t2
            II(n+(k-1)*N)=sin(100*pi*t+(k-1)/180*pi);   原函数
            III(n+(k-1)*N)=(t-t1)*(t-t2)/(t0-t1)/(t0-t2)*sin(100*pi*t0+(k-1)/180*pi)+(t-t0)*(t-t2)/(t1-t0)/(t1-t2)*sin(100*pi*t1+(k-1)/180*pi)+(t-t0)*(t-t1)/(t2-t0)/(t2-t1)*sin(100*pi*t2+(k-1)/180*pi);二次插值函数
            Rmax(n+(k-1)*N)=abs(II(n+(k-1)*N)-III(n+(k-1)*N));  求解两者的误差   
            if Rmax(n+(k-1)*N)>a
                a=Rmax(n+(k-1)*N);
            end            
        end
        Rmax(n+(k-1)*N)=a;
    end
end

plot(Rmax(1:250))

我是了了

悲剧，被百度搜索都到了前十几页了，但是没人能解答！

我是了了

有人懂这一块么？给点提示！谢谢！