From: <cd...@us...> - 2011-10-25 15:18:41
|
Revision: 8858 http://octave.svn.sourceforge.net/octave/?rev=8858&view=rev Author: cdf Date: 2011-10-25 15:18:35 +0000 (Tue, 25 Oct 2011) Log Message: ----------- bugfix Modified Paths: -------------- trunk/octave-forge/extra/secs3d/inst/Utilities/Uphilipsmobiltymodeln.m Modified: trunk/octave-forge/extra/secs3d/inst/Utilities/Uphilipsmobiltymodeln.m =================================================================== --- trunk/octave-forge/extra/secs3d/inst/Utilities/Uphilipsmobiltymodeln.m 2011-10-25 15:17:47 UTC (rev 8857) +++ trunk/octave-forge/extra/secs3d/inst/Utilities/Uphilipsmobiltymodeln.m 2011-10-25 15:18:35 UTC (rev 8858) @@ -54,12 +54,18 @@ N_sc = N_D_star + N_A_star + p; - F = @(P) fF (P, m_n_star, m_p_star); - G = @(P) fG (P, T, T_300, m_n_star, m_0, a_g, b_g, c_g, alpha_g, alpha_prime_g, beta_g, gamma_g); P = (T ./ T_300).^2 ./ (f_CW ./ (3.97e13 .* N_sc .^ (-2/3)) + f_BH ./ ( 1.36e20 .* m_n_star ./ (m_0 .* (n + p)))); - N_sc_eff = N_D_star + G(P) .* N_A_star + p ./ F(P); + %% F = @(P) fF (P, m_n_star, m_p_star); + F = (0.7643 .* P .^ (0.6478) + 2.2999 + 6.5502 .* (m_n_star ./ m_p_star)) ./ ... + (P .^ (0.6478) + 2.3670 - 0.8552 .* (m_n_star ./ m_p_star)); + %% G = @(P) fG (P, T, T_300, m_n_star, m_0, a_g, b_g, c_g, alpha_g, alpha_prime_g, beta_g, gamma_g); + G = 1 - a_g ./ (b_g + P .* (m_0 .* T ./ (m_n_star .* T_300)) .^ alpha_g) .^ beta_g + ... + c_g ./ (P .* (m_0 .* T ./ (m_n_star .* T_300)) .^ alpha_prime_g) .^ gamma_g; + + N_sc_eff = N_D_star + G .* N_A_star + p ./ F; + mu_c = ((mu_max .* mu_min) ./ (mu_max - mu_min)) .* (T./T_300).^(1/2); mu_N = (mu_max.^2 ./ (mu_max - mu_min)) .* (T./T_300).^(3.*(alpha - 1/2)); mu_DAeh = mu_N .* (N_sc ./ N_sc_eff) .* (N_ref ./ N_sc).^alpha + mu_c .* ((n + p) ./ N_sc_eff); @@ -70,16 +76,16 @@ endfunction -function F = fF (P, m_n_star, m_p_star) - F = (0.7643 .* P .^ (0.6478) + 2.2999 + 6.5502 .* (m_n_star ./ m_p_star)) ./ ... - (P .^ (0.6478) + 2.3670 - 0.8552 .* (m_n_star ./ m_p_star)); -endfunction +%%function F = fF (P, m_n_star, m_p_star) +%% F = (0.7643 .* P .^ (0.6478) + 2.2999 + 6.5502 .* (m_n_star ./ m_p_star)) ./ ... +%% (P .^ (0.6478) + 2.3670 - 0.8552 .* (m_n_star ./ m_p_star)); +%%endfunction +%% +%%function G = fG (P, T, T_300, m_n_star, m_0, a_g, b_g, c_g, alpha_g, alpha_prime_g, beta_g, gamma_g) +%% G = 1 - a_g ./ (b_g + P .* (m_0 .* T ./ (m_n_star .* T_300)) .^ alpha_g) .^ beta_g + ... +%% c_g ./ (P .* (m_0 .* T ./ (m_n_star .* T_300)) .^ alpha_prime_g) .^ gamma_g; +%%endfunction -function G = fG (P, T, T_300, m_n_star, m_0, a_g, b_g, c_g, alpha_g, alpha_prime_g, beta_g, gamma_g) - G = 1 - a_g ./ (b_g + P .* (m_0 .* T ./ (m_n_star .* T_300)) .^ alpha_g) .^ beta_g + ... - c_g ./ (P .* (m_0 .* T ./ (m_n_star .* T_300)) .^ alpha_prime_g) .^ gamma_g; -endfunction - %!demo %! N_D = logspace (14, 20, 30); %! N_A = 1e12; This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site. |