A Example of Mathematical Modelling - Rocket Wreckage Ten Hour Blitz

“你微微地笑着，不同我说什么话，而我觉得，为了这个，我已等待得很久了。”——泰戈尔《飞鸟集》。

You smiled and talked to me of nothing and I felt that for this I had been waiting long. -Tagore, The Book of Birds.

目录

0. Abstract:

1. Introduction

2.Data analysis and visualisation

3. Summary

4. Acknowledgements

5. References

---

Citation format:HI Engie. an Example of Mathematical Modelling - Rocket Wreckage Ten Hour Blitz [EB/OL]. (2024-05-18)[2024-x-x]. http://t.csdnimg.cn/xxx.

0. Abstract:

This paper is about the 2024 Shenzhen Cup Mathematical Modelling (Eastern Provinces Mathematical Modelling League) A Question. A draft solution to the problem ‘Accurate localisation of multiple rocket debris’. At the same time, some codes and formulas are given.
Keywords: mathematical modelling; MATLAB 2023b; LaTeX

1. Introduction

Four ground-based detection devices completely track the trajectory of the wreckage.

2.Data analysis and visualisation

%Fig.1 Sonic boom source in polar coordinate system
deg1=7;
deg2=57;
deg3=97;
deg4=137;
deg5=169;
deg6=204;
deg7=253;
deg8=324;
deg21=27;
deg22=77;
deg23=127;
deg24=157;
deg25=191;
deg26=225;
deg27=277;
deg28=344;
deg91=45;
deg92=300;
degVector=[deg1,deg2,deg3,deg4,deg5,deg6,deg7,deg8,deg21,deg22,deg23,deg24,deg25,deg26,deg27,deg28,deg91,deg92];
theta=deg2rad(degVector);
rho=[0.95,0.9,0.85,0.8,0.95,0.9,0.85,0.97,0.9,0.85,0.8,0.95,0.9,0.85,0.95,1,0.9,0.9];
z=rho.*exp(1i*theta);
compass(z)
%set(gca,'FontName','Times New Roman','FontSize',7,'LineWidth',1.5);
xlabel('\it longitudes','FontSize',10,'FontName','Arial');
title('图1 极坐标系中的音爆源','FontSize',14,'FontName','宋体');
legend('超音速冲击波','FontSize',10);

%Fig.2 Sonic boom source in a 3D Cartesian coordinate system
[X,Y,Z] = sphere(100);
[U,V,W] = surfnorm(X,Y,Z);
quiver3(X,Y,Z,U,V,W,0)
axis equal
set(gca,'FontName','Times New Roman');
xlabel('\it X','FontSize',10,'FontName','Arial');
ylabel('\it Y','FontSize',10,'FontName','Arial');
zlabel('\it Z','FontSize',10,'FontName','Arial');
title('图2 三维直角坐标系中的音爆源','FontSize',14,'FontName','宋体', 'HorizontalAlignment', 'center','verticalAlignment','bottom');
%title(['图2 三维直角坐标系中的音爆源' sprintf('\n') 'Fig.2 Sonic Boom Source in 3D Cartesian Coordinate System.'],'FontSize',20,'FontName','宋体');
legend('超音速冲击波','FontSize',10); 

%Fig.3 3D Scatterplot
cla
load accidents hwydata                             % load data
long = -hwydata(:,2);                              % longitude data
lat = hwydata(:,3);                                % latitude data
rural = 100 - hwydata(:,17);                       % percent rural data
fatalities = hwydata(:,11);                        % fatalities data
scatter3(long,lat,rural,40,fatalities,'filled')    % draw the scatter plot
ax = gca;
ax.XDir = 'reverse';
view(-31,14)
xlabel('W. Longitude')
ylabel('N. Latitude')
zlabel('% Rural Population')
cb = colorbar;                                     % create and label the colorbar
cb.Label.String = 'Fatalities per 100M vehicle-miles';

%Figu.4 Contour map
Z = peaks(100);
zmin = floor(min(Z(:))); 
zmax = ceil(max(Z(:)));
zinc = (zmax - zmin) / 40;
zlevs = zmin:zinc:zmax;
figure
contour(Z,zlevs)
zindex = zmin:2:zmax;
hold on
contour(Z,zindex,'LineWidth',2)
hold off

%fig.5 3D Scatterplot 2
tbl = readtable('patients.xls');
s = scatter3(tbl,'Systolic','Diastolic','Weight','filled', ...
    'ColorVariable','Diastolic');

%Fig.6 3D coordinate plot
load seamount
tiledlayout(2,1) 
ax1 = nexttile.
ax2 = nexttile.
scatter3(ax1,x,y,z,‘MarkerFaceColor’,[0 .75 .75])
scatter3(ax2,x,y,z,‘*’)

%Fig.7 Visualisation of the detection device layout
clc
clear all
close all
x=[
10726890.26
10779337.12
10772720.45
10727863.3
10754427.3
10748880.97
10708013.29
];
y=[
3026798.652
3054836.928
3091442.455
3095892.975
3072750.271
3106574.223
3017563.823
];
z=[
824
727
742
850
786
678
575];
scatter3(x,y,z)
set(gca,'FontName','Times New Roman');
xlabel('\it X','FontSize',10,'FontName','Arial');
ylabel('\it Y','FontSize',10,'FontName','Arial');
zlabel('\it Z','FontSize',10,'FontName','Arial');
title('火箭残骸坠落的准确定位中可视化设备布局 ','FontSize',14,'FontName','宋体', 'HorizontalAlignment', 'center','verticalAlignment','bottom');
legend('设备','FontSize',10); 

%Fig. 8 Visualisation of the rocket debris fallout point
[X,Y,Z] = peaks;
contourf(X,Y,Z,LineColor="#4F4F4F")
set(gca,'FontName','Times New Roman');
xlabel('\it X','FontSize',10,'FontName','Arial'); 
ylabel('\it Y','FontSize',10,'FontName','Arial');
title('智能优化算法粒子群优化模型求解可视化','FontSize',14,'FontName','宋体', 'HorizontalAlignment', 'center','verticalAlignment','bottom');
colorbar
annotation("textarrow",[.53 .41],[.65 .47],String="火箭残骸坠落点")
annotation("textarrow",[.53 .59],[.65 .55])

%Fig.9 Visualisation of the coordinates of the position of the wreckage where 
lat = 27.52;
lon = 110.42;
A=lon;
geolimits([0 90],[0 180])
geoscatter(lat,lon,A,'^')
geobasemap streets
set(gca,'FontName','Times New Roman');
title('残骸发生音爆的位置坐标可视化','FontSize',14,'FontName','宋体', 'HorizontalAlignment', 'center','verticalAlignment','bottom');
colorbar
annotation("textarrow",[.53 .48],[.65 .56],String="音爆点")

//LaTeX
\left\{
\begin{aligned}
m_i=\sum_{j=1}^{3}{a_{ij}k_j}\\
n_i=\sum_{j=1}^{3}{a_{ij}∆r_j}
\end{aligned}
\right.
//others
//Word
m_i=\sum_{j=1}^{3}{a_{ij}k_j}
n_i=\sum_{j=1}^{3}{a_{ij}∆r_j}
//CSDN
$m_i=\sum_{j=1}^{3}{a_{ij}k_j},n_i=\sum_{j=1}^{3}{a_{ij}{\Delta r}_j}$

3. Summary

Also available.

4. Acknowledgements

Thank you for your time.

5. References

[1] Wang Pei, Xu Shan, Xiong Wei, et al. Research on trajectory localisation and monitoring technology of bullet and arrow wreckage based on TDOA [J]. Fire and Command Control, 2023, 48 (11): 139-144+151.
[2] ZHANG Snap, LIU Bingjie, WANG Ruichen. Calculation model of rocket debris fallout based on separation point information [J]. Computer Measurement and Control, 2022, 30 (11): 161-167. DOI:10.16526/j.cnki.11-4762/tp.2022.11.024.
[3] Wang Q, Li W, Gong JZ, et al. Drop point calculation model based on real-time positioning information of rocket debris [J]. Computer Measurement and Control, 2021, 29 (05): 154-158. DOI:10.16526/j.cnki.11-4762/tp.2021.05.031.
[4] Feng Song, Zeng Xiangjin, Huang Yuhao. Three-dimensional spatial gas source localisation based on improved fireworks optimisation algorithm[J]. China Safety Production Science and Technology,2024,20(03):69-76.
[5] https://www.bilibili.com/video/BV1wP411J7H3/?share_source=copy_web&vd_source=d3c3f249411283fe45758626abac08f1
[6] https://www.bilibili.com/video/BV18h4y187hC/?share_source=copy_web&vd_source=d3c3f249411283fe45758626abac08f1
[7] https://neuralwriter.com/zh/summarize-tool/
[8] https://www.deepl.com/zh/translator
[9] https://developer.baidu.com/article/details/2828883
[10] https://image.cnki.net/
[11] https://www.latexlive.com/
[12] https://blog.csdn.net/m0_37864814/article/details/102854015
[13] https://www.4kbizhi.com/wallpaper/6785.html
[14] https://www.toopic.cn/dnbz/12637.html
 

This article is also published on the following platforms: Baidu Bai Jiahao, bilibili, Tik Tok, Jin RItoutiao, WeChat Gong Zhonghao, Weibo, Zhihu, and China's professional IT community(CSDN).