关于超越方程求解的问题（完整题目已经跟帖，希望能得到帮助）

wanglu

x、y、fai、误差曲线：

wanglu

z、R曲线数值较大，单独绘制：

wanglu

怎么与其他软件结果不一样呢？楼主看看我的公式有没有错误。

wanglu

又验证了一下，在不同的初始值处求解，所得结果、精度都不一样，以上是在初始值fai=100处的结果。
在初始值fai=50处求解精度更高一些，约－5e－9～5e－9。
在某些初始值如fai＝1处求解，误差较大，故曲线不规则。

wanglu

以下绘制t＝1，初始值fai取0～100所求得的fai的图形及精度。
可以看出，取不同初始值时，求得的fai不一样。在某些初始值处，变化太大。

1. 
   
2. FCC文件
   
3. 缺省设置：
   
4. (XNote=请修改为X轴单位) (YNote=请修改为Y轴单位)
   
5. (AutoY=1) (XMin=0) (XMax=1) (YMin=0) (YMax=1)
   
6. (BorderPixels=60) (MultiplyX=1) (MultiplyY=1) (Grid=0) (DivideXY=10) (XYNumWidth=3) (DataMax=2)
   
7. (ForMax=50) (LoadDll=)
   
8. [CODE]
   
9. // 通用设置：
   
10. // (XNote=误差：红色；fai：蓝色) (YNote=请修改为Y轴单位)
    
11. // (AutoY=0) (XMin=0) (XMax=100) (YMin=-10) (YMax=100)
    
12. // (BorderPixels=80) (MultiplyX=1) (MultiplyY=1) (Grid=0) (DivideXY=10)  (XYNumWidth=3) (DataMax=2)
    
13. // (ForMax=50) (LoadDll="dll\FcData32W" "dll\XSLSF32W" ) (DotColor=0) (DotSize=10)
    
14. /*[LINE]
    
15. (_fai=1,1000,0,0,16711680)
    
16. (_误差=1,1000,0,0,255)
    
17. [LEND]*/
    
18. // [BODY]
    
19. !const["pi",acos(0)*2,1];
    
20. !const["r",pi/2-0.6685,1];
    
21. f(fai::t)=1/16*(-2932650*cos(1/2*t-fai-r)*pi^2+2932650*cos(-fai-r+3/2*t)*pi^2+1857345*cos(-fai-r+t)*pi^2-1857345*cos(2*t-fai-r)*pi^2-3870720*fai*sin(-fai-r+t)+3677184*fai*sin(-fai-r+3/2*t))/pi^2;
    
22. _fai(fai:s,i:t)=
    
23. {
    
24.   t=1,s=fai,
    
25.   i=XSLSF::pqrt(HFor("f"),&s,1e-6),
    
26.   if{i==10,printff{"已迭代{1,i}次，但未满足精度要求，返回值仅作参考。fai={2,r}\r\n",i,s}},
    
27.   s
    
28. };
    
29. 误差(fai,t)=1/16*(-2932650*cos(1/2*t-fai-r)*pi^2+2932650*cos(-fai-r+3/2*t)*pi^2+1857345*cos(-fai-r+t)*pi^2-1857345*cos(2*t-fai-r)*pi^2-3870720*fai*sin(-fai-r+t)+3677184*fai*sin(-fai-r+3/2*t))/pi^2;
    
30. _误差(fai:x,s:t)=
    
31. {
    
32.   t=1,s=fai,
    
33.   XSLSF::pqrt(HFor("f"),&s,1e-6),
    
34.   x=误差(s,t),
    
35.   if{abs(x)>1e-6,printff{"t={1,r}时，fai={2,r}，误差={3,r}\r\n",t,s,x}},
    
36.   x
    
37. };
    

复制代码

wanglu

以下是初始值fai取0～10的图形。

wanglu

发现初始值fai取0时精度较高，图形如下：
注意初始值fai取不同值时，所得结果不一样。

wanglu

这是初始值fai取2时的图形：

laihai718

非常感谢rocwoods 、TBE_Legend 、messager版主及 wanglu 同志，感谢你们的帮助，让我解决了自己的问题，本人是学机械的，研究方向是产品开发，目前在做螺杆泵设计及搅拌车的流体运动仿真，如果需要我帮忙的，尽管开口，我定尽全力。再一次感谢你们的帮助。QQ409693958.

flzs

楼主的问题得到解决了。。真有诚意