2017-11-11 Sa Oct 消参
2017-11-11 Sa Oct 消参
Prior versions:
- 2017-11-04 Sa Oct 消参
- 2017-11-10 Fr Oct 消参
2017-11-04 Sa
$ P(-3, 0) $ 在圆C $ (x-3)^2 + y^2 = 8^2 $ 内，动圆M与圆相切且过P点，求M点轨迹。
设切点 $ A(a, b) $，圆心 \(M(x, y)\)，则有 \(R_M = MA = MP = R_C-MC\)：
(a-3)^2 + b^2 &= 8^2 \\
\sqrt{(x-a)^2+(y-b)^2} &= \sqrt{(x+3)^2+y^2} \\
&= 8-\sqrt{(x-3)^2+y^2}
\end{aligned} \right. \]
UPD 2017-11-10 Fr 10:39PM
I thought about it at school on Monday. We can solve the relationship between \(a\) and \(b\) in $ (a-3)^2 + b^2 = 8^2 $, and \(x\) and \(y\) in $ \sqrt{(x+3)2+y2} = 8-\sqrt{(x-3)2+y2} $. Now there are just two unknowns left, and another equation \(\sqrt{(x-a)^2+(y-b)^2} = \sqrt{(x+3)^2+y^2}\). Then we can solve it.
On Monday I simplfied it by hand but the last step involved a polynomial with too much terms, thus I decided to go to Maxima at the weekend.
UPD 2017-11-11 Sa 2:50PM
Yesterday night, I have made some attempts on Maxima but failed. Let's do it again.
张惠妹01年专辑里的《记得》声音真细啊，可惜听了几遍就没有一开始听到的那种惊艳的感觉了。
第一次听到这首歌是JJ在一百天Live里唱记得+心墙+她说，挑的三首写给女生的作品。
Maxima 5.25.0 http://maxima.sourceforge.net
using Lisp Clozure Common Lisp Version 1.7-r14925M (WindowsX8632)
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
The function bug_report() provides bug reporting information.
(%i1) expand((x+3)^2+y^2-(8-sqrt((x-3)^2+y^2))^2)
;
2 2
(%o1) 16 sqrt(y + x - 6 x + 9) + 12 x - 64
(%i2) expand(16*16*(y^2+x^2-6*x+9)-(64-12*x)^2)
;
2 2
(%o2) 256 y + 112 x - 1792
(%i3) plot2d(%o2, [x,-10,10], [y,-10,10]);
plot2d: expression evaluates to non-numeric value everywhere in plotting range.
plot2d: nothing to plot.
(%o3) false
(%i4) plot2d(%o2=0, [x,-10,10], [y,-10,10]);
plot2d: expression evaluates to non-numeric value everywhere in plotting range.
plot2d: nothing to plot.
(%o4) false
(%i5) contour_plot(%o2=0, [x,-10,10], [y,-10,10]);contour_plot(%o2=0, [x,-10,10], [y,-10,10]);
(%o5)
(%i6)
(%o6)
(%i7) contour_plot(%o2, [x,-10,10], [y,-10,10]);
(%o7)
(%i8) expand((a-3)^2 + b^2 = 8^2)
;
2 2
(%o8) b + a - 6 a + 9 = 64
(%i9) expand((x-a)^2+(y-b)^2 = (x+3)^2+y^2)
;
2 2 2 2 2 2
(%o9) y - 2 b y + x - 2 a x + b + a = y + x + 6 x + 9
(%i10) expand((x-a)^2+(y-b)^2 - (x+3)^2+y^2);
2 2 2
(%o10) 2 y - 2 b y - 2 a x - 6 x + b + a - 9
(%i11) solve([(a-3)^2 + b^2 = 8^2], [b]);
2 2
(%o11) [b = - sqrt(- a + 6 a + 55), b = sqrt(- a + 6 a + 55)]
(%i12) solve([(a-3)^2 + b^2 = 8^2], [a]);
2 2
(%o12) [a = 3 - sqrt(64 - b ), a = sqrt(64 - b ) + 3]
(%i13) solve([%o2], [y]);
2 2
sqrt(7) sqrt(16 - x ) sqrt(7) sqrt(16 - x )
(%o13) [y = - ---------------------, y = ---------------------]
4 4
(%i14) solve([(x-3)^2 + y^2 = 8^2], [y]);
2 2
(%o14) [y = - sqrt(- x + 6 x + 55), y = sqrt(- x + 6 x + 55)]
(%i15) tentex(%o9);
2 2 2 2 2 2
(%o15) tentex(y - 2 b y + x - 2 a x + b + a = y + x + 6 x + 9)
(%i16) load(tentex);
;Compiler warnings for "C:/PROGRA~2/MAXIMA~1.0/share/maxima/5.25.0/share/tensor/tentex.lisp" :
; In TENREFORM: Undefined function CONTI
; In TENREFORM: Undefined function CONTI
; In TENREFORM: Undefined function NAME
; In TENREFORM: Undefined function NAME
; In TENREFORM: Undefined function DERI
; In TENREFORM: Undefined function DERI
; In TENREFORM: Undefined function COVI
; In TENREFORM: Undefined function DERI
; In TENREFORM: Undefined function COVI
; In TENREFORM: Undefined function COVI
; In TENREFORM: Undefined function RPOBJ
(%o16) C:/PROGRA~2/MAXIMA~1.0/share/maxima/5.25.0/share/tensor/tentex.lisp
(%i17) tentex(%o9);
Maxima encountered a Lisp error:
Undefined function RPOBJ called with arguments (((MEQUAL SIMP) ((MPLUS SIMP) ((MEXPT SIMP) $A 2) ((MEXPT SIMP) $B 2) ((MTIMES SIMP) -2 $A $X) ((MEXPT SIMP) $X 2) ((MTIMES SIMP) -2 $B $Y) ((MEXPT SIMP) $Y 2)) ((MPLUS SIMP) 9 ((MTIMES SIMP) 6 $X) ((MEXPT SIMP) $X 2) ((MEXPT SIMP) $Y 2)))) .
Automatically continuing.
To enable the Lisp debugger set *debugger-hook* to nil.
(%i18) tentex(a+b);
Maxima encountered a Lisp error:
Undefined function RPOBJ called with arguments (((MPLUS SIMP) $A $B)) .
Automatically continuing.
To enable the Lisp debugger set *debugger-hook* to nil.
(%i19) load(itensor);
;Compiler warnings for "C:/PROGRA~2/MAXIMA~1.0/share/maxima/5.25.0/share/tensor/itensor.lisp" :
; In $LC_L: Unused lexical variable L2
;Compiler warnings for "C:/PROGRA~2/MAXIMA~1.0/share/maxima/5.25.0/share/tensor/itensor.lisp" :
; In $LC_U: Unused lexical variable L1
(%o19) C:/PROGRA~2/MAXIMA~1.0/share/maxima/5.25.0/share/tensor/itensor.lisp
(%i20) load(tentex);
(%o20) C:/PROGRA~2/MAXIMA~1.0/share/maxima/5.25.0/share/tensor/tentex.lisp
(%i21) tentex(%o9);
$$y^2-2\,b\,y+x^2-2\,a\,x+b^2+a^2=y^2+x^2+6\,x+9$$
(%o21) false
(%i22) tentetx(%o8);
2 2
(%o22) tentetx(b + a - 6 a + 9 = 64)
(%i23) tentex(%o8);
$$b^2+a^2-6\,a+9=64$$
(%o23) false
(%i24)
\]
\]
(%i24) solve([%o8], [a^2+b^2]);
(%o24) []
(%i25) %o9
;
2 2 2 2 2 2
(%o25) y - 2 b y + x - 2 a x + b + a = y + x + 6 x + 9
(%i26) expand(2 a);
incorrect syntax: A is not an infix operator
(%i26) incorrect syntax: Premature termination of input at ;.
(%i26) expand(2*a);
(%o26) 2 a
(%i27) expand(y^2-2*b*y+x^2-2*a*x+64-9+6*a-y^2-x^2-6*x-9);
(%o27) - 2 b y - 2 a x - 6 x + 6 a + 46
(%i28) tentex(%);
$$-2\,b\,y-2\,a\,x-6\,x+6\,a+46$$
(%o28) false
(%i29)
\]
3:14 PM
(%i29) expend(2*b*y = - 2*a*x - 6*x + 6*a + 46);
(%o29) expend(2 b y = - 2 a x - 6 x + 6 a + 46)
(%i30) tentex(%);
$${\it expend}\left(2\,b\,y=-2\,a\,x-6\,x+6\,a+46\right)$$
(%o30) false
(%i31) expand(2*b*y = - 2*a*x - 6*x + 6*a + 46);
(%o31) 2 b y = - 2 a x - 6 x + 6 a + 46
(%i32) tentex(%);
$$2\,b\,y=-2\,a\,x-6\,x+6\,a+46$$
(%o32) false
\]
(%i33) %i8
;
2 2
(%o33) expand(b + (a - 3) = 64)
(%i34) %o8;
2 2
(%o34) b + a - 6 a + 9 = 64
(%i35) expand((2*b*y)^2 = (- 2*a*x - 6*x + 6*a + 46)^2);
2 2 2 2 2 2 2 2
(%o35) 4 b y = 4 a x + 24 a x + 36 x - 24 a x - 256 a x - 552 x + 36 a
+ 552 a + 2116
(%i36) b^2 := 64 - a^2 + 6*a - 9;
define: argument cannot be an atom or a subscripted memoizing function; found:
2
-- an error. To debug this try: debugmode(true);
(%i37) expand(4*(64 - a^2 + 6*a - 9)*y*y - (- 2*a*x - 6*x + 6*a + 46)^2);
2 2 2 2 2 2 2 2 2
(%o37) - 4 a y + 24 a y + 220 y - 4 a x - 24 a x - 36 x + 24 a x
2
+ 256 a x + 552 x - 36 a - 552 a - 2116
(%i38) %o2;
2 2
(%o38) 256 y + 112 x - 1792
(%i39) solve([%o2], [y]);
2 2
sqrt(7) sqrt(16 - x ) sqrt(7) sqrt(16 - x )
(%o39) [y = - ---------------------, y = ---------------------]
4 4
(%i40) solve([%o2], [y*y]);
2
2 7 x - 112
(%o40) [y = - ----------]
16
(%i41) solve([%o2], [y^2]);
2
2 7 x - 112
(%o41) [y = - ----------]
16
(%i42) %o41;
2
2 7 x - 112
(%o42) [y = - ----------]
16
(%i43) expand(4*(64 - a^2 + 6*a - 9)*%o41 - (- 2*a*x - 6*x + 6*a + 46)^2);
2 2 2 2 2 2 2 2 2
(%o43) [- 4 a y + 24 a y + 220 y - 4 a x - 24 a x - 36 x + 24 a x
2
+ 256 a x + 552 x - 36 a - 552 a - 2116 =
2 2 2 2
9 a x 69 a x 529 x 2 2
- ------- - ------- - ------ + 24 a x + 256 a x + 552 x - 64 a - 384 a - 576]
4 2 4
(%i44) expand(4*(64 - a^2 + 6*a - 9)*(-1/16)*(7*x^2-112) - (- 2*a*x - 6*x + 6*a + 46)^2);
2 2 2 2
9 a x 69 a x 529 x 2 2
(%o44) - ------- - ------- - ------ + 24 a x + 256 a x + 552 x - 64 a - 384 a
4 2 4
- 576
(%i45) solve([%o44], [x]);
16 a + 48
(%o45) [x = ---------]
3 a + 23
(%i46) solve([y^2=(-1/16)*(7*((16*a+48)/(3*a+23))^2-112)], [y^2]);
2
2 49 a - 294 a - 2695
(%o46) [y = - --------------------]
2
9 a + 138 a + 529
(%i47) x:=%o44;
define: argument cannot be an atom or a subscripted memoizing function; found:
x
-- an error. To debug this try: debugmode(true);
(%i48) x:=(16*a+48)/(3*a+23);
define: argument cannot be an atom or a subscripted memoizing function; found:
x
-- an error. To debug this try: debugmode(true);
(%i49) x
;
(%o49) x
(%i50) x(a):=(16*a+48)/(3*a+23);
16 a + 48
(%o50) x(a) := ---------
3 a + 23
(%i51) x^2(a);
incorrect syntax: Syntax error
(%i51) incorrect syntax: Too many )'s
(%i51) incorrect syntax: Premature termination of input at ;.
(%i51) (x(a))^2;
2
(16 a + 48)
(%o51) ------------
2
(3 a + 23)
(%i52) y2(a):=(-1/16)*(7*((16*a+48)/(3*a+23))^2-112);
- 1 16 a + 48 2
(%o52) y2(a) := --- (7 (---------) - 112)
16 3 a + 23
(%i53) x2(a):=(x(a))^2;
2
(%o53) x2(a) := x (a)
(%i54) plot2d([parametric, sqrt(x2), sqrt(y2), [a,-5,11], [nticks,2000]]);
set_plot_option: unknown plot option: a
-- an error. To debug this try: debugmode(true);
(%i55) plot2d([parametric, sqrt(x2(a)), sqrt(y2(a)), [a,-5,11], [nticks,2000]]);
set_plot_option: unknown plot option: a
-- an error. To debug this try: debugmode(true);
(%i56) x2(a)/16+y2(a)/7;
2
7 (16 a + 48)
-------------- - 112
2 2
(16 a + 48) (3 a + 23)
(%o56) -------------- - --------------------
2 112
16 (3 a + 23)
(%i57) expand(%o56);
2
256 a 1536 a 2304
(%o57) ---------------------- + ---------------------- + ----------------------
2 2 2
144 a + 2208 a + 8464 144 a + 2208 a + 8464 144 a + 2208 a + 8464
2
16 a 96 a 144
- ------------------ - ------------------ - ------------------ + 1
2 2 2
9 a + 138 a + 529 9 a + 138 a + 529 9 a + 138 a + 529
(%i58)
about that magic number:
圆心m到p(-3,0) 的距离为r
圆心m到圆心c距离为:R-r.
MP+MC=R=8
所以M的轨迹为 以（-3,0）和（3,0）为焦点的椭圆.
x^2/16 + y^2/7 =1
3:32 PM
I thought I have found the answer at %i50 but it's not. Let's see what's the wrong.
(a-3)^2 + b^2 &= 8^2 \\
\sqrt{(x-a)^2+(y-b)^2} &= \sqrt{(x+3)^2+y^2} \\
&= 8-\sqrt{(x-3)^2+y^2}
\end{aligned} \right. \]
(%i8) expand((a-3)^2 + b^2 = 8^2)
;
2 2
(%o8) b + a - 6 a + 9 = 64
\]
(%i1) expand((x+3)^2+y^2-(8-sqrt((x-3)^2+y^2))^2)
;
2 2
(%o1) 16 sqrt(y + x - 6 x + 9) + 12 x - 64
(%i2) expand(16*16*(y^2+x^2-6*x+9)-(64-12*x)^2)
;
2 2
(%o2) 256 y + 112 x - 1792
(%i41) solve([%o2], [y^2]);
2
2 7 x - 112
(%o41) [y = - ----------]
16
\sqrt{(x+3)^2+y^2} &= 8-\sqrt{(x-3)^2+y^2} \\
\Longrightarrow 16\,\sqrt{y^2+x^2-6\,x+9}+12\,x-64 &= 0 \\
\Longrightarrow 256\,y^2+112\,x^2-1792 &= 0 \\
\Longrightarrow y^2 &= -{{7\,x^2-112}\over{16}}
\end{aligned} \]
(%i10) expand((x-a)^2+(y-b)^2 - (x+3)^2+y^2);
2 2 2
(%o10) 2 y - 2 b y - 2 a x - 6 x + b + a - 9
\sqrt{(x-a)^2+(y-b)^2} &= \sqrt{(x+3)^2+y^2} \\
\Longrightarrow (x-a)^2+(y-b)^2 &= (x+3)^2+y^2 \\
\Longrightarrow 2\,y^2-2\,b\,y-2\,a\,x-6\,x+b^2+a^2-9 &= 0
\end{aligned} \]
4:04 PM
I'll do it next week. Do course-selection spider first.
2017-11-11 Sa Oct 消参的更多相关文章
- 2017-11-04 Sa Oct 消参
2017-11-04 Sa $ P(-3, 0) $ 在圆C $ (x-3)^2 + y^2 = 8^2 $ 内,动圆M与圆相切且过P点,求M点轨迹. 设切点 $ A(a, b) $,圆心 \(M(x ...
- 2017-11-10 Fr Oct 消参
2017-11-04 Sa $ P(-3, 0) $ 在圆C $ (x-3)^2 + y^2 = 8^2 $ 内,动圆M与圆相切且过P点,求M点轨迹. 设切点 $ A(a, b) $,圆心 \(M(x ...
- 2017年11月Dyn365/CRM用户社区活动报名
UG是全球最大Dynamics的用户组织,由最终用户自发组织,由行业有经验的专家自愿贡献知识和经验的非营利机构,与会人员本着务实中立的态度,不进行推介产品,服务以及其他营销行为.在美国,微软Dynam ...
- WPS 表格筛选两列相同数据-完美-2017年11月1日更新
应用: 1.选出A列中的数据是否在B列中出现过: 2.筛选出某一批序号在一个表格里面的位置(整批找出) 3.其实还有其他很多应用,难描述出来... ... A列中有几百的名字,本人想帅选出B列中的名字 ...
- 2017年11月GitHub上最热门的Java项目出炉
2017年11月GitHub上最热门的Java项目出炉~ 一起来看看这些项目你使用过哪些呢? 1分布式 RPC 服务框架 dubbohttps://github.com/alibaba/dubbo S ...
- 2017.11.11 B201 练习题思路及解题方法
2017.11.11 B201 练习题思路及解题方法 题目类型及涵盖知识点 本次总共有6道题目,都属于MISC分类的题目,涵盖的知识点有 信息隐藏 暴力破解 音轨,摩斯电码 gif修改,base64原 ...
- 【主席树维护mex】 【SG函数递推】 Problem H. Cups and Beans 2017.8.11
Problem H. Cups and Beans 2017.8.11 原题: There are N cups numbered 0 through N − 1. For each i(1 ≤ i ...
- NOIp 11.11/12
最后一场比较正式的NOIp模拟赛,写一发小总结.题目没什么好说的,大部分很简单,先贴一下代码. 1111 T1 //string //by Cydiater //2016.11.11 #include ...
- 11.11光棍节工作心得——github/MVP
11.11光棍节工作心得 1.根据scrum meeting thirdday中前辈的指导进行学习 我在博客中贴了链接,竟然TrackBack引来了原博主,
随机推荐
- 第一章 XAML概览
1.1 XAML是什么? XAML是eXtensible Application Markup Language的英文缩写,相应的中文名称为可扩展应用程序标记语言.也就是说在开发一个应用程序时,我们可 ...
- 剑指offer 7. 递归和循环 斐波那契数列
题目描述 大家都知道斐波那契数列,现在要求输入一个整数n,请你输出斐波那契数列的第n项(从0开始,第0项为0). n<=39 简简单单 废话不多说,直接上代码: public class Sol ...
- 【java】异常
异常体系:Throwable ---Error 通常出现重大问题:如类不存在或者内存溢出.停止运行,不对代码进行处理 ---Exception 在运行时出现的情况,可以进行针对性处理.可以通过try ...
- 筛选最小值---verilog
筛选最小值---verilog `timescale 1ns / 1ps /////////////////////////////////////////////////////////////// ...
- sequelize 中文文档
https://demopark.github.io/sequelize-docs-Zh-CN/
- 7.6.2 break 语句
7.6.2 break 语句 程序执行到循环中的break语句时,会种植包含它的循环,并继续执行下一阶段. 如果break语句位于嵌套循环内,它只会影响包含它的当前循环. break还可用于因其他原因 ...
- 7.2.4 else与if配对
规则是,如果没有花括号,else与离它最近的if匹配,除非最近的if被花括号括起来. 注意:要缩进"语句","语句"可以是一条简单语句或复合语句. 记住,编译器 ...
- TableVie优化方法和优化机制
1,UITableview的优化方法(缓存高度,异步绘制,减少层级,hide,避免离屏渲染) (1)缓存高度:当我们创建frame模型的时候,计算出来cell的高度的时候,我们可以将cell的高度缓存 ...
- 提问：MicrosoftUnderlying input stream returned zero bytes
报错信息:MicrosoftUnderlying input stream returned zero bytes 报错截图: 查阅资料后,提示 jdbc的bug,不能将一个"NULL&qu ...
- 适用于nodercms的打包构建脚本
背景 最近自己用nodercms搭建了一个简单的博客系统,用户发布一些自己谁便谢谢的文章.感谢nodercms团队,这个cms轻量易用,用于做个人博客太方便了.开发了博客系统,肯定设计到部署到AWS或 ...