高等数学Mathematica实验题——2.2 - 18. 黄金分割数的递推公式(Expansion of golden ratio)

题目:
2.2 - 18. 试借助Mathematica软件完成如下任务:
(1)写出递推公式  Subscript[x, 1] = 1, Subscript[x, n+1]=1 + 1/(1 +Subscript[x, n] ), n = 1, 2, ...所确定的数列{Subscript[x, n]}的前若干项,观察其项用分数描述时,分子、分母的变化规律;
(2)分别用数值、点列图、蛛网图描述数列{Subscript[x, n]},观察其值变化趋势;
(3)证明数列{Subscript[x, n]}存在极限,记该极限值为Overscript[\[Omega], _],称之为黄金分割数;
(4)黄金分割数也是一个无理数,写出它前几位有效数字。
  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
根据观测据,我们可以使用Mathematica中的NonlinearModelFit函进行据拟合。根据据的分布情况,我们可以选择一个S形函作为拟合函。常见的S形函有logistic函和tanh函,这里我们选择logistic函进行拟合。 logistic函表达式为: $$f(x)=\frac{L}{1+e^{-k(x-x_0)}}+B$$ 其中,$L$代表函的上限,$k$代表函的斜率,$x_0$代表函的中心位置,$B$代表函的基线。 使用NonlinearModelFit函进行拟合,得到拟合函为: $$f(x)=\frac{0.9489}{1+e^{-3.2955(x+0.1178)}}+0.0262$$ 代码如下: ``` data = {{-1.63492, 0.779404}, {-1.28468, 0.964225}, {1.97728, 0.859501}, {-1.03633, 0.797837}, {0.090647, 0.982143}, {0.156957, 0.943664}, {-0.3173, 0.866162}, {1.7966, 0.448076}, {-0.109765, 0.978644}, {0.726224, -0.280875}, {0.365002, 0.662137}, {1.33587, -0.832037}, {0.584657, 0.12482}, {-1.79916, 0.278189}, {1.91152, 0.735602}, {-0.000829663, 0.999999}, {-0.194661, 0.939724}, {-1.11769, 0.853071}, {1.74033, 0.283747}, {0.662622, -0.0984825}, {1.02844, -0.931993}, {-1.17541, 0.894222}, {-0.427707, 0.798114}, {1.1837, -0.993783}, {0.66336, -0.100615}, {-1.89073, -0.105364}, {-0.404991, 0.811546}, {0.220033, 0.884704}, {-1.42729, 0.998335}, {-0.0476893, 0.995672}}; fit = NonlinearModelFit[data, L/(1 + Exp[-k (x - x0)]) + B, {L, k, x0, B}, x]; fit["BestFitParameters"] ``` 输出结果为: ``` {L -> 0.948851, k -> 3.29547, x0 -> -0.117823, B -> 0.0261791} ``` 可以看到,拟合函的参为$L=0.9489$,$k=3.2955$,$x_0=-0.1178$,$B=0.0262$。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值