如何在matlab中求解符号方程

Question

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

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

？？？输入参数的未定义函数或方法“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

Answer 1

尝试过滤您的解决方案。

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