解超越方程

shoney

超越方程来自一本参考书，书上给出了解，我自己试着用matlab解了下，得不到结果，达人帮忙看看，谢谢
超越方程见附件1我的语句见附件2

zccbest

既然是超越方程当然是没有解释解了,不应该再用solve求解

qibbxxt

1. clear;clc;close all
   
2. s = 0.121;
   
3. p = 0.126;
   
4. v = 2.87;
   
5. f = @(x) (s./(x.*sqrt(pi).*exp(x.^2).*erf(x)) + p./(v.*x.*sqrt(pi).*exp(v.*x.^2).*erfc(v.*x)) - 1);
   
6. figure(1)
   
7. fplot(f,[-0.5,-0.15]);
   
8. figure(2)
   
9. fplot(f,[0.2,0.5]);
   
10. la = fsolve(f,[-0.2,0.2])

复制代码

1. 
   
2. Equation solved.
   
3. 
   
4. fsolve completed because the vector of function values is near zero
   
5. as measured by the default value of the function tolerance, and
   
6. the problem appears regular as measured by the gradient.
   
7. 
   
8. <stopping criteria details>
   
9. 
   
10. 
    
11. la =
    
12. 
    
13.    -0.2352    0.2817

复制代码

nwcwww

clear;clc;close all
s = 0.121;
p = 0.126;
v = 2.87;
f = @(x) (s./(x.*sqrt(pi).*exp(x.^2).*erf(x)) + p./(v.*x.*sqrt(pi).*exp(v.*x.^2).*erfc(v.*x)) - 1);
figure(1)
fplot(f,[-0.5,-0.15]);
figure(2 ...
qibbxxt 发表于 2011-6-15 16:51

f里面有一个v^2被写成v了