数学中国

用户名  找回密码
 注册
帖子
热搜: 活动 交友 discuz
查看: 5511|回复: 6

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

[复制链接]
发表于 2023-3-24 12:39 | 显示全部楼层 |阅读模式
本帖最后由 永远 于 2023-3-24 12:51 编辑

关于公式的编程:

C=π(a+b)(1+3λ210+43λ2+3217λ10(1+μλ(1+(1λu)v)w))

μ=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]
复制代码

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?注册

x
 楼主| 发表于 2023-3-24 14:04 | 显示全部楼层
来个老师帮我改一下吧
回复 支持 反对

使用道具 举报

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


本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?注册

x
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-3-24 17:49 | 显示全部楼层
本帖最后由 永远 于 2023-3-24 22:18 编辑
elim 发表于 2023-3-24 15:42
这是正在运行的版本。来个老师帮一下永远同学把。


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

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

σ(x)=1+(μ1)xa(1+(1xb)c)d
回复 支持 反对

使用道具 举报

发表于 2023-3-24 21:18 | 显示全部楼层
永远 发表于 2023-3-24 02:49
谢谢e老师指导,终于搞定了!其提供的程序还是有小瑕疵,不过还好已解决。

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

你的程序的错误是怎么产生,又是怎么搞定的?你的‘事实上’的事实是什么?弄懂我的不依赖目测确定参量的梯度法了吗?可否依样画葫芦弄个四参量的版本看看?
我的建议是 σ(x)=1+(μ1)x2(1+a(1xb)c)d,μ=2173(4π1411)
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-3-24 21:27 | 显示全部楼层
elim 发表于 2023-3-24 21:18
你的程序的错误是怎么产生,又是怎么搞定的?你的‘事实上’的事实是什么?弄懂我的不依赖目测确定参量的 ...

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

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

LaTEX预览输入 教程 符号库 加行内标签 加行间标签 
对应的 LaTEX 效果:

Archiver|手机版|小黑屋|数学中国 ( 京ICP备05040119号 )

GMT+8, 2025-7-12 15:48 , Processed in 0.093919 second(s), 16 queries .

Powered by Discuz! X3.4

Copyright © 2001-2020, Tencent Cloud.

快速回复 返回顶部 返回列表
\frac{\square}{\square}\sqrt{\square}\square_{\baguet}^{\baguet}\overarc{\square}\ \dot{\baguet}\left(\square\right)\binom{\square}{\square}\begin{cases}\square\\\square\end{cases}\ \begin{bmatrix}\square&\square\\\square&\square\end{bmatrix}\to\Rightarrow\mapsto\alpha\ \theta\ \pi\times\div\pm\because\angle\ \infty
\frac{\square}{\square}\sqrt{\square}\sqrt[\baguet]{\square}\square_{\baguet}\square^{\baguet}\square_{\baguet}^{\baguet}\sum_{\baguet}^{\baguet}\prod_{\baguet}^{\baguet}\coprod_{\baguet}^{\baguet}\int_{\baguet}^{\baguet}\lim_{\baguet}\lim_{\baguet}^{\baguet}\bigcup_{\baguet}^{\baguet}\bigcap_{\baguet}^{\baguet}\bigwedge_{\baguet}^{\baguet}\bigvee_{\baguet}^{\baguet}
\underline{\square}\overline{\square}\overrightarrow{\square}\overleftarrow{\square}\overleftrightarrow{\square}\underrightarrow{\square}\underleftarrow{\square}\underleftrightarrow{\square}\dot{\baguet}\hat{\baguet}\vec{\baguet}\tilde{\baguet}
\left(\square\right)\left[\square\right]\left\{\square\right\}\left|\square\right|\left\langle\square\right\rangle\left\lVert\square\right\rVert\left\lfloor\square\right\rfloor\left\lceil\square\right\rceil\binom{\square}{\square}\boxed{\square}
\begin{cases}\square\\\square\end{cases}\begin{matrix}\square&\square\\\square&\square\end{matrix}\begin{pmatrix}\square&\square\\\square&\square\end{pmatrix}\begin{bmatrix}\square&\square\\\square&\square\end{bmatrix}\begin{Bmatrix}\square&\square\\\square&\square\end{Bmatrix}\begin{vmatrix}\square&\square\\\square&\square\end{vmatrix}\begin{Vmatrix}\square&\square\\\square&\square\end{Vmatrix}\begin{array}{l|l}\square&\square\\\hline\square&\square\end{array}
\to\gets\leftrightarrow\nearrow\searrow\downarrow\uparrow\updownarrow\swarrow\nwarrow\Leftarrow\Rightarrow\Leftrightarrow\rightharpoonup\rightharpoondown\impliedby\implies\Longleftrightarrow\leftharpoonup\leftharpoondown\longleftarrow\longrightarrow\longleftrightarrow\Uparrow\Downarrow\Updownarrow\hookleftarrow\hookrightarrow\mapsto
\alpha\beta\gamma\Gamma\delta\Delta\epsilon\varepsilon\zeta\eta\theta\Theta\iota\kappa\varkappa\lambda\Lambda\mu\nu\xi\Xi\pi\Pi\varpi\rho\varrho\sigma\Sigma\tau\upsilon\Upsilon\phi\Phi\varphi\chi\psi\Psi\omega\Omega\digamma\vartheta\varsigma\mathbb{C}\mathbb{H}\mathbb{N}\mathbb{P}\mathbb{Q}\mathbb{R}\mathbb{Z}\Re\Im\aleph\partial\nabla
\times\cdot\ast\div\pm\mp\circ\backslash\oplus\ominus\otimes\odot\bullet\varnothing\neq\equiv\not\equiv\sim\approx\simeq\cong\geq\leq\ll\gg\succ\prec\in\ni\cup\cap\subset\supset\not\subset\not\supset\notin\not\ni\subseteq\supseteq\nsubseteq\nsupseteq\sqsubset\sqsupset\sqsubseteq\sqsupseteq\sqcap\sqcup\wedge\vee\neg\forall\exists\nexists\uplus\bigsqcup\bigodot\bigotimes\bigoplus\biguplus\bigcap\bigcup\bigvee\bigwedge
\because\therefore\angle\parallel\perp\top\nparallel\measuredangle\sphericalangle\diamond\diamondsuit\doteq\propto\infty\bowtie\square\smile\frown\bigtriangledown\triangle\triangleleft\triangleright\bigcirc \wr\amalg\models\preceq\mid\nmid\vdash\dashv\nless\ngtr\ldots\cdots\vdots\ddots\surd\ell\flat\sharp\natural\wp\clubsuit\heartsuit\spadesuit\oint\lfloor\rfloor\lceil\rceil\lbrace\rbrace\lbrack\rbrack\vert\hbar\aleph\dagger\ddagger

MathQuill输入:

Latex代码输入: