LPF
一阶低通滤波
同样拿那个炸机气压计日志分析,可以更好对比
- 新建变量,把高度存进去BaroAlt.mat
- 写获取高度值函数
- 写LPF函数
function output = LPF(x)
persistent lastx Init
if isempty(Init)
lastx = x;
Init = 1;
end
alpha = 0.9;
output = alpha*lastx + (1 - alpha)*x;
lastx = output;
- 写个调用脚本RunLPF.m
clear all
Nsamples = 6800;
Xsaved = zeros(Nsamples, 1);
Xmsaved = zeros(Nsamples, 1);
for k=1:Nsamples
xm = GetAlt();
x = LPF(xm);
Xsaved(k) = x;
Xmsaved(k) = xm;
end
dt = 0.02;
t = 0:dt:Nsamples*dt-dt;
figure
hold on
plot(t, Xmsaved, 'r.');
plot(t, Xsaved, 'b');
legend('Measured', 'LPF')
- 实验结果如下: