2

根据 X 参数，我在 Matlab 中有一个方程。我想为 F(x) 的随机数量找到 X 的数量。我尝试了下面的代码。但它给了我两个不同的结果，而我的方程应该只有一个结果。

即使我尝试了roots(f)，solve(f)但它给了我一个错误：

？？？'sym' 类型的输入参数的未定义函数或方法 'isfinite'。

任何人都可以帮助我吗？我该怎么办 ？即使我对解决这个问题有错误的想法，请告诉我。谢谢

函数 betaDistribution_2(a,b)

    syms  x ;
    y=inline((x^(a-1))*((1-x)^(b-1)));
    beta=quad(y,0,1);
    g=(1/beta)*(x^(a-1))*((1-x)^(b-1));
    % I have this equation and I want to find the amount of x for the random 
    %amounts of p 
    p=int(g,x,0,x); 

    for i=0:50
       fxi=rand(1);
       f=p-fxi;
       xi=solve(f);
       result=eval(xi);
       disp(result)
    end

    end
4

1 回答 1

0

尝试过滤您的解决方案。

a=1;
b=2;
% a threshold for imagery part
SMALL=1e-9;

syms x real;
y=inline((x^(a-1))*((1-x)^(b-1)));
beta=quad(y,0,1);
g=(1/beta)*(x^(a-1))*((1-x)^(b-1));
p=int(g,x,0,x);

% return true for physically meaningfull results
filter = @(xc) abs(imag(xc))<SMALL && real(xc)>0 && subs(g, x, xc) > 0 && subs(p, x, xc)>0;

for m=0:50
  fxi=rand(1);
  f=p-fxi;
  xi=solve(f, x);
  result=eval(xi);

  idx = arrayfun (filter, result);
  result = result(idx);
  % make sure it is OK
  assert(length(result)==1);
  disp(result)
end
于 2013-01-02T11:18:29.780 回答