MATLAB实现数据延拓

一、前项与后向预测

1.函数调用：

2.范例：

滤波后的数据一般会在输出信号的初始端有瞬态效应，利用对信号的两端延拓，消除瞬态效应，此处使用的是零相位滤波器。

clear all; clc; close all
N = 1000;                    % 数据长度
Fs = 1000;                   % 采样频率
t = (0:N-1)/Fs;              % 时间刻度
% 滤波器设计
fp=[3 15];                   % 滤波器通带阻带参数设定
fs = [0.5 30];
rp = 1.5;                    % 通带波纹
rs = 20;                     % 阻带衰减
wp = fp*2/Fs;                % 归一化频率
ws = fs*2/Fs;
[n,wn]=buttord(wp,ws,rp,rs); % 计算滤波器阶数
[b,a] = butter(n,wn);        % 计算滤波器系数
[h,w] = freqz(b,a,1000,Fs);  % 求滤波器响应
h = 20*log10(abs(h));        % 计算滤波器幅值响应 
%信号的产生
f1 = 0.001;                  % 分量1,准直流
f2 = 5;                      % 分量2,有用信号
f3 = 50;                     % 分量3,工频干扰
x1 = 100+10*sin(2*pi*f1*t);  % 产生3个分量的信号
x2 = 10*sin(2*pi*f2*t);
x3 = 10*sin(2*pi*f3*t);
xn = x1+x2+x3;               % 合并为信号xn
y1