标题: 血氧饱和度的c程序和matlab程序 [打印本页]

---

作者: 布安于室    时间: 2018-4-22 12:47
标题: 血氧饱和度的c程序和matlab程序
血氧饱和度的c程序和matlab程序

1. clear,clc;
   
2. t=0:1/500:5;
   
3. s_sin=sin(2*pi*1.5*t);
   
4. s_940nm=s_sin*0.2+0.75;
   
5. s_660nm=s_sin*0.1+0.8;
   
6. plot(t,s_sin);axis([0 5 0 1]);
   
7. figure,plot(t,s_940nm);axis([0 5 0 1]);
   
8. figure,plot(t,s_660nm);axis([0 5 0 1]);
   
9. 
   
10. max_total=0;
    
11. max_num=0;
    
12. min_total=0;
    
13. min_num=0;
    
14. b_value=0;
    
15. b_f=0;
    
16. n_value=0;
    
17. n_f=0;
    
18. 
    
19. for i=1:2500
    
20.     b_value=n_value;
    
21.     b_f=n_f;
    
22.     n_value=s_940nm(i);
    
23.     if n_value>b_value
    
24.         n_f=1;
    
25.     elseif n_value<b_value
    
26.         n_f=-1;
    
27.     end
    
28.    
    
29.     if(n_f==-1)&&(b_f==1)
    
30.         max_total=max_total+s_940nm(i);
    
31.         max_num=max_num+1;
    
32.     end
    
33.    
    
34.     if(n_f==1)&&(b_f==-1)
    
35.         min_total=min_total+s_940nm(i);
    
36.         min_num=min_num+1;
    
37.     end
    
38. end
    
39. 
    
40. max_940nm=max_total/max_num;
    
41. min_940nm=min_total/min_num;
    
42. 
    
43. max_total=0;
    
44. max_num=0;
    
45. min_total=0;
    
46. min_num=0;
    
47. b_value=0;
    
48. b_f=0;
    
49. n_value=0;
    
50. n_f=0;
    
51. 
    
52. first_p=0;
    
53. last_p=0;
    
54. f_p=0;
    
55. 
    
56. for i=1:2500
    
57.     b_value=n_value;
    
58.     b_f=n_f;
    
59.     n_value=s_660nm(i);
    
60.     if n_value>b_value
    
61.         n_f=1;
    
62.     elseif n_value<b_value
    
63.         n_f=-1;
    
64.     end
    
65.    
    
66.     if(n_f==-1)&&(b_f==1)
    
67.         max_total=max_total+s_660nm(i);
    
68.         max_num=max_num+1;
    
69.         if f_p==0
    
70.             f_p=1;
    
71.             first_p=i;
    
72.         else
    
73.             last_p=i;
    
74.         end
    
75.     end
    
76.    
    
77.     if(n_f==1)&&(b_f==-1)
    
78.         min_total=min_total+s_660nm(i);
    
79.         min_num=min_num+1;
    
80.     end
    
81. end
    
82. 
    
83. max_660nm=max_total/max_num;
    
84. min_660nm=min_total/min_num;
    
85. heart_beat=60/((last_p-first_p)/(max_num-1)/500);
    
86. 
    
87. spo2=110-25*(max_660nm-min_660nm)/max_660nm/((max_940nm-min_940nm)/max_940nm);
    
88. 
    

复制代码

SPO2（matlab）.doc

2018-4-22 12:47 上传

点击文件名下载附件

下载积分: 黑币 -5

18.5 KB, 下载次数: 22, 下载积分: 黑币 -5

spo2（c).doc

2018-4-22 12:47 上传

点击文件名下载附件

下载积分: 黑币 -5

23.5 KB, 下载次数: 27, 下载积分: 黑币 -5

---

作者: 无亟待    时间: 2019-8-6 10:59
这个可以研究一下

---

作者: 蒜毫炒肉    时间: 2020-7-9 16:22
老哥这也没个注释啊，萌新看不懂

---

作者: xmdhfyz    时间: 2020-7-13 19:00
这个可以研究一下

---

欢迎光临 (http://www.51hei.com/bbs/)

Powered by Discuz! X3.1