求助于uk702、天山草老师等网友。

永远

关于公式的编程：

\(\small C=\pi(a+b)\big(1+\dfrac{3\lambda^2}{10+\sqrt{4-3\lambda^2}}+\dfrac{3}{2^{17}}\lambda^{10}\big(1+\dfrac{\mu\lambda}{(1+(1-\lambda^u)^v)^w}\big)\big)\)

\(\mu=21.381505539,\;u=2.02777,\;v=0.823239,\;w=4.84343\)

我按照e老师19楼http://www.mathchina.com/bbs/for ... 4&fromuid=80637做的，怎么出现一大堆无聊的东西，怎么修改。

1. a = 0.996; b = 0.723; c = 6.365;
   
2. p[n_] := 1 - (1 - n/20)^2
   
3. H[x_] := (
   
4. Hypergeometric2F1[-1/2, -1/2, 1,
   
5.    x^2] - (1 + (3 x^2)/(10 + Sqrt[4 - 3 x^2])))/(3/2^17 x^10)
   
6. W[x_, a_, b_, c_] := (21.381505539 x)/(1 + (1 - x^a)^b)^c
   
7. U[a_, b_, c_] := 1/2 \!\(
   
8. \*UnderoverscriptBox[\(\[Sum]\), \(k =
   
9.      1\), \(19\)]\(\((1 + W[P[k], a, b, c] - H[P[k]])\)^2\)\)
   
10. u[a_, b_, c_] := \!\(
    
11. \*UnderoverscriptBox[\(\[Sum]\), \(k =
    
12.     1\), \(19\)]\(\((1 + W[P[k], a, b, c] - H[P[k]])\)
    
13. \*FractionBox[\(b\ c \((P[k])\)^a
    
14. \*SuperscriptBox[\((1 - \((P[k])\)^a)\), \(b - 1\)] Log[
    
15.       P[k]]\), \(1 + \((1 - \((P[k])\)^a)\)^b\)] W[P[k], a, b, c]\)\)
    
16. v[a_, b_, c_] := \!\(
    
17. \*UnderoverscriptBox[\(\[Sum]\), \(k =
    
18.     1\), \(19\)]\(\((1 + W[P[k], a, b, c] - H[P[k]])\)
    
19. \*FractionBox[\(\(-c\)
    
20. \*SuperscriptBox[\((1 - \((P[k])\)^a)\), \(b\)] Log[
    
21.       1 - \((P[k])\)^a]\), \(1 + \((1 - \((P[k])\)^a)\)^b\)] W[P[k],
    
22.     a, b, c]\)\)
    
23. w[a_, b_, c_] := -\!\(
    
24. \*UnderoverscriptBox[\(\[Sum]\), \(k =
    
25.      1\), \(19\)]\(\((1 + W[P[k], a, b, c] - H[P[k]])\) W[P[k], a, b,
    
26.      c] Log[1 + \((1 - \((P[k])\)^a)\)^b]\)\)
    
27. a = a - t u[a, b, c]
    
28. b = b - t v[a, b, c]
    
29. c = c - t w[a, b, c]

复制代码

永远

来个老师帮我改一下吧

Nicolas2050

巨婴

elim

这是正在运行的版本。来个老师帮一下永远同学把。

永远

elim 发表于 2023-3-24 15:42
这是正在运行的版本。来个老师帮一下永远同学把。

谢谢e老师指导，终于搞定了！其提供的程序还是有小瑕疵，不过还好已解决。

事实上，楼上老师你提供的3个参数没有4个的好。可否再编个程序看看。形如：

\(\Large\sigma \left( x \right) = 1 + \left( {\mu  - 1} \right)\frac{{{x^a}}}{{{{\left( {1 + {{\left( {1 - {x^b}} \right)}^c}} \right)}^d}}}\)

elim

永远 发表于 2023-3-24 02:49
谢谢e老师指导，终于搞定了！其提供的程序还是有小瑕疵，不过还好已解决。

事实上，楼上老师你提供的3 ...

你的程序的错误是怎么产生，又是怎么搞定的？你的‘事实上’的事实是什么？弄懂我的不依赖目测确定参量的梯度法了吗？可否依样画葫芦弄个四参量的版本看看？
我的建议是 \(\small\sigma(x)=1+{\large\frac{(\mu-1) x^2}{(1+a(1-x^b)^c)^d}},\;\mu=\frac{2^{17}}{3}\big(\frac{4}{\large\pi}-\frac{14}{11}\big)\)

永远

elim 发表于 2023-3-24 21:18
你的程序的错误是怎么产生，又是怎么搞定的？你的‘事实上’的事实是什么？弄懂我的不依赖目测确定参量的 ...

我只是想说目前我按照你的方法的程序能正常运行！其中程序里部分定义不明所以然，显然搞4个参数还是老师你来吧，我能力有限。