n="justify"> xlabel (Координата X, дам)
ylabel (Координата Y, дам)
% Зміна швидкостей за фронтом відбитої хвилі
for zet=0.1: 0.1: 0.9x=xwol: dx: dly=1: dy: shirDP0 (x, y) lt;=3 (x, y)=(- 1.33) * DP0 (x, y) DP0 (x, y) gt; 10 (x, y)=0 (x, y)=0.33 * DP0 (x, y) - 5.6DP0 (x, y) lt;=1 (x, y)=3.2 * DP0 (x, y) - 0.75 (x, y)=de (x, y) + 6.4 * DP0 (x, y)/(1 + 0.725 * DP0 (x, y) * zet) ( x, y)=PSK (x, y) * (1-zet) ^ 2 * exp ((- 1) * be (x, y) * zet) (x, y)=sqrt (2 * PSKNOV (x, y)/R0 (x, y))
End
% Побудова графіка поля швидкостей за фронтом відбитої хвилі
[X, Y]=meshgrid ([1: dx: dl]) (zet * 10 + 10) tall (1,1,1) grayfill (VNOV) (Поле швидкостей за фронтом відбитої хвилі , zet =)
xlabel (Координата X, дам)
ylabel (Координата Y, дам)
end
ДОДАТОК 2
ПРОГРАМА РОЗРАХУНКУ ПОЛІВ КОНЦЕНТРАЦІЙ в ударній хвилі
% Дана програма розроблена як продовження програми розрахунку
% полів швидкостей і використовує наявні в ній дані. Для практи -
% чеського використання тіло даної програми повинно бути
% розміщено безпосередньо після закінчення програми розрахунку
% полів швидкостей. Графічне зображення результатів здійснюється
% так само, як у першій частині розрахунку.
% Інтервал розбиття по відстані
dr=10
% Початкові концентрації
C (1,0)=1693000
For i=2: 1: 10
C (i, 0)=0
end
% Перехід до сферичної системі
for x=1: dx: dly=1: dy: shirrr=1: 1: 10x=y (rr)=varargin (x)
end
end
end
% Перехід до кінцевого кроці за часом
for zet=0: 0.1: 0.9=10 * zetk=0: 1: 9=T0 (rr) * zet
End
% Проведення розрахунку
F=dr*((r(1))^2+r(1)*dr/2+dr^2/12)*RO(1,k+1)=dr((r(1))^2+r(1)*dr/2+dr^2/12)*Ro(1,k)=RO(1,k+1)*VNOV(1,k+1)*dt*((- 1)*r(1)^2+r(1)*dr+dr^2/4)=RO(2,k+1)*VNOV(2,k+1)*dt*((r(1))^2+r(1)*dr+dr^2/4)(1,k+1)=P*C(2,k+1)/(F-H)+G*C(1,k)/(F-H)(1)=P/(F-H)(1)=G*C(1,0)/(F-H)i=2:1:9k=0:1:9(i)=(2/3)*(3*(r(i))^2+dr^3/4)*RO(i,k+1)(i)=(2/3)*( 3*(r(i))^2+dr^3/4)*RO(i,k)(i)=2*r(i)*dr*dt*RO(i,k+1)*VNOV(i,k+1)(i)=RO(i+1,k+1)*VNOV(i+1,k+1)*dt*((r(i)^2+r(i)*dr+dr^2/4)(i)=RO(i- 1, k + 1) * VNOV (i - 1, k + 1) * dt * ((r (i) ^ 2 + r (i) * dr + dr ^ 2/4) (i)=D/( A-C + E * al (i - 1)) (i)=(B (i) * C (i, k) -E * be (i - 1))/(A-C + E * al (i- 1))(i,k+1)=al(i)*C(i+1,k+1)+be(i)k=0:1:9=RO(10,k+1)*dr*((r(10))^2-r(10)*dr/2+dr^2/12)=C(10,k)*RO(10,k)*dr*((r(10))^2-r(10)*dr/2+dr^2/12)=RO(9,k+1)*VNOV(9,k+1*((r(10))^2-r(10)*dr+dr^2/4)*dt=RO(9,k+1)*VNOV(10,k+1) *(-(r(10))^2-r(10)*dr+dr^2/4)*dt(10,k+1)=(R*C(10,k)-S*(be(9))/(Q+S*al(9)+T)i=9:1:1k=0:1:9(i,k+1)=al(i)*C(i+1,k+1)+be(i)
