%打开反馈仪数据
[fname1,pathname1]=uigetfile('*.dat','"C:\Documents and Settings\wff\桌面\反馈仪临床\新data','*.dat');
clear LL1
clear LL2
FIDS1=fopen([pathname1 fname1],'r');
LL1=fread(FIDS1,20001,'uint8');
LL2=fread(FIDS1,20001,'uint8');
fclose(FIDS1);
num=0;
for i=1:10001
if LL1(i)>0 & LL1(i)<17 & LL2(i)>0 & LL2(i)<17
num=num+1 ;
Light1(num)=LL1(i);
Light2(num)=LL2(i);
end
end
cha=Light1-Light2;
subplot(2,1,1),plot((1:num)/50,Light1(1:num),'k')
AXIS([0 30 1 16])
TITLE([fname1 ' 同类上市肌电生物反馈仪'])
XLABEL('时间/分钟')
YLABEL('反馈灯显示序号')
subplot(2,1,2),plot((1:num)/50,Light2(1:num),'k')
AXIS([0 30 1 16])
%subplot(3,1,3),plot((1:num)/50,cha(1:num),'k')
%AXIS([0 30 -5 5])
crr= corr2(Light1(1:num),Light2(1:num))
TITLE([fname1 ' 本肌电生物反馈仪'])
XLABEL('时间/分钟')
YLABEL('反馈灯显示序号')
set(1,'Color','white')
%subplot(3,1,3),plot(1:num,Light2(1:num)-Light1(1:num),'k')
m1=mean(Light1(1:num))
m2=mean(Light2(1:num))
std1=std(Light1(1:num))
std2=std(Light2(1:num))
评论0