load('shimi/bujiangweiduan')
load('shimi/bujiangweiduan1')
load('shimi/bujiangweishimi')
load('shimi/bujiangweishimi1')
load('allpositive')
xishu=positive4*shimi
allmo1=xishu*shimi'
allmo=positive4-allmo1
xishu1=positive4*shimi2
allmo2=xishu1*shimi2'
allmo2=positive4-allmo1
mo1=sqrt(sum(allmo.*allmo))
mo2=sqrt(sum(allmo2.*allmo2))
w=isequal(shimi,shimi2)