作者captainVince (purple rain N )
看板comm_and_RF
標題[問題] Lagrange Interpolation的Freq spectrum怎麼畫?
時間Tue May 8 17:46:00 2007
請問一下喔 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