医学信号数字处理技术作业.docxVIP

作业 1：一、获取生理信号1）从 MIT 生物医学工程数据库 /获取生理信号并保存为文本文件。获取生理信号的步骤登录数据库网站如图鼠标置于PHYSIONET，即出现下拉菜单如图，鼠标置于PHYSIOBANK，出现下拉菜单如图单击PhysioBank ATM显示如下选择Database，Record等如下图将右上方Toolbox选为show sampples as text 如下图鼠标置于samples.txt单击鼠标右键，在出现的下拉菜单中单击目标另存为，将数据文件存为文本文件备用如图显示生理信号波形的所有程序语句及运行结果：运行结果：Person_01/rec_2_ECG Ir10.edf_ECGs001/s001R02.edf_Fc5slp14_ECGslp14_BP程序clear allload Person_01rec_2load r10_edfload s001_s001R02_edfload slp14%%Person_01rec_2%%t = Person01rec2(:,1);%??data?Dμ?μúò?áDêy?Y?3??tecg1 = Person01rec2(:,2); %??data?Dμ?μú?táDêy?Y?3??x1iltered = Person01rec2(:,3); %??data?Dμ?μúèyáDêy?Y?3??x2figure; plot(t,ecg1)%ò?t?aoá×?±ê￡?x1?a×Y×?±ê×÷í?figure; plot(t,iltered)%%r10_edf%%t = r10edf(:,1);x1 = r10edf(:,2); %??data?Dμ?μú?táDêy?Y?3??x1x2 = r10edf(:,3); %??data?Dμ?μúèyáDêy?Y?3??x2x3 = r10edf(:,4);x4 = r10edf(:,5);x5 = r10edf(:,6);x6 = r10edf(:,7);figure; plot(t,x1)%ò?t?aoá×?±ê￡?x1?a×Y×?±ê×÷í?figure; plot(t,x2)figure; plot(t,x3)figure; plot(t,x4)figure;plot(t,x5)figure; plot(t,x6)%%s001_s001R02_edf%%t = s001s001R02edf(:,1);Fc5 = s001s001R02edf(:,2); %??data?Dμ?μú?táDêy?Y?3??x1figure; plot(t,Fc5)%ò?t?aoá×?±ê￡?x1?a×Y×?±ê×÷í?%%slp14%%t = slp14(:,1);ecg = slp14(:,2); %??data?Dμ?μú?táDêy?Y?3??x1bp = slp14(:,3); %??data?Dμ?μúèyáDêy?Y?3??x2C3 = slp14(:,4);nasal = slp14(:,5);figure; plot(t,ecg)%ò?t?aoá×?±ê￡?x1?a×Y×?±ê×÷í?figure; plot(t,bp)figure;plot(t,C3)figure;plot(t,nasal)2）你认为有哪几种方法可以将 physionet 中的生理信号转换为 matlab 可调用的数据。请写出具体实现的步骤。答：三种可以在command window中输入load test.txt ，然后就会产生一个test的数据文件，内容跟test.txt中的数据一样；另一种方法是在file/import data....../next/finish 也可产生一个叫test的数据文件。[a1,a2,a3,a4]=textread(test1.txt,%s%s%s%s,headerlines,4)说明：%s可以是其他形式，跟读入的数据类型有关，比如这里也可以用%n,%f等。这里%s的个数和[a1,a2,a3,a4]对应。fidin=fopen(test.txt);?????????????????????????????? % 打开test2.txt文件fidout=fopen(mkmatlab.txt,w);?????????????????????? % 创建MKMATLAB.txt文件while ~feof(fidin)????????????????????????????????????? % 判断是否为文件末尾tline=fgetl(fidin);???????????????????????????????? % 从文件读行if double(tline(1))=48double(tline(1))=57??????