東京大学 情報理工学系研究科 コンピュータ科学専攻 2016年8月実施 専門科目II 問題1

标签：

• 东京大学
• 舍入误差与截断误差

Author

祭音Myyura (co-authored with GPT 5.6 SOL)

Description

Answer the following questions.

(1) Express each of $\cos(\alpha+\beta)$ and $\sin(\alpha+\beta)$ with $\cos(\alpha)$, $\sin(\alpha)$, $\cos(\beta)$, and $\sin(\beta)$.

(2) Suppose $0\le\alpha<\pi/2$. Express each of $\cos(\alpha/2)$ and $\sin(\alpha/2)$ with $\cos(\alpha)$ and $\sin(\alpha)$.

(3) Suppose $0\le\alpha<\pi/2$. Consider computing each of $\cos(\alpha/2)$ and $\sin(\alpha/2)$ using $\cos(\alpha)$ and $\sin(\alpha)$ with finite precision, by the expressions obtained in Question (2). Is there a risk of losing significant digits when $\alpha$ is close to zero? If there is, improve the expressions to avoid this loss.

(4) Suppose that $\alpha$ is expressed as $\frac{\pi}{2}\sum_{i=1}^{N}b_i2^{-i}$, where $N$ is a natural number and $b_i$ $(i=1,2,\ldots,N)$ is either 0 or 1.

(a) For $n=1,2,\ldots,N$, we define $\alpha_n=\frac{\pi}{2}\sum_{i=1}^{n}b_i2^{-i}$. Give expressions to compute $\cos(\alpha_{n+1})$ and $\sin(\alpha_{n+1})$ using $\cos(\alpha_n)$ and $\sin(\alpha_n)$, for $n=1,2,\ldots,N-1$.

(b) Give an algorithm to compute $\cos(\alpha)$ and $\sin(\alpha)$ using only the four arithmetic operations $(+,-,\times,/)$ and the extraction of square root $(\sqrt{\ })$.

(5) We define $\cosh(\alpha)=\frac{e^\alpha+e^{-\alpha}}{2}$ and $\sinh(\alpha)=\frac{e^\alpha-e^{-\alpha}}{2}$. Suppose that $0\le x<1$ and $x$ is expressed in a binary representation with finite precision. Design an algorithm to compute $\cosh(x\log_e2)$ and $\sinh(x\log_e2)$ using only the four arithmetic operations $(+,-,\times,/)$ and the extraction of square root $(\sqrt{\ })$, and describe an outline of the algorithm.

题目描述

回答下列问题。

（1）用 $\cos\alpha,\sin\alpha,\cos\beta,\sin\beta$ 表示 $\cos(\alpha+\beta)$ 和 $\sin(\alpha+\beta)$。

（2）设 $0\le\alpha<\pi/2$，用 $\cos\alpha,\sin\alpha$ 表示 $\cos(\alpha/2)$ 和 $\sin(\alpha/2)$。

（3）设 $0\le\alpha<\pi/2$。有限精度下，直接使用第（2）问公式在 $\alpha$ 接近 $0$ 时是否会丢失有效数字？若会，改写公式以避免该问题。

（4）设

$$\alpha=\frac\pi2\sum_{i=1}^{N}b_i2^{-i},\qquad b_i\in\{0,1\},$$

其中 $N$ 为自然数，$i=1,\ldots,N$。对 $n=1,\ldots,N$，定义 $\alpha_n=\frac\pi2\sum_{i=1}^{n}b_i2^{-i}$。

1. 对 $n=1,\ldots,N-1$，用 $\cos\alpha_n,\sin\alpha_n$ 表示 $\cos\alpha_{n+1},\sin\alpha_{n+1}$；
2. 设计只使用四则运算和开平方计算 $\cos\alpha,\sin\alpha$ 的算法。

（5）定义 $\cosh a=(e^a+e^{-a})/2$、$\sinh a=(e^a-e^{-a})/2$。设 $0\le x<1$ 有有限位二进制表示，设计只使用四则运算和开平方计算 $\cosh(x\log_e2)$、$\sinh(x\log_e2)$ 的算法，并简述其过程。

Kai

(1)

$$\cos(\alpha+\beta)=\cos\alpha\cos\beta-\sin\alpha\sin\beta,$$

$$\sin(\alpha+\beta)=\sin\alpha\cos\beta+\cos\alpha\sin\beta.$$

(2)

因 $0\le\alpha/2<\pi/4$，两根均取正号：

$$\cos\frac\alpha2=\sqrt{\frac{1+\cos\alpha}{2}},\qquad \sin\frac\alpha2=\sqrt{\frac{1-\cos\alpha}{2}}.$$

(3)

$\alpha\approx0$ 时 $\cos\alpha\approx1$，计算 $1-\cos\alpha$ 会发生相近数相减，导致 $\sin(\alpha/2)$ 丢失有效数字。利用

$$1-\cos\alpha=\frac{\sin^2\alpha}{1+\cos\alpha}$$

改写为

$$\sin\frac\alpha2=\frac{\sin\alpha}{\sqrt{2(1+\cos\alpha)}}.$$

$\cos(\alpha/2)$ 的原公式没有该消减问题。

(4)

令

$$\theta_{n+1}=\frac{\pi}{2^{n+2}},\qquad \alpha_{n+1}=\alpha_n+b_{n+1}\theta_{n+1}.$$

记 $C_n=\cos\alpha_n,S_n=\sin\alpha_n,c=\cos\theta_{n+1},s=\sin\theta_{n+1}$，则

$$(C_{n+1},S_{n+1})= \begin{cases} (C_n,S_n),&b_{n+1}=0,\\ (C_nc-S_ns,\ S_nc+C_ns),&b_{n+1}=1. \end{cases}$$

算法从 $(C_0,S_0)=(1,0)$ 开始。先置

$$c_1=s_1=\frac1{\sqrt2}$$

（对应 $\theta_1=\pi/4$），并递推

$$c_{i+1}=\sqrt{\frac{1+c_i}{2}},\qquad s_{i+1}=\frac{s_i}{\sqrt{2(1+c_i)}}.$$

按 $i=1,\ldots,N$ 依次用上面的加法公式：若 $b_i=1$，就把 $(C,S)$ 更新为 $(Cc_i-Ss_i,Sc_i+Cs_i)$；若 $b_i=0$ 则不变。最终 $(C,S)=(\cos\alpha,\sin\alpha)$。全过程只含四则运算与开平方。

(5)

写成 $x=\sum_{i=1}^{N}b_i2^{-i}$。从 $r_0=2$ 开始反复开平方，得到

$$r_i=\sqrt{r_{i-1}}=2^{2^{-i}}.$$

令 $q=1$，对每个 $b_i=1$ 乘上 $r_i$，则

$$q=\prod_{i=1}^{N}r_i^{b_i}=2^x=e^{x\log_e2}.$$

因此

$$\cosh(x\log_e2)=\frac{q+q^{-1}}2,\qquad \sinh(x\log_e2)=\frac{q-q^{-1}}2.$$

若还要避免 $x$ 很小时计算 $q-q^{-1}$ 的相消，可直接递推双曲函数。置

$$c_0=\cosh(\log 2)=\frac54,\qquad s_0=\sinh(\log 2)=\frac34,$$

并用

$$c_i=\sqrt{\frac{c_{i-1}+1}{2}},\qquad s_i=\frac{s_{i-1}}{2c_i}$$

得到 $c_i=\cosh(2^{-i}\log 2)$、$s_i=\sinh(2^{-i}\log 2)$。从 $(C,S)=(1,0)$ 出发，在 $b_i=1$ 时同步更新

$$(C,S)\leftarrow(Cc_i+Ss_i,\ Sc_i+Cs_i).$$

最终即得所求的双曲余弦和双曲正弦。该算法仍只用四则运算与开平方，且各步不需要相近数相减。