|
本帖最后由 永远 于 2023-3-24 12:51 编辑
关于公式的编程:
\(\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做的,怎么出现一大堆无聊的东西,怎么修改。
- a = 0.996; b = 0.723; c = 6.365;
- p[n_] := 1 - (1 - n/20)^2
- H[x_] := (
- Hypergeometric2F1[-1/2, -1/2, 1,
- x^2] - (1 + (3 x^2)/(10 + Sqrt[4 - 3 x^2])))/(3/2^17 x^10)
- W[x_, a_, b_, c_] := (21.381505539 x)/(1 + (1 - x^a)^b)^c
- U[a_, b_, c_] := 1/2 \!\(
- \*UnderoverscriptBox[\(\[Sum]\), \(k =
- 1\), \(19\)]\(\((1 + W[P[k], a, b, c] - H[P[k]])\)^2\)\)
- u[a_, b_, c_] := \!\(
- \*UnderoverscriptBox[\(\[Sum]\), \(k =
- 1\), \(19\)]\(\((1 + W[P[k], a, b, c] - H[P[k]])\)
- \*FractionBox[\(b\ c \((P[k])\)^a
- \*SuperscriptBox[\((1 - \((P[k])\)^a)\), \(b - 1\)] Log[
- P[k]]\), \(1 + \((1 - \((P[k])\)^a)\)^b\)] W[P[k], a, b, c]\)\)
- v[a_, b_, c_] := \!\(
- \*UnderoverscriptBox[\(\[Sum]\), \(k =
- 1\), \(19\)]\(\((1 + W[P[k], a, b, c] - H[P[k]])\)
- \*FractionBox[\(\(-c\)
- \*SuperscriptBox[\((1 - \((P[k])\)^a)\), \(b\)] Log[
- 1 - \((P[k])\)^a]\), \(1 + \((1 - \((P[k])\)^a)\)^b\)] W[P[k],
- a, b, c]\)\)
- w[a_, b_, c_] := -\!\(
- \*UnderoverscriptBox[\(\[Sum]\), \(k =
- 1\), \(19\)]\(\((1 + W[P[k], a, b, c] - H[P[k]])\) W[P[k], a, b,
- c] Log[1 + \((1 - \((P[k])\)^a)\)^b]\)\)
- a = a - t u[a, b, c]
- b = b - t v[a, b, c]
- c = c - t w[a, b, c]
复制代码 |
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有帐号?注册
x
|