✅作者简介热爱科研的Matlab仿真开发者擅长毕业设计辅导、数学建模、数据处理、程序设计科研仿真。完整代码获取 定制创新 论文复现点击Matlab科研工作室 关注我领取海量matlab电子书和数学建模资料个人信条做科研博学之、审问之、慎思之、明辨之、笃行之是为博学慎思明辨笃行。 内容介绍一、引言在现代军事战略领域弹道导弹作为一种具有强大威慑力和打击能力的武器其轨迹的精确控制与跟踪至关重要。基于卫星星座的红外跟踪技术为可配置弹道导弹轨迹提供了全新的可能性使得从地球上任意起点到目的地的精确打击成为可能。这种技术的融合不仅提升了导弹的作战效能也深刻影响着现代战争的战略布局。二、卫星星座与红外跟踪技术基础一卫星星座架构星座组成卫星星座通常由多颗不同轨道高度、轨道倾角的卫星构成。例如低地球轨道LEO卫星具有较高的分辨率和快速的数据更新能力能够提供详细的目标信息中地球轨道MEO卫星则在覆盖范围和数据传输稳定性上具有优势地球静止轨道GEO卫星可实现对特定区域的持续监测。这些卫星相互协作形成一个全方位、多层次的观测网络。功能互补不同轨道卫星各司其职又相互补充。LEO 卫星凭借其靠近地球的优势能够快速获取高分辨率的红外图像对目标进行精确识别和定位MEO 卫星则负责在较大范围内进行目标搜索和跟踪将初步信息传递给其他卫星GEO 卫星保持对重点区域的持续凝视确保目标不脱离监测范围。通过这种功能互补卫星星座实现了对地球表面的全面覆盖和对目标的连续跟踪。二红外跟踪原理红外辐射特性任何物体只要温度高于绝对零度都会向外辐射红外线。弹道导弹在飞行过程中其发动机尾焰、弹体与空气摩擦产生的热量等都会使其成为一个强大的红外辐射源。不同部件的温度不同辐射的红外线波长和强度也有所差异这为红外跟踪提供了物理基础。跟踪技术实现卫星搭载的红外传感器通过探测目标的红外辐射获取目标的位置、速度、姿态等信息。红外传感器中的探测器将红外辐射转化为电信号经过信号处理和分析确定目标在图像中的位置。通过连续监测目标在多帧图像中的位置变化结合卫星自身的轨道参数和姿态信息就可以精确计算出目标的运动轨迹实现对弹道导弹的跟踪。三、可配置弹道导弹轨迹设计一轨迹规划原则目标导向根据打击目标的位置、性质和防御情况确定导弹的飞行轨迹。例如对于高价值且防御严密的目标可能需要设计复杂的变轨轨迹以避开敌方的防空系统对于一般目标则可以采用较为直接的弹道以提高打击效率。环境适应性考虑地球的重力场、大气环境以及可能存在的空间碎片等因素。在不同的高度和区域大气密度、重力加速度等参数会发生变化这会影响导弹的飞行性能。因此轨迹规划需要充分考虑这些因素确保导弹能够按照预定的轨迹飞行。二实现方式多阶段变轨弹道导弹通常采用多阶段飞行模式在不同阶段进行变轨操作。例如在助推段导弹依靠发动机推力快速上升到一定高度并获得足够的速度在自由飞行段可以通过微调发动机推力或利用姿态控制发动机进行变轨改变飞行方向和轨道高度在再入段根据目标位置和大气条件进行最后的轨道修正确保精确命中目标。智能算法应用利用先进的智能算法如遗传算法、粒子群算法等对导弹轨迹进行优化。这些算法可以在复杂的约束条件下快速搜索到最优的轨迹方案。例如遗传算法通过模拟生物进化过程中的选择、交叉和变异操作不断迭代优化轨迹参数使导弹在满足各种约束条件的前提下以最小的能量消耗或最短的飞行时间到达目标。四、基于卫星星座的红外跟踪与可配置弹道导弹轨迹的融合一信息交互机制实时数据传输卫星星座将实时获取的弹道导弹红外图像和跟踪数据通过高速通信链路传输到地面控制中心或导弹自身的控制系统。地面控制中心对这些数据进行分析和处理生成轨迹调整指令并及时反馈给导弹。这种实时的数据交互确保了对导弹轨迹的精确控制和调整。数据融合处理地面控制中心将卫星传来的红外跟踪数据与其他来源的信息如雷达数据、情报信息等进行融合处理。通过数据融合可以更全面、准确地了解导弹的飞行状态和目标信息为轨迹调整提供更可靠的依据。例如结合雷达提供的目标距离信息和红外图像中的目标姿态信息可以更精确地计算出导弹的命中点从而及时调整轨迹。二协同工作流程目标探测与识别卫星星座中的卫星利用红外传感器对地球表面进行大面积扫描探测潜在的导弹发射迹象。一旦发现疑似目标通过对红外图像的分析和特征提取识别目标是否为弹道导弹并确定其大致的飞行方向和初始轨道参数。轨迹跟踪与预测多颗卫星协同对已识别的导弹进行跟踪实时获取其位置和速度信息。利用这些数据结合导弹的动力学模型和轨道预测算法对导弹的未来轨迹进行预测。预测结果用于提前规划拦截方案或为导弹自身的轨迹调整提供参考。轨迹调整与优化根据卫星跟踪数据和预测结果地面控制中心或导弹自身的控制系统判断是否需要对轨迹进行调整。如果需要调整将计算出最优的调整参数并通过指令传输给导弹。导弹根据指令进行变轨操作实现对轨迹的精确控制确保准确命中目标。⛳️ 运行结果 部分代码classdef Satellite handle% Class for forming objects of type satellite, managing their orbital% elements/state, determining missile observability, and forming% measurements and simulated measurements.propertiesorbElems % [a; e; i; O; w; f] initially, after that the state x is used for everythingx % [px; py; pz; vx; vy; vz] in ECIconstellationName % string denoting which constellation satellite belongs tosatID % int denoting the satellite IDz % IR sensor measurementy % predicted IR sensor measurementcanObserve false; % boolean denoting whether the satellite can observe the ballistic missile at any given timeH_funcHhasLOS false;meetsIntensity false;IR_sensor_gain 1; % make IR sensors more or less sensitive to detectionI_threshold 1e-3; % TODO editmu 3.986e14; % Earths gravitational constant, m^3/s^2endmethodsfunction obj Satellite(satID, constellationName, a, e, i, O, w, f)obj.satID satID;obj.constellationName constellationName;obj.orbElems [a;e;i;O;w;f];[r, v] obj.elements_to_rv(a, e, i, O, w, f, obj.mu);obj.x [r;v];obj.formHfunction();endfunction [xdot] satelliteDynamics(obj, x)r x(1:3);v x(4:6);xdot zeros(6,1);xdot(1:3) v;xdot(4:6) (-obj.mu/(norm(r)^3))*r;endfunction [obj] propagateSatellite(obj, dt)% Updates object state obj.x with satellite dynamics function% and RK4 (faster in a large scale simulation than ode45())k1 obj.satelliteDynamics(obj.x);k2 obj.satelliteDynamics(obj.x (dt/2)*k1);k3 obj.satelliteDynamics(obj.x (dt/2)*k2);k4 obj.satelliteDynamics(obj.x dt*k3);obj.x obj.x (dt/6)*(k1 2*k2 2*k3 k4);endfunction [obj] checkLOS(obj, r_missile_ECI)% Function to check if missile has line of sight to satelliter_satellite_ECI obj.x(1:3);lambda 0:0.01:1;line_to_missile r_satellite_ECI (r_missile_ECI-r_satellite_ECI)*lambda;norm_line vecnorm(line_to_missile, 2);min_altitude_visible 6371000 20000; % Radius of Earth plus where atmosphere is thickif any(norm_line min_altitude_visible)obj.hasLOS false;elseobj.hasLOS true;endendfunction [obj] checkIntensity(obj, r_missile_ECI, beta)% Function that checks if the visible intensity of the% satellite to an IR sensor is high enough to be detected% beta: current infrared output of missile based on plume (if% in boost), area observed, emissivity, etc.r_satellite_ECI obj.x(1:3);r_sat_to_missile r_missile_ECI - r_satellite_ECI;I (obj.IR_sensor_gain*beta)/(norm(r_sat_to_missile)^2);if I obj.I_thresholdobj.meetsIntensity true;elseobj.meetsIntensity false;endendfunction [obj] checkObservability(obj, r_missile_ECI, beta)% Overall function to check if a satellite can observe/detect% the missile, conditions are having line of sight and meeting% the intensity thresholdobj.checkLOS(r_missile_ECI);obj.checkIntensity(r_missile_ECI, beta);if obj.hasLOS % meetsIntensityobj.canObserve true;elseobj.canObserve false;endendfunction [y] getMeasurement(obj, x_missile_ECI, R)% Function to collect either a predicted measurement (R 0)% or an actual measurement (R \neq 0)% Measurement vector: [az; el; az rate; el rate]x_satellite_ECI obj.x;r_sat_to_missile x_missile_ECI(1:3) - x_satellite_ECI(1:3);v_sat_to_missile x_missile_ECI(4:6) - x_satellite_ECI(4:6);az atan2(r_sat_to_missile(2), r_sat_to_missile(1)) normrnd(0, sqrt(R(1,1)));el atan2(-r_sat_to_missile(3), sqrt(r_sat_to_missile(1)^2 r_sat_to_missile(2)^2)) normrnd(0, sqrt(R(2,2)));az_rate (r_sat_to_missile(1)*v_sat_to_missile(2) - r_sat_to_missile(2)*v_sat_to_missile(1))/(r_sat_to_missile(1)^2 r_sat_to_missile(2)^2) normrnd(0, sqrt(R(3,3)));el_rate (r_sat_to_missile(3)*(r_sat_to_missile(1)*v_sat_to_missile(1) r_sat_to_missile(2)*v_sat_to_missile(2)) - v_sat_to_missile(3)*(r_sat_to_missile(1)^2 r_sat_to_missile(2)^2))/ ...(r_sat_to_missile(1)^2 r_sat_to_missile(2)^2 (r_sat_to_missile(3)^2)*sqrt(r_sat_to_missile(1)^2 r_sat_to_missile(2)^2)) normrnd(0, sqrt(R(4,4)));y [az; el; az_rate; el_rate];endfunction [obj] formHfunction(obj)% Function to form a matrix-valued function which, when% evaluated, forms the H matrix corresponding to the current% satellites measurementsyms x y z xdot ydot zdot realh1 atan2(y, x);h2 atan2(-z, sqrt(x^2 y^2));h3 (x*ydot - y*xdot)/(x^2 y^2);h4 (z*(x*xdot y*ydot) - zdot*(x^2 y^2))/ ...(x^2 y^2 (z^2)*sqrt(x^2 y^2));H_sym jacobian([h1;h2;h3;h4], [x; y; z; xdot; ydot; zdot]);obj.H_func matlabFunction(H_sym, Vars, {x, y, z, xdot, ydot, zdot});endfunction [H] evaluateH(obj, x_missile)obj.H obj.H_func(x_missile(1), x_missile(2), x_missile(3), x_missile(4), x_missile(5), x_missile(6));H obj.H;endfunction[r, v] elements_to_rv(obj, a, e, i, O, w, f, mu)%Input all element angles in degrees%Input mu in SI units, no km stuffi i*pi/180;O O*pi/180;w w*pi/180;f f*pi/180;r_norm (a*(1-e^2))/(1e*cos(f));r_peri [r_norm*cos(f); r_norm*sin(f); 0];p a*(1-e^2);v_peri [sqrt(mu/p)*-1*sin(f); sqrt(mu/p)*(ecos(f)); 0];% Rotation matrix termsR_11 cos(O)*cos(w) - sin(O)*sin(w)*cos(i);R_12 -cos(O)*sin(w) - sin(O)*cos(w)*cos(i);R_13 sin(O)*sin(i);R_21 sin(O)*cos(w) cos(O)*sin(w)*cos(i);R_22 -sin(O)*sin(w) cos(O)*cos(w)*cos(i);R_23 -cos(O)*sin(i);R_31 sin(w)*sin(i);R_32 cos(w)*sin(i);R_33 cos(i);R [R_11, R_12, R_13;R_21, R_22, R_23;R_31, R_32, R_33];r R*r_peri;v R*v_peri;endendend 参考文献[1]都欣欣,张其善.GPS信号模拟器卫星状态参数的算法研究[J].电子测量技术, 2007, 30(7):3.DOI:10.3969/j.issn.1002-7300.2007.07.034.更多免费数学建模和仿真教程关注领取
【导弹】基于卫星星座的红外跟踪可配置弹道导弹轨迹（从地球上任何起点和目的地）附Matlab代码
✅作者简介热爱科研的Matlab仿真开发者擅长毕业设计辅导、数学建模、数据处理、程序设计科研仿真。完整代码获取 定制创新 论文复现点击Matlab科研工作室 关注我领取海量matlab电子书和数学建模资料个人信条做科研博学之、审问之、慎思之、明辨之、笃行之是为博学慎思明辨笃行。 内容介绍一、引言在现代军事战略领域弹道导弹作为一种具有强大威慑力和打击能力的武器其轨迹的精确控制与跟踪至关重要。基于卫星星座的红外跟踪技术为可配置弹道导弹轨迹提供了全新的可能性使得从地球上任意起点到目的地的精确打击成为可能。这种技术的融合不仅提升了导弹的作战效能也深刻影响着现代战争的战略布局。二、卫星星座与红外跟踪技术基础一卫星星座架构星座组成卫星星座通常由多颗不同轨道高度、轨道倾角的卫星构成。例如低地球轨道LEO卫星具有较高的分辨率和快速的数据更新能力能够提供详细的目标信息中地球轨道MEO卫星则在覆盖范围和数据传输稳定性上具有优势地球静止轨道GEO卫星可实现对特定区域的持续监测。这些卫星相互协作形成一个全方位、多层次的观测网络。功能互补不同轨道卫星各司其职又相互补充。LEO 卫星凭借其靠近地球的优势能够快速获取高分辨率的红外图像对目标进行精确识别和定位MEO 卫星则负责在较大范围内进行目标搜索和跟踪将初步信息传递给其他卫星GEO 卫星保持对重点区域的持续凝视确保目标不脱离监测范围。通过这种功能互补卫星星座实现了对地球表面的全面覆盖和对目标的连续跟踪。二红外跟踪原理红外辐射特性任何物体只要温度高于绝对零度都会向外辐射红外线。弹道导弹在飞行过程中其发动机尾焰、弹体与空气摩擦产生的热量等都会使其成为一个强大的红外辐射源。不同部件的温度不同辐射的红外线波长和强度也有所差异这为红外跟踪提供了物理基础。跟踪技术实现卫星搭载的红外传感器通过探测目标的红外辐射获取目标的位置、速度、姿态等信息。红外传感器中的探测器将红外辐射转化为电信号经过信号处理和分析确定目标在图像中的位置。通过连续监测目标在多帧图像中的位置变化结合卫星自身的轨道参数和姿态信息就可以精确计算出目标的运动轨迹实现对弹道导弹的跟踪。三、可配置弹道导弹轨迹设计一轨迹规划原则目标导向根据打击目标的位置、性质和防御情况确定导弹的飞行轨迹。例如对于高价值且防御严密的目标可能需要设计复杂的变轨轨迹以避开敌方的防空系统对于一般目标则可以采用较为直接的弹道以提高打击效率。环境适应性考虑地球的重力场、大气环境以及可能存在的空间碎片等因素。在不同的高度和区域大气密度、重力加速度等参数会发生变化这会影响导弹的飞行性能。因此轨迹规划需要充分考虑这些因素确保导弹能够按照预定的轨迹飞行。二实现方式多阶段变轨弹道导弹通常采用多阶段飞行模式在不同阶段进行变轨操作。例如在助推段导弹依靠发动机推力快速上升到一定高度并获得足够的速度在自由飞行段可以通过微调发动机推力或利用姿态控制发动机进行变轨改变飞行方向和轨道高度在再入段根据目标位置和大气条件进行最后的轨道修正确保精确命中目标。智能算法应用利用先进的智能算法如遗传算法、粒子群算法等对导弹轨迹进行优化。这些算法可以在复杂的约束条件下快速搜索到最优的轨迹方案。例如遗传算法通过模拟生物进化过程中的选择、交叉和变异操作不断迭代优化轨迹参数使导弹在满足各种约束条件的前提下以最小的能量消耗或最短的飞行时间到达目标。四、基于卫星星座的红外跟踪与可配置弹道导弹轨迹的融合一信息交互机制实时数据传输卫星星座将实时获取的弹道导弹红外图像和跟踪数据通过高速通信链路传输到地面控制中心或导弹自身的控制系统。地面控制中心对这些数据进行分析和处理生成轨迹调整指令并及时反馈给导弹。这种实时的数据交互确保了对导弹轨迹的精确控制和调整。数据融合处理地面控制中心将卫星传来的红外跟踪数据与其他来源的信息如雷达数据、情报信息等进行融合处理。通过数据融合可以更全面、准确地了解导弹的飞行状态和目标信息为轨迹调整提供更可靠的依据。例如结合雷达提供的目标距离信息和红外图像中的目标姿态信息可以更精确地计算出导弹的命中点从而及时调整轨迹。二协同工作流程目标探测与识别卫星星座中的卫星利用红外传感器对地球表面进行大面积扫描探测潜在的导弹发射迹象。一旦发现疑似目标通过对红外图像的分析和特征提取识别目标是否为弹道导弹并确定其大致的飞行方向和初始轨道参数。轨迹跟踪与预测多颗卫星协同对已识别的导弹进行跟踪实时获取其位置和速度信息。利用这些数据结合导弹的动力学模型和轨道预测算法对导弹的未来轨迹进行预测。预测结果用于提前规划拦截方案或为导弹自身的轨迹调整提供参考。轨迹调整与优化根据卫星跟踪数据和预测结果地面控制中心或导弹自身的控制系统判断是否需要对轨迹进行调整。如果需要调整将计算出最优的调整参数并通过指令传输给导弹。导弹根据指令进行变轨操作实现对轨迹的精确控制确保准确命中目标。⛳️ 运行结果 部分代码classdef Satellite handle% Class for forming objects of type satellite, managing their orbital% elements/state, determining missile observability, and forming% measurements and simulated measurements.propertiesorbElems % [a; e; i; O; w; f] initially, after that the state x is used for everythingx % [px; py; pz; vx; vy; vz] in ECIconstellationName % string denoting which constellation satellite belongs tosatID % int denoting the satellite IDz % IR sensor measurementy % predicted IR sensor measurementcanObserve false; % boolean denoting whether the satellite can observe the ballistic missile at any given timeH_funcHhasLOS false;meetsIntensity false;IR_sensor_gain 1; % make IR sensors more or less sensitive to detectionI_threshold 1e-3; % TODO editmu 3.986e14; % Earths gravitational constant, m^3/s^2endmethodsfunction obj Satellite(satID, constellationName, a, e, i, O, w, f)obj.satID satID;obj.constellationName constellationName;obj.orbElems [a;e;i;O;w;f];[r, v] obj.elements_to_rv(a, e, i, O, w, f, obj.mu);obj.x [r;v];obj.formHfunction();endfunction [xdot] satelliteDynamics(obj, x)r x(1:3);v x(4:6);xdot zeros(6,1);xdot(1:3) v;xdot(4:6) (-obj.mu/(norm(r)^3))*r;endfunction [obj] propagateSatellite(obj, dt)% Updates object state obj.x with satellite dynamics function% and RK4 (faster in a large scale simulation than ode45())k1 obj.satelliteDynamics(obj.x);k2 obj.satelliteDynamics(obj.x (dt/2)*k1);k3 obj.satelliteDynamics(obj.x (dt/2)*k2);k4 obj.satelliteDynamics(obj.x dt*k3);obj.x obj.x (dt/6)*(k1 2*k2 2*k3 k4);endfunction [obj] checkLOS(obj, r_missile_ECI)% Function to check if missile has line of sight to satelliter_satellite_ECI obj.x(1:3);lambda 0:0.01:1;line_to_missile r_satellite_ECI (r_missile_ECI-r_satellite_ECI)*lambda;norm_line vecnorm(line_to_missile, 2);min_altitude_visible 6371000 20000; % Radius of Earth plus where atmosphere is thickif any(norm_line min_altitude_visible)obj.hasLOS false;elseobj.hasLOS true;endendfunction [obj] checkIntensity(obj, r_missile_ECI, beta)% Function that checks if the visible intensity of the% satellite to an IR sensor is high enough to be detected% beta: current infrared output of missile based on plume (if% in boost), area observed, emissivity, etc.r_satellite_ECI obj.x(1:3);r_sat_to_missile r_missile_ECI - r_satellite_ECI;I (obj.IR_sensor_gain*beta)/(norm(r_sat_to_missile)^2);if I obj.I_thresholdobj.meetsIntensity true;elseobj.meetsIntensity false;endendfunction [obj] checkObservability(obj, r_missile_ECI, beta)% Overall function to check if a satellite can observe/detect% the missile, conditions are having line of sight and meeting% the intensity thresholdobj.checkLOS(r_missile_ECI);obj.checkIntensity(r_missile_ECI, beta);if obj.hasLOS % meetsIntensityobj.canObserve true;elseobj.canObserve false;endendfunction [y] getMeasurement(obj, x_missile_ECI, R)% Function to collect either a predicted measurement (R 0)% or an actual measurement (R \neq 0)% Measurement vector: [az; el; az rate; el rate]x_satellite_ECI obj.x;r_sat_to_missile x_missile_ECI(1:3) - x_satellite_ECI(1:3);v_sat_to_missile x_missile_ECI(4:6) - x_satellite_ECI(4:6);az atan2(r_sat_to_missile(2), r_sat_to_missile(1)) normrnd(0, sqrt(R(1,1)));el atan2(-r_sat_to_missile(3), sqrt(r_sat_to_missile(1)^2 r_sat_to_missile(2)^2)) normrnd(0, sqrt(R(2,2)));az_rate (r_sat_to_missile(1)*v_sat_to_missile(2) - r_sat_to_missile(2)*v_sat_to_missile(1))/(r_sat_to_missile(1)^2 r_sat_to_missile(2)^2) normrnd(0, sqrt(R(3,3)));el_rate (r_sat_to_missile(3)*(r_sat_to_missile(1)*v_sat_to_missile(1) r_sat_to_missile(2)*v_sat_to_missile(2)) - v_sat_to_missile(3)*(r_sat_to_missile(1)^2 r_sat_to_missile(2)^2))/ ...(r_sat_to_missile(1)^2 r_sat_to_missile(2)^2 (r_sat_to_missile(3)^2)*sqrt(r_sat_to_missile(1)^2 r_sat_to_missile(2)^2)) normrnd(0, sqrt(R(4,4)));y [az; el; az_rate; el_rate];endfunction [obj] formHfunction(obj)% Function to form a matrix-valued function which, when% evaluated, forms the H matrix corresponding to the current% satellites measurementsyms x y z xdot ydot zdot realh1 atan2(y, x);h2 atan2(-z, sqrt(x^2 y^2));h3 (x*ydot - y*xdot)/(x^2 y^2);h4 (z*(x*xdot y*ydot) - zdot*(x^2 y^2))/ ...(x^2 y^2 (z^2)*sqrt(x^2 y^2));H_sym jacobian([h1;h2;h3;h4], [x; y; z; xdot; ydot; zdot]);obj.H_func matlabFunction(H_sym, Vars, {x, y, z, xdot, ydot, zdot});endfunction [H] evaluateH(obj, x_missile)obj.H obj.H_func(x_missile(1), x_missile(2), x_missile(3), x_missile(4), x_missile(5), x_missile(6));H obj.H;endfunction[r, v] elements_to_rv(obj, a, e, i, O, w, f, mu)%Input all element angles in degrees%Input mu in SI units, no km stuffi i*pi/180;O O*pi/180;w w*pi/180;f f*pi/180;r_norm (a*(1-e^2))/(1e*cos(f));r_peri [r_norm*cos(f); r_norm*sin(f); 0];p a*(1-e^2);v_peri [sqrt(mu/p)*-1*sin(f); sqrt(mu/p)*(ecos(f)); 0];% Rotation matrix termsR_11 cos(O)*cos(w) - sin(O)*sin(w)*cos(i);R_12 -cos(O)*sin(w) - sin(O)*cos(w)*cos(i);R_13 sin(O)*sin(i);R_21 sin(O)*cos(w) cos(O)*sin(w)*cos(i);R_22 -sin(O)*sin(w) cos(O)*cos(w)*cos(i);R_23 -cos(O)*sin(i);R_31 sin(w)*sin(i);R_32 cos(w)*sin(i);R_33 cos(i);R [R_11, R_12, R_13;R_21, R_22, R_23;R_31, R_32, R_33];r R*r_peri;v R*v_peri;endendend 参考文献[1]都欣欣,张其善.GPS信号模拟器卫星状态参数的算法研究[J].电子测量技术, 2007, 30(7):3.DOI:10.3969/j.issn.1002-7300.2007.07.034.更多免费数学建模和仿真教程关注领取
