画图，matlab，
clear;close all;clc;tic;dirOutput = dir('*.dat'); % 罗列所有后缀-1.dat的文件列表，罗列BDDATA的数据
filenames = string({dirOutput.name}'); % 提取文件名%% 丢包统计
FILENAMES = [""];
LOSS_YTJ = [ ];
LOSS_RAD = [ ];
LOSS_ETH = [ ];
LOSS_YK = [ ];
j=1;%% 批处理
for i = 1:length(filenames)%% 文件名filename = char(filenames(i)); filename = string(filename([1:end-16])); % 例如：20090111-190130-2%% 加载数据文件TELEMETRY_fileID = fopen(filename+'TELEMETRY-v1.dat'); % 获取 TELEMETRY 文件ID DATA_TELEMETRY = fread(TELEMETRY_fileID); % 读取 TELEMETRY 文件数据fclose(TELEMETRY_fileID); % 关闭文件%% 解析数据文件TELEMETRY_RAW = reshape(DATA_TELEMETRY,[ 239 length(DATA_TELEMETRY)/239])'; % 帧格式转换cnt_yaoce_link = TELEMETRY_RAW(:,17); % 遥测链路cnt （239字节的第17字节）time_yitiji = (TELEMETRY_RAW(:,234)*2^8 + TELEMETRY_RAW(:,233)) ; % 一体机时间 (239字节的第234、233字节)cnt_gnd_eth = TELEMETRY_RAW(:,13); % 地面接口网口cnt (239字节的第13字节)ip_config = TELEMETRY_RAW(:,30); % 电台反馈IP (239字节的第30字节)cfk_ml = TELEMETRY_RAW(:,31); % 最近测发控指令 (239字节的第31字节)signal_power = TELEMETRY_RAW(:,32); % 遥控信号能量值单字节 (239字节的第32字节)zijian = TELEMETRY_RAW(:,33); % 自检完成 (239字节的第33字节)state_yaokong = TELEMETRY_RAW(:,32); % 遥控链路通断状态 (239字节的第32字节)control_yaokong = TELEMETRY_RAW(:,36)*2^8+TELEMETRY_RAW(:,35); % 地面遥控指令是否连续 (239字节的第36、35字节)angle_cabrage = (TELEMETRY_RAW(:,50)*2^8 + TELEMETRY_RAW(:,49)) * 0.01; % 飞行姿态 俯仰 （239字节第50、49）angle_roll = (TELEMETRY_RAW(:,52)*2^8 + TELEMETRY_RAW(:,51)) * 0.01; % 飞行姿态 滚动 （239字节第52、51）angle_Heading = (TELEMETRY_RAW(:,54)*2^8 + TELEMETRY_RAW(:,53)) * 0.01; % 飞行姿态 航向 （239字节第54、53）speed_bd = TELEMETRY_RAW(:,149)*0.005; % 靶弹飞行速度 （239字节第149字节）high_bd = (TELEMETRY_RAW(:,88)*2^8+TELEMETRY_RAW(:,87))*0.2-500; % 靶弹飞行高度 （239字节第88、87字节）len_tele = length(cnt_yaoce_link); % TELEMETRY 数据长度%% 数据处理% 遥测链路是否丢包，前后相减是否为1diubao_yaoce_link = diff(cnt_yaoce_link);diubao_yaoce_link(diubao_yaoce_link <= 0) = diubao_yaoce_link(diubao_yaoce_link <= 0) + 256; % 处理负数flag_diubao_yaoce_link = diubao_yaoce_link ~= 1; % 遥测链路丢包时刻count_diubao_yaoce_link = sum(diubao_yaoce_link-1); % 遥测链路丢包数% 遥测是否丢包，基于一体机，前后相减是否为1 diubao_yaoce_yitiji = diff(time_yitiji);
% diubao_yaoce_yitiji(diubao_yaoce_yitiji == -65535) = 1; % 处理负数diubao_yaoce_yitiji(diubao_yaoce_yitiji < 0) = diubao_yaoce_yitiji(diubao_yaoce_yitiji < 0) + 65536; % 处理负数flag_diubao_yaoce_yitiji = diubao_yaoce_yitiji ~= 1; % 遥测丢包时刻count_diubao_yaoce_yitiji = sum(diubao_yaoce_yitiji-1); % 遥测总丢包数% 地面网口接口是否丢包，前后相减是否为1diubao_gnd_eth = diff(cnt_gnd_eth);diubao_gnd_eth(diubao_gnd_eth <= 0) = diubao_gnd_eth(diubao_gnd_eth <= 0) + 256; % 处理负数flag_diubao_gnd_eth = diubao_gnd_eth ~= 1; % 地面网口丢包时刻count_diubao_gnd_eth = sum(diubao_gnd_eth-1); % 地面网口丢包数% 链路能量情况，整个字节显示% 自检完成，加载成功，判断是否等于FFstate_zijian = zijian == 0xff;% 遥控链路通断状态（大于0即可）state_yk_link = state_yaokong > 0;% 地面遥控指令是否连续（是否存在00）state_yk_zl = control_yaokong == 0;k=1;while state_yk_zl(k) == 1 && k < 500 % 清除错误的数值,容错前50个。state_yk_zl(k) = 0;k = k+1; endcount_state_yk_zl = sum(state_yk_zl); % 遥控总丢包数flag_state_yk_zl = state_yk_zl ~= 0;% 遥控指令切换状态，切换计数器% 角度负数处理angle_cabrage(angle_cabrage>327.67) = angle_cabrage(angle_cabrage>327.67) - 655.36; angle_roll(angle_roll>327.67) = angle_roll(angle_roll>327.67) - 655.36; angle_Heading(angle_Heading>327.67) = angle_Heading(angle_Heading>327.67) - 655.36; % 统计丢包数（网络、链路、遥测、遥控）FILENAMES(j) = filename; % 文件名称
% LOSS_YTJ(j) = count_diubao_yaoce_yitiji - count_diubao_yaoce_link; % 一体机丢包
% LOSS_RAD(j) = count_diubao_yaoce_link - count_diubao_gnd_eth; % 遥测链路丢包LOSS_YTJ(j) = count_diubao_yaoce_yitiji ; % 一体机丢包LOSS_RAD(j) = count_diubao_yaoce_link ; % 遥测链路丢包LOSS_ETH(j) = count_diubao_gnd_eth; % 以太网丢包LOSS_YK(j) = count_state_yk_zl; % 遥控丢包disp('如上图所示，网络丢'+string(LOSS_ETH(j))+'包，遥测链路丢'+string(LOSS_RAD(j))+'包，一体机丢'+string(LOSS_YTJ(j))+'包。'); % 快速报告j=j+1;%% 创建对应文件名的文件夹mkdir(filename);%% 绘制第一类图 ： 丢包 帧计数（网络、链路、遥测、指令） （ BD 指令和指令切换时刻）img = gcf;clf(img);hold on;plot(flag_diubao_gnd_eth*150,'--');plot(flag_diubao_yaoce_link*140,'.');plot(flag_diubao_yaoce_yitiji*130,'s');plot(cfk_ml,'r');plot([0,len_tele],[97,97],'linestyle','--');plot([0,len_tele],[113,113],'linestyle','--');legend('网络丢包','链路丢包','一体机丢包','指令','0x61大小功率','0x71数据装订','Location','NorthEastOutside');title('丢包与帧计数（网络、链路、遥测、指令）');ylim([-0.3 155]);hold off;print(img,'-dpng','-r800',filename+'/1_1丢包帧计数与指令.png');savefig(img,filename+'/1_1丢包帧计数与指令.fig');% clf(img);
% hold on;
% plot(state_yk_qie(:,2),'r');
% plot(yk_zl_qiehuan*150,'g--');
% plot([0,len_bd],[97,97],'linestyle','--');
% plot([0,len_bd],[113,113],'linestyle','--');
% legend('指令','指令切换时刻','0x61大小功率','0x71数据装订','Location','NorthEastOutside');
% ylim([-0.3 155]);
% hold off;
% print(img,'-dpng','-r800',filename+'/1_2靶弹指令和指令切换时刻.png');
% savefig(img,filename+'/1_2靶弹指令和指令切换时刻.fig');%% 绘制第二类图 遥控信号能量 & 遥控丢包 & 自检完成clf(img);hold on;plot(signal_power,'r');% 通断、其余未断链plot([0,len_tele],[0,0],'linestyle','--'); % 大信号plot([0,len_tele],[48,48],'linestyle','--');plot([0,len_tele],[95,95],'linestyle','--');% 中信号plot([0,len_tele],[96,96],'linestyle','--');plot([0,len_tele],[143,143],'linestyle','--');% 小信号plot([0,len_tele],[144,144],'linestyle','--');plot([0,len_tele],[191,191],'linestyle','--');% 遥控丢包plot(flag_state_yk_zl*150,'g');% 自检plot(state_zijian * 10); % 自检title('遥控信号能量 & 遥控丢包 & 自检完成');legend('遥控信号能量','通断','大信号上限','大信号下限', ...'中信号下限','中信号上限','小信号下限','小信号上限',...'遥控丢包','自检完成','Location','NorthEastOutside');hold off;print(img,'-dpng','-r800',filename+'/2_遥控丢包_遥控能量.png');savefig(img,filename+'/2_遥控丢包_遥控能量.fig');%% 绘制第三类图 飞行 （俯仰、滚动、航向）
% clf(img);
% hold on;
% plot(angle_cabrage);
% plot(angle_roll);
% plot(angle_Heading);
% title('俯仰、滚动、航向');
% ylim([-180 180]);
% hold off;
% print(img,'-dpng','-r800',filename+'/3_飞行角度.png');
% savefig(img,filename+'/3_飞行角度.fig'); %% 绘制第四类图 IP分析图
% clf(img);
% plot(ip_config);
% title('IP反馈图');
% print(img,'-dpng','-r800',filename+'/4_IP反馈图.png');
% savefig(img,filename+'/4_IP反馈图.fig'); %% 纯链路丢包clf(img);hold on;plot((flag_diubao_yaoce_link-flag_diubao_gnd_eth)*150,'--');plot(cfk_ml,'r');plot([0,len_tele],[97,97],'linestyle','--');plot([0,len_tele],[113,113],'linestyle','--');legend('链路丢包','指令','0x61大小功率','0x71数据装订','Location','NorthEastOutside');title('链路丢包与帧计数（链路、指令）');ylim([-0.3 155]);hold off;print(img,'-dpng','-r800',filename+'/5丢包帧计数与指令.png');savefig(img,filename+'/5纯链路丢包帧计数与指令.fig');%% 纯一体机丢包clf(img);hold on;plot((flag_diubao_yaoce_yitiji-flag_diubao_yaoce_link)*150,'--');plot(cfk_ml,'r');plot([0,len_tele],[97,97],'linestyle','--');plot([0,len_tele],[113,113],'linestyle','--');legend('一体机丢包','指令','0x61大小功率','0x71数据装订','Location','NorthEastOutside');title('一体机丢包与帧计数（链路、指令）');ylim([-0.3 155]);hold off;print(img,'-dpng','-r800',filename+'/6丢包帧计数与指令.png');savefig(img,filename+'/6纯一体丢包帧计数与指令.fig');%% 纯网络丢包clf(img);hold on;plot(flag_diubao_gnd_eth*150,'--');plot(cfk_ml,'r');plot([0,len_tele],[97,97],'linestyle','--');plot([0,len_tele],[113,113],'linestyle','--');legend('网络链路丢包','指令','0x61大小功率','0x71数据装订','Location','NorthEastOutside');title('网络丢包与帧计数（链路、指令）');ylim([-0.3 155]);hold off;print(img,'-dpng','-r800',filename+'/7丢包帧计数与指令.png');savefig(img,filename+'/7纯网络丢包帧计数与指令.fig');save(filename+'/8_data.mat', ...'cnt_yaoce_link','diubao_yaoce_link','flag_diubao_yaoce_link', ...'time_yitiji','diubao_yaoce_yitiji','flag_diubao_yaoce_yitiji', ...'cnt_gnd_eth','diubao_gnd_eth','flag_diubao_gnd_eth');
end%% 数据统计与导出
FILENAMES=FILENAMES';
LOSS_YTJ = LOSS_YTJ';
LOSS_RAD=LOSS_RAD';
LOSS_ETH=LOSS_ETH';
LOSS_YK=LOSS_YK';writetable(table(FILENAMES,LOSS_YTJ,LOSS_RAD,LOSS_ETH,LOSS_YK,'VariableNames',["FILENAMES","LOSS_YTJ","LOSS_RAD","LOSS_ETH","LOSS_YK"]),'丢包数据统计（一体机丢包-遥测链路丢包-地面网口丢包-遥控丢包）.csv');toc;
%% 清空
clear all;close all;clc;
tic; % 计时dirOutput = dir('*-1.dat'); % 罗列所有后缀-1.dat的文件列表，罗列BDDATA的数据
filenames = string({dirOutput.name}'); % 提取文件名%% 丢包统计
FILENAMES = [""];
LOSS_YTJ = [ ];
LOSS_RAD = [ ];
LOSS_ETH = [ ];
LOSS_YK = [ ];
j=1;%% 批处理
for i = 1:length(filenames)%% 文件名filename = char(filenames(i)); filename = string(filename([1:end-16])); % 例如：20090111-190130-2%% 加载数据文件BDDATA_fileID = fopen(filename+'-A2-BDDATA-1.dat'); % 获取 BDDATA 文件ID TELEMETRY_fileID = fopen(filename+'-C2-TELEMETRY-v1.dat'); % 获取 TELEMETRY 文件ID DATA_BDDATA = fread(BDDATA_fileID); % 读取 BDDATA 文件数据DATA_TELEMETRY = fread(TELEMETRY_fileID); % 读取 TELEMETRY 文件数据fclose(BDDATA_fileID); fclose(TELEMETRY_fileID); % 关闭文件%% 解析数据文件TELEMETRY_RAW = reshape(DATA_TELEMETRY,[ 239 length(DATA_TELEMETRY)/239])'; % 帧格式转换BDDATA_RAW = reshape(DATA_BDDATA ,[ 178 length(DATA_BDDATA )/178])'; % 帧格式转换cnt_yaoce_link = TELEMETRY_RAW(:,17); % 遥测链路cnt （239字节的第17字节）time_yitiji = (TELEMETRY_RAW(:,234)*2^8 + TELEMETRY_RAW(:,233)) ; % 一体机时间 (239字节的第234、233字节)cnt_gnd_eth = TELEMETRY_RAW(:,13); % 地面接口网口cnt (239字节的第13字节)ip_config = TELEMETRY_RAW(:,30); % 电台反馈IP (239字节的第30字节)cfk_ml = TELEMETRY_RAW(:,31); % 最近测发控指令 (239字节的第31字节)signal_power = TELEMETRY_RAW(:,32); % 遥控信号能量值单字节 (239字节的第32字节)zijian = TELEMETRY_RAW(:,33); % 自检完成 (239字节的第33字节)state_yaokong = TELEMETRY_RAW(:,32); % 遥控链路通断状态 (239字节的第32字节)control_yaokong = TELEMETRY_RAW(:,36)*2^8+TELEMETRY_RAW(:,35); % 地面遥控指令是否连续 (239字节的第36、35字节)angle_cabrage = (TELEMETRY_RAW(:,50)*2^8 + TELEMETRY_RAW(:,49)) * 0.01; % 飞行姿态 俯仰 （239字节第50、49）angle_roll = (TELEMETRY_RAW(:,52)*2^8 + TELEMETRY_RAW(:,51)) * 0.01; % 飞行姿态 滚动 （239字节第52、51）angle_Heading = (TELEMETRY_RAW(:,54)*2^8 + TELEMETRY_RAW(:,53)) * 0.01; % 飞行姿态 航向 （239字节第54、53）speed_bd = TELEMETRY_RAW(:,149)*0.005; % 靶弹飞行速度 （239字节第149字节）high_bd = (TELEMETRY_RAW(:,88)*2^8+TELEMETRY_RAW(:,87))*0.2-500; % 靶弹飞行高度 （239字节第88、87字节）state_yk_qie = BDDATA_RAW(:,[11 10]); % 遥控指令切换状态（阶梯状），BD_DATA里面第11（计数）、10（指令）字节。time_bd_yitiji = (BDDATA_RAW(:,24)*2^24+BDDATA_RAW(:,23)*2^16+BDDATA_RAW(:,22)*2^8+BDDATA_RAW(:,21))*0.00001;len_tele = length(cnt_yaoce_link); % TELEMETRY 数据长度len_bd = length(state_yk_qie); % BDDATA 数据长度%% 数据处理% 遥测链路是否丢包，前后相减是否为1diubao_yaoce_link = diff(cnt_yaoce_link);diubao_yaoce_link(diubao_yaoce_link <= 0) = diubao_yaoce_link(diubao_yaoce_link <= 0) + 256; % 处理负数flag_diubao_yaoce_link = diubao_yaoce_link ~= 1; % 遥测链路丢包时刻count_diubao_yaoce_link = sum(diubao_yaoce_link-1); % 遥测链路丢包数% 遥测是否丢包，基于一体机，前后相减是否为1 diubao_yaoce_yitiji = diff(time_yitiji);diubao_yaoce_yitiji(diubao_yaoce_yitiji < 0) = diubao_yaoce_yitiji(diubao_yaoce_yitiji < 0) + 65536; % 处理负数flag_diubao_yaoce_yitiji = diubao_yaoce_yitiji ~= 1; % 遥测丢包时刻count_diubao_yaoce_yitiji = sum(diubao_yaoce_yitiji-1); % 遥测总丢包数% 地面网口接口是否丢包，前后相减是否为1diubao_gnd_eth = diff(cnt_gnd_eth);diubao_gnd_eth(diubao_gnd_eth <= 0) = diubao_gnd_eth(diubao_gnd_eth <= 0) + 256; % 处理负数flag_diubao_gnd_eth = diubao_gnd_eth ~= 1; % 地面网口丢包时刻count_diubao_gnd_eth = sum(diubao_gnd_eth-1); % 地面网口丢包数% 链路能量情况，整个字节显示% 自检完成，加载成功，判断是否等于FFstate_zijian = zijian == 0xff;% 遥控链路通断状态（大于0即可）state_yk_link = state_yaokong > 0;% 地面遥控指令是否连续（是否存在00）state_yk_zl = control_yaokong == 0;k=1;while state_yk_zl(k) == 1 && k < 500 % 清除错误的数值,容错前50个。state_yk_zl(k) = 0;k = k+1; endcount_state_yk_zl = sum(state_yk_zl); % 遥控总丢包数flag_state_yk_zl = state_yk_zl ~= 0;% 遥控指令切换状态，切换计数器yk_zl_qiehuan = diff(state_yk_qie(:,1));% 角度负数处理angle_cabrage(angle_cabrage>327.67) = angle_cabrage(angle_cabrage>327.67) - 655.36; angle_roll(angle_roll>327.67) = angle_roll(angle_roll>327.67) - 655.36; angle_Heading(angle_Heading>327.67) = angle_Heading(angle_Heading>327.67) - 655.36; % 统计丢包数（网络、链路、遥测、遥控）FILENAMES(j) = filename; % 文件名称
% LOSS_YTJ(j) = count_diubao_yaoce_yitiji - count_diubao_yaoce_link; % 一体机丢包
% LOSS_RAD(j) = count_diubao_yaoce_link - count_diubao_gnd_eth; % 遥测链路丢包LOSS_YTJ(j) = count_diubao_yaoce_yitiji ; % 一体机丢包LOSS_RAD(j) = count_diubao_yaoce_link ; % 遥测链路丢包LOSS_ETH(j) = count_diubao_gnd_eth; % 以太网丢包LOSS_YK(j) = count_state_yk_zl; % 遥控丢包disp('如上图所示，网络丢'+string(LOSS_ETH(j))+'包，遥测链路丢'+string(LOSS_RAD(j))+'包，一体机丢'+string(LOSS_YTJ(j))+'包。'); % 快速报告j=j+1;%% 创建对应文件名的文件夹mkdir(filename);%% 绘制第一类图 ： 丢包 帧计数（网络、链路、遥测、指令） （ BD 指令和指令切换时刻）img = gcf;clf(img);hold on;plot(flag_diubao_gnd_eth*150,'--');plot(flag_diubao_yaoce_link*140,'.');plot(flag_diubao_yaoce_yitiji*130,'s');plot(cfk_ml,'r');plot([0,len_tele],[97,97],'linestyle','--');plot([0,len_tele],[113,113],'linestyle','--');legend('网络丢包','链路丢包','一体机丢包','指令','0x61大小功率','0x71数据装订','Location','NorthEastOutside');title('丢包与帧计数（网络、链路、遥测、指令）');ylim([-0.3 155]);hold off;print(img,'-dpng','-r800',filename+'/1_1丢包帧计数与指令.png');savefig(img,filename+'/1_1丢包帧计数与指令.fig');clf(img);hold on;plot(state_yk_qie(:,2),'r');plot(yk_zl_qiehuan*150,'g--');plot([0,len_bd],[97,97],'linestyle','--');plot([0,len_bd],[113,113],'linestyle','--');legend('指令','指令切换时刻','0x61大小功率','0x71数据装订','Location','NorthEastOutside');ylim([-0.3 155]);hold off;print(img,'-dpng','-r800',filename+'/1_2靶弹指令和指令切换时刻.png');savefig(img,filename+'/1_2靶弹指令和指令切换时刻.fig');%% 绘制第二类图 遥控信号能量 & 遥控丢包 & 自检完成clf(img);hold on;plot(signal_power,'r');% 通断、其余未断链plot([0,len_tele],[0,0],'linestyle','--'); % 大信号plot([0,len_tele],[48,48],'linestyle','--');plot([0,len_tele],[95,95],'linestyle','--');% 中信号plot([0,len_tele],[96,96],'linestyle','--');plot([0,len_tele],[143,143],'linestyle','--');% 小信号plot([0,len_tele],[144,144],'linestyle','--');plot([0,len_tele],[191,191],'linestyle','--');% 遥控丢包plot(flag_state_yk_zl*150,'g');% 自检plot(state_zijian * 10); % 自检title('遥控信号能量 & 遥控丢包 & 自检完成');legend('遥控信号能量','通断','大信号上限','大信号下限', ...'中信号下限','中信号上限','小信号下限','小信号上限',...'遥控丢包','自检完成','Location','NorthEastOutside');hold off;print(img,'-dpng','-r800',filename+'/2_遥控丢包_遥控能量.png');savefig(img,filename+'/2_遥控丢包_遥控能量.fig');%% 绘制第三类图 飞行 （俯仰、滚动、航向）
% clf(img);
% hold on;
% plot(angle_cabrage);
% plot(angle_roll);
% plot(angle_Heading);
% title('俯仰、滚动、航向');
% ylim([-180 180]);
% hold off;
% print(img,'-dpng','-r800',filename+'/3_飞行角度.png');
% savefig(img,filename+'/3_飞行角度.fig'); %% 绘制第四类图 IP分析图clf(img);plot(ip_config);title('IP反馈图');print(img,'-dpng','-r800',filename+'/4_IP反馈图.png');savefig(img,filename+'/4_IP反馈图.fig'); %% 纯链路丢包clf(img);hold on;plot((flag_diubao_yaoce_link-flag_diubao_gnd_eth)*150,'--');plot(cfk_ml,'r');plot([0,len_tele],[97,97],'linestyle','--');plot([0,len_tele],[113,113],'linestyle','--');legend('链路丢包','指令','0x61大小功率','0x71数据装订','Location','NorthEastOutside');title('链路丢包与帧计数（链路、指令）');ylim([-0.3 155]);hold off;print(img,'-dpng','-r800',filename+'/5丢包帧计数与指令.png');savefig(img,filename+'/5纯链路丢包帧计数与指令.fig');%% 纯一体机丢包clf(img);hold on;plot((flag_diubao_yaoce_yitiji-flag_diubao_yaoce_link)*150,'--');plot(cfk_ml,'r');plot([0,len_tele],[97,97],'linestyle','--');plot([0,len_tele],[113,113],'linestyle','--');legend('一体机丢包','指令','0x61大小功率','0x71数据装订','Location','NorthEastOutside');title('一体机丢包与帧计数（链路、指令）');ylim([-0.3 155]);hold off;print(img,'-dpng','-r800',filename+'/6丢包帧计数与指令.png');savefig(img,filename+'/6纯一体丢包帧计数与指令.fig');%% 纯网络丢包clf(img);hold on;plot(flag_diubao_gnd_eth*150,'--');plot(cfk_ml,'r');plot([0,len_tele],[97,97],'linestyle','--');plot([0,len_tele],[113,113],'linestyle','--');legend('网络链路丢包','指令','0x61大小功率','0x71数据装订','Location','NorthEastOutside');title('网络丢包与帧计数（链路、指令）');ylim([-0.3 155]);hold off;print(img,'-dpng','-r800',filename+'/7丢包帧计数与指令.png');savefig(img,filename+'/7纯网络丢包帧计数与指令.fig');save(filename+'/8_data.mat', ...'cnt_yaoce_link','diubao_yaoce_link','flag_diubao_yaoce_link', ...'time_yitiji','diubao_yaoce_yitiji','flag_diubao_yaoce_yitiji', ...'cnt_gnd_eth','diubao_gnd_eth','flag_diubao_gnd_eth');
end%% 数据统计与导出
FILENAMES=FILENAMES';
LOSS_YTJ = LOSS_YTJ';
LOSS_RAD=LOSS_RAD';
LOSS_ETH=LOSS_ETH';
LOSS_YK=LOSS_YK';writetable(table(FILENAMES,LOSS_YTJ,LOSS_RAD,LOSS_ETH,LOSS_YK,'VariableNames',["FILENAMES","LOSS_YTJ","LOSS_RAD","LOSS_ETH","LOSS_YK"]),'丢包数据统计（一体机丢包-遥测链路丢包-地面网口丢包-遥控丢包）.csv');toc;
ALL
% 计时
tic; %% 丢包统计
FILENAMES = [""];
LOSS_YTJ = [ ];
LOSS_RAD = [ ];
LOSS_ETH = [ ];
LOSS_YK = [ ];%% 加载数据文件
BD_fileID = fopen('bd.dat'); % 获取 靶弹 文件ID
ZHU_fileID = fopen('zhu.dat'); % 获取 主站 文件ID
FU_fileID = fopen('fu.dat'); % 获取 副站 文件ID DATA_BD = fread(BD_fileID); % 读取 靶弹 文件数据
DATA_ZHU = fread(ZHU_fileID); % 读取 主站 文件数据
DATA_FU = fread(FU_fileID); % 读取 副站 文件数据fclose(BD_fileID); fclose(ZHU_fileID); fclose(FU_fileID); % 关闭文件%% 解析数据文件
ZHU_RAW = reshape(DATA_ZHU,[239 length(DATA_ZHU)/239])'; % 主站帧格式转换
FU_RAW = reshape(DATA_FU,[ 239 length(DATA_FU )/239])'; % 副站帧格式转换
BDDATA_RAW = reshape(DATA_BD,[ 178 length(DATA_BD )/178])'; % 靶弹帧格式转换%% 靶弹数据解析
state_yk_qie = BDDATA_RAW(:,[11 10]); % 遥控指令切换状态（阶梯状），BD_DATA里面第11（计数）、10（指令）字节。
time_bd_yitiji = (BDDATA_RAW(:,24)*2^24 + BDDATA_RAW(:,23)*2^16 + ...BDDATA_RAW(:,22)*2^8 + BDDATA_RAW(:,21))*0.00001;
BD_UTC_raw = BDDATA_RAW(:,28)*2^24 + BDDATA_RAW(:,27)*2^16 + ...BDDATA_RAW(:,26)*2^8 + BDDATA_RAW(:,25); % UTC 时间原码 (BD_DATA里面第25、26、27、28字节)
BD_UTC = ConvertNum2UTC(BD_UTC_raw); % BD 的 UTC 时间
len_bd = length(state_yk_qie); % BDDATA 数据长度%% 主站数据解析
zhu_cnt_yaoce_link = ZHU_RAW(:,17); % 遥测链路cnt （239字节的第17字节）
zhu_time_yitiji = (ZHU_RAW(:,234)*2^8 + ZHU_RAW(:,233)) ; % 一体机时间 (239字节的第234、233字节)
zhu_sky_control_time= (ZHU_RAW(:,110)*2^24 + ZHU_RAW(:,109)*2^16 + ...ZHU_RAW(:,108)*2^8 + ZHU_RAW(:,107)); % 飞控时间 （239字节第110、109、108、107）
zhu_cnt_gnd_eth = ZHU_RAW(:,13); % 地面接口网口cnt (239字节的第13字节)
zhu_ip_config = ZHU_RAW(:,30); % 电台反馈IP (239字节的第30字节)
zhu_cfk_ml = ZHU_RAW(:,31); % 最近测发控指令 (239字节的第31字节)
zhu_signal_power = ZHU_RAW(:,32); % 遥控信号能量值单字节 (239字节的第32字节)
zhu_zijian = ZHU_RAW(:,33); % 自检完成 (239字节的第33字节)
zhu_state_yaokong = ZHU_RAW(:,32); % 遥控链路通断状态 (239字节的第32字节)
zhu_control_yaokong = ZHU_RAW(:,36)*2^8 + ZHU_RAW(:,35); % 地面遥控指令是否连续 (239字节的第36、35字节)
zhu_angle_cabrage = (ZHU_RAW(:,50)*2^8 + ZHU_RAW(:,49)) * 0.01; % 飞行姿态 俯仰 （239字节第50、49）
zhu_angle_roll = (ZHU_RAW(:,52)*2^8 + ZHU_RAW(:,51)) * 0.01; % 飞行姿态 滚动 （239字节第52、51）
zhu_angle_Heading = (ZHU_RAW(:,54)*2^8 + ZHU_RAW(:,53)) * 0.01; % 飞行姿态 航向 （239字节第54、53）
zhu_longitude = (ZHU_RAW(:,83)*2^16 + ZHU_RAW(:,82)*2^8 + ...ZHU_RAW(:,81)); % 经度 （239字节第83、82、81）原始数据
zhu_latitude = (ZHU_RAW(:,86)*2^16 + ZHU_RAW(:,85)*2^8 + ...ZHU_RAW(:,84)); % 纬度 （239字节第86、85、84）原始数据
zhu_positiong_high = (ZHU_RAW(:,88)*2^8 + ZHU_RAW(:,87)) * 0.2 - 500; % 高度 (239字节的第88、87字节)
zhu_speed_bd = ZHU_RAW(:,149)*0.005; % 靶弹飞行速度 （239字节第149字节）
zhu_high_bd = (ZHU_RAW(:,88)*2^8 + ZHU_RAW(:,87))*0.2 - 500; % 靶弹飞行高度 （239字节第88、87字节）
zhu_UTC_raw =(ZHU_RAW(:,104)*2^24 + ZHU_RAW(:,103)*2^16 + ...ZHU_RAW(:,102)*2^8 + ZHU_RAW(:,101)); % UTC 时间原码 （239字节第104、103、102、101）
zhu_UTC = ConvertNum2UTC(zhu_UTC_raw); % UTC 时间zhu_len_tele = length(zhu_cnt_yaoce_link); % 主站数据包长度
%% 副站数据解析
fu_cnt_yaoce_link = FU_RAW(:,17); % 遥测链路cnt （239字节的第17字节）
fu_time_yitiji = (FU_RAW(:,234)*2^8 + FU_RAW(:,233)) ; % 一体机时间 (239字节的第234、233字节)
fu_sky_control_time= (FU_RAW(:,110)*2^24 + FU_RAW(:,109)*2^16 + ...FU_RAW(:,108)*2^8 + FU_RAW(:,107)); % 飞控时间 （239字节第110、109、108、107）
fu_cnt_gnd_eth = FU_RAW(:,13); % 地面接口网口cnt (239字节的第13字节)
fu_ip_config = FU_RAW(:,30); % 电台反馈IP (239字节的第30字节)
fu_cfk_ml = FU_RAW(:,31); % 最近测发控指令 (239字节的第31字节)
fu_signal_power = FU_RAW(:,32); % 遥控信号能量值单字节 (239字节的第32字节)
fu_zijian = FU_RAW(:,33); % 自检完成 (239字节的第33字节)
fu_state_yaokong = FU_RAW(:,32); % 遥控链路通断状态 (239字节的第32字节)
fu_control_yaokong = FU_RAW(:,36)*2^8 + FU_RAW(:,35); % 地面遥控指令是否连续 (239字节的第36、35字节)
fu_angle_cabrage = (FU_RAW(:,50)*2^8 + FU_RAW(:,49)) * 0.01; % 飞行姿态 俯仰 （239字节第50、49）
fu_angle_roll = (FU_RAW(:,52)*2^8 + FU_RAW(:,51)) * 0.01; % 飞行姿态 滚动 （239字节第52、51）
fu_angle_Heading = (FU_RAW(:,54)*2^8 + FU_RAW(:,53)) * 0.01; % 飞行姿态 航向 （239字节第54、53）
fu_longitude = (FU_RAW(:,83)*2^16 + FU_RAW(:,82)*2^8 + ...FU_RAW(:,81)); % 经度 （239字节第83、82、81）原始数据
fu_latitude = (FU_RAW(:,86)*2^16 + FU_RAW(:,85)*2^8 + ...FU_RAW(:,84)); % 纬度 （239字节第86、85、84）原始数据
fu_positiong_high = (FU_RAW(:,88)*2^8 + FU_RAW(:,87)) * 0.2 - 500; % 高度 (239字节的第88、87字节)
fu_speed_bd = FU_RAW(:,149)*0.005; % 靶弹飞行速度 （239字节第149字节）
fu_high_bd = (FU_RAW(:,88)*2^8 + FU_RAW(:,87))*0.2-500; % 靶弹飞行高度 （239字节第88、87字节）
fu_UTC_raw =(FU_RAW(:,104)*2^24 + FU_RAW(:,103)*2^16 + ...FU_RAW(:,102)*2^8 + FU_RAW(:,101)); % UTC 时间原码 （239字节第104、103、102、101）
fu_UTC = ConvertNum2UTC(fu_UTC_raw); % UTC 时间
fu_len_tele = length(fu_cnt_yaoce_link); % 副站数据包长度
%% 主站数据处理% 遥测链路是否丢包，前后相减是否为1zhu_diubao_yaoce_link = diff(zhu_cnt_yaoce_link);zhu_diubao_yaoce_link(zhu_diubao_yaoce_link <= 0) = zhu_diubao_yaoce_link(zhu_diubao_yaoce_link <= 0) + 256; % 处理负数zhu_flag_diubao_yaoce_link = zhu_diubao_yaoce_link ~= 1; % 遥测链路丢包时刻zhu_count_diubao_yaoce_link = sum(zhu_diubao_yaoce_link-1); % 遥测链路丢包数% 遥测是否丢包，基于一体机，前后相减是否为1 zhu_diubao_yaoce_yitiji = diff(zhu_time_yitiji);zhu_diubao_yaoce_yitiji(zhu_diubao_yaoce_yitiji < 0) = zhu_diubao_yaoce_yitiji(zhu_diubao_yaoce_yitiji < 0) + 65536; % 处理负数zhu_flag_diubao_yaoce_yitiji = zhu_diubao_yaoce_yitiji ~= 1; % 遥测丢包时刻zhu_count_diubao_yaoce_yitiji = sum(zhu_diubao_yaoce_yitiji-1); % 遥测总丢包数% 地面网口接口是否丢包，前后相减是否为1zhu_diubao_gnd_eth = diff(zhu_cnt_gnd_eth);zhu_diubao_gnd_eth(zhu_diubao_gnd_eth <= 0) = zhu_diubao_gnd_eth(zhu_diubao_gnd_eth <= 0) + 256; % 处理负数zhu_flag_diubao_gnd_eth = zhu_diubao_gnd_eth ~= 1; % 地面网口丢包时刻zhu_count_diubao_gnd_eth = sum(zhu_diubao_gnd_eth-1); % 地面网口丢包数% 遥控链路通断状态（大于0即可）zhu_state_yk_link = zhu_state_yaokong > 0;% 地面遥控指令是否连续（是否存在00）zhu_state_yk_zl = zhu_control_yaokong == 0;k=1;while zhu_state_yk_zl(k) == 1 && k < 500 % 清除错误的数值,容错前50个。zhu_state_yk_zl(k) = 0;k = k+1; endzhu_count_state_yk_zl = sum(zhu_state_yk_zl); % 遥控总丢包数zhu_flag_state_yk_zl = zhu_state_yk_zl ~= 0;% 遥控指令切换状态，切换计数器yk_zl_qiehuan = diff(state_yk_qie(:,1));% 角度负数处理zhu_angle_cabrage(zhu_angle_cabrage>327.67) = zhu_angle_cabrage(zhu_angle_cabrage>327.67) - 655.36; zhu_angle_roll(zhu_angle_roll>327.67) = zhu_angle_roll(zhu_angle_roll>327.67) - 655.36; zhu_angle_Heading(zhu_angle_Heading>327.67) = zhu_angle_Heading(zhu_angle_Heading>327.67) - 655.36; % 经纬度负数处理 比例尺换算zhu_longitude(zhu_longitude>8388607) = zhu_longitude(zhu_longitude>8388607) - 16777216;zhu_latitude(zhu_latitude>8388607) = zhu_latitude(zhu_latitude>8388607) - 16777216;zhu_longitude = zhu_longitude * 0.00005;zhu_latitude = zhu_latitude * 0.00005;% 轨迹与丢包图zhu_remove_eth = zhu_flag_diubao_yaoce_link-zhu_flag_diubao_gnd_eth;zhu_longitude_err = zhu_longitude(zhu_remove_eth==1);zhu_latitude_err = zhu_latitude(zhu_remove_eth==1);zhu_high_err = zhu_positiong_high(zhu_remove_eth==1);zhu_longitude_right = zhu_longitude(zhu_remove_eth==0);zhu_latitude_right = zhu_latitude(zhu_remove_eth==0);zhu_high_right = zhu_positiong_high(zhu_remove_eth==0);% 统计丢包数（网络、链路、遥测、遥控）FILENAMES(1) = 'main'; % 文件名称LOSS_YTJ(1) = zhu_count_diubao_yaoce_yitiji ; % 一体机丢包LOSS_RAD(1) = zhu_count_diubao_yaoce_link ; % 遥测链路丢包LOSS_ETH(1) = zhu_count_diubao_gnd_eth; % 以太网丢包LOSS_YK(1) = zhu_count_state_yk_zl; % 遥控丢包disp('如上图所示，网络丢'+string(LOSS_ETH(1))+'包，遥测链路丢'+ ...string(LOSS_RAD(1))+'包，一体机丢'+string(LOSS_YTJ(1))+ ...'包,遥控丢'+string(LOSS_YK(1))+'包。'); % 快速报告%% 副站数据处理% 遥测链路是否丢包，前后相减是否为1fu_diubao_yaoce_link = diff(fu_cnt_yaoce_link);fu_diubao_yaoce_link(fu_diubao_yaoce_link <= 0) = fu_diubao_yaoce_link(fu_diubao_yaoce_link <= 0) + 256; % 处理负数fu_flag_diubao_yaoce_link = fu_diubao_yaoce_link ~= 1; % 遥测链路丢包时刻fu_count_diubao_yaoce_link = sum(fu_diubao_yaoce_link-1); % 遥测链路丢包数% 遥测是否丢包，基于一体机，前后相减是否为1 fu_diubao_yaoce_yitiji = diff(fu_time_yitiji);fu_diubao_yaoce_yitiji(fu_diubao_yaoce_yitiji < 0) = fu_diubao_yaoce_yitiji(fu_diubao_yaoce_yitiji < 0) + 65536; % 处理负数fu_flag_diubao_yaoce_yitiji = fu_diubao_yaoce_yitiji ~= 1; % 遥测丢包时刻fu_count_diubao_yaoce_yitiji = sum(fu_diubao_yaoce_yitiji-1); % 遥测总丢包数% 地面网口接口是否丢包，前后相减是否为1fu_diubao_gnd_eth = diff(fu_cnt_gnd_eth);fu_diubao_gnd_eth(fu_diubao_gnd_eth <= 0) = fu_diubao_gnd_eth(fu_diubao_gnd_eth <= 0) + 256; % 处理负数fu_flag_diubao_gnd_eth = fu_diubao_gnd_eth ~= 1; % 地面网口丢包时刻fu_count_diubao_gnd_eth = sum(fu_diubao_gnd_eth-1); % 地面网口丢包数% 遥控链路通断状态（大于0即可）fu_state_yk_link = fu_state_yaokong > 0;% 地面遥控指令是否连续（是否存在00）fu_state_yk_zl = fu_control_yaokong == 0;k=1;while fu_state_yk_zl(k) == 1 && k < 500 % 清除错误的数值,容错前50个。fu_state_yk_zl(k) = 0;k = k+1; endfu_count_state_yk_zl = sum(fu_state_yk_zl); % 遥控总丢包数fu_flag_state_yk_zl = fu_state_yk_zl ~= 0;% 角度负数处理fu_angle_cabrage(fu_angle_cabrage>327.67) = fu_angle_cabrage(fu_angle_cabrage>327.67) - 655.36; fu_angle_roll(fu_angle_roll>327.67) = fu_angle_roll(fu_angle_roll>327.67) - 655.36; fu_angle_Heading(fu_angle_Heading>327.67) = fu_angle_Heading(fu_angle_Heading>327.67) - 655.36; % 经纬度负数处理 比例尺换算fu_longitude(fu_longitude>8388607) = fu_longitude(fu_longitude>8388607) - 16777216;fu_latitude(fu_latitude>8388607) = fu_latitude(fu_latitude>8388607) - 16777216;fu_longitude = fu_longitude * 0.00005;fu_latitude = fu_latitude * 0.00005;% 轨迹与丢包图fu_remove_eth = fu_flag_diubao_yaoce_link-fu_flag_diubao_gnd_eth;fu_longitude_err = fu_longitude(fu_remove_eth==1);fu_latitude_err = fu_latitude(fu_remove_eth==1);fu_high_err = fu_positiong_high(fu_remove_eth==1);fu_longitude_right = fu_longitude(fu_remove_eth==0);fu_latitude_right = fu_latitude(fu_remove_eth==0);fu_high_right = fu_positiong_high(fu_remove_eth==0);% 统计丢包数（网络、链路、遥测、遥控）FILENAMES(2) = 'slave'; % 文件名称LOSS_YTJ(2) = fu_count_diubao_yaoce_yitiji ; % 一体机丢包LOSS_RAD(2) = fu_count_diubao_yaoce_link ; % 遥测链路丢包LOSS_ETH(2) = fu_count_diubao_gnd_eth; % 以太网丢包LOSS_YK(2) = fu_count_state_yk_zl; % 遥控丢包disp('如上图所示，网络丢'+string(LOSS_ETH(2))+'包，遥测链路丢'+ ...string(LOSS_RAD(2))+'包，一体机丢'+string(LOSS_YTJ(2))+ ...'包,遥控丢'+string(LOSS_YK(2))+'包。'); % 快速报告[Angle_Antenna,Length,fuyang_test] = ...Antenna_direction(109.621197,40.847384,fu_longitude,fu_latitude,...fu_positiong_high,fu_angle_Heading); % 副站 仰角，距离，[Angle_Antenna_zhu,Length_zhu,fuyang_test_zhu] = ...Antenna_direction(109.621197,40.847384,zhu_longitude,zhu_latitude,...zhu_positiong_high,zhu_angle_Heading); % 主站 仰角，距离
%% 联合数据分析
% 遥控信号能量 & 主站遥控丢包
img = gcf;
clf(img);
hold on;
plot(zhu_signal_power,'r');plot([0,zhu_len_tele],[48,48],'linestyle','--'); % 大信号
plot([0,zhu_len_tele],[95,95],'linestyle','--'); % 大信号
plot([0,zhu_len_tele],[96,96],'linestyle','--'); % 中信号
plot([0,zhu_len_tele],[143,143],'linestyle','--'); % 中信号
plot([0,zhu_len_tele],[144,144],'linestyle','--'); % 小信号
plot([0,zhu_len_tele],[191,191],'linestyle','--'); % 小信号
plot(zhu_flag_state_yk_zl*150,'g'); % 遥控丢包title('遥控信号能量 & 主站遥控丢包');
legend('遥控信号能量','大信号上限','大信号下限', ...'中信号下限','中信号上限','小信号下限','小信号上限',...'遥控丢包','Location','NorthEastOutside');
hold off;
savefig(img,'1_遥控丢包_遥控能量zhu.fig');% 遥控信号能量 & 副站遥控丢包
img = gcf;
clf(img);
hold on;
plot(fu_signal_power,'r');plot([0,fu_len_tele],[48,48],'linestyle','--'); % 大信号
plot([0,fu_len_tele],[95,95],'linestyle','--'); % 大信号
plot([0,fu_len_tele],[96,96],'linestyle','--'); % 中信号
plot([0,fu_len_tele],[143,143],'linestyle','--'); % 中信号
plot([0,fu_len_tele],[144,144],'linestyle','--'); % 小信号
plot([0,fu_len_tele],[191,191],'linestyle','--'); % 小信号
plot(fu_flag_state_yk_zl*150,'g'); % 遥控丢包title('遥控信号能量 & 副站遥控丢包');
legend('遥控信号能量','大信号上限','大信号下限', ...'中信号下限','中信号上限','小信号下限','小信号上限',...'遥控丢包','Location','NorthEastOutside');
hold off;
savefig(img,'1_遥控丢包_遥控能量fu.fig');% 主站遥控丢包与飞行轨迹
clf(img);
subplot(1,2,1);
scatter3(zhu_longitude,zhu_latitude,zhu_positiong_high,10,"blue");
hold on;
scatter3( zhu_longitude(zhu_flag_state_yk_zl==1), ...zhu_latitude(zhu_flag_state_yk_zl==1), ...zhu_positiong_high(zhu_flag_state_yk_zl==1),100,"red");
ylim([40.75 40.95]);
xlim([109 109.7]);
zlim([1200 2000]);
xlabel("longitude");ylabel("latitude");zlabel("high");
title('靶弹飞行轨迹图 & 主站遥控丢包');
legend('飞行轨迹','遥控丢包','Location','NorthEastOutside');
% savefig(img,'1_2_轨迹与遥控丢包主站.fig');
% hold off;% 副站遥控丢包 与 飞行轨迹
% clf(img);
subplot(1,2,2);
scatter3(fu_longitude,fu_latitude,fu_positiong_high,10,"blue");
hold on;
scatter3( fu_longitude(fu_flag_state_yk_zl==1), ...fu_latitude(fu_flag_state_yk_zl==1), ...fu_positiong_high(fu_flag_state_yk_zl==1),100,"red");
ylim([40.75 40.95]);
xlim([109 109.7]);
zlim([1200 2000]);
xlabel("longitude");ylabel("latitude");zlabel("high");
title('靶弹飞行轨迹图 & 副站遥控丢包');
legend('飞行轨迹','遥控丢包','Location','NorthEastOutside');
savefig(img,'1_2_轨迹与遥控丢包主副站比较图.fig');
hold off;% 飞行（俯仰、滚动、航向）& 丢包时刻
% 飞行（经纬度、高度）& 丢包时刻
clf(img);
hold on;
plot(zhu_angle_cabrage);
plot(zhu_angle_roll);
plot(zhu_angle_Heading);
plot(zhu_flag_diubao_yaoce_link*150,'g--');
legend('俯仰角','滚动角','航向角','遥测丢包','Location','NorthEastOutside');
xlabel('时刻');ylabel("角度值");
title('主站：俯仰、滚动、航向 & 丢包图');
ylim([-180 180]);
hold off;
% print(img,'-dpng','-r800','2_1_飞行角度zhu.png');
savefig(img,'2_1_飞行角度zhu.fig'); clf(img);
hold off;
scatter3(zhu_longitude_right,zhu_latitude_right,zhu_high_right,10,"blue");
hold on;
scatter3(zhu_longitude_err,zhu_latitude_err,zhu_high_err,50,"red");
xlabel("longitude");ylabel("latitude");zlabel("high");
ylim([40.75 40.95]);
xlim([109 109.7]);
zlim([1200 2000]);
title('主站：飞行轨迹 & 丢包图');
hold off;
savefig(img,'2_2_轨迹与丢包zhu.fig'); clf(img);
hold on;
plot(fu_angle_cabrage);
plot(fu_angle_roll);
plot(fu_angle_Heading);
plot(fu_flag_diubao_yaoce_link*150,'g--');
legend('俯仰角','滚动角','航向角','遥测丢包','Location','NorthEastOutside');
title('副站：俯仰、滚动、航向 & 丢包图');
ylim([-180 180]);
hold off;
% print(img,'-dpng','-r800','2_1_飞行角度fu.png');
savefig(img,'2_1_飞行角度fu.fig'); clf(img);
hold off;
scatter3(fu_longitude_right,fu_latitude_right,fu_high_right,10,"blue");
hold on;
scatter3(fu_longitude_err,fu_latitude_err,fu_high_err,50,"red");
ylim([40.75 40.95]);
xlim([109 109.7]);
zlim([1200 2000]);
xlabel("longitude");ylabel("latitude");zlabel("high");
title('副站：飞行轨迹 & 丢包图');
hold off;
savefig(img,'2_2_轨迹与丢包fu.fig'); clf(img);
hold off;
scatter3(fu_longitude_right,fu_latitude_right,fu_high_right,1,"black");
hold on;
scatter3(fu_longitude_err,fu_latitude_err,fu_high_err,200,"red");
scatter3(zhu_longitude_err,zhu_latitude_err,zhu_high_err,50,"blue");
ylim([40.75 40.95]);
xlim([109 109.7]);
zlim([1200 2000]);
xlabel("longitude");ylabel("latitude");zlabel("high");
title('主副站：俯仰、滚动、航向 & 丢包图');
legend('轨迹图','副站丢包','主站丢包','Location','NorthEastOutside');
hold off;
savefig(img,'2_3_主副站轨迹与丢包.fig'); % 副站 丢包时刻的三个角度信息
fu_diubao_index = find(fu_remove_eth==1);
fu_diubao_pitch = [];
fu_diubao_yaw = [];
fu_diubao_roll = [];
fu_diubao_angle = [];
fu_diubao_jin = [];
fu_diubao_wei = [];
fu_diubao_gao = [];
j=1;
for i = fu_diubao_index'fu_diubao_pitch(j) = fu_angle_cabrage(i);fu_diubao_yaw(j) = fu_angle_Heading(i);fu_diubao_roll(j) = fu_angle_roll(i);fu_diubao_angle(j) = Angle_Antenna(i); fu_diubao_jin(j) = fu_longitude(i);fu_diubao_wei(j) = fu_latitude(i);fu_diubao_gao(j) = fu_positiong_high(i);j=j+1;
end% 主站 丢包时刻的三个角度信息
zhu_diubao_index = find(zhu_remove_eth==1);
zhu_diubao_pitch = [];
zhu_diubao_yaw = [];
zhu_diubao_roll = [];
zhu_diubao_angle = [];
j=1;
for i = zhu_diubao_index'zhu_diubao_pitch(j) = zhu_angle_cabrage(i);zhu_diubao_yaw(j) = zhu_angle_Heading(i);zhu_diubao_roll(j) = zhu_angle_roll(i);zhu_diubao_angle(j) = Angle_Antenna_zhu(i); j=j+1;
end% 副站 全数据点分析（包含起飞和降落阶段）
fu_all_analysis_pitch = fu_angle_cabrage(fu_remove_eth==1);
fu_all_analysis_yaw = fu_angle_Heading(fu_remove_eth==1);
fu_all_analysis_roll = fu_angle_roll(fu_remove_eth==1);
scatter3(fu_all_analysis_pitch,fu_all_analysis_yaw,fu_all_analysis_roll,100,"red");
xlabel("俯仰角 Pitch");ylabel("航向角 Yaw");zlabel("滚动角 Roll");
title('副站：全过程的俯仰、滚动、航向 & 丢包图');
savefig(img,'3_1_全数据姿态与丢包.fig'); % 副站 飞行数据点分析
fu_fly_pitch = fu_angle_cabrage(32727:148411);
fu_fly_yaw = fu_angle_Heading(32727:148411);
fu_fly_roll = fu_angle_roll(32727:148411);
fu_fly_remove= fu_remove_eth(32727:148411);
fu_fly_analysis_pitch = fu_fly_pitch(fu_fly_remove==1);
fu_fly_analysis_yaw = fu_fly_yaw(fu_fly_remove==1);
fu_fly_analysis_roll = fu_fly_roll(fu_fly_remove==1);
scatter3(fu_fly_analysis_pitch,fu_fly_analysis_yaw,fu_fly_analysis_roll,200,"blue");
xlabel("俯仰角 Pitch");ylabel("航向角 Yaw");zlabel("滚动角 Roll");
title('副站：飞行过程的俯仰、滚动、航向 & 丢包图');
savefig(img,'3_2_飞行数据与丢包.fig'); % 主站 全数据点分析（包含起飞和降落阶段）
zhu_all_analysis_pitch = zhu_angle_cabrage(zhu_remove_eth==1);
zhu_all_analysis_yaw = zhu_angle_Heading(zhu_remove_eth==1);
zhu_all_analysis_roll = zhu_angle_roll(zhu_remove_eth==1);
scatter3(zhu_all_analysis_pitch,zhu_all_analysis_yaw,zhu_all_analysis_roll,100,"red");
xlabel("俯仰角 Pitch");ylabel("航向角 Yaw");zlabel("滚动角 Roll");
title('主站：全过程的俯仰、滚动、航向 & 丢包图');
savefig(img,'3_3_全数据姿态与丢包.fig'); % 主站 飞行数据点分析
zhu_fly_pitch = zhu_angle_cabrage(32727:148411);
zhu_fly_yaw = zhu_angle_Heading(32727:148411);
zhu_fly_roll = zhu_angle_roll(32727:148411);
zhu_fly_remove= zhu_remove_eth(32727:148411);
zhu_fly_analysis_pitch = zhu_fly_pitch(zhu_fly_remove==1);
zhu_fly_analysis_yaw = zhu_fly_yaw(zhu_fly_remove==1);
zhu_fly_analysis_roll = zhu_fly_roll(zhu_fly_remove==1);
scatter3(zhu_fly_analysis_pitch,zhu_fly_analysis_yaw,zhu_fly_analysis_roll,200,"blue");
xlabel("俯仰角 Pitch");ylabel("航向角 Yaw");zlabel("滚动角 Roll");
title('主站：飞行过程的俯仰、滚动、航向 & 丢包图');
savefig(img,'3_4_飞行数据与丢包.fig'); % 主副站 全数据点分析（包含起飞和降落阶段）
scatter3(fu_all_analysis_pitch,fu_all_analysis_yaw,fu_all_analysis_roll,300,"red");
hold on;
scatter3(zhu_all_analysis_pitch,zhu_all_analysis_yaw,zhu_all_analysis_roll,100,"blue");
xlabel("俯仰角 Pitch");ylabel("航向角 Yaw");zlabel("滚动角 Roll");
legend('副站丢包','主站丢包','Location','NorthEastOutside');
title('主副站：全过程的俯仰、滚动、航向 & 丢包图');
savefig(img,'3_5_全数据姿态与丢包.fig');
hold off;% 主副站 飞行数据点分析
scatter3(fu_fly_analysis_pitch,fu_fly_analysis_yaw,fu_fly_analysis_roll,300,"red");
hold on;
scatter3(zhu_fly_analysis_pitch,zhu_fly_analysis_yaw,zhu_fly_analysis_roll,100,"blue");
xlabel("俯仰角 Pitch");ylabel("航向角 Yaw");zlabel("滚动角 Roll");
legend('副站丢包','主站丢包','Location','NorthEastOutside');
title('主副站：飞行过程的俯仰、滚动、航向 & 丢包图');
savefig(img,'3_6_飞行数据与丢包.fig');
hold off;% 副站 天线方向图
plot(fu_diubao_roll,fu_diubao_angle,'.');
xlabel("滚动角 Roll");ylabel("天线方向图角");
title('副站 天线方向图');
savefig(img,'4_1_副站天线方向图.fig'); % 副站 天线仰角与丢包 ｜｜ 副站天线 俯仰角20度，10度
clf(img);
hold on;
plot(find(fu_remove_eth==1),fuyang_test(fu_remove_eth==1),'rs','LineWidth',2,'MarkerSize',5);
plot(find(fu_remove_eth==0),fuyang_test(fu_remove_eth==0),'b');
plot([0,length(fu_remove_eth)],[10,10],'linestyle','--','LineWidth',2); % 天线波束范围
legend('丢包','不丢包','天线范围','Location','NorthEastOutside');
xlabel("飞行过程");ylabel("角度");
title('副站：地面与靶弹的俯仰角 丢包示意图');
savefig(img,'4_2_靶弹与副站天线仰角丢包示意图.fig');
hold off;% 主站 天线方向图
plot(zhu_diubao_roll,zhu_diubao_angle,'.');
xlabel("滚动角 Roll");ylabel("天线方向图角");
title('主站 天线方向图');
savefig(img,'4_3_主站天线方向图.fig');
hold off;% 主站 天线仰角与丢包 ｜｜ 主站天线 H波束宽度33度，一半有16.5度 , 上翘9度，和起来有 25.5度 相对副站好
clf(img);
hold on;
plot(find(zhu_remove_eth==1),fuyang_test_zhu(zhu_remove_eth==1),'rs','LineWidth',2,'MarkerSize',5);
plot(find(zhu_remove_eth==0),fuyang_test_zhu(zhu_remove_eth==0),'b');
plot([0,length(zhu_remove_eth)],[25.5,25.5],'linestyle','--','LineWidth',2); % 天线波束范围
legend('丢包','不丢包','天线范围','Location','NorthEastOutside');
xlabel("飞行过程");ylabel("角度");
title('主站：地面与靶弹的俯仰角 丢包示意图');
savefig(img,'4_4_靶弹与主站天线仰角丢包示意图.fig'); %% 数据统计与导出
FILENAMES = FILENAMES';
LOSS_YTJ = LOSS_YTJ';
LOSS_RAD = LOSS_RAD';
LOSS_ETH = LOSS_ETH';
LOSS_YK = LOSS_YK';writetable(table(FILENAMES,LOSS_YTJ,LOSS_RAD,LOSS_ETH,LOSS_YK, ...'VariableNames',["FILENAMES","LOSS_YTJ","LOSS_RAD","LOSS_ETH", ..."LOSS_YK"]),'丢包数据统计（一体机丢包-遥测链路丢包-地面网口丢包-遥控丢包）.csv');fu_diubao_pitch = fu_diubao_pitch';
fu_diubao_yaw = fu_diubao_yaw';
fu_diubao_roll = fu_diubao_roll';
fu_diubao_angle = fu_diubao_angle';
fu_diubao_jin = fu_diubao_jin';
fu_diubao_wei = fu_diubao_wei';
fu_diubao_gao = fu_diubao_gao';writetable(table(fu_diubao_index,fu_diubao_pitch,fu_diubao_yaw,fu_diubao_roll, ...fu_diubao_angle, fu_diubao_jin,fu_diubao_wei, fu_diubao_gao, ...'VariableNames',["index","pitch","yaw","roll","angle","jin","wei","gao" ...]),'副站飞行过程中丢包与角度分析.csv');toc;%% 计算所需的函数
function UTC = ConvertNum2UTC(x)len = length(x);for i = 1:lenh = fix(x(i)/1000/3600);m = fix((x(i)/1000-h*3600)/60);s = fix(x(i)/1000-h*3600-m*60);UTC(i) = string(h)+"h"+string(m)+"m"+string(s)+"s";endUTC=UTC';
endfunction Write_UTC_control(control,change,UTC,filename)len = length(change);j=1;LIST_CONTROL = [""];LIST_UTC = [""];LIST_INDEX = [];for i=1:lenif(change(i)) == 1LIST_CONTROL(j) = dec2hex(control(i+1));LIST_UTC(j) = UTC(i+1);LIST_INDEX(j) = i;j=j+1;endendLIST_CONTROL = LIST_CONTROL';LIST_UTC = LIST_UTC';LIST_INDEX = LIST_INDEX';writetable(table(LIST_INDEX,LIST_CONTROL,LIST_UTC,'VariableNames',["index","control","utc"]),filename+'/指令切换与UTC时间.csv');
endfunction [Angle_Antenna,Length,fuyang_test] = Antenna_direction(local_j,local_w,aim_j,aim_w,high,yaw)%Angle_Antenna = Antenna_direction(109.621197,40.847384,fu_longitude,fu_latitude,fu_positiong_high,fu_angle_Heading)earthR = 6371000000 ;%单位厘米local_j_hudu= local_j*pi/180;%角度与弧度local_w_hudu= local_w*pi/180;%角度与弧度aim_j_hudu = aim_j.*pi/180 ;%角度与弧度aim_w_hudu = aim_w.*pi/180 ;%角度与弧度V1=sin(aim_j_hudu-local_j_hudu).*cos(aim_w_hudu);V2=cos(local_w_hudu).*sin(aim_w_hudu) -sin(local_w_hudu).*cos(aim_w_hudu).*cos(aim_j_hudu-local_j_hudu);AB2=sin((local_w_hudu-aim_w_hudu)/2).*sin((local_w_hudu-aim_w_hudu)/2) + cos(local_w_hudu).*cos(aim_w_hudu).*sin((local_j_hudu-aim_j_hudu)/2).*sin((local_j_hudu-aim_j_hudu)/2);Length=2*earthR.*asin(sqrt(AB2))/1000; %单位厘米 直线距离 Angle=atan2(V1,V2).*180./pi; %该单位是-pi至pi，正值不变，需要将负值加上360°。而UB482是0至360°与真北的夹角fuyang_test =asin( abs((high - 1284))./sqrt((high - 1284).^2 + (Length).^2))./pi.*180;%假设地面天线架设5米% 改变航向角的真北角定义 原始定义：顺时针为负。 % 计算天线方向图Angles = - yaw - Angle;% 转化为0到360Angles_360 = Angles;for i = 1:length(Angles_360)if(Angles_360(i) < 0)Angles_360(i) = Angles_360(i) + 360;elseif(Angles_360(i) > 360)Angles_360(i) = Angles_360(i) - 360;elseAngles_360(i) = Angles_360(i);endend% 转化为 -180～+180Angles_zf_180 = Angles_360;for i = 1:length(Angles_zf_180)if(Angles_zf_180(i) > 180)Angles_zf_180(i) = Angles_zf_180(i) - 360;elseAngles_zf_180(i) = Angles_zf_180(i);endend% 转化为天线方向图Angle_Antenna = Angles_zf_180;for i = 1:length(Angle_Antenna)if(Angle_Antenna(i) >= 0)Angle_Antenna(i) = Angle_Antenna(i) - 180;elseAngle_Antenna(i) = Angle_Antenna(i) + 180;endend
end
function plot_k_means(pitch,yaw,roll,n)[idx,C] = kmeans([pitch;yaw;roll],n);disp(size(idx))for i = 1:nscatter3(pitch(idx==i),yaw(idx==i),roll(idx==i),".");hold on;endhold off;
end
相关文章：
画图，matlab，
clear;close all;clc;tic;dirOutput dir(*.dat); % 罗列所有后缀-1.dat的文件列表，罗列BDDATA的数据 filenames string({dirOutput.name}); % 提取文件名%% 丢包统计 FILENAMES [""]; LOSS_YTJ [ ]; LOSS_RAD [ ]; LOSS_ETH [ ]…...
any/all 子查询优化规则的原理与解析 | OceanBase查询优化
背景 在通常情况下，当遇到包含any/all子查询的语句时，往往需要遵循嵌套执行的方式，因此其查询效率较低。Oceanbase中制定了相应的any/all子查询优化规则，能够能够识别并优化符合条件的any/all子查询，从而有效提升查询…...
Visio——导出的PDF文件缺乏嵌入字体的解决办法 / 设置导出的PDF文件添加嵌入字体的方法
导出PDF时，勾选 “符合PDF/A” 选项 这样就导出的PDF文件添加了嵌入字体了。...
python：用 sklearn SVM 构建分类模型，并评价
编写 test_sklearn_5.py 如下 # -*- coding: utf-8 -*- """ 使用 sklearn 估计器构建分类模型，并评价 """ import numpy as np import matplotlib.pyplot as plt from sklearn import datasets from sklearn.svm import SVC from sk…...
【Python】制作函数，并且实现【注册】【登录】功能
这段代码是一个简单的命令行论坛模拟系统，包含了用户注册、登录和退出的功能。让我们逐行分析代码，并解释每个部分的功能与逻辑： ### 1. 引入 hashlib 模块 python import hashlib - **功能**：引入 Python 内置的 hashlib 模块…...
【大模型】LLaMA-2:Open Foundation and Fine-Tuned Chat Models, July. 2023.
论文：LLaMA-2：Open Foundation and Fine-Tuned Chat Models, July. 2023. 链接：https://arxiv.org/abs/2307.09288 Introduction 创新点 7B - 70B 预训练 微调 开源Llama 2 和Llama 2-Chat，针对对话用例进行了优化Motivation A…...
 RTX NACK概述)
WebRTC服务质量（04）- 重传机制（01) RTX NACK概述
WebRTC服务质量（01）- Qos概述 WebRTC服务质量（02）- RTP协议 WebRTC服务质量（03）- RTCP协议 WebRTC服务质量（04）- 重传机制（01) RTX NACK概述 WebRTC服务质量（…...
Qt之connectSlotsByName分析
简介 用于界面设置信号槽自动生成，要求槽函数定义形式为on_< objectName >_< signal > connectSlotsByName 其定义在QMetaObject类中，为静态方法，在文件qobjectdefs.h struct Q_CORE_EXPORT QMetaObject {....static void c…...
【Unity3D】无限循环列表（扩展版）
基础版：【Unity技术分享】UGUI之ScrollRect优化_ugui scrollrect 优化-CSDN博客 using UnityEngine; using UnityEngine.UI; using System.Collections.Generic;public delegate void OnBaseLoopListItemCallback(GameObject cell, int index); public class BaseLo…...
校园点餐订餐外卖跑腿Java源码
简介： 一个非常实用的校园外卖系统，基于 SpringBoot 和 Vue 的开发。这一系统源于黑马的外卖案例项目 经过站长的进一步改进和优化，提供了更丰富的功能和更高的可用性。 这个项目的架构设计非常有趣。虽然它采用了SpringBoot和Vue的组合&am…...
如何使用 Python 连接 PostgreSQL 数据库？
在Python开发中，连接PostgreSQL数据库是一个常见的需求。 我们可以使用多种库来实现这一功能，其中最常用的是psycopg2。 下面我将详细介绍如何使用psycopg2来连接PostgreSQL数据库，并提供一些实际开发中的建议和注意事项。 1. 使用 psycop…...
音频开发中常见的知识体系
在 Linux 系统中，/dev/snd 目录包含与声音设备相关的文件。每个文件代表系统中的一部分音频硬件或音频控制接口。以下是你列出的文件及其含义： 一.基本术语 样本长度(sample)：样本是记录音频数据最基本的单位，计算机对每个通道采…...
国内外人工智能AI工具网站大全（一键收藏，应有尽有）
本文由 大侠(AhcaoZhu)原创，转载请声明。 链接: https://blog.csdn.net/Ahcao2008 国内外人工智能AI工具网站大全（一键收藏，应有尽有） 摘要一、AI写作工具二、AI图像工具2.1、常用AI图像工具2.2、AI图片插画生成2.3、AI图片背景移…...
HB1910数字IP程控交换机generate.php存在RCE漏洞
免责声明: 本文旨在提供有关特定漏洞的深入信息，帮助用户充分了解潜在的安全风险。发布此信息的目的在于提升网络安全意识和推动技术进步，未经授权访问系统、网络或应用程序，可能会导致法律责任或严重后果。因此，作者不对读者基于本文内容所采取的任何行为承担责任。读者在…...
基于Spring Boot的社区药房系统
一、系统背景与目的 随着医疗改革的深入和社区医疗服务的不断完善，社区药房在居民健康保障中扮演着越来越重要的角色。然而，传统的药房管理方式存在着库存管理混乱、药品销售不透明、客户信息管理不规范等问题。为了解决这些问题，基于Spring…...
level2逐笔委托查询接口
沪深逐笔委托队列查询 前置步骤 分配数据库服务器 查询模板 以下是沪深委托队列查询的请求模板： http://<数据库服务器>/sql?modeorder_book&code<股票代码>&offset<offset>&token<token>查询参数说明 参数名类型说明mo…...
固定资产分类，提升资产盘活效益
固定资产是企业长期使用的重要资源，涵盖范围广、种类多，不同的资产需要针对性管理。通过科学的分类与高效的盘活策略，不仅可以优化资源配置，还能提升企业资产的利用效率和经济效益。以下将详细解析固定资产的分类方式和盘活效益的…...
图像根据mask拼接时，边缘有色差 解决
目录 渐变融合（Feathering） 沿着轮廓线模糊： 代码： 泊松融合 效果比较好： 效果图： 源代码： 泊松融合，mask不扩大试验 效果图： 源代码： 两个图像根据mask拼接时，边缘有色差 渐变融合（Feathering） import numpy as np import cv2# 假设 img1, img2 是两个…...
练习题：一维数组
练习题 第一题 键盘录入一组数列，利用冒泡排序将数据由大到小排序 代码 #include <stdio.h>int arr_home01() {int arr[10];int i,j,temp;printf("请输入10个测试整数：\n");int len sizeof(arr) / sizeof(arr[0]);for(i 0;i < …...
Linux系统安装node.js
一、node官网下载想要的node版本 https://nodejs.org/en/download/package-manager 二、将tar.xz文件解压 tar -xvf node-vxxx.tar.xz 三、改文件夹的名字，改成nodejs mv node-xxx nodejs 四、复制nodejs文件，并上传到linux 服务器 /usr/local 目录下…...
lettuce 默认情况下连接池参数不生效，源码分析
先说结论： 1.LettuceConnectionFactory 属性 shareNativeConnection 默认为true，要想连接池生效，该参数设置为false; 2.使用redisTemplate模版封装的pipeline没有意义，autoFlashCommands 默认为true;spring2.0开始默认使用lettuc…...
三维无人机航迹算法的目标函数如何确定
一、定义目标函数 在三维无人机航迹算法中，目标函数的确定通常基于具体的任务需求和飞行约束。以下是一个简单的例子，展示了如何为三维无人机航迹规划定义一个目标函数。 例子：最小化飞行时间和避障的三维无人机航迹规划 1.任务描述：无人机需要从起点飞到终点，同时避开一些…...
Linux docker离线部署
1. Docker下载 Docker下载地址：https://mirrors.dahuatech.com/docker-ce/。本文下载当前最新版本，链接如下：https://mirrors.dahuatech.com/docker-ce/linux/static/stable/aarch64/docker-27.4.0.tgz。 2. 安装Docker 将压缩包上传到服务器…...
短剧系统开发教程概要
引言 随着移动互联网的快速发展，短剧内容因其简短、精炼、情节紧凑的特点，吸引了大量观众的喜爱。为了满足市场需求，开发一款功能完善、体验优良的短剧平台显得尤为重要。本文将详细介绍短剧源码的开发搭建过程，包括需求分析、技…...
MySQL分区建表例子
以下为MySQL分区建表的例子 CREATE TABLE mz_mjzcfxxmx (SERIALNUM_ID varchar(96) NOT NULL COMMENT 业务角度唯一性ID,DATAGENERATE_DATE datetime NOT NULL COMMENT 业务数据产生时间,ROW_ID varchar(300) DEFAULT NULL COMMENT 数据角度唯一性ID,TASK_ID varchar(96) DEFAU…...
nginx变量
一、Nginx 变量概述 Nginx 变量是一种在 Nginx 配置中用于存储和操作数据的机制，它们可以在不同的配置块（如 http、server、location 等）中使用，以实现动态配置和灵活的请求处理。变量的值可以根据各种条件（如请求头信…...
网络编程 02：IP 地址，IP 地址的作用、分类，通过 Java 实现 IP 地址的信息获取
一、概述 记录时间 [2024-12-18] 前置文章：网络编程 01：计算机网络概述，网络的作用，网络通信的要素，以及网络通信协议与分层模型 本文讲述网络编程相关知识——IP 地址，包括 IP 地址的作用、分类ÿ…...
小红书笔记详情API接口：解锁社交媒体商业价值的钥匙
在当今数字化时代，社交媒体平台已成为企业营销和品牌推广的重要渠道。小红书，作为一个以内容分享为核心的社交媒体平台，汇聚了大量用户和创作者，他们在这里分享生活心得、购物体验、美妆技巧等多元化内容。小红书笔记详情API接口&…...
达梦查询表字段详细信息脚本（字段名称、描述、类型、长度及是否为空）
达梦查询表字段详细信息脚本（字段名称、描述、类型、长度及是否为空） 该SQL 脚本，用于查询表中字段的基本信息，包括字段名称、描述、数据类型、数据长度、是否为空及是否为主键等属性。 SQL 脚本 -- 输入变量 DECLAREp_owner VA…...
群落生态学研究进展】Hmsc包开展单物种和多物种分析的技术细节及Hmsc包的实际应用
联合物种分布模型（Joint Species Distribution Modelling，JSDM）在生态学领域，特别是群落生态学中发展最为迅速，它在分析和解读群落生态数据的革命性和独特视角使其受到广大国内外学者的关注。在学界不同研究团队研发出…...
CSS 语法
CSS 语法 CSS（层叠样式表）是一种用于描述HTML或XML文档样式的样式表语言。它允许您将样式信息与文档内容分离，从而更有效地控制网页的布局和外观。本文将详细介绍CSS的基本语法和结构，帮助您更好地理解和应用CSS。 CSS的基本结构 CSS由一系列的规则组成，每个规则包含一…...
【区块链】区块链密码学基础
🌈个人主页: 鑫宝Code 🔥热门专栏: 闲话杂谈｜ 炫酷HTML | JavaScript基础 💫个人格言: "如无必要，勿增实体" 文章目录 区块链密码学基础引言一、哈希函数1.1 基本概念1.2 数学表达 二、非对称加密2.1…...
ik分词器了解 和 通过zip安装包的方式 将ik分词器安装到elasticsearch中
目录 1. ik分词器的作用（效果） （1）标准分析器效果 （2）ik_smart分词 （3）ik_max_word分词 2. 首先根据自己的elasticsearch的版本下载对应的ik分词器版本 3. 将下载好的ik分词器…...
【数字化】华为数字化转型架构蓝图
导读：华为的数字化转型规划团队在2016年年底基于对愿景的系统诠释，整合出了数字化转型架构蓝图。该蓝图共分为5层，旨在通过数字化转型实现客户交互方式的转变、作战方式的转变、公司各平台业务能力的数字化、服务化以及运营模式的转变。 目录…...
K8s中 statefulset 和deployment的区别
在 Kubernetes 中，StatefulSet 和 Deployment 是两种管理 Pod 的控制器，它们的主要区别在于 状态管理 和 Pod 的标识。以下是详细对比： 1. 功能定位 Deployment 用途：用于 无状态应用 的部署，例如 Web 服务、API 服务…...
C#多线程系列章节五
线程创建的几种写法 thread类 threadpool 线程池 task类 Thread t = new Thread(func); t.Start();//线程开始执行threadPool.QueueUserWorkItem(test)创建var t1 = new Task(() => test(“Task 1”));t.start();Task.Run(() => test(“Task 2”));Task.Factory.StartNew…...
了解 JVM 运行原理，掌握常见的内存模型以及性能调优的基本方法
JVM运行原理概述 Java Virtual Machine (JVM) 是 Java 程序运行的核心组件，它的运行机制包括加载字节码、字节码解释或编译成本地机器代码执行。下面是 JVM 的运行流程和核心部分： JVM 的运行机制 类加载（Class Loading）…...
游戏AI实现-寻路算法（DFS）
深度优先搜索算法（英语：Depth-First-Search，缩写为DFS）是一种用于遍历或搜索树或图的算法。 寻路地图搭建： 游戏AI实现-寻路地图搭建-CSDN博客 算法过程：遍历方向为从竖直向上沿顺时针方向 1.首先将开…...
Qlib专为AI而生的量化投资平台
1、三层 Infrastructure层 该层为量化研究提供了基础支持。Data Server模块为用户管理和检索原始数据提供了高性能的基础架构。Trainer模块提供了灵活的接口来定义模型的训练过程，同时也为分布式训练提供了接口。Model Manager模块可以让用户更好地管理繁多的AI模型…...
圣乔ERP系统downloadFile.action存在任意文件读取漏洞
免责声明: 本文旨在提供有关特定漏洞的深入信息，帮助用户充分了解潜在的安全风险。发布此信息的目的在于提升网络安全意识和推动技术进步，未经授权访问系统、网络或应用程序，可能会导致法律责任或严重后果。因此，作者不对读者基于本文内容所采取的任何行为承担责任。读者在…...
YOLOv11改进，YOLOv11添加DLKA-Attention可变形大核注意力，WACV2024 ，二次创新C3k2结构
摘要 作者引入了一种称为可变形大核注意力 （D-LKA Attention） 的新方法来增强医学图像分割。这种方法使用大型卷积内核有效地捕获体积上下文，避免了过多的计算需求。D-LKA Attention 还受益于可变形卷积，以适应不同的数据模式。 理论介绍 大核卷积（Large Kernel Convolu…...
深入解析 `DataFrame.groupby` 和 `agg` 的用法及使用场景
深入解析 DataFrame.groupby 和 agg 的用法及使用场景 1. groupby 的基本用法语法：示例： 2. agg 的基本用法语法：示例： 3. first、sum、lambda 的用法3.1 first示例： 3.2 sum示例： 3.3 lambda示例ÿ…...
用.Net Core框架创建一个Web API接口服务器
我们选择一个Web Api类型的项目创建一个解决方案为解决方案取一个名称我们这里选择的是。Net 8.0框架 注意，需要勾选的项。 我们找到appsetting.json配置文件 appsettings.json配置文件内容如下 {"Logging": {"LogLevel": {"Default&quo…...
[论文阅读]Universal and transferable adversarial attacks on aligned language models
Universal and transferable adversarial attacks on aligned language models http://arxiv.org/abs/2307.15043 图 1：Aligned LLMs 不是对抗性 Aligned。我们的攻击构建了一个单一的对抗性提示，该提示始终绕过最先进的商业模式（包括 ChatG…...
【自动化】Python SeleniumUtil 油猴 工具 自动安装用户脚本
【自动化】Python SeleniumUtil 油猴 工具 【自动化】Python SeleniumUtil 工具-CSDN博客【自动化】Python SeleniumUtil 工具。https://blog.csdn.net/G971005287W/article/details/144565691 油猴工具 import timefrom selenium.webdriver.support.wait import WebDriverW…...
[LeetCode-Python版]21. 合并两个有序链表（迭代+递归两种解法）
题目 将两个升序链表合并为一个新的 升序 链表并返回。新链表是通过拼接给定的两个链表的所有节点组成的。 示例 1： 输入：l1 [1,2,4], l2 [1,3,4] 输出：[1,1,2,3,4,4] 示例 2： 输入：l1 [], l2 [] 输出&#x…...
WPF 相比 winform 的优势
wpf 相比 winform 的一些优点，网上也是众说纷纭，总的来说包括下面几点： 丰富的视觉效果：能够创建更具吸引力和现代化的用户界面，支持更复杂的图形和动画效果。不需要像 winform 一样，稍微做一点效果&#x…...
机器学习预处理-表格数据的分析与可视化
机器学习预处理-表格数据的分析与可视化 最近在做一些模型部署的工作，但是发现对于数据的处理、分析、训练方面还是缺少一些系统的学习，因此抽空余时间分析总结一些使用python进行数据处理的实用案例，希望能够方便自己已经其他人的Ctrl CV。…...
Linux——Shell
if 语句 格式：if list; then list; [ elif list; then list; ] ... [ else list; ] fi 单分支 if 条件表达式; then 命令 fi 示例： #!/bin/bash N10 if [ $N -gt 5 ]; then echo yes fi # bash test.sh yes 双分支 if 条件表达式; then 命令 else 命令…...
《深入探究：C++ 在多方面对 C 语言实现的优化》
目录 一、C 在 C 上进行的优化二、C 关键字（C 98）三、C 的输入输出1. cin 和 cout 的使用2. cin、cout 和 scanf()、printf() 的区别 三、命名空间1. 命名空间的使用2. 嵌套命名空间3. 在多个头文件中使用相同的命名空间 四、函数缺省值1. 缺省值的使用2…...