sicp每日一题[2.13-2.16]

Exercise 2.13

Show that under the assumption of small percentage tolerances there is a simple formula for the approximate percentage tolerance of the product of two intervals in terms of the tolerances of the factors. You may simplify the problem by assumingthat all numbers are positive.

After considerable work, Alyssa P. Hacker delivers her finished system. Several years later, after she has forgotten all about it, she gets a frenzied call from an irate user, Lem E. Tweakit. It seems that Lem has noticed that the formula for parallel resistors can be written in two algebraically equivalent ways:

He has written the following two programs, each of which computes the parallel-resistors formula differently:

复制代码
(define (par1 r1 r2)
  (div-interval (mul-interval r1 r2)
                (add-interval r1 r2)))

(define (par2 r1 r2)
  (let ((one (make-interval 1 1)))
    (div-interval
     one (add-interval (div-interval one r1)
                       (div-interval one r2)))))

Lem complains that Alyssa's program gives different answers for the two ways of computing. This is a serious complaint.


这道题一是要修改乘法函数,二是要演示两种计算并联电阻的方法计算的结果不一致。

复制代码
(define (mul-interval x y)
  (let ((cx (center x))
        (cy (center y))
        (px (/ (percent x) 100))
        (py (/ (percent y) 100)))
    (make-interval (* (- cx (* cx px)) (- cy (* cy py)))
                   (* (+ cx (* cx px)) (+ cy (* cy py))))))

(define r1 (make-interval 6.12 7.48))
(define r2 (make-interval 4.465 4.935))

(display-interval-tolerance (par1 r1 r2))
(display-interval-tolerance (par2 r1 r2))

; 执行结果
2.844199964577264 ± 22.613352145193332%
2.777440701636504 ± 7.05260392723452%

Exercise 2.14

Demonstrate that Lem is right. Investigate the behavior of the system on a variety of arithmetic expressions. Make some intervals A and B, and use them in computing the expressions A/A and A/B. You will get the most insight by using intervals whose width is a small percentage of the center value. Examine the results of the computation in center-percent form (see Exercise 2.12).


复制代码
(define r1 (make-interval 6.12 7.48))
(define r2 (make-interval 4.465 4.935))

(display-interval (div-interval r1 r1))
(display-interval (div-interval r1 r2))
(newline)
(display-interval-tolerance (div-interval r1 r1))
(display-interval-tolerance (div-interval r1 r2))

; 执行结果
[0.8181818181818182, 1.222222222222222]
[1.2401215805471126, 1.6752519596864504]

1.02020202020202 ± 19.801980198019795%
1.4576867701167815 ± 14.925373134328357%

根据执行结果来看,Lem显然是对的,对于 A/A,答案显然应该是 1 才对,但是程序计算出来的却是一个相对很大的范围,这是因为计算区间除法时,我们并没有定义"相等"这个概念,所以计算 A/A 时,它也按照通用的方法去计算了。

Exercise 2.15

Eva Lu Ator, another user, has also noticed the different intervals computed by different but algebraically equivalent expressions. She says that a formula to compute with intervals using Alyssa's system will produce tighter error bounds if it can be written in such a form that no variable that represents an uncertain number is repeated. Thus, she says, par2 is a "better" program for parallel resistances than par1. Is she right? Why?


Eva 说的是对的,对于练习 2.13 计算出的两个答案,par2 更准确,跟书上给出的答案是一致的,至于原因我觉得应该是区间之间的运算每次都会增加不确定性,par2 其实只有相加那一步是两个区间进行运算,其他的都可以看作确定数字与区间的运算,没有增加额外的不确定性;而 par1 进行了3次区间之间的运算------------乘法、加法和除法,乘法和除法都会使结果变得不确定。

Exercise 2.16

Explain, in general, why equivalent algebraic expressions may lead to different answers. Can you devise an interval-arithmetic package that does not have this shortcoming, or is this task impossible? (Warning: This problem is very difficult.)


复制代码
(define a (make-interval 2 4))
 
(define b (make-interval -2 0))
 
(define c (make-interval 3 8))
 
(define x (mul-interval a (add-interval b c)))
 
(define y (add-interval (mul-interval a b)
                        (mul-interval a c)))


(display-interval x)
(display-interval y)

; 执行结果
[2, 32]
[-2, 32]

根据上面的结果可以看到,区间运算连乘法分配律都无法保证。至于有没有无缺陷的区间运算规则,我大概想了一下,没有什么头绪,就没继续想了。。

相关推荐
AI大模型顾潇4 天前
[特殊字符] Prompt如何驱动大模型对本地文件实现自主变更:Cline技术深度解析
前端·人工智能·llm·微调·prompt·编程·ai大模型
Java知识库4 天前
Java基础知识面试题(已整理Java面试宝典pdf版)
java·面试·程序员·编程·面试题
zidea15 天前
Rust 闭包:捕获环境的魔法函数
rust·ai编程·函数式编程
jimin_callon19 天前
VBA第三十八期 VBA自贡分把表格图表生成PPT
开发语言·python·powerpoint·编程·vba·deepseek
chem411122 天前
Conmon lisp Demo
服务器·数据库·lisp
莲动渔舟23 天前
Nyquist插件基础:打印格式化字符串(LISP语言)
开发语言·lisp·音频处理·audacity
莲动渔舟23 天前
Nyquist插件基础:LISP语法-自定义函数
服务器·开发语言·lisp·音频处理·audacity
莲动渔舟23 天前
Nyquist插件基础:LISP语法-运算符
开发语言·lisp·音频处理·audacity
xinxiyinhe1 个月前
github免费编程类工具汇总与评估(二)
前端·后端·编程
独泪了无痕1 个月前
Optional 使用指南:彻底告别 NPE
后端·函数式编程