tag 标签: 喇叭图

相关博文
  • 热度 12
    2018-7-6 16:54
    6129 次阅读|
    0 个评论
    输液泵注射泵喇叭图的理解和生成
    以下是 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 = 1min T0 = 120min P 观测窗宽 min m = (T0 - P)/S + 1 观测窗数量 1 120 2 119 5 116 11 110 19 102 31 90 间隔 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 = ; % 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 = j:(j+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 流速下的喇叭图形曲线