作者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