本帖最后由 elim 于 2023-8-23 02:56 编辑
相信这里的探讨的理论阐述会相当繁琐，或许可以从优化参数序列的(包络)极限来理解？总之很妙很有趣！
现在可设法用初等超越函数$\varphi_2$拟合楼上的 $Er$, 使 $\|E_r-\varphi_2\|_{\infty}\le 10^{-12}$
作为阶段性的有理论合实际意义的结果.
补充内容 (2023-9-15 05:35):
应为非初等函数。
本帖最后由 笨笨 于 2023-8-23 07:09 编辑
elim 发表于 2023-8-23 02:54
相信这里的探讨的理论阐述会相当繁琐，或许可以从优化参数序列的(包络)极限来理解？总之很妙很有趣！
现 ...
10^(-12)是楼主一直以来的目标，这样的初等函数族或初等超越函数族估计没人能找出来！
光搞个10^（-7）都这么费劲，10^（-12）想都别想！
【被拟合函数】
有一个比较复杂的函数：
\[\begin{matrix} F(x)=\sum_{n=0}^ \infty \begin{pmatrix}1/2\\n\end{pmatrix}^2x^{2n},&x\in\end{matrix}\]
设想找一个由初等函数构成的函数 \(G(x)\) 来拟合 \(F(x)\)，
使得两者之间的误差函数 \(Er(x)=F(x)-G(x)\) 的绝对值在定义域范围内越小越好。
目前希望能做到\(|Er(x)|_{MAX}<10^{-8}\) 的水平 。
在题目中给出了这样一个等式：
\[\begin{pmatrix}1/2\\k\end{pmatrix}=\frac{1}{k!}\prod_{j=0}^{k-1}\bigg(\frac{1}{2}-j\bigg)\]
这个等式只有在 k 取正整数时成立。
现在主要关注等式的左边，因为它就出现在 \(F(x)\) 的定义中。
在 Mathematica 中它对应 \(Binomial[\frac{1}{2},k]\) 函数。
用下列命令简化 \(F(x)\) 函数：
FullSimplify[\!\(TraditionalForm\`
\*UnderoverscriptBox[\(\\), \(n = 0\), \(\\)]\((
\*SuperscriptBox[\(Binomial[
\*FractionBox[\(1\), \(2\)], n]\), \(2\)]
\*SuperscriptBox[\(x\), \(2n\)])\)\)]
\[\text{FullSimplify}\left[\sum _{n=0}^{\infty } \binom{\frac{1}{2}}{n}^2 x^{2 n}\right]\]
结果得到 \(Hypergeometric2F1[-\frac{1}{2},-\frac{1}{2},1,x^2]\)。
所以定义的这个比较复杂的函数实际上是一个超几何函数：
\\tag{EQ:1}\]
当然，如果用上述等式右边来定义 \(F(x)\)
可以得到另外一个等价的函数。这里不作进一步讨论。
在 Mathematica 中用下列命令定义 \(F(x)\) ：
f := Hypergeometric2F1[-(1/2), -(1/2), 1, x^2]
本帖最后由 Jack315 于 2023-8-24 16:19 编辑
【拟合函数】
在题目中给出了这样一个拟合函数：
\
在 Mathematica 中用下列命令定义 \(G(x)\) 函数：
g := 1 + (3 x^2)/(10 + Sqrt)
接下来用下列命令画出图来看一下这两个函数之间的误差曲线：
LogPlot[(f - g), {x, 0, 1},
PlotRange -> All,
GridLines -> Automatic,
Frame -> True,
PlotLabel -> "待拟合函数与拟合函数的误差"]
说明：图能不上不上了，新人权限太低，一直被限制。
大家可以自己在软件里画出来看。
结果发现，随着 x 的增加，误差快速增大，使得拟合精度变差。
所以就需要再用一个修正函数 \(H(x)\) 来修正这个误差。
本帖最后由 Jack315 于 2023-8-24 16:48 编辑
【修正函数】
在题目中给出了一个含参数的修正函数：
\[\begin{matrix}H(x,a,b,c,d)=\frac{3x^{10}}{2^{17}}(1+\frac{[\frac{79}{48}+(\mu-1-\frac{79}{48})x]x^2}{^d}),&\mu=\frac{2^{17}}{3}(\frac{4}{\pi}-\frac{14}{11})\end{matrix}\tag{EQ:3}\]
用下列命令定义 \(H(x)\) 函数：
h := ((3 x^10)/2^17 (1 + ((79/48 + (\ - 1 - 79/48) x) x^2)/(1 + x^a (1 - x^b)^c)^d)) /. \ -> 2^17/3 (4/\ - 14/11)
题目要求能找出这样一组 \((a,b,c,d)\) 的值，
使得拟合的糖度能达到要求。
本帖最后由 Jack315 于 2023-8-24 17:49 编辑
【误差函数】
定义下列拟合的误差函数：
\
用下列命令定义误差函数：
err := f - g - h
题目中给了个例子：
例：当 \((x,a,b,c,d)=(0.997,0.53,2.75,0.82,3.38)\) 时，\(Er\approx6\times10^{-6}\)。
可以使用下列命令来查看误差函数的情况：
(*定义参数*)
x0 = x -> 0.997;
param = {a -> 0.53, b -> 2.75, c -> 0.82, d -> 3.38};
(*查看误差函数值*)
err /. x0 /. param
(*查看定义域内误差曲线*)
Table /. param, {x, (i - 1)/4, i/4},
PlotRange -> All,
GridLines -> Automatic,
PlotLegends -> "Expressions",
Frame -> True,
{i, 4}]
本帖最后由 笨笨 于 2023-8-24 17:42 编辑
Jack315 发表于 2023-8-24 17:02
【误差函数】
定义下列拟合的误差函数：
\
感谢先生闲暇时间回帖.有以下内容可供参考。
【解题思路】
假设给一组 (a,b,c,d) 的值，利用 (EQ:4) 定义的
Err 函数就能求出 x 在其定义域 范围内的最大值。
将这个过程包装到一个新的函数中：
maxErr :=
Module[{a = aVal, b = bVal, c = cVal, d = dVal},
MaxValue[{Abs], 0 <= x <= 1}, x \ Reals]
]
然后用下列命令进行求解：
sol = NMinimize[{maxErr, {a, b, c, d} \ Reals}, {a, b, c, d}]
运行一段时间后，给出了结果……
结果不贴出来了，因为答案是错误的。
原因是 \(G(x)+H(x)\) 函数不是在整个参数空间都能计算出结果，
因而导致（自动的）优化过程步入了歧途。
解决的方法大致有两种：
1. 修改代码，使得代码能规避上述导致错误的原因……期待高人出手。
2. 利用相关的方法（如 DOE、RSM：即梯度优化），一步步的进行手工优化。
这个方法效率特别低，所以不建议，但最终应该能得到一个结果。
至少可以得到一个局部解。
Jack315 发表于 2023-8-24 17:02
【误差函数】
定义下列拟合的误差函数：
\
先生你好，代码少一个中括号：
本帖最后由 Jack315 于 2023-8-24 18:21 编辑
笨笨 发表于 2023-8-24 18:05
先生你好，代码少一个中括号：
复制的时候漏掉了。重新来过：
(*定义参数*)
x0 = x -> 0.997;
param = {a -> 0.53, b -> 2.75, c -> 0.82, d -> 3.38};
(*查看误差函数值*)
err /. x0 /. param
(*查看定义域内误差曲线*)
Table /. param, {x, (i - 1)/4, i/4},
PlotRange -> All,
GridLines -> Automatic,
PlotLegends -> "Expressions",
Frame -> True,
PlotLabel -> "超几何函数与（修正后）拟合函数的误差"],
{i, 4}]