利用求均值方法实现4倍降采样的MATLAB代码

function [y] = mean_downsample(x,N)
% x:原始信号，格式为1*n的行向量
% y：平均降采样后的信号
% N：降采样的倍数，即降采样后频率=降采样前频率/N
% 作者：张亦弛 2023/5/31
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
downdata = zeros(1,size(x,2)/N); % 建立空数组存放降采样后的和
for phase = 0:N-1 % 偏置从0取到N-1
temp = downsample(x,N,phase); % N倍降采样，偏移phase个点进行
downdata = downdata + temp; % 累加N次降采样后的结果
end
y = downdata/N; % 平均累加结果，得到平均降采样的值
end