Which should I use for this integration.

1 次查看（过去 30 天）

显示更早的评论

Ashwani Kumar MAlviya 2020-5-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/524219-which-should-i-use-for-this-integration

评论： Ashwani Kumar MAlviya 2020-5-16

在 MATLAB Online 中打开

Hello everyone

I have to make a function to calculate significance of given arrays (lets say FC = C12 which is 47X1 double)

The fucnction should look like

function alpha = significance (FC)
% Variable FC is consistant values ​​coloum. For example FC = C12 will give the significane of C12 coloum.
n = length (FC); % Calculates the Number of consistent respondents (Nt)
x = sum (FC == mode (FC)); % Calculates the Number of respondents answered mode
p = 1/17; % Probability of selecting mode (Pi)
mu = n * p; 
mu1 = n * (1-p);
sigma = sqrt (n. * p. * (1 - p)); % Standerd deviation
if n <= 20
    i = 0;
    syms i
    alpha = symsum ((factorial (n) / (factorial (i) * factorial (ni)) * (1/17). ^ (i) * (1-1 / 17). ^ (ni)), i, 0 , x);
elseif n> 20
    syms y
    f = (- (y-mu). ^ 2) / (2 * sigma. ^ 2); % 5.05 is 2 * sigma ^ 2, & 2.7 = mu
    fun = (1 / (sigma * 2.5)) * (exp (f)); % 0.2516 = 1 / (sigma * 2.5), e = 2.71
    a = -inf;
    b = x; % x; % number of respondent selected mode
    alpha = int (fun, a, b);
  %%% ^ Which Integration function should I use here?
end

The problem here is, when I am running the script in commond window seperatly with values given as comment, it give me the currect value.

>> syms y
    f = (- (y-2.7). ^ 2) / (5.05); % 5.05 is 2 * sigma ^ 2, & 2.7 = mu
    fun = (0.2516) * (exp (f)); % 0.2516 = 1 / (sigma * 2.5), e = 2.71
    a = -inf;
    b = 19; % x; % number of respondent selected mode
    alpha = int (fun, a, b);
>> alpha
 
alpha =
 
(629 * 2020 ^ (1/2) * pi ^ (1/2) * (erf ((163 * 2020 ^ (1/2)) / 1010) + 1)) / 100000
 
>> ans = (629 * 2020 ^ (1/2) * pi ^ (1/2) * (erf ((163 * 2020 ^ (1/2)) / 1010) + 1)) / 100000
ans =
    1.0021