interps.m

来自「B3g_phase2_C语言_Matlab程序及说明」· M 代码 · 共 34 行

function y = InterPS(x,h)

% Lx=length(x); Lh=length(h);Lz=(Lx-1)*3+Lh;
% z=zeros(1,Lz);
% l=1;
% for n=1:Lx
%     z(l:l+Lh-1)=z(l:l+Lh-1)+x(n)*h;
%     l=l+N;
% end
% 
% y=z((Lh+1)/2:Lz-(Lh-1)/2);


h1=h(3:3:end);
h2=h(2:3:end);

Lx=length(x); Lh=length(h1);Ly=Lx-Lh+1;

for n=0:Ly-1
    y(3*n+1)=x(n+6);
    y(3*n+2)=x(n+1:n+12)*h1';
    y(3*n+3)=x(n+1:n+12)*h2';
end

% a=y(16:45)-yy
% 
% plot(real(y))
% hold on 
% plot(real(yy))
% 
% hold off
% pause