外弹道设计

eata) / VA2;РHA3;Рsin (theata) / V;РРРРРРРРРРРРРРРРРРРРРРРРРРРРРРРРРРРРРРРРdx_ds = cos (theata);Рdy_ds = sin(theata);Рdt = h / V;Рx = x + RK(dx_ds,h);Рy = y + RK(dy_ds,h);РV = V + RK(dV_ds,h);Рtheata = theata + RK(dtheata_ds,h); arc = arc 4- h;Рn = n + 1;РiteraTime = iteraTime + dt;Рend time = iteraTime;Рfunction temp = tao(y)Рtemp = 288.9;РAA = 230;BB = -6.32 8e-3;CC = 1.172e-6;Rd =2 87.05;Рif (y <= 9300)Рtemp = 288.9 - y * 0.006328; end if (y > 9300 && v <12000)Рtemp = AA + (y - 9300) * BB + pow((y - 9300)r2) * CC; end if (y >= 12000 && y < 30000)Рtemp = 221.5;Рend function value = pi_y(y)Рvalue = 1;РRd = 287.05;Рif (y <= 9300)Рvalue = (1 - 2.1904e-5 夫 y)A5.4;РendРif (y > 9300 && v <12000)Рvalue = 0.2922575 * exp(-2.1206426 * (atan( (2.344 * (y - 9300)- 6328) / 32221.057) + 0.19392520));