原创 输液泵注射泵喇叭图的理解和生成

以下是9706.27上的定义：

流速 r 由操作者选择的输液速度
流量 flow 单位时间内测得的输出容量
丸剂 bolus 作为一次注人但不属于灌注程序的一部分，在短时间内传输的液体的离散量。
取样间隔 5 连续的质量读数或液滴计数之间的时间
试验周期 T 试验开始到结束的总的持续时间
分析周期 T o 试验周期的第一个 2 h
分析周期 T 试验周期的第 2 h
分析周期 T 。 试验周期的最后 1 h
分析周期 T * 作为 Ta, T ， 或 几 的分析周期
W 总 质 量
W 规定的分析周期中第2 次取样质量
W ， 规定的分析周期或试验周期结束时的取样质量
W ， 规定的分析周期开始时的取样质量
A 分析周期 T : 中测得的总的平均百分比流量误差
B 分析周期 T 中测得的总的平均百分比流量误差
尸 观 测 窗 期 间
E , (m ax) 规定期间内在观测窗中测得的最大测量误差

分析周期是人为定义的，前120分钟为T0。第1小时到2小时期间为T1 最后一小时为T2

S为采样间隔，不同的测试在标准上有明确的规定，对于计滴数的，采样间隔为1min.

P 观测窗 是一个固定的时间，用这个固定的时间在采样间隔上滑动。

观测窗的数量m = (Tx - P)/S + 1  标准中的最大观测窗数量容易引起误解，实际上对于固定的Tx , P和S值，m也是确定的。

比如T0 = 120 S = 1 如果P = 1 观测窗为 120个。

比如T0 = 120 S = 1 如果P = 2 观测窗为 119个。

比如T0 = 120 S = 1 如果P = 31 观测窗为 90个。

计算Ep

这个标准上给的两个公式，实际上是一个。

S 采样间隔是固定的，对于滴速控制器是 1min，其它泵类要求会不一样，有0.5min的。

P 观测窗分别为 1m in,2m in,5 min,11m in,19m in和31m in时 对应有不同的m值。如下表：

间隔S  这里是1分钟有一个滴速误差Qi

对于P = 1 Ep 就是找各个Qi的最大值和最小值。

对于P = 2 Ep 就是移动找相邻2个Qi的平均值的最大值和最小值。

(Q1+Q2)/2  (Q2+Q3)/2...... (Q119+Q120)/2

对于P = 5 Ep 就是移动找相邻5个Qi的平均值的最大值和最小值。

对于P = 31 Ep 就是移动找相邻31个Qi的平均值的最大值和最小值。

以下是matlab 做的摸拟

clear;

T0 = 120; % totle test time   120min

S  = 1;   % measture interval 1min   

P = [1 2 5 11 19 31]; % obsive window

r = 20;   % wanted rate 20 drop/min

a = 0:r/S:T0*r;       %这个是模拟数据  

n = random('norm',1,0.1,1,length(a));  % generator some errors

a = a + n;  %这个是模拟数据加上误差。可以用实际测到的数据代入，得到真的结果

Q = a(2:length(a))-a(1:length(a)-1);

%加入误差 刚开始的几个数据不能太好看了，不然就假了。

Q(1) = Q(1)-r/15;

Q(2) = Q(2)+r/20;

Q(3) = Q(3)-r/25;

B = r* ones(1,length(Q));

x = 1:length(Q);

plot(x,Q,'b',x,B,'r--');

xlabel('开始后的时间/min');

ylabel('滴速');

title('设定滴速 = 20滴/min');

for i = 1:length(P)      % 分别求每个观测窗的误差最大值与最小值

    m = (T0 - P(i))/S + 1;

   for j = 1:m

       temp = 0;

       for k = jj+P(i)/S-1)

           temp = temp + 100* (Q(k)-r)/r; % 计算误差和

       end

       Ep(j) =  temp * S / P(i);    % 求平均

   end

   Epmax(i) = max(Ep);

   Epmin(i) = min(Ep);

   clear Ep

end

A = 100 * ones(1,length(P))*(a(length(a))/T0 - r);

plot(P,Epmax,'bx-',P,Epmin,'rx-',P,A,'g')

ylabel('百分比流量误差 %');

xlabel('观测窗 /min');

legend('Ep(max)','Ep(min)','总的百分比误差');

这是运行上面这段程序，得到的一个运行结果

20滴

20滴/min 流速下的喇叭图形曲线