Scribd
Upload
11 views3 pages

Matlab Assignment Group 4

The document describes a process to calculate the activation energy and frequency factor for a reaction using temperature and reaction rate data. It performs linear regression on the natural logarithm of the reaction rates against the inverse of temperature to extract these values. The results indicate an activation energy of approximately 98674.96 J/mol, a frequency factor of 5.133e+14 s^-1, and a reaction constant at 305 K of 0.00647 s^-1.

Uploaded by

gdpsarmad
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
11 views3 pages

Matlab Assignment Group 4

The document describes a process to calculate the activation energy and frequency factor for a reaction using temperature and reaction rate data. It performs linear regression on the natural logarithm of the reaction rates against the inverse of temperature to extract these values. The results indicate an activation energy of approximately 98674.96 J/mol, a frequency factor of 5.133e+14 s^-1, and a reaction constant at 305 K of 0.00647 s^-1.

Uploaded by

gdpsarmad
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd

% Given data

>>T = [303.0 310.0 314.0 319.0 325.0]; % Temperature in Kelvin


>>k = [0.0053 0.0113 0.0190 0.0365 0.0727]; % Reaction rate
constant in s^-1

% Constants
>>R = 8.314; % Gas constant in J/mol.K

% Linearize data
>>x = 1./T; % 1/T
>>y = log(k); % ln(k)

% Perform linear regression (fit line)


>>p = polyfit(x, y, 1); % p(1) is slope, p(2) is intercept

% Extract activation energy E


>>slope = p(1);
>>E = -slope * R; % Activation energy in J/mol

% Calculate A (frequency factor)


>>A = exp(p(2));

% Calculate k at T = 305 K
>>T_new = 305; % K
>>k_305 = A * exp(-E/(R*T_new));

% Display results
>>fprintf('Activation Energy (E): %.2f J/mol\n', E);
>>fprintf('Frequency Factor (A): %.3e s^-1\n', A);
>>fprintf('Reaction constant at 305 K (k_305): %.5f s^-1\n',
k_305);

% Optional: Plot ln(k) vs 1/T


>>figure;
>>plot(x, y, 'o', 'MarkerFaceColor', 'b'); hold on;
>>plot(x, polyval(p, x), '-r');
>>xlabel('1/T (1/K)');
>>ylabel('ln(k)');
>>title('Arrhenius Plot');
>>legend('Data', 'Linear Fit');
>>grid on;
Activation Energy (E): 98674.96 J/mol
Frequency Factor (A): 5.133e+14 s^-1
Reaction constant at 305 K (k_305): 0.00647 s^-1

You might also like