能谱（EDS）五大问答解析

%clear all %close all %clc function ballfitting(filelocation) %六路汇聚-硬X负载实验炮号设置 filestart=6039; fileend=filestart; %filelocation='D:\实验数据档案\硬X射线负载第二轮\'; %电缆长度矩阵,28*6矩阵，依次表示从101示波器到140示波器电缆长度 dianlanchangdu=zeros(28,6);%定义电缆长度矩阵 for i=1:28 dianlanchangdu(:,1)=59.9+10; dianlanchangdu(:,2)=56.7+10; dianlanchangdu(:,3)=59.8+10; dianlanchangdu(:,4)=62.9+10; dianlanchangdu(:,5)=65.9+10; dianlanchangdu(:,6)=62.9+10; end %1#电缆修正 dianlanchangdu(15,1)=59.9+20; dianlanchangdu(16,1)=59.9+20; dianlanchangdu(22,1)=59.9+11.7; dianlanchangdu(23,1)=59.9+11.6; dianlanchangdu(24,1)=59.9+11.6; dianlanchangdu(25,1)=59.9+11.6; dianlanchangdu(26,1)=59.9+11.6; dianlanchangdu(27,1)=59.9+11.6; %2#电缆修正 dianlanchangdu(15,2)=56.7+20; dianlanchangdu(16,2)=56.7+20; dianlanchangdu(22,2)=56.7+11.6+0.2;%经巴伦后，经0.2m电缆至示波器 dianlanchangdu(23,2)=56.7+11.6; dianlanchangdu(24,2)=56.7+11.6+0.2;%经巴伦后，经0.2m电缆至示波器 dianlanchangdu(25,2)=56.7+11.6; dianlanchangdu(26,2)=56.7+11.6+0.2;%经巴伦后，经0.2m电缆至示波器 dianlanchangdu(27,2)=56.7+11.6; %3#电缆修正 dianlanchangdu(1,3)=59.9+10; dianlanchangdu(3,3)=59.9+10; dianlanchangdu(4,3)=59.9+10; dianlanchangdu(7,3)=59.9+10; dianlanchangdu(15,3)=59.8+20; dianlanchangdu(16,3)=59.9+20; dianlanchangdu(19,3)=59.9+10; dianlanchangdu(22,3)=59.8+11.7; dianlanchangdu(23,3)=59.8+11.7; dianlanchangdu(24,3)=59.8+11.7; dianlanchangdu(25,3)=59.8+11.2; dianlanchangdu(28,3)=59.9+10; %4#电缆修正 dianlanchangdu(1,4)=62.8+10; dianlanchangdu(3,4)=62.8+10; dianlanchangdu(6,4)=62.7+10; dianlanchangdu(7,4)=62.7+10; dianlanchangdu(8,4)=62.8+10; dianlanchangdu(9,4)=62.8+10; dianlanchangdu(10,4)=62.7+10; dianlanchangdu(14,4)=62.7+10; dianlanchangdu(15,4)=62.9+20; dianlanchangdu(16,4)=62.7+20; dianlanchangdu(17,4)=62.7+10; dianlanchangdu(19,4)=62.8+10; dianlanchangdu(22,4)=62.7+12.1; dianlanchangdu(23,4)=62.7+12.1; dianlanchangdu(24,4)=62.9+12.1; dianlanchangdu(25,4)=62.9+12.1; dianlanchangdu(26,4)=62.8+12.1; dianlanchangdu(27,4)=62.6+12.1; %5#电缆修正 dianlanchangdu(4,5)=65.2+10; dianlanchangdu(9,5)=65.8+10; dianlanchangdu(10,5)=65+10; dianlanchangdu(15,5)=65.9+20; dianlanchangdu(16,5)=65.9+20; dianlanchangdu(18,5)=66+10; dianlanchangdu(19,5)=66+10; dianlanchangdu(22,5)=65.9+10.4; dianlanchangdu(23,5)=65.9+10.4; dianlanchangdu(24,5)=65.9+10.4; dianlanchangdu(25,5)=66+11.1; dianlanchangdu(26,5)=65.8+11.1; dianlanchangdu(27,5)=66+11.1; %6#电缆修正 dianlanchangdu(1,6)=62.8+10; dianlanchangdu(2,6)=62.8+10; dianlanchangdu(3,6)=62.8+10; dianlanchangdu(4,6)=62.7+10; dianlanchangdu(5,6)=62.8+10; dianlanchangdu(6,6)=62.7+10; dianlanchangdu(7,6)=62.8+10; dianlanchangdu(9,6)=62.7+10; dianlanchangdu(11,6)=62.8+10; dianlanchangdu(12,6)=62.8+10; dianlanchangdu(13,6)=62.8+10; dianlanchangdu(15,6)=62.9+20; dianlanchangdu(16,6)=62.7+20; dianlanchangdu(18,6)=62.8+10; dianlanchangdu(22,6)=62.8+10+1.5; dianlanchangdu(23,6)=62.8+10+2+0.5; dianlanchangdu(24,6)=62.8+10+1.5+0.5; dianlanchangdu(25,6)=62.8+10+1.5+0.5; dianlanchangdu(27,6)=62.7+10; A=10;%触发支路波形光滑系数 B=20; C=1;%六路汇聚系统测得波形光滑点数 %读取数据 for p=filestart:1:fileend filename=filelocation;%文件位置 file=csvread(filename,3,0); time=file(:,1)*10^9; for i=1:28 for j=1:1:6 time_xiuzheng(:,i,j)=time-dianlanchangdu(i,j)*5; end end %101 upperMITL_5(p-filestart+1,:)=smooth(162925.9/1000*(file(:,2)),C);%上MITL 5#，即上MITL 240°，单位kA time_upperMITL_5=time_xiuzheng(:,26,1); upperMITL_6(p-filestart+1,:)=smooth(166000/1000*(file(:,3)),C);%上MITL 6#，即上MITL 300°，单位kA time_upperMITL_6=time_xiuzheng(:,27,1); upperstackD_3(p-filestart+1,:)=smooth(3661290.3*(file(:,4)),C)/1000;%上堆D 3#，即上堆D 240°，采用巴伦，单位kV time_upperstackD_3=time_xiuzheng(:,26,2); upperstackD_3gong(p-filestart+1,:)=smooth(file(:,5),C);%上堆D 3#的共模信号，未经过巴伦，为微分信号 time_upperstackD_3gong=time_xiuzheng(:,27,2); %102 %LTD_3_LRLcurrent1(p-filestart+1,:)=-smooth(23.3*file(:,6),C);%未使用此通道 %time_LTD_3_LRLcurrent1=time_xiuzheng(:,26,3); %LTD_3_LRLcurrent2(p-filestart+1,:)=-smooth(21.5*file(:,7),C);%未使用此通道 %time_LTD_3_LRLcurrent2=time_xiuzheng(:,27,3); lowerMITL_5(p-filestart+1,:)=smooth(144264.1/1000*(file(:,8)),C);%下MITL 5#，即下MITL 240°，单位kA time_lowerMITL_5=time_xiuzheng(:,26,4); lowerMITL_6(p-filestart+1,:)=smooth(141959.8/1000*(file(:,9)),C);%下MITL 6#，即下MITL 300°，单位kA time_lowerMITL_6=time_xiuzheng(:,27,4); %103 lowerstackB_2(p-filestart+1,:)=smooth(142740/1000*(file(:,10)),C);%下堆B 2#，即下堆B 150°，单位kA time_lowerstackB_2=time_xiuzheng(:,26,5); lowerstackB_3(p-filestart+1,:)=smooth(142185.9/1000*(file(:,11)),C);%下堆B 3#，即下堆B 270°，单位kA time_lowerstackB_3=time_xiuzheng(:,27,5); %LTD_6_LRLcurrent1(p-filestart+1,:)=smooth(24.7*file(:,12),C);%未使用此通道 %time_LTD_6_LRLcurrent1=time_xiuzheng(:,26,6); PIN(p-filestart+1,:)=file(:,13);%PIN输出信号 time_PIN=time_xiuzheng(:,27,6);%PIN时间修正 %104 LTD_1_triggerbrick(:,1)=smooth(-file(:,14),A); time_LTD_1_triggerbrick(:,1)=time_xiuzheng(:,1,1); LTD_1_triggerbrick(:,2)=smooth(-file(:,15),A); time_LTD_1_triggerbrick(:,2)=time_xiuzheng(:,2,1); LTD_1_triggerbrick(:,3)=smooth(file(:,16),A); time_LTD_1_triggerbrick(:,3)=time_xiuzheng(:,3,1); LTD_1_triggerbrick(:,4)=smooth(-file(:,17),A); time_LTD_1_triggerbrick(:,4)=time_xiuzheng(:,4,1); %105 LTD_2_triggerbrick(:,1)=smooth(file(:,18),A); time_LTD_2_triggerbrick(:,1)=time_xiuzheng(:,1,2); LTD_2_triggerbrick(:,2)=smooth(-file(:,19),A); time_LTD_2_triggerbrick(:,2)=time_xiuzheng(:,2,2); LTD_2_triggerbrick(:,3)=smooth(-file(:,20),A); time_LTD_2_triggerbrick(:,3)=time_xiuzheng(:,3,2); LTD_2_triggerbrick(:,4)=smooth(file(:,21),A); time_LTD_2_triggerbrick(:,4)=time_xiuzheng(:,4,2); %106 LTD_2_triggerbrick(:,5)=smooth(file(:,22),A); time_LTD_2_triggerbrick(:,5)=time_xiuzheng(:,5,2); LTD_2_triggerbrick(:,6)=smooth(file(:,23),A); time_LTD_2_triggerbrick(:,6)=time_xiuzheng(:,6,2); LTD_2_triggerbrick(:,7)=smooth(file(:,24),A); time_LTD_2_triggerbrick(:,7)=time_xiuzheng(:,7,2); LTD_2_triggerbrick(:,8)=smooth(file(:,25),A); time_LTD_2_triggerbrick(:,8)=time_xiuzheng(:,8,2); %107 LTD_2_triggerbrick(:,9)=smooth(file(:,26),A); time_LTD_2_triggerbrick(:,9)=time_xiuzheng(:,9,2); LTD_2_triggerbrick(:,10)=smooth(file(:,27),A); time_LTD_2_triggerbrick(:,10)=time_xiuzheng(:,10,2); LTD_2_shuixianshouduanvoltage1(p-filestart+1,:)=smooth(71*file(:,28),B); time_LTD_2_shuixianshouduanvoltage1=time_xiuzheng(:,11,2); LTD_2_shuixianshouduanvoltage2(p-filestart+1,:)=smooth(75.2*file(:,29),B); time_LTD_2_shuixianshouduanvoltage2=time_xiuzheng(:,12,2); %108 LTD_2_shuixianmoduancurrent(p-filestart+1,:)=smooth(-214.2*file(:,30),B); time_LTD_2_shuixianmoduancurrent=time_xiuzheng(:,13,2); LTD_2_shuixianmoduanvoltage(p-filestart+1,:)=smooth(65.8*file(:,31),B); time_LTD_2_shuixianmoduanvoltage=time_xiuzheng(:,14,2); LTD_2_triggermachine1=smooth(file(:,32),10); time_LTD_2_triggermachine1=time_xiuzheng(:,15,2); LTD_2_triggermachine2=smooth(file(:,33),10); time_LTD_2_triggermachine2=time_xiuzheng(:,16,2); %109 LTD_3_triggerbrick(:,1)=smooth(-file(:,34),A); time_LTD_3_triggerbrick(:,1)=time_xiuzheng(:,1,3); LTD_3_triggerbrick(:,2)=smooth(file(:,35),A); time_LTD_3_triggerbrick(:,2)=time_xiuzheng(:,2,3); LTD_3_triggerbrick(:,3)=smooth(file(:,36),A); time_LTD_3_triggerbrick(:,3)=time_xiuzheng(:,3,3); LTD_3_triggerbrick(:,4)=smooth(-file(:,37),A); time_LTD_3_triggerbrick(:,4)=time_xiuzheng(:,4,3); %110 LTD_3_triggerbrick(:,5)=smooth(-file(:,38),A); time_LTD_3_triggerbrick(:,5)=time_xiuzheng(:,5,3); LTD_3_triggerbrick(:,6)=smooth(-file(:,39),A); time_LTD_3_triggerbrick(:,6)=time_xiuzheng(:,6,3); LTD_3_triggerbrick(:,7)=smooth(file(:,40),A); time_LTD_3_triggerbrick(:,7)=time_xiuzheng(:,7,3); LTD_3_triggerbrick(:,8)=smooth(file(:,41),A); time_LTD_3_triggerbrick(:,8)=time_xiuzheng(:,8,3); %111 PCD(p-filestart+1,:)=file(:,42);%PCD输出信号 time_PCD=time_xiuzheng(:,18,3);%PCD时间修正 %未使用此通道 %LTD_3_cijidianya(:,2)=smooth(-file(:,43),A); guangdianguan1(p-filestart+1,:)=file(:,44);%光电管1输出信号 time_guangdianguan1=time_xiuzheng(:,20,3);%光电管1时间修正 guangdianguan2(p-filestart+1,:)=file(:,45);%光电管2输出信号 time_guangdianguan2=time_xiuzheng(:,21,3);%光电管2时间修正 %112 upperstackB_1(p-filestart+1,:)=smooth(142277.2/1000*(file(:,46)),C);%上堆B 1#，即上堆B 60°，单位kA time_upperstackB_1=time_xiuzheng(:,22,3); upperstackB_2(p-filestart+1,:)=smooth(136930.2/1000*(file(:,47)),C);%上堆B 2#，即上堆B 180°，单位kA time_upperstackB_2=time_xiuzheng(:,23,3); upperstackB_3(p-filestart+1,:)=smooth(139473.2/1000*(file(:,48)),C);%上堆B 3#，即上堆B 300°，单位kA time_upperstackB_3=time_xiuzheng(:,24,3); LoadCurrent(p-filestart+1,:)=smooth(797486.7/1000*(file(:,49)),C);%短路负载罗氏线圈，单位kA time_LoadCurrent=time_xiuzheng(:,25,3); %113 LTD_4_triggerbrick(:,9)=smooth(-file(:,50),A); time_LTD_4_triggerbrick(:,9)=time_xiuzheng(:,9,4); LTD_4_triggerbrick(:,10)=smooth(-file(:,51),A); time_LTD_4_triggerbrick(:,10)=time_xiuzheng(:,10,4); LTD_4_shuixianshouduanvoltage1(p-filestart+1,:)=smooth(57.3*file(:,52),B); time_LTD_4_shuixianshouduanvoltage1=time_xiuzheng(:,11,4); LTD_4_shuixianshouduanvoltage2(p-filestart+1,:)=smooth(59.3*file(:,53),B); time_LTD_4_shuixianshouduanvoltage2=time_xiuzheng(:,12,4); %114 LTD_4_shuixianmoduancurrent(p-filestart+1,:)=smooth(345*file(:,54),B); time_LTD_4_shuixianmoduancurrent=time_xiuzheng(:,13,4); LTD_4_shuixianmoduanvoltage(p-filestart+1,:)=smooth(67.1*file(:,55),B); time_LTD_4_shuixianmoduanvoltage=time_xiuzheng(:,14,4); LTD_4_triggermachine1=smooth(file(:,56),10); time_LTD_4_triggermachine1=time_xiuzheng(:,15,4); LTD_4_triggermachine2=smooth(file(:,57),10); time_LTD_4_triggermachine2=time_xiuzheng(:,16,4); %115 LTD_5_triggerbrick(:,1)=smooth(-file(:,58),A); time_LTD_5_triggerbrick(:,1)=time_xiuzheng(:,1,5); LTD_5_triggerbrick(:,2)=smooth(-file(:,59),A); time_LTD_5_triggerbrick(:,2)=time_xiuzheng(:,2,5); LTD_5_triggerbrick(:,3)=smooth(-file(:,60),A); time_LTD_5_triggerbrick(:,3)=time_xiuzheng(:,3,5); LTD_5_triggerbrick(:,4)=smooth(-file(:,61),A); time_LTD_5_triggerbrick(:,4)=time_xiuzheng(:,4,5); %116 LTD_5_triggerbrick(:,5)=smooth(-file(:,62),A); time_LTD_5_triggerbrick(:,5)=time_xiuzheng(:,5,5); LTD_5_triggerbrick(:,6)=smooth(-file(:,63),A); time_LTD_5_triggerbrick(:,6)=time_xiuzheng(:,6,5); LTD_5_triggerbrick(:,7)=smooth(-file(:,64),A); time_LTD_5_triggerbrick(:,7)=time_xiuzheng(:,7,5); LTD_5_triggerbrick(:,8)=smooth(-file(:,65),A); time_LTD_5_triggerbrick(:,8)=time_xiuzheng(:,8,5); %117 LTD_5_cijidianya(:,1)=smooth(-file(:,66),A); LTD_5_cijidianya(:,2)=smooth(-file(:,67),A); LTD_5_cijidianya(:,3)=smooth(-file(:,68),A); LTD_5_cijidianya(:,4)=smooth(-file(:,69),A); %118 lowerstackD_1(p-filestart+1,:)=smooth(1730769.2/1000*(file(:,70)),C);%下堆D 1#，即下堆D 90°，单位kV time_lowerstackD_1=time_xiuzheng(:,22,5); lowerstackD_2(p-filestart+1,:)=smooth(1666666.7/1000*(file(:,71)),C);%下堆D 2#，即下堆D 210°，单位kV time_lowerstackD_2=time_xiuzheng(:,23,5); lowerstackD_3(p-filestart+1,:)=smooth(1682243/1000*(file(:,72)),C);%下堆D 3#，即下堆D 330°，单位kV time_lowerstackD_3=time_xiuzheng(:,24,5); lowerstackB_1(p-filestart+1,:)=smooth(146433.8/1000*(file(:,73)),C);%下堆B 1#，即下堆B 30°，单位kA time_lowerstackB_1=time_xiuzheng(:,25,5); %119 LTD_6_triggerbrick(:,9)=smooth(-file(:,74),A); time_LTD_6_triggerbrick(:,9)=time_xiuzheng(:,9,6); LTD_6_triggerbrick(:,10)=smooth(-file(:,75),A); time_LTD_6_triggerbrick(:,10)=time_xiuzheng(:,10,6); LTD_6_shuixianshouduanvoltage1(p-filestart+1,:)=smooth(71.6*file(:,76),B); time_LTD_6_shuixianshouduanvoltage1=time_xiuzheng(:,11,6); LTD_6_shuixianshouduanvoltage2(p-filestart+1,:)=smooth(86.3*file(:,77),B); time_LTD_6_shuixianshouduanvoltage2=time_xiuzheng(:,12,6); %120 LTD_6_shuixianmoduancurrent(p-filestart+1,:)=smooth(259.0*file(:,78),B); time_LTD_6_shuixianmoduancurrent=time_xiuzheng(:,13,6); LTD_6_shuixianmoduanvoltage(p-filestart+1,:)=smooth(80.3*file(:,79),B); time_LTD_6_shuixianmoduanvoltage=time_xiuzheng(:,14,6); LTD_6_triggermachine1=smooth(file(:,80),10); time_LTD_6_triggermachine1=time_xiuzheng(:,15,6); LTD_6_triggermachine2=smooth(file(:,81),10); time_LTD_6_triggermachine2=time_xiuzheng(:,16,6); %122 LTD_1_triggerbrick(:,5)=smooth(file(:,82),A); time_LTD_1_triggerbrick(:,5)=time_xiuzheng(:,5,1); LTD_1_triggerbrick(:,6)=smooth(file(:,83),A); time_LTD_1_triggerbrick(:,6)=time_xiuzheng(:,6,1); LTD_1_triggerbrick(:,7)=smooth(-file(:,84),A); time_LTD_1_triggerbrick(:,7)=time_xiuzheng(:,7,1); LTD_1_triggerbrick(:,8)=smooth(file(:,85),A); time_LTD_1_triggerbrick(:,8)=time_xiuzheng(:,8,1); %123 LTD_1_triggerbrick(:,9)=smooth(-file(:,86),A); time_LTD_1_triggerbrick(:,9)=time_xiuzheng(:,9,1); LTD_1_triggerbrick(:,10)=smooth(file(:,87),A); time_LTD_1_triggerbrick(:,10)=time_xiuzheng(:,10,1); LTD_1_shuixianshouduanvoltage1(p-filestart+1,:)=smooth(66.3*file(:,88),B); time_LTD_1_shuixianshouduanvoltage1=time_xiuzheng(:,11,1); LTD_1_shuixianshouduanvoltage2(p-filestart+1,:)=smooth(68.7*file(:,89),B); time_LTD_1_shuixianshouduanvoltage2=time_xiuzheng(:,12,1); %124 LTD_1_shuixianmoduancurrent(p-filestart+1,:)=smooth(-226*file(:,90),B); time_LTD_1_shuixianmoduancurrent=time_xiuzheng(:,13,1); LTD_1_shuixianmoduanvoltage(p-filestart+1,:)=smooth(66.8*file(:,91),B); time_LTD_1_shuixianmoduanvoltage=time_xiuzheng(:,14,1); LTD_1_triggermachine1=smooth(file(:,92),10); time_LTD_1_triggermachine1=time_xiuzheng(:,15,1); LTD_1_triggermachine2=smooth(file(:,93),10); time_LTD_1_triggermachine2=time_xiuzheng(:,16,1); %125 EDS1(p-filestart+1,:)=file(:,94);%能谱仪1输出信号 time_EDS1=time_xiuzheng(:,18,1);%能谱仪1时间修正 EDS2(p-filestart+1,:)=file(:,95);%能谱仪2输出信号 time_EDS2=time_xiuzheng(:,19,1);%能谱仪2时间修正 EDS3(p-filestart+1,:)=file(:,96);%能谱仪3输出信号 time_EDS3=time_xiuzheng(:,20,1);%能谱仪3时间修正 EDS4(p-filestart+1,:)=file(:,97);%能谱仪4输出信号 time_EDS4=time_xiuzheng(:,21,1);%能谱仪4时间修正 %126 upperMITL_1(p-filestart+1,:)=smooth(168252.6/1000*(file(:,98)),C);%上MITL 1#，即上MITL 0°，单位kA time_upperMITL_1=time_xiuzheng(:,22,1); upperMITL_2(p-filestart+1,:)=smooth(166000/1000*(file(:,99)),C);%上MITL 2#，即上MITL 60°，单位kA time_upperMITL_2=time_xiuzheng(:,23,1); upperMITL_3(p-filestart+1,:)=smooth(165108/1000*(file(:,100)),C);%上MITL 3#，即上MITL 120°，单位kA time_upperMITL_3=time_xiuzheng(:,24,1); upperMITL_4(p-filestart+1,:)=smooth(166666/1000*(file(:,101)),C);%上MITL 4#，即上MITL 180°，单位kA time_upperMITL_4=time_xiuzheng(:,25,1); %127 EDS5(p-filestart+1,:)=file(:,102);%能谱仪5输出信号 time_EDS5=time_xiuzheng(:,18,2);%能谱仪5时间修正 EDS6(p-filestart+1,:)=file(:,103);%能谱仪6输出信号 time_EDS6=time_xiuzheng(:,19,2);%能谱仪6时间修正 EDS7(p-filestart+1,:)=file(:,104);%能谱仪7输出信号 time_EDS7=time_xiuzheng(:,20,2);%能谱仪7时间修正 EDS8(p-filestart+1,:)=file(:,105);%能谱仪8输出信号 time_EDS8=time_xiuzheng(:,21,2);%能谱仪8时间修正 %128 upperstackD_1(p-filestart+1,:)=smooth(3681818.2/1000*(file(:,106)),C);%上堆D 1#，即上堆D 0°，采用巴伦，单位kV time_upperstackD_1=time_xiuzheng(:,22,2); %LTD_2_zhuixingduancurrent2(p-filestart+1,:)=smooth(26.2*file(:,107),C); %time_LTD_2_zhuixingduancurrent2=time_xiuzheng(:,23,2); upperstackD_2(p-filestart+1,:)=smooth(3645038.1/1000*(file(:,108)),C);%上堆D 2#，即上堆D 120°，采用巴伦，单位kV time_upperstackD_2=time_xiuzheng(:,24,2); %LTD_2_zhuixingduanvoltage2(p-filestart+1,:)=smooth(file(:,109),C); %129 LTD_3_triggerbrick(:,9)=smooth(-file(:,110),A); time_LTD_3_triggerbrick(:,9)=time_xiuzheng(:,9,3); LTD_3_triggerbrick(:,10)=smooth(file(:,111),A); time_LTD_3_triggerbrick(:,10)=time_xiuzheng(:,10,3); LTD_3_shuixianshouduanvoltage1(p-filestart+1,:)=smooth(67.1*file(:,112),B); time_LTD_3_shuixianshouduanvoltage1=time_xiuzheng(:,11,3); LTD_3_shuixianshouduanvoltage2(p-filestart+1,:)=smooth(84.5*file(:,113),B); time_LTD_3_shuixianshouduanvoltage2=time_xiuzheng(:,12,3); %130 LTD_3_shuixianmoduancurrent(p-filestart+1,:)=smooth(-279.9*file(:,114),B); time_LTD_3_shuixianmoduancurrent=time_xiuzheng(:,13,3); LTD_3_shuixianmoduanvoltage(p-filestart+1,:)=smooth(76.2*file(:,115),B); time_LTD_3_shuixianmoduanvoltage=time_xiuzheng(:,14,3); LTD_3_triggermachine1=smooth(file(:,116),10); time_LTD_3_triggermachine1=time_xiuzheng(:,15,3); LTD_3_triggermachine2=smooth(file(:,117),10); time_LTD_3_triggermachine2=time_xiuzheng(:,16,3); %131 LTD_4_triggerbrick(:,1)=smooth(-file(:,118),A); time_LTD_4_triggerbrick(:,1)=time_xiuzheng(:,1,4); LTD_4_triggerbrick(:,2)=smooth(-file(:,119),A); time_LTD_4_triggerbrick(:,2)=time_xiuzheng(:,2,4); LTD_4_triggerbrick(:,3)=smooth(-file(:,120),A); time_LTD_4_triggerbrick(:,3)=time_xiuzheng(:,3,4); LTD_4_triggerbrick(:,4)=smooth(-file(:,121),A); time_LTD_4_triggerbrick(:,4)=time_xiuzheng(:,4,4); %132 LTD_4_triggerbrick(:,5)=smooth(-file(:,122),A); time_LTD_4_triggerbrick(:,5)=time_xiuzheng(:,5,4); LTD_4_triggerbrick(:,6)=smooth(-file(:,123),A); time_LTD_4_triggerbrick(:,6)=time_xiuzheng(:,6,4); LTD_4_triggerbrick(:,7)=smooth(-file(:,124),A); time_LTD_4_triggerbrick(:,7)=time_xiuzheng(:,7,4); LTD_4_triggerbrick(:,8)=smooth(-file(:,125),A); time_LTD_4_triggerbrick(:,8)=time_xiuzheng(:,8,4); %133 EDS9(p-filestart+1,:)=file(:,126);%能谱仪9输出信号 time_EDS9=time_xiuzheng(:,18,4);%能谱仪9时间修正 %LTD_4_cijidianya(:,1)=smooth(-file(:,126),A); %LTD_4_cijidianya(:,2)=smooth(-file(:,127),A); %LTD_4_cijidianya(:,3)=smooth(-file(:,128),A); %LTD_4_cijidianya(:,4)=smooth(-file(:,129),A); %134 lowerMITL_1(p-filestart+1,:)=smooth(141174.4/1000*(file(:,130)),C);%下MITL 1#，即下MITL 0°，单位kA time_lowerMITL_1=time_xiuzheng(:,22,4); lowerMITL_2(p-filestart+1,:)=smooth(147767.6/1000*(file(:,131)),C);%下MITL 2#，即下MITL 60°，单位kA time_lowerMITL_2=time_xiuzheng(:,23,4); lowerMITL_3(p-filestart+1,:)=smooth(143505.1/1000*(file(:,132)),C);%下MITL 3#，即下MITL 120°，单位kA time_lowerMITL_3=time_xiuzheng(:,24,4); lowerMITL_4(p-filestart+1,:)=smooth(145618.6/1000*(file(:,133)),C);%下MITL 4#，即下MITL 180°，单位kA time_lowerMITL_4=time_xiuzheng(:,25,4); %135 LTD_5_triggerbrick(:,9)=smooth(-file(:,134),A); time_LTD_5_triggerbrick(:,9)=time_xiuzheng(:,9,5); LTD_5_triggerbrick(:,10)=smooth(-file(:,135),A); time_LTD_5_triggerbrick(:,10)=time_xiuzheng(:,10,5); LTD_5_shuixianshouduanvoltage1(p-filestart+1,:)=smooth(79.8*file(:,136),B); time_LTD_5_shuixianshouduanvoltage1=time_xiuzheng(:,11,5); LTD_5_shuixianshouduanvoltage2(p-filestart+1,:)=smooth(75.3*file(:,137),B); time_LTD_5_shuixianshouduanvoltage2=time_xiuzheng(:,12,5); %136 LTD_5_shuixianmoduancurrent(p-filestart+1,:)=smooth(244.3*file(:,138),B); time_LTD_5_shuixianmoduancurrent=time_xiuzheng(:,13,5); LTD_5_shuixianmoduanvoltage(p-filestart+1,:)=smooth(78.8*file(:,139),B); time_LTD_5_shuixianmoduanvoltage=time_xiuzheng(:,14,5); LTD_5_triggermachine1=smooth(file(:,140),10); time_LTD_5_triggermachine1=time_xiuzheng(:,15,5); LTD_5_triggermachine2=smooth(file(:,141),10); time_LTD_5_triggermachine2=time_xiuzheng(:,16,5); %137 LTD_6_triggerbrick(:,1)=smooth(-file(:,142),A); time_LTD_6_triggerbrick(:,1)=time_xiuzheng(:,1,6); LTD_6_triggerbrick(:,2)=smooth(-file(:,143),A); time_LTD_6_triggerbrick(:,2)=time_xiuzheng(:,2,6); LTD_6_triggerbrick(:,3)=smooth(-file(:,144),A); time_LTD_6_triggerbrick(:,3)=time_xiuzheng(:,3,6); LTD_6_triggerbrick(:,4)=smooth(-file(:,145),A); time_LTD_6_triggerbrick(:,4)=time_xiuzheng(:,4,6); %138 LTD_6_triggerbrick(:,5)=smooth(-file(:,146),A); time_LTD_6_triggerbrick(:,5)=time_xiuzheng(:,5,6); LTD_6_triggerbrick(:,6)=smooth(-file(:,147),A); time_LTD_6_triggerbrick(:,6)=time_xiuzheng(:,6,6); LTD_6_triggerbrick(:,7)=smooth(-file(:,148),A); time_LTD_6_triggerbrick(:,7)=time_xiuzheng(:,7,6); LTD_6_triggerbrick(:,8)=smooth(-file(:,149),A); time_LTD_6_triggerbrick(:,8)=time_xiuzheng(:,8,6); %139 lowerID1=file(:,150)*360;%下ID1，B-dot接积分器，系数360kA/V，输出信号单位kA time_lowerID1=time_xiuzheng(:,18,6);%下ID1时间修正 lowerID3=file(:,151)*334;%下ID3，B-dot接积分器，系数334kA/V，输出信号单位kA time_lowerID3=time_xiuzheng(:,19,6);%下ID3时间修正 upperID1=file(:,152)*767;%上ID1，B-dot接积分器，系数767kA/V，输出信号单位kA time_upperID1=time_xiuzheng(:,20,6);%上ID1时间修正 upperID3=file(:,153)*763;%上ID3，B-dot接积分器，系数763kA/V，输出信号单位kA time_upperID3=time_xiuzheng(:,21,6);%上ID3时间修正 %140 lowerID2=file(:,154);%下ID2，B-dot输出微分信号，输出信号单位V time_lowerID2=time_xiuzheng(:,22,6);%下ID2时间修正 lowerID4=file(:,155);%下ID4，B-dot输出微分信号，输出信号单位V time_lowerID4=time_xiuzheng(:,23,6);%下ID4时间修正 upperID2=file(:,156);%上ID2，B-dot输出微分信号，输出信号单位V time_upperID2=time_xiuzheng(:,24,6);%上ID2时间修正 upperID4=file(:,157);%上ID4，B-dot输出微分信号，输出信号单位V time_upperID4=time_xiuzheng(:,25,6);%上ID4时间修正 %---------------------------------------- %读取 LTD_1触发支路起点 figure(1) for j=1:1:10 %去零漂，取前1000个点的平均值 for i=1:1:1000 LTD_1_triggerbrick_qian1000(i,j)=LTD_1_triggerbrick(i,j); end LTD_1_triggerbrick_lingpiao(j)=mean(LTD_1_triggerbrick_qian1000(:,j)); LTD_1_triggerbrick(:,j)=LTD_1_triggerbrick(:,j)-LTD_1_triggerbrick_lingpiao(j); %归一化，除以最大值 LTD_1_triggerbrickmax(j)=max(LTD_1_triggerbrick(:,j)); LTD_1_triggerbrick(:,j)=LTD_1_triggerbrick(:,j)/LTD_1_triggerbrickmax(j); end LTD_1_triggerbrick_qidian=zeros(10,1); for j=1:1:10 for i=1:1:10000 if LTD_1_triggerbrick(i,j)>0.2 LTD_1_triggerbrick_qidian(j)=time_LTD_1_triggerbrick(i); break end end end %读取触发支路开关延时数据用于计算抖动 for i=1:1:10000 if LTD_1_triggermachine1(i)<-0.1 LTD_1_triggermachine1_qidian=time_LTD_1_triggermachine1(i); break end end for i=1:1:10000 if LTD_1_triggermachine2(i)<-0.1 LTD_1_triggermachine2_qidian=time_LTD_1_triggermachine1(i); break end end for x=1:1:5 LTD_1_chufazhiluyanshi(p-filestart+1,x)=LTD_1_triggerbrick_qidian(x)-LTD_1_triggermachine1_qidian+50; end for x=6:1:10 LTD_1_chufazhiluyanshi(p-filestart+1,x)=LTD_1_triggerbrick_qidian(x)-LTD_1_triggermachine2_qidian+50; end %LTD_1触发支路时序画图 XX_1=1:1:10; subplot(2,3,1) plot(XX_1,LTD_1_triggerbrick_qidian(XX_1)-LTD_1_triggerbrick_qidian(1),'-o','MarkerSize',10,'linewidth',1) hold on xx_1=linspace(0,11,100); yy_1=7*(xx_1-1); plot(xx_1,yy_1,'linewidth',2) grid on axis([0 11 -10 80]) title('LTD-1#触发时序') xlabel('级数') ylabel('触发支路电流起点') %---------------------------------------- %读取 LTD_2触发支路起点 for j=1:1:10 %去零漂，取前1000个点的平均值 for i=1:1:1000 LTD_2_triggerbrick_qian1000(i,j)=LTD_2_triggerbrick(i,j); end LTD_2_triggerbrick_lingpiao(j)=mean(LTD_2_triggerbrick_qian1000(:,j)); LTD_2_triggerbrick(:,j)=LTD_2_triggerbrick(:,j)-LTD_2_triggerbrick_lingpiao(j); %归一化，除以最大值 LTD_2_triggerbrickmax(j)=max(LTD_2_triggerbrick(:,j)); LTD_2_triggerbrick(:,j)=LTD_2_triggerbrick(:,j)/LTD_2_triggerbrickmax(j); end LTD_2_triggerbrick_qidian=zeros(10,1); for j=1:1:10 for i=1:1:10000 if LTD_2_triggerbrick(i,j)>0.2 LTD_2_triggerbrick_qidian(j)=time_LTD_2_triggerbrick(i); break end end end %读取触发支路开关延时数据用于计算抖动 for i=1:1:10000 if LTD_2_triggermachine1(i)<-0.1 LTD_2_triggermachine1_qidian=time_LTD_2_triggermachine1(i); break end end for i=1:1:10000 if LTD_2_triggermachine2(i)<-0.1 LTD_2_triggermachine2_qidian=time_LTD_2_triggermachine1(i); break end end for x=1:1:5 LTD_2_chufazhiluyanshi(p-filestart+1,x)=LTD_2_triggerbrick_qidian(x)-LTD_2_triggermachine1_qidian+50; end for x=6:1:10 LTD_2_chufazhiluyanshi(p-filestart+1,x)=LTD_2_triggerbrick_qidian(x)-LTD_2_triggermachine2_qidian+50; end %LTD_2触发支路时序画图 XX_2=1:1:10; subplot(2,3,2) plot(XX_2,LTD_2_triggerbrick_qidian(XX_2)-LTD_2_triggerbrick_qidian(1),'-o','MarkerSize',10,'linewidth',1) hold on xx_2=linspace(0,11,100); yy_2=7*(xx_2-1); plot(xx_2,yy_2,'linewidth',2) grid on axis([0 11 -10 80]) title('LTD-2#触发时序') xlabel('级数') ylabel('触发支路电流起点') %---------------------------------------- %读取 LTD_3触发支路起点 for j=1:1:10 %去零漂，取前1000个点的平均值 for i=1:1:1000 LTD_3_triggerbrick_qian1000(i,j)=LTD_3_triggerbrick(i,j); end LTD_3_triggerbrick_lingpiao(j)=mean(LTD_3_triggerbrick_qian1000(:,j)); LTD_3_triggerbrick(:,j)=LTD_3_triggerbrick(:,j)-LTD_3_triggerbrick_lingpiao(j); %归一化，除以最大值 LTD_3_triggerbrickmax(j)=max(LTD_3_triggerbrick(:,j)); LTD_3_triggerbrick(:,j)=LTD_3_triggerbrick(:,j)/LTD_3_triggerbrickmax(j); end LTD_3_triggerbrick_qidian=zeros(10,1); for j=1:1:10 for i=1:1:10000 if LTD_3_triggerbrick(i,j)>0.2 LTD_3_triggerbrick_qidian(j)=time_LTD_3_triggerbrick(i); break end end end %读取触发支路开关延时数据用于计算抖动 for i=1:1:10000 if LTD_3_triggermachine1(i)<-0.1 LTD_3_triggermachine1_qidian=time_LTD_3_triggermachine1(i); break end end for i=1:1:10000 if LTD_3_triggermachine2(i)<-0.1 LTD_3_triggermachine2_qidian=time_LTD_3_triggermachine1(i); break end end for x=1:1:5 LTD_3_chufazhiluyanshi(p-filestart+1,x)=LTD_3_triggerbrick_qidian(x)-LTD_3_triggermachine1_qidian+50; end for x=6:1:10 LTD_3_chufazhiluyanshi(p-filestart+1,x)=LTD_3_triggerbrick_qidian(x)-LTD_3_triggermachine2_qidian+50; end %LTD_3触发支路时序画图 XX_3=1:1:10; subplot(2,3,3) plot(XX_3,LTD_3_triggerbrick_qidian(XX_3)-LTD_3_triggerbrick_qidian(1),'-o','MarkerSize',10,'linewidth',1) hold on xx_3=linspace(0,11,100); yy_3=7*(xx_3-1); plot(xx_3,yy_3,'linewidth',2) grid on axis([0 11 -10 80]) title('LTD-3#触发时序') xlabel('级数') ylabel('触发支路电流起点') %---------------------------------------- %读取 LTD_4触发支路起点 for j=1:1:10 %去零漂，取前1000个点的平均值 for i=1:1:1000 LTD_4_triggerbrick_qian1000(i,j)=LTD_4_triggerbrick(i,j); end LTD_4_triggerbrick_lingpiao(j)=mean(LTD_4_triggerbrick_qian1000(:,j)); LTD_4_triggerbrick(:,j)=LTD_4_triggerbrick(:,j)-LTD_4_triggerbrick_lingpiao(j); %归一化，除以最大值 LTD_4_triggerbrickmax(j)=max(LTD_4_triggerbrick(:,j)); LTD_4_triggerbrick(:,j)=LTD_4_triggerbrick(:,j)/LTD_4_triggerbrickmax(j); end LTD_4_triggerbrick_qidian=zeros(10,1); for j=1:1:10 for i=1:1:10000 if LTD_4_triggerbrick(i,j)>0.2 LTD_4_triggerbrick_qidian(j)=time_LTD_4_triggerbrick(i); break end end end %读取触发支路开关延时数据用于计算抖动 for i=1:1:10000 if LTD_4_triggermachine1(i)<-0.1 LTD_4_triggermachine1_qidian=time_LTD_4_triggermachine1(i); break end end for i=1:1:10000 if LTD_4_triggermachine2(i)<-0.1 LTD_4_triggermachine2_qidian=time_LTD_4_triggermachine1(i); break end end for x=1:1:5 LTD_4_chufazhiluyanshi(p-filestart+1,x)=LTD_4_triggerbrick_qidian(x)-LTD_4_triggermachine1_qidian+50; end for x=6:1:10 LTD_4_chufazhiluyanshi(p-filestart+1,x)=LTD_4_triggerbrick_qidian(x)-LTD_4_triggermachine2_qidian+50; end %LTD_4触发支路时序画图 XX_4=1:1:10; subplot(2,3,4) plot(XX_4,LTD_4_triggerbrick_qidian(XX_4)-LTD_4_triggerbrick_qidian(1),'-o','MarkerSize',10,'linewidth',1) hold on xx_4=linspace(0,11,100); yy_4=7*(xx_4-1); plot(xx_4,yy_4,'linewidth',2) grid on axis([0 11 -10 80]) title('LTD-4#触发时序') xlabel('级数') ylabel('触发支路电流起点') %---------------------------------------- %读取 LTD_5触发支路起点 for j=1:1:10 %去零漂，取前1000个点的平均值 for i=1:1:1000 LTD_5_triggerbrick_qian1000(i,j)=LTD_5_triggerbrick(i,j); end LTD_5_triggerbrick_lingpiao(j)=mean(LTD_5_triggerbrick_qian1000(:,j)); LTD_5_triggerbrick(:,j)=LTD_5_triggerbrick(:,j)-LTD_5_triggerbrick_lingpiao(j); %归一化，除以最大值 LTD_5_triggerbrickmax(j)=max(LTD_5_triggerbrick(:,j)); LTD_5_triggerbrick(:,j)=LTD_5_triggerbrick(:,j)/LTD_5_triggerbrickmax(j); end LTD_5_triggerbrick_qidian=zeros(10,1); for j=1:1:10 for i=1:1:10000 if LTD_5_triggerbrick(i,j)>0.2 LTD_5_triggerbrick_qidian(j)=time_LTD_5_triggerbrick(i); break end end end %读取触发支路开关延时数据用于计算抖动 for i=1:1:10000 if LTD_5_triggermachine1(i)<-0.1 LTD_5_triggermachine1_qidian=time_LTD_5_triggermachine1(i); break end end for i=1:1:10000 if LTD_5_triggermachine2(i)<-0.1 LTD_5_triggermachine2_qidian=time_LTD_5_triggermachine1(i); break end end for x=1:1:5 LTD_5_chufazhiluyanshi(p-filestart+1,x)=LTD_5_triggerbrick_qidian(x)-LTD_5_triggermachine1_qidian+50; end for x=6:1:10 LTD_5_chufazhiluyanshi(p-filestart+1,x)=LTD_5_triggerbrick_qidian(x)-LTD_5_triggermachine2_qidian+50; end %LTD_5触发支路时序画图 XX_5=1:1:10; subplot(2,3,5) plot(XX_5,LTD_5_triggerbrick_qidian(XX_5)-LTD_5_triggerbrick_qidian(1),'-o','MarkerSize',10,'linewidth',1) hold on xx_5=linspace(0,11,100); yy_5=7*(xx_4-1); plot(xx_5,yy_5,'linewidth',2) grid on axis([0 11 -10 80]) title('LTD-5#触发时序') xlabel('级数') ylabel('触发支路电流起点') %---------------------------------------- %读取 LTD_6触发支路起点 for j=1:1:10 %去零漂，取前1000个点的平均值 for i=1:1:1000 LTD_6_triggerbrick_qian1000(i,j)=LTD_6_triggerbrick(i,j); end LTD_6_triggerbrick_lingpiao(j)=mean(LTD_6_triggerbrick_qian1000(:,j)); LTD_6_triggerbrick(:,j)=LTD_6_triggerbrick(:,j)-LTD_6_triggerbrick_lingpiao(j); %归一化，除以最大值 LTD_6_triggerbrickmax(j)=max(LTD_6_triggerbrick(:,j)); LTD_6_triggerbrick(:,j)=LTD_6_triggerbrick(:,j)/LTD_6_triggerbrickmax(j); end LTD_6_triggerbrick_qidian=zeros(10,1); for j=1:1:10 for i=1:1:10000 if LTD_6_triggerbrick(i,j)>0.2 LTD_6_triggerbrick_qidian(j)=time_LTD_6_triggerbrick(i); break end end end %读取触发支路开关延时数据用于计算抖动 for i=1:1:10000 if LTD_6_triggermachine1(i)<-0.1 LTD_6_triggermachine1_qidian=time_LTD_6_triggermachine1(i); break end end for i=1:1:10000 if LTD_6_triggermachine2(i)<-0.1 LTD_6_triggermachine2_qidian=time_LTD_6_triggermachine1(i); break end end for x=1:1:5 LTD_6_chufazhiluyanshi(p-filestart+1,x)=LTD_6_triggerbrick_qidian(x)-LTD_6_triggermachine1_qidian+50; end for x=6:1:10 LTD_6_chufazhiluyanshi(p-filestart+1,x)=LTD_6_triggerbrick_qidian(x)-LTD_6_triggermachine2_qidian+50; end %LTD_6触发支路时序画图 XX_6=1:1:10; subplot(2,3,6) plot(XX_6,LTD_6_triggerbrick_qidian(XX_6)-LTD_6_triggerbrick_qidian(1),'-o','MarkerSize',10,'linewidth',1) hold on xx_6=linspace(0,11,100); yy_6=7*(xx_4-1); plot(xx_6,yy_6,'linewidth',2) grid on axis([0 11 -10 80]) title('LTD-6#触发时序') xlabel('级数') ylabel('触发支路电流起点') end set(gcf,'Position',get(0,'Screensize')); %计算触发器时序 triggermachine_qidian=[LTD_1_triggermachine1_qidian,LTD_2_triggermachine1_qidian,LTD_3_triggermachine1_qidian,LTD_4_triggermachine1_qidian,LTD_5_triggermachine1_qidian,LTD_6_triggermachine1_qidian; LTD_1_triggermachine2_qidian,LTD_2_triggermachine2_qidian,LTD_3_triggermachine2_qidian,LTD_4_triggermachine2_qidian,LTD_5_triggermachine2_qidian,LTD_6_triggermachine2_qidian; LTD_1_triggermachine2_qidian-LTD_1_triggermachine1_qidian,LTD_2_triggermachine2_qidian-LTD_2_triggermachine1_qidian,LTD_3_triggermachine2_qidian-LTD_3_triggermachine1_qidian,LTD_4_triggermachine2_qidian-LTD_4_triggermachine1_qidian,LTD_5_triggermachine2_qidian-LTD_5_triggermachine1_qidian,LTD_6_triggermachine2_qidian-LTD_6_triggermachine1_qidian]; %其他各级1#触发器与第1路1#触发器的时延差 triggermachine_differ=[LTD_1_triggermachine1_qidian+230-28 LTD_1_triggermachine2_qidian+202-28;LTD_2_triggermachine1_qidian+246-28 LTD_2_triggermachine2_qidian+220-28;LTD_3_triggermachine1_qidian+234-28 LTD_3_triggermachine2_qidian+201-27;LTD_4_triggermachine1_qidian+267-28 LTD_4_triggermachine2_qidian+237-28;LTD_5_triggermachine1_qidian+257-28 LTD_5_triggermachine2_qidian+233-30;LTD_6_triggermachine1_qidian+263-28 LTD_6_triggermachine2_qidian+233-27]; figure(2) bar(triggermachine_differ) title('各级触发器时序') axis([0 7 -10 10]) grid on hold on set(gcf,'Position',get(0,'Screensize')); %---------------------------------------- %六路LTD测量结果显示 %六路触发测量结果绘图 figure(3) subplot(2,3,1) plot(time_LTD_1_triggermachine1,LTD_1_triggermachine1,time_LTD_1_triggermachine2,LTD_1_triggermachine2,'LineWidth',1.5); title('第一路触发器时序') axis([-400 600 -2 0.5]) legend('1-1#','1-2#') xlabel('时间/ns') ylabel('信号幅值/V') grid on hold on subplot(2,3,2) plot(time_LTD_2_triggermachine1,LTD_2_triggermachine1,time_LTD_2_triggermachine2,LTD_2_triggermachine2,'LineWidth',1.5) title('第二路触发器时序') axis([-400 600 -2 0.5]) legend('2-1#','2-2#') xlabel('时间/ns') ylabel('信号幅值/V') grid on hold on subplot(2,3,3) plot(time_LTD_3_triggermachine1,LTD_3_triggermachine1,time_LTD_3_triggermachine2,LTD_3_triggermachine2,'LineWidth',1.5) title('第三路触发器时序') axis([-400 600 -2 0.5]) legend('3-1#','3-2#') xlabel('时间/ns') ylabel('信号幅值/V') grid on hold on subplot(2,3,4) plot(time_LTD_4_triggermachine1,LTD_4_triggermachine1,time_LTD_4_triggermachine2,LTD_4_triggermachine2,'LineWidth',1.5) title('第四路触发器时序') axis([-400 600 -2 0.5]) legend('4-1#','4-2#') xlabel('时间/ns') ylabel('信号幅值/V') grid on hold on subplot(2,3,5) plot(time_LTD_5_triggermachine1,LTD_5_triggermachine1,time_LTD_5_triggermachine2,LTD_5_triggermachine2,'LineWidth',1.5) title('第五路触发器时序') axis([-400 600 -2 0.5]) legend('5-1#','5-2#') xlabel('时间/ns') ylabel('信号幅值/V') grid on hold on subplot(2,3,6) plot(time_LTD_6_triggermachine1,LTD_6_triggermachine1,time_LTD_6_triggermachine2,LTD_6_triggermachine2,'LineWidth',1.5) title('第六路触发器时序') axis([-400 600 -2 0.5]) legend('6-1#','6-2#') xlabel('时间/ns') ylabel('信号幅值/V') grid on hold on set(gcf,'Position',get(0,'Screensize')); %---------------------------------------- %六路触发支路测量结果绘图 figure(4) subplot(2,3,1) for trbri=1:1:10 plot(time_LTD_1_triggerbrick(:,trbri),LTD_1_triggerbrick(:,trbri),'LineWidth',1.5) hold on end title('第一路触发支路测量结果') legend('1-1#','1-2#','1-3#','1-4#','1-5#','1-6#','1-7#','1-8#','1-9#','1-10#') axis([-200 300 -1 1.5]) xlabel('时间/ns') ylabel('信号幅值/V') grid on hold on subplot(2,3,2) for trbri=1:1:10 plot(time_LTD_2_triggerbrick(:,trbri),LTD_2_triggerbrick(:,trbri),'LineWidth',1.5) hold on end title('第二路触发支路测量结果') legend('2-1#','2-2#','2-3#','2-4#','2-5#','2-6#','2-7#','2-8#','2-9#','2-10#') axis([-200 300 -1 1.5]) xlabel('时间/ns') ylabel('信号幅值/V') grid on hold on subplot(2,3,3) for trbri=1:1:10 plot(time_LTD_3_triggerbrick(:,trbri),LTD_3_triggerbrick(:,trbri),'LineWidth',1.5) hold on end title('第三路触发支路测量结果') legend('3-1#','3-2#','3-3#','3-4#','3-5#','3-6#','3-7#','3-8#','3-9#','3-10#') axis([-200 300 -1 1.5]) xlabel('时间/ns') ylabel('信号幅值/V') grid on hold on subplot(2,3,4) for trbri=1:1:10 plot(time_LTD_4_triggerbrick(:,trbri),LTD_4_triggerbrick(:,trbri),'LineWidth',1.5) hold on end title('第四路触发支路测量结果') legend('4-1#','4-2#','4-3#','4-4#','4-5#','4-6#','4-7#','4-8#','4-9#','4-10#') axis([-200 300 -1 1.5]) xlabel('时间/ns') ylabel('信号幅值/V') grid on hold on subplot(2,3,5) for trbri=1:1:10 plot(time_LTD_5_triggerbrick(:,trbri),LTD_5_triggerbrick(:,trbri),'LineWidth',1.5) hold on end title('第五路触发支路测量结果') legend('5-1#','5-2#','5-3#','5-4#','5-5#','5-6#','5-7#','5-8#','5-9#','5-10#') axis([-200 300 -1 1.5]) xlabel('时间/ns') ylabel('信号幅值/V') grid on hold on subplot(2,3,6) for trbri=1:1:10 plot(time_LTD_6_triggerbrick(:,trbri),LTD_6_triggerbrick(:,trbri),'LineWidth',1.5) hold on end title('第六路触发支路测量结果') legend('6-1#','6-2#','6-3#','6-4#','6-5#','6-6#','6-7#','6-8#','6-9#','6-10#') axis([-200 300 -1 1.5]) xlabel('时间/ns') ylabel('信号幅值/V') grid on hold on set(gcf,'Position',get(0,'Screensize')); %---------------------------------------- %六路水线测量结果绘图 figure(5) subplot(2,2,1) plot(time_LTD_1_shuixianshouduanvoltage1,LTD_1_shuixianshouduanvoltage1,time_LTD_2_shuixianshouduanvoltage1,LTD_2_shuixianshouduanvoltage1,time_LTD_3_shuixianshouduanvoltage1,LTD_3_shuixianshouduanvoltage1,time_LTD_4_shuixianshouduanvoltage1,LTD_4_shuixianshouduanvoltage1,time_LTD_5_shuixianshouduanvoltage1,LTD_5_shuixianshouduanvoltage1,time_LTD_6_shuixianshouduanvoltage1,LTD_6_shuixianshouduanvoltage1,'LineWidth',1.5) title('水线首端电压1') axis([-200 1000 -1000 1000]) legend('LTD-1#','LTD-2#','LTD-3#','LTD-4#','LTD-5#','LTD-6#') xlabel('时间/ns') ylabel('电压/kV') grid on hold on subplot(2,2,2) plot(time_LTD_1_shuixianshouduanvoltage2,LTD_1_shuixianshouduanvoltage2,time_LTD_2_shuixianshouduanvoltage2,LTD_2_shuixianshouduanvoltage2,time_LTD_3_shuixianshouduanvoltage2,LTD_3_shuixianshouduanvoltage2,time_LTD_4_shuixianshouduanvoltage2,LTD_4_shuixianshouduanvoltage2,time_LTD_5_shuixianshouduanvoltage2,LTD_5_shuixianshouduanvoltage2,time_LTD_6_shuixianshouduanvoltage2,LTD_6_shuixianshouduanvoltage2,'LineWidth',1.5) title('水线首端电压2') axis([-200 1000 -1000 1000]) legend('LTD-1#','LTD-2#','LTD-3#','LTD-4#','LTD-5#','LTD-6#') xlabel('时间/ns') ylabel('电压/kV') grid on hold on subplot(2,2,3) plot(time_LTD_1_shuixianmoduanvoltage,LTD_1_shuixianmoduanvoltage,time_LTD_2_shuixianmoduanvoltage,LTD_2_shuixianmoduanvoltage,time_LTD_3_shuixianmoduanvoltage,LTD_3_shuixianmoduanvoltage,time_LTD_4_shuixianmoduanvoltage,LTD_4_shuixianmoduanvoltage,time_LTD_5_shuixianmoduanvoltage,LTD_5_shuixianmoduanvoltage,time_LTD_6_shuixianmoduanvoltage,LTD_6_shuixianmoduanvoltage,'LineWidth',1.5) title('水线末端电压') axis([-200 1000 -1000 1000]) legend('LTD-1#','LTD-2#','LTD-3#','LTD-4#','LTD-5#','LTD-6#') xlabel('时间/ns') ylabel('电压/kV') grid on hold on subplot(2,2,4) plot(time_LTD_1_shuixianmoduancurrent,LTD_1_shuixianmoduancurrent,time_LTD_2_shuixianmoduancurrent,LTD_2_shuixianmoduancurrent,time_LTD_3_shuixianmoduancurrent,LTD_3_shuixianmoduancurrent,time_LTD_4_shuixianmoduancurrent,LTD_4_shuixianmoduancurrent,time_LTD_5_shuixianmoduancurrent,LTD_5_shuixianmoduancurrent,time_LTD_6_shuixianmoduancurrent,LTD_6_shuixianmoduancurrent,'LineWidth',1.5) title('水线末端电流') axis([-200 1000 -600 600]) legend('LTD-1#','LTD-2#','LTD-3#','LTD-4#','LTD-5#','LTD-6#') xlabel('时间/ns') ylabel('电压/kV') grid on hold on set(gcf,'Position',get(0,'Screensize')); %---------------------------------------- %六路Stack-MITL测量结果绘图 figure(6) subplot(2,2,1) plot(time_upperstackB_1,upperstackB_1,time_upperstackB_2,upperstackB_2,time_upperstackB_3,upperstackB_3,time_lowerstackB_1,lowerstackB_1,time_lowerstackB_2,lowerstackB_2,time_lowerstackB_3,lowerstackB_3,'LineWidth',1.5) axis([100 800 -500 1600]) title('堆电流（正北方向为0°，顺时针编号）') legend('上堆B 60°','上堆B 180°','上堆B 300°','下堆B 30°','下堆B 150°','下堆B 270°') xlabel('时间/ns') ylabel('电流/kA') grid on hold on subplot(2,2,2) plot(time_upperstackD_1,upperstackD_1,time_upperstackD_2,upperstackD_2,time_upperstackD_3,upperstackD_3,time_lowerstackD_1,lowerstackD_1,time_lowerstackD_2,lowerstackD_2,time_lowerstackD_3,lowerstackD_3,'LineWidth',1.5) axis([100 800 -1000 1200]) title('堆电压（正北方向为0°，顺时针编号）') legend('上堆D 0°','上堆D 120°','上堆D 240°','下堆D 90°','下堆D 210°','下堆D 330°') xlabel('时间/ns') ylabel('电压/kV') grid on hold on subplot(2,2,3) plot(time_upperMITL_1,upperMITL_1,time_upperMITL_2,upperMITL_2,time_upperMITL_3,upperMITL_3,time_upperMITL_4,upperMITL_4,time_upperMITL_5,upperMITL_5,time_upperMITL_6,upperMITL_6,time_lowerMITL_1,lowerMITL_1,time_lowerMITL_2,lowerMITL_2,time_lowerMITL_3,lowerMITL_3,time_lowerMITL_4,lowerMITL_4,time_lowerMITL_5,lowerMITL_5,time_lowerMITL_6,lowerMITL_6,'LineWidth',1.5) axis([100 800 -500 1600]) title('MITL电流（正北方向为0°，顺时针编号）') legend('上MITL 0°','上MITL 60°','上MITL 120°','上MITL 180°','上MITL 240°','上MITL 300°','下MITL 0°','下MITL 60°','下MITL 120°','下MITL 180°','下MITL 240°','下MITL 300°') xlabel('时间/ns') ylabel('电流/kA') grid on hold on set(gcf,'Position',get(0,'Screensize')); %---------------------------------------- %硬X负载测量结果绘图 figure(7) subplot(2,2,1) plot(time_lowerID1,lowerID1,time_lowerID3,lowerID3,time_upperID1,upperID1,time_upperID3,upperID3,'LineWidth',1.5) title('负载电流积分信号') legend('下ID1','下ID3','上ID1','上ID3') xlabel('时间/ns') ylabel('电流/kA') axis([0 800 -2000 2000]) grid on hold on subplot(2,2,2) %plot(time_lowerID2,lowerID2,time_lowerID4,lowerID4,time_upperID2,upperID2,time_upperID4,upperID4,'LineWidth',1.5) %对sumlowerID2和sumlowerID4的计算过程进行修改，实现去零飘；不再构造sumupperID2和sumupperID4两个变量。2025.04.11，张梦凡 % sumlowerID2=zeros(10000,1); % sumlowerID4=zeros(10000,1); % sumupperID2=zeros(10000,1); % sumupperID4=zeros(10000,1); % for sumi=2:10000 % sumlowerID2(sumi,1)=sumlowerID2(sumi-1)+(lowerID2(sumi)+lowerID2(sumi-1))*(0.4E-9)/2;%积分信号，lowerID2 % sumlowerID4(sumi,1)=sumlowerID4(sumi-1)+(lowerID4(sumi)+lowerID4(sumi-1))*(0.4E-9)/2;%积分信号，lowerID4 % sumupperID2(sumi,1)=sumupperID2(sumi-1)+(upperID2(sumi)+upperID2(sumi-1))*(0.4E-9)/2;%积分信号，upperID2 % sumupperID4(sumi,1)=sumupperID4(sumi-1)+(upperID4(sumi)+upperID4(sumi-1))*(0.4E-9)/2;%积分信号，upperID4 % end % sumlowerID2=sumlowerID2; % sumlowerID4=sumlowerID4; sumlowerID2=cumtrapz(time_lowerID2,lowerID2-sum(lowerID2(1:1000))/1000); sumlowerID4=cumtrapz(time_lowerID4,lowerID4-sum(lowerID4(1:1000))/1000); lowerID1guiyi=lowerID1/min(lowerID1); lowerID2guiyi=(sumlowerID2)/min(sumlowerID2); lowerID3guiyi=lowerID3/min(lowerID3); lowerID4guiyi=(sumlowerID4)/min(sumlowerID4); plot(time_lowerID1,lowerID1guiyi,time_lowerID2,lowerID2guiyi,time_lowerID3,lowerID3guiyi,time_lowerID4,lowerID4guiyi,'LineWidth',1.5) title('归一化负载电流信号对比') legend('下ID1','下ID2','下ID3','下ID4') xlabel('时间/ns') ylabel('电压/V') axis([0 800 -1.5 1.5]) grid on hold on subplot(2,2,3) plot(time_EDS1,EDS1,time_EDS2,EDS2,time_EDS3,EDS3,time_EDS4,EDS4,time_EDS5,EDS5,time_EDS6,EDS6,time_EDS7,EDS7,time_EDS8,EDS8,time_EDS9,EDS9,'LineWidth',1.5) title('能谱仪测量结果') axis([-200 800 -1 30]) legend('能谱仪1','能谱仪2','能谱仪3','能谱仪4','能谱仪5','能谱仪6','能谱仪7','能谱仪8','能谱仪9') xlabel('时间/ns') ylabel('电压/V') grid on hold on set(gcf,'Position',get(0,'Screensize')); subplot(2,2,4) plot(time_PCD,PCD,time_PIN,PIN,time_guangdianguan1,guangdianguan1,time_guangdianguan2,guangdianguan2,'LineWidth',1.5) title('PCD、PIN及光电管测试结果') axis([-200 800 -40 40]) legend('PCD1','PCD2','光电管1','光电管2') xlabel('时间/ns') ylabel('电压/V') grid on for ii=1:1:10000 upperstackDaverage(ii)=[(upperstackD_1(ii)+upperstackD_2(ii)+upperstackD_3(ii))/3];% lowerstackDaverage(ii)=[(lowerstackD_1(ii)+lowerstackD_2(ii)+lowerstackD_3(ii))/3]; upperstackBaverage(ii)=[(upperstackB_1(ii)+upperstackB_2(ii)+upperstackB_3(ii))/3]; lowerstackBaverage(ii)=[(lowerstackB_1(ii)+lowerstackB_2(ii)+lowerstackB_3(ii))/3]; upperMITLaverage(ii)=(upperMITL_1(ii)+upperMITL_2(ii)+upperMITL_3(ii)+upperMITL_4(ii)+upperMITL_5(ii)+upperMITL_6(ii))/6; lowerMITLaverage(ii)=(lowerMITL_1(ii)+lowerMITL_2(ii)+lowerMITL_3(ii)+lowerMITL_4(ii)+lowerMITL_5(ii)+lowerMITL_6(ii))/6; lowerIDaverage(ii)=(lowerID1(ii)+lowerID3(ii))/2; end shangduigonglvliu=upperstackDaverage.*upperstackBaverage; xiaduigonglvliu=lowerstackDaverage.*lowerstackBaverage; figure(8) plot(time_upperstackD_1,shangduigonglvliu,time_lowerstackD_1,xiaduigonglvliu,'LineWidth',1.5) hold on title('堆栈功率') legend('上堆功率流','下堆功率流') grid on set(gcf,'Position',get(0,'Screensize')); axis([100 800 -400000 600000]) %---------------------------------------- %下层电流结果绘图 figure(9) plot(time_lowerstackB_1,lowerstackBaverage,time_lowerMITL_1,lowerMITLaverage,time_lowerID1,-lowerIDaverage); hold on title('下层电流') legend('下层绝缘堆栈','下层MITL','下层负载') grid on set(gcf,'Position',get(0,'Screensize')); axis([100 800 -600 1500]) %---------------------------------------- %检查数据是否存在全为0的列 all_zero_cols=all(file==0,1); if any(all_zero_cols) msg='存在全为零的列'; uiwait(warndlg(msg,'警告')); end %clear all %close all %clc end 将这个matlab函数 转化为c#

eds能谱图的元素选择