脉搏波的数学分解

clc
clear
close
% 定义时间轴
t = 0:0.01:2*pi;

% 前向压力波
a1 = 12;  % 振幅
mu1 = 1.0*pi/3;  % 水平位置
sigma1 =0.35;  % 波宽
forward_pressure_wave = a1 * exp(-(t - mu1-sigma1).^2 / (2 * sigma1^2));

% 反射波形
a2 = 5;  % 振幅
mu2 = 1.7*pi/2;  % 水平位置
sigma2 = 0.5;  % 波宽
reflection_wave = a2 * exp(-(t - mu2).^2 / (2 * sigma2^2));

% 重搏波形
a3 = 3;  % 振幅
mu3 = 5.5 * pi / 4; % 水平位置
sigma3 = 0.55;   % 波宽
rebounding_wave = a3 * exp(-(t - mu3).^2 / (2 * sigma3^2));

% 合成桡动脉脉搏波（由前三个波形合成的结果）
combined_wave = forward_pressure_wave + reflection_wave + rebounding_wave;

% 画图
figure;
hold on;
plot(t, forward_pressure_wave, '--', 'LineWidth', 1.5); 
plot(t, reflection_wave, '--', 'LineWidth', 1.5); 
plot(t, rebounding_wave, '--', 'LineWidth', 1.5); 
plot(t, combined_wave, 'k-', 'LineWidth', 2); 
hold off;

title('波形合成模型');
legend('前向压力波', '反射波', '重搏波', '桡动脉脉搏波');
xlabel('时间');
ylabel('幅度');
grid on;