RCC.zip_matlab 核反应_互功率谱密度_材料属性反演_核 谱_核反应

function cc=Rcc(data1,data2,N) %打开CH1.dat文件，并读取段索引表内容 fid1=fopen(data1,'r'); status=fseek(fid1,52,'bof'); segpulsenum1=fread(fid1,[1,24],'uint32'); status=fseek(fid1,16440,'bof'); %打开CH2.dat文件，并读取段索引表内容 fid2=fopen(data2,'r'); status=fseek(fid2,52,'bof'); segpulsenum2=fread(fid2,[1,24],'uint32'); status=fseek(fid2,16440,'bof'); i=1; j=1; b=0; b1=0; b2=0; k1=0; k2=0; p1=0; p2=0; a=0; over1=0; over2=0; C=zeros(2,2*N-1); t=0; %参数初始 %设置循环 while (k1<=23)|(k2<=23) if a==0 pos1=fread(fid1,1,'uint32'); over1=feof(fid1); if over1==1 break end k1=k1+(pos1<=p1); p1=pos1; b1=floor(p1/N); i=i+(t==2); elseif a==1 pos2=fread(fid2,1,'uint32'); over2=feof(fid2); if over2==1 break end k2=k2+(pos2<=p2); p2=pos2; b2=floor(p2/N); j=j+(t==1); end %由a确定取值文件。从ch1开始每次从两个文件中一个取一个值， %与前一个值比较，确定下一次取值文件为值小的那个。 if ((b2>b)&(t==1)) t=0; blo(2,11)=j-1; blo(2,1:j-1)=block(2,1:j-1); C(2,:)=Cab(blo,N); C(1,:)=sum(C); j=1; i=1; end if b1==b2 t=2; b=b1; elseif (b1~=b2)&(t==2) t=1; blo(1,11)=i-1; blo(1,1:i-1)=block(1,1:i-1); i=1; end if (t==1)&(b2<b) t=0; blo(2,11)=j-1; blo(2,1:j-1)=block(2,1:j-1); C(2,:)=Cab(blo,N); C(1,:)=sum(C); j=1; i=1; end block(1,i)=p1; block(2,j)=p2; if k1~=k2 a=(k1>k2); else a=~(b1<=b2); end end if t==2 blo(1,1)=block(1,1); blo(2,1)=block(2,1); blo(1,11)=1; blo(2,11)=1; C(2,:)=Cab(blo,N); C(1,:)=sum(C); end cc=C(1,:); status=fclose(fid1); status=fclose(fid2);