sysrep.m

来自「里面囊括了基于matlab滤波器设计的各种.m文件」· M 代码 · 共 20 行

% y = sysrep(h0,h1,g0,g1,M,N,x)
% System input and output are x, y respectively.
% The output y is a column vector. M and N are
% decimation, interpolation factors respectively.
% h0, h1 and v0, v1 are system impulse response.

function y = sysrep(h0,h1,g0,g1,M,N,x);

v0 = filter(h0,1,x);
v1 = filter(h1,1,x);
v0 = dsample(v0,M);
v1 = dsample(v1,M);
v0 = usample(v0,N);
v1 = usample(v1,N);
v0 = filter(g0,1,v0);
v1 = filter(g1,1,v1);
v0 = v0(:);
v1 = v1(:);
y = v0 + v1;
y = y(:);