x?=?xlsread（‘data1.xlsx‘‘B6?:?D8‘）;
P?=?xlsread（‘data1.xlsx‘‘C3?:?D3‘）;
P?=?P‘;
%?idx=?x==0;?
%?x（idx）?=?inf;?
idx?=?find（isnan（x））;?
x（idx）?=?inf;?
b?=?1?./?x;
[rc]?=?size（b）;
B?=?zeros（r-1）;
for?i?=?1?:?（r-1）
????for?j?=?1?:?（r-1）
????????if?i?==j
????????????bi?=?b（i+1:）;
????????????B（ij）?=?sum（bi）;??
????????else
????????????B（ij）?=?-b（i+1j+1）;
????????end
????end
end