符号方程求解（提示找不到解答）

liuyang5299

syms x a
solve('2*x*3.14-sin(2*x*3.14)-2*a=0',x)
请网友帮忙解答一下，是不是该方程无解析解

qibbxxt

1. Warning: Explicit solution could not be found.

复制代码

这不就是找不到解析解的意思吗？

liuyang5299

2# qibbxxt
请问能利用数值解来解吗 谢谢

zhouyang664

显然可以，但不知你的a是变量还是常数呢？
help - fsolve

qibbxxt

3# liuyang5299

1. >> f=@(x,a)2*x*3.14-sin(2*x*3.14)-2*a
   
2. 
   
3. f =
   
4. 
   
5. @(x,a)2*x*3.14-sin(2*x*3.14)-2*a
   
6. 
   
7. >> fplot(@(x)f(x,1),[-3,3])
   
8. >> fzero(@(x)f(x,1),0)
   
9. 
   
10. ans =
    
11. 
    
12. 0.4067
    
13. 
    
14. >> grid on
    
15. >> grid minor

复制代码

Walker0728

我觉得你的方程好像有问题，用solve或是fsolve都能不能求解，你可以看看帮助文档
SOLVE  Symbolic solution of algebraic equations.
    SOLVE('eqn1','eqn2',...,'eqnN')
    SOLVE('eqn1','eqn2',...,'eqnN','var1,var2,...,varN')
    SOLVE('eqn1','eqn2',...,'eqnN','var1','var2',...'varN')

    The eqns are symbolic expressions or strings specifying equations.  The
    vars are symbolic variables or strings specifying the unknown variables.
    SOLVE seeks zeros of the expressions or solutions of the equations.
    If not specified, the unknowns in the system are determined by FINDSYM.
    If no analytical solution is found and the number of equations equals
    the number of dependent variables, a numeric solution is attempted.

    Three different types of output are possible.  For one equation and one
    output, the resulting solution is returned, with multiple solutions to
    a nonlinear equation in a symbolic vector.  For several equations and
    an equal number of outputs, the results are sorted in lexicographic
    order and assigned to the outputs.  For several equations and a single
    output, a structure containing the solutions is returned.

    Examples:

       solve('p*sin(x) = r') chooses 'x' as the unknown and returns

         ans =
         asin(r/p)

       [x,y] = solve('x^2 + x*y + y = 3','x^2 - 4*x + 3 = 0') returns

         x =
         [ 1]
         [ 3]

         y =
         [    1]
         [ -3/2]

       S = solve('x^2*y^2 - 2*x - 1 = 0','x^2 - y^2 - 1 = 0') returns
       the solutions in a structure.

         S =
           x: [8x1 sym]
           y: [8x1 sym]

       [u,v] = solve('a*u^2 + v^2 = 0','u - v = 1') regards 'a' as a
       parameter and solves the two equations for u and v.

       S = solve('a*u^2 + v^2','u - v = 1','a,u') regards 'v' as a
       parameter, solves the two equations, and returns S.a and S.u.

       [a,u,v] = solve('a*u^2 + v^2','u - v = 1','a^2 - 5*a + 6') solves
       the three equations for a, u and v.
我看上面的例子就有符合你的方程的形式，你可以换用其他类似的方程试试

liuyang5299

4# zhouyang664
a是常数，用a来表示x

lin2009

可以参考：http://forum.simwe.com/thread-902236-2-1.html  27#的方法求数值解。

1. 
   
2. f = @(x, a) 2*x*3.14 - sin(2*x*3.14) - 2*a;
   
3. g = @(a) fzero(@(x)f(x, a), 0);
   
4. 
   
5. subplot(2, 1, 1);
   
6. fplot(@(x)f(x, 1), [ 0, 5 ]);  % 5# qibbxxt 的程序
   
7. 
   
8. a = linspace(0, 5, 100);   
   
9. y = arrayfun(g, a);             % 对于一组a，求出对应的根y
   
10. subplot(2, 1, 2);
    
11. plot(a, y);
    

复制代码

zhouyang664

楼上几位正解！
但不是很理解，楼主说：a是常数，用a来表示x，那岂不是说必须用符号求解，不能用数值方法，但符号解又没有闭式表达式，我看楼主要悲剧了...
楼主将就一下，参考楼上用数值解法吧...

liuyang5299

9# zhouyang664
谢谢以上朋友的回复和解答。
如果能给出一个x(a)表达式，近似的能够逼近真实解也可以。数值解，能否给出定义域范围的关系式？

lin2009

用Maple来解有一个结果  0.1592356688*RootOf(_Z-2*a-sin(_Z)) ，但却不是“闭式表达式”。
实际上判断根的情况，可以结合图形来分析：
可以将 2*x*3.14-sin(2*x*3.14)-2*a = 0的实根
看成是曲线  y1(x) =  2*x*3.14 - 2*a 与曲线 y2(x) =  sin(2*x*3.14)的交点。
因此画出二者的曲线，就可以判断根的情况和定义域。
从图可知，定义域 -inf < a < inf。