Revision 0ab53102c02cd5c97fc58c71489748c2546cd368 (click the page title to view the current version)
function praktika(shot,tlak)
objemKomory = 0.079; % [m^3]
objem = 0.057; % [m^3]
k_B = 1.38e-23; % [J/K]
k_Rogowski = 1.1e7;
load([int2str(shot) '/CH1.txt']);
load([int2str(shot) '/CH2.txt']);
load([int2str(shot) '/CH3.txt']);
load([int2str(shot) '/CH4.txt']);
casy = CH1(:,1); % [s]
dt = casy(2)-casy(1); % [s]
derivaceProudu = k_Rogowski * CH1(:,2); % [A/s]
proud = cumsum(derivaceProudu)*dt; % [s]
% keyboard
B_t = cumsum(CH4(:,2));
H_alfa = CH2(:,2);
H_alfa_c = conv(H_alfa,ones(101,1)/101,'same');
H_alfa_sorted = sort(H_alfa_c(1:50:end));
H_alfa_min = H_alfa_sorted(round(length(H_alfa_sorted)*0.1));
H_alfa_max = H_alfa_sorted(round(length(H_alfa_sorted)*0.95));
plazmaInd = H_alfa_c>H_alfa_min+(H_alfa_max-H_alfa_min)*0.1;
delkaVyboje = sum(plazmaInd)*dt;
fprintf('delka vyboje: %f ms\n',delkaVyboje*1000);
U_loop = CH3(:,2);
I_p = proud - U_loop/0.01029;
I_p = I_p * 1.14;
fprintf('maximalni I_p: %f kA\n',max(conv(I_p,ones(101,1)/101,'same'))/1000);
P_ohm = U_loop .* I_p;
fprintf('maximalni P_ohm: %f kW\n',max(conv(P_ohm,ones(101,1)/101,'same'))/1000);
fprintf('prumerny P_ohm: %f kW\n',mean(P_ohm(plazmaInd))/1000);
T_e = (0.7*(I_p./U_loop)).^0.6666666667;
T_e(imag(T_e)~=0) = 0;
T_e(T_e > 50) = 0;
[T_e_max,T_e_max_ind] = max(conv(T_e,ones(101,1)/101,'same'));
fprintf('maximalni T_e: %f eV\n',T_e_max);
hustota = tlak/1000/k_B/300;
E_th_max = 3*k_B*(T_e_max*11600)*objem*hustota;
fprintf('maximalni E_th: %f J\n',E_th_max);
I_p_conv = conv(I_p,ones(101,1)/101,'same');
U_loop_conv = conv(U_loop,ones(101,1)/101,'same');
fprintf('tau_E: %f ms\n',E_th_max/I_p_conv(T_e_max_ind)/U_loop_conv(T_e_max_ind)*1000);
close all
figure;
plot(casy*1e3,U_loop,'g');
xlabel('t [ms]');
ylabel('U_loop [V]');
figure;
plot(casy*1e3,B_t,'b');
xlabel('t [ms]');
ylabel('B_t [??]');
figure;
plot(casy*1e3,I_p,'r');
xlabel('t [ms]');
ylabel('I_p [A]');
end