[問題] Lagrange Interpolation的Freq spectrum怎麼畫?

[captainVince (purple rain￾ ￾N￾ )]

請問一下喔 lagrange interpolation要怎麼畫出他在

frequency domain上的response勒?

我只能畫出整數倍的interpolation impulse response 像是3倍 4倍 10倍

然後再轉成Freq domain畫出spectrum

不知道該怎麼畫出友小數點的倍數 像是1.5 5.2 7.4這種的 想請問一下

該怎麼畫呢?
下面是我的Matlab code

%這個function會產生係數
function [C,B]=Lagrange(N,mu)
C=ones(1,N);
D=ones(1,N);
for m=1:N,
    for n=1:N,
        if n~=m,
            C(m)=C(m)*(N/2+mu-n)/(m-n);
            A(m,n)=(-N/2+n);
            D(m)=D(m)*(m-n);
        end;
    end;
end;


%plot freq
clear;
order=6;
ratio=10;
k=1;
u=0;
C=[];


for u=0:1/ratio:1
        [C1(k,:)]=Lagrange(order,u);
        k=k+1;
end
C1=transpose(C1);
[m,n]=size(C1);
C1=C1(: ,1: n-1 );
for l=1:m
        C=[ C1(l,:) C ];
end
C=C/sum(C);
fvtool(C)

有人知道該怎麼做嗎?


--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 60.248.150.178