已知加速度信号，如何基于高通滤波和数值积分的方法求位移和速度，给出具体的matlab代码，速出结果用图像表示

以下是基于高通滤波和数值积分的方法求位移和速度的MATLAB代码示例，并使用图像表示结果：

% 加载加速度信号数据
load('acceleration_signal.mat');

% 定义采样频率和采样时间间隔
fs = 100;  % 采样频率为100Hz
dt = 1/fs; % 采样时间间隔为0.01s

% 高通滤波器参数设置
fc = 0.5; % 截止频率为0.5Hz
[b, a] = butter(1, fc*dt*2, 'high'); % 一阶Butterworth高通滤波器

% 高通滤波
acceleration_filtered = filtfilt(b, a, acceleration_signal);

% 数值积分得到速度
velocity_signal = cumtrapz(acceleration_filtered) * dt;

% 数值积分得到位移
displacement_signal = cumtrapz(velocity_signal) * dt;

% 绘制加速度、速度和位移的时域图像
t = 0:dt:(length(acceleration_signal)-1)*dt;
figure;
subplot(3,1,1);
plot(t, acceleration_signal);
title('Acceleration Signal');
xlabel('Time (s)');
ylabel('Acceleration (m/s^2)');

subplot(3,1,2);
plot(t, velocity_signal);
title('Velocity Signal');
xlabel('Time (s)');
ylabel('Velocity (m/s)');

subplot(3,1,3);
plot(t, displacement_signal);
title('Displacement Signal');
xlabel('Time (s)');
ylabel('Displacement (m)');

运行以上代码，可以得到加速度、速度和位移的时域图像，如下图所示：