5楼: Originally posted by study3040 at 2014-06-20 16:50:24
谢谢。请问，您说的数值近似公式指什么？请指教，谢谢。...

2楼: Originally posted by dingd at 2014-06-20 14:13:29
楼主图方便自己贴个图片省事，难道还要帮助的人照着图片一个个输入？应先学会懂得如何求助。

3楼: Originally posted by arising2010 at 2014-06-20 15:49:57
可以试试将积分项用数值近似公式表示，把整个方程组转化成一个近似的代数方程组，然后利用相关求解非线性代数的方法求解

Constant s0=0.055, r=0.04, c=1000;
Parameter theta1,alfa,p,q;
ConstStr g=-p*cos(alfa)-q*sin(alfa),
         fun_cal1=1/sqrt(abs(2*(p*cos(aaa)+q*sin(aaa)+g)/c)),
         fun_cal2=cos(bbb)/sqrt(abs(2*(p*cos(bbb)+q*sin(bbb)+g)/c)),
         fun_cal3=sin(ccc)/sqrt(abs(2*(p*cos(ccc)+q*sin(ccc)+g)/c));
Function
    s0-theta1*r-Integrate(fun_cal1,aaa=theta1,alfa)+Integrate(fun_cal1,aaa=alfa,pi/2);
    r*(1-sin(theta1))-Integrate(fun_cal2,bbb=theta1,alfa)+Integrate(fun_cal2,bbb=alfa,pi/2);
    r*cos(theta1)-Integrate(fun_cal3,ccc=theta1,alfa)+Integrate(fun_cal3,ccc=alfa,pi/2);
    r-1/sqrt(abs(2*(p*cos(theta1)+q*sin(theta1)+g)/c));

7楼: Originally posted by dingd at 2014-06-20 23:00:16
1stOpt求解：

Constant s0=0.055, r=0.04, c=1000;
Parameter theta1,alfa,p,q;
ConstStr g=-p*cos(alfa)-q*sin(alfa),
         fun_cal1=1/sqrt(abs(2*(p*cos(aaa)+q*sin(aaa)+g)/c)),
         fun_cal2= ...

6楼: Originally posted by arising2010 at 2014-06-20 18:26:57
求定积分有一些近似公式，简单的如“梯形公式”，复杂点的有龙贝格公式
你看看数值分析的书就明白了
...

7楼: Originally posted by dingd at 2014-06-20 23:00:16
1stOpt求解：

Constant s0=0.055, r=0.04, c=1000;
Parameter theta1,alfa,p,q;
ConstStr g=-p*cos(alfa)-q*sin(alfa),
         fun_cal1=1/sqrt(abs(2*(p*cos(aaa)+q*sin(aaa)+g)/c)),
         fun_cal2= ...