Loading [MathJax]/jax/output/HTML-CSS/jax.js

2022年5月31日 星期二

用矩陣方法處理橢圓旋轉後的方程式

目錄

1  矩陣的轉置
2  矩陣轉置的乘法
3  內積運算用矩陣乘法表達
4  橢圓的標準式
5  旋轉矩陣及其逆
6  旋轉橢圓的方程式

1  矩陣的轉置

        考慮域K 上的m×n矩陣A:

A=[a11a1nam1amn].

其轉置AT是域K 上的n×m矩陣:

AT=[a11am1a1namn].

若命AT=B,於是有

bij=aji,(1in,1jm.)

其中bij意味著矩陣B中的(i,j)位置元素,而aji是矩陣A中的(j,i)位置元素。

例1

[123456]T=[142536].

(例題終了)

2  矩陣轉置的乘法

定理1在域K上,設Am×n矩陣、Bn×k矩陣,於是有(AB)T=BTAT

[]. 命AB=C,(AB)T=D,BT=E,AT=F,於是

dij=cji=Ar(j)Bc(i)=ns=1ajsbsi=ns=1fsjeis=ns=1eisfsj=Er(i)Fc(j).

其中Ar(i)=[ai1ain]Bc(j)=[b1jbnj],餘類推。

(證明終了)

例2:取

A=[123456],B=[101].

於是

AT=[142536],BT=[101],

然後

AB=[123456][101]=[22]

因此

(AB)T=[22]T=[22].

BTAT=[101][142536]=[22]

所以有

(AB)T=BTAT.

(例題終了)

3  內積運算用矩陣乘法表達

        考慮域K上的n維內積空間Kn,對於其中任意兩個向量vw

v=[v1vn],w=[w1wn],

定義其內積(inner product)為

vw=v1w1++vnwn.

如果我們將域K視作1×1矩陣空間,於是有

vw=[v1w1++vnwn]=[v1vn][w1wn]=vTw.

定理2對於域K上的n維內積空間Kn中的任意兩個向量vw,有

vw=vTw.

4  橢圓標準式的矩陣表述

        在R2中,二元二次方程式

x2a2+y2b2=1

稱為「中心在原點的橢圓的標準式」。更精細地說,用方程式x2a2+y2b2=1定義之橢圓Γ的集合描述為

Γ={P=(x,y)R2x2a2+y2b2=1}.

現在,我們重新改寫橢圓標準式,改用矩陣語言來表示:

x2a2+y2b2=xxa2+yyb2=[xy][xa2yb2]=[xy][1a2001b2][xy]=[xy]T[1a2001b2][xy].

OP=[xy]=vD=[1a2001b2],便可將橢圓標準式寫為

vTDv=1.

(其中矩陣取名為D之理由為該矩陣是對角矩陣,diagonal matrix)於是橢圓Γ的集合描述可改寫為

Γ={vR2vTDv=1,D=[1a2001b2]}.

例3:橢圓x24+y29=1的矩陣方程式為

[xy][140019][xy]=1.

(例題終了)

例4:橢圓3x2+7y2=1的矩陣方程式為

[xy][3007][xy]=1.

(例題終了)

5  旋轉矩陣及其逆

        在R2上,以原點O為旋轉中心,將向量v循逆時針方向旋轉角度θ的變換矩陣為

Rθ=[cosθsinθsinθcosθ].

而此變換的逆的幾何意義就是將向量v循逆時針方向旋轉角度θ,因此可得到

R1θ=Rθ=[cos(θ)sin(θ)sin(θ)cos(θ)]=[cosθsinθsinθcosθ].

再根據轉置的定義,可發現其實[cosθsinθsinθcosθ]=[cosθsinθsinθcosθ]T,亦即有

R1θ={RθRTθ.

所以在計算旋轉矩陣的逆的時候,不需要引用逆矩陣公式,直接取其轉置即可。

例5R30=[cos30sin30sin30cos30]=[32121232],於是R130=RT30=[32121232]T=[32121232]

(例題終了)

6  旋轉橢圓的方程式

        現在考慮橢圓

Γ={(x,y)R2x2a2+y2b2=1},

其矩陣方程式為vTDv=1v=[xy]D=[1a2001b2]。然後以原點O為中心,逆時針方向旋轉角度θ,命所得到的新橢圓為Γ,於是對於任意(x,y)Γ,必有

(Rθv)TD(Rθv)=1.

然後開始化簡左式,

(Rθv)TD(Rθv)=vT(Rθ)TDR1θv=vT(RTθ)TDR1θv=vTRθDR1θv.

因此Γ上的點必滿足方程式

vTRθDR1θv=1.

        但請注意,我們現在只是證明Γ上的點會滿足此方程式,但不表示滿足此方程式的所有點所構成的集合就是Γ。例如以原點O為圓心、半徑為1,座落於第一、二象限的上半圓,其上每個點的坐標都滿足方程式x2+y2=1,但方程式x2+y2=1所代表的圖形卻是整個圓!

        取方程式vTRθDR1θv=1的任一解v=[xy],命w=Rθv,於是

wTDw=(Rθv)TD(Rθv)=vT(Rθ)TDRθv=vT(RTθ)TDR1θv=vTRθDR1θv=1.

這意味著w=Rθv滿足方程式vTDv=1,也就是說w向量的終點在橢圓Γ上,而wvRθ變換而來,所以v必可由wRθ變換而得,從而可確定v的終點在橢圓Γ上。亦即方程式vTRθDR1θv=1的任一解的終點都會是橢圓Γ上的一點。

        因此我們就能完整確定方程式vTRθDR1θv=1代表旋轉後的新橢圓!

定理3:橢圓Γ:vTDv以原點O為旋轉中心、逆時針旋轉角度θ後的新橢圓的方程式為Γ:vTRDR1v

例6:橢圓x24+y29=1以原點O為中心、旋轉30後,新橢圓的方程式為何?

[] 取v=[xy]D=[140019]R=[cos30sin30sin30cos30]=[32121232],於是

RDR1=[32121232][140019][32121232]=[31144531445314421144].

所以新橢圓方程式為

[xy][31144531445314421144][xy]=1.

展開得

31144x2+103144xy+21144y2=1.

(例題終了)

2022年5月30日 星期一

一道四球相切疊橘子求高度的問題

==問題== 

三個半徑1的球,和一個半徑為31的球疊成兩層在桌上,較小的球放上層,四球兩兩相切,則上層小球的最高點到桌面的距離為何?

==答案== 

53+3

==解答== 

        設大球三顆的球心分別為A,B,C,而設小球的球心為D,於是四點A,B,C,D構成一個三角錐,其中¯AB=¯AC=¯BC=1+1=2,而¯DA=¯DB=¯DC=(31)+1=3

        所求高度即為

小球半徑+D與平面ABC距離+大球半徑.

        現在開始計算D與平面ABC距離。取D在平面ABC上的投影點為G,則GΔABC的重心。設AG¯BC交於M¯BM=¯CM=1¯AM=3¯AG=233,因此

¯DG=32(233)2=53.

所以

所求高度=(31)+53+1=53+3.

(解答終了)

一道空間幾何的三角形面積問題

==問題== 

在正四面體ABCD中,於¯AB,¯AC,¯AD上分別取P,Q,R,已知¯AD垂直於平面PQR,且¯AP=6,試求出ΔPQR的面積。

==答案==

92

==解答==

題目略圖如下:

分析三角錐APQR,其中BAC=60=PAR=QAR,再利用三角形邊長關係得下圖:


於是ΔPQR是邊長為33,33,6的三角形,其中

cosR=32+3222233=13,

sinR=+1cos2R=223.

因此

ΔPQR=123333223=92.

(解答終了)

2022年5月26日 星期四

一道空間三平面關係的問題

==問題==

若點P(1,2,3)為三相異平面E1:a1x+b1y+c1z=0,E2:a2x+b2y+c2z=0,E3:a3x+b3y+c3z=0的共同點,且點Q(2,3,1)為另外三平面F1:a1x+b1y+c1z=d1,F2:a2x+b2y+c2z=d2,F3:a3x+b3y+c3z=d3的共同點,則F1,F2,F3三平面的交線參數式為何?

==解答==

        首先注意題目條件稱E1,E2,E3為三相異平面,所以必定不發生重合。

        接著注意由E1,E2,E3構成的線性方程組{a1x+b1y+c1z=0a2x+b2y+c2z=0a3x+b3y+c3z=0為齊次方程組(Homogeneous system),必有全零解(x,y,z)=(0,0,0),所以三平面的交點至少一個。

        但又由題目條件「P(1,2,3)為三相異平面E1,E2,E3的共同點」可知三平面交點並不唯一。於是根據下圖中所揭示之空間中三平面的關係,可知E1,E2,E3必交於同一直線!

(欲觀看更清晰原圖,請用滑鼠左鍵點擊一下以開啟)

命該直線為L,於是我們可寫出L的方程式為

L:{x=0+1ty=0+2tz=0+3t,tR.

注意其中起點為O(0,0,0),方向向量v=OP=(1,2,3),並且OP同於垂直於三平面E1,E2,E3的法向量n1=(a1,b1,c1),n2=(a2,b2,c2),n3=(a3,b3,c3)

        接著假設F1,F2,F3三平面的交線為K,其方向向量v當然也會同時與三平面F1,F2,F3的法向量垂直,而F1E1的法向量同樣是n1,其餘F2E2F3E3皆然。因此我們可以取v=v=(1,2,3),從而直線K的方程式為

K:{x=2+1ty=3+2tz=1+3t,tR.

其中起點為Q(2,3,1)

(解答終了)

2022年5月24日 星期二

107指考數甲的複數問題

==問題== 

O為複數平面上的原點,並令點A,B分別代表複數z1,z2,且滿足|z1|=2,|z2|=3,|z2z1|=5。若z2z1=a+bi,其中a,b為實數,i=1。試選出正確的選項。

(1)  cosAOB=23

(2)  |z2+z1|=23

(3)  a>0

(4)  b>0

(5)  設點C代表z2z1,則BOC可能等於π2

[107,指考,數甲]

==答案==

(1), (3), (5)

==解答==

        由|z1|=2可假設z1=2(cosθ1+isinθ1)

        由|z2|=3可假設z2=3(cosθ2+isinθ2)

        於是

z2z2=(3cosθ22cosθ1)+i(3sinθ22sinθ1).

|z2z1|=5

(3cosθ22cosθ1)2+(3sinθ22sinθ1)2=5.

左右平方得

(3cosθ22cosθ1)2+(3sinθ22sinθ1)2=5.

對左式展開得

9cos2θ212cosθ2cosθ1+4cos2θ1+9sin2θ212sinθ2sinθ1+4sin2θ1=5.

整理為

9(cos2θ2+sin2θ)12(cosθ2cosθ1+sinθ2sinθ1)+4(cos2θ1+sin2θ1)=5,

912cos(θ2θ1)+4=5,

就得到

cos(θ2θ1)=23.

由平方恆等式cos2φ+sin2φ=1

sin(θ2θ1)=±1(23)2=±53.

z2z1=a+bi

3(cosθ2+isinθ2)2(cosθ1+isinθ1)=a+bi.

利用複數極式的除法關係得

32[cos(θ2θ1)+isin(θ2θ1)]=a+bi.

因此

a=32cos(θ2θ1)=3223=1,

b=32sin(θ2θ1)=32(±53)=±52.

        下面開始討論每一個選項。

        (1) 正確。

z2=z1z2z1=z1(1±52i)=z132(23±53)可知點B是由點A以原點O為旋轉中心、逆/順時針旋轉arccos23後,再伸縮32倍而得。注意此時點B的位置有2種可能。無論如何,必有cosAOB=23

        (2) 錯誤。

因為

|z2+z1|=|3(cosθ2+isinθ2)+2(cosθ1+isinθ1)|=|(3cosθ2+2cosθ1)+i(3sinθ2+2sinθ1)|=(3cosθ2+2cosθ1)2+(3sinθ2+2sinθ1)2=9+12cos(θ2θ1)+4=21.

       (3) 正確。

因為a=32

       (4) 錯誤。

因為b=±52

       (5) 正確。

因為z2=z2z1z1=z2z12(cosθ1+isinθ1),所以可知點B是由點C以原點O為旋轉中心、逆時針旋轉θ1後,再伸縮2倍而得。設θ1的最小正同界角為θ+1,則BOC=θ+12πθ+1。特別當我們取θ1=π2時,就有BOC=π2。因此BOC確實可能等於π2

(解答終了)

==附註==

雖然這是一道幾何問題,但我硬是把解答寫的像是代數問題一樣嚴謹而仔細,沒有放示意圖。我當然可以很幾何地處理這道問題,不過倘若太依賴圖形,選項(5)很容易出問題。我大部分都是幾何學取向,偶爾會突然切換為代數取向,追求極致的邏輯嚴謹。

2022年5月17日 星期二

一道三角形三邊長關係的線性方程組問題

==問題== 

a,b,cΔABC的三邊長,且使方程組{ax+by=cbx+cy=a有無限多組解,則ΔABC是何種三角形?

==解答==

由於已知線性方程組為無限多組解,而由二元一次線性方程組的幾何意義可知,此時兩條方程式可視為同一條直線,因此係數與常數之比例相同,亦即有

a:b=b:c=c:a,

其中依序是x係數比例、y係數比例與常數比例。將比例式改為分數式,

ab=bc=ca.

假定比值為k,得

ab=bc=ca=k.

(注意k為正數相除的結果,所以k必為正數)於是可得

a=bk,b=ck,c=ak.

注意到這三個關係式中,a,b,c各在等號兩邊恰出現1次,所以考慮三數乘積得

abc=(bk)(ck)(ak),

整理得

abc=abck3,

於是解出k=1,故a=b=c,即ΔABC為正三角形。

(解答終了)

2022年5月12日 星期四

一道身高排列問題

==問題==

設有身高不等的10人排成一列,若想讓任一較矮者不夾在兩較高者之間,則排列方法有多少種?

==解答==

        假定身高分別為:0, 1, 2, 3, 4, 5, 6, 7, 8, 9。先寫出幾個可能的情況:

0, 1, 2, 3, 4, 5, 6, 7, 8, 9

9, 8, 7, 6, 5, 4, 3, 2, 1, 0

1, 3, 6, 9, 8, 7, 5, 4, 2, 0

2, 5, 7, 8, 9, 6, 4, 3, 1, 0

(排列組合最重要的就是要舉例!)

這些排列有什麼特徵?

        仔細觀察可以發現,一旦最高的9確定後,接著比較矮的8就要選擇在9的左或右,然後是7選擇在左或右,依此類推。

=方法1=

        依照上面的觀察,我們可以從9開始去「生成」所有可行的排法。例如:

(i) 9 → 89 → 897 → 8976 → 89765 → 489765 → 4897653 → 48976532 → 148976532 → 1489765320

(ii) 9 → 98 → 987 → 6987 → 69875 → 698754 → 3698754 → 36987542 → 136987542 → 0136987542

也就是說,一開始就有一個9,然後放8,而8可以選擇在9的左邊或是右邊;接著放7,而7可以選擇左邊或是右邊;...依此類推。從8開始到0,每個數字都可選擇左邊或右邊,於是算式為

2×2××29=29=512.

=方法2=

        另外一個方法是,首先想像有10個空位:

⬜⬜⬜⬜⬜⬜⬜⬜⬜⬜

接著填入9,有以下可能情況:

(i) 9⬜⬜⬜⬜⬜⬜⬜⬜⬜

(ii) ⬜9⬜⬜⬜⬜⬜⬜⬜⬜

...

(xi) ⬜⬜⬜⬜⬜⬜⬜⬜9⬜

(x) ⬜⬜⬜⬜⬜⬜⬜⬜⬜9

        對於情形(i),我們在9的左邊沒有任何空格可以填數字,而在右方有9個位置可以填數字,於是方法為C90×C99

        對於情形(ii),我們在9的左邊有1個空格可以填數字,而在右方有8個位置可以填數字,於是方法為C91×C88

...

        對於情形(xi),我們在9的左邊有8個空格可以填數字,而在右方有1個位置可以填數字,於是方法為C98×C11

        對於情形(x),我們在9的左邊有9個空格可以填數字,而在右方有0個位置可以填數字,於是方法為C99×C00

        綜上所述,所求方法數為
C90×C99+C91×C88++C98×C11+C99×C00=C90×1+C91×1++C98×1+C99×1=C90+C91++C98+C99=29=512

(這裡引用了二項式定理:(A+B)n=Cn0AnB0+Cn1An1B1+Cn2An2B2++CnnA0Bn,其中代入A=1,B=1,可得Cn0+Cn1+Cn2++Cnn=2n。)

==註記==

        我個人喜歡用方法2,覺得比較自然,也比較好玩(可以跟組合級數扯上關係)。對於方法1那種慢慢生長的方式感到不太適應(雖然是自己想出來的)。

一道求高次多項式除法餘式的問題

        W小姐下午發來一道求多項式除法餘式的問題(不知道是她在得勝者(👎)那邊算到的題目還是啥的),我起初以為用假設餘式的方式、再代入特殊值就可以解出,然而進展不順利。後來我改用二項式定理處理,但是計算量十分大,跟直接寫直式除法沒兩樣,對於考場解題幾乎沒任何用處。我想了三個小時左右,總算勉強搞出來一個考試中可行的方法,這題算是一種特殊的求餘式問題吧!記錄下來,以後或許再碰到就有素材可以用。

==問題==

f(t)=t15+3t10t52,若f(t)除以t4t2之餘式為at3+bt2+ct+d,求abc+d

==解答==

        假設除法商式為q(t),於是可寫出恆等式

t15+3t10t52=(t4t2)q(t)+at3+bt2+ct+d.

整理為

t15+3t10t5at3bt2ct(2+d)=(t4t2)q(t),

然後

t15+3t10t5at3bt2ct(2+d)=t2(t1)(t+1)q(t).

注意到右式中有個因子為t2,所以左式必為t2的倍式,從而可判定左式沒有1次項與常數項,故得到

c=0,2+d=0,

c=0,d=2

        重新再寫一遍式子,

t15+3t10t5at3bt2=t2(t1)(t+1)q(t).

根據多項式乘法的消去律,左右兩端同時消去t2,得到

t13+3t8t3atb=(t1)(t+1)q(t).

接著分別代入t=1t=1,得

{1+31ab=01+3+1+ab=0

解出a=0,b=3。因此所求

abc+d=030+(2)=5.

(解答終了)

==註記==

題目被除式中的未知數的次數15, 10, 5實在很有迷惑性,一直讓我想要使用變數變換z=t5,但發現這樣變換後沒什麼鳥用。也許是我沒看出題目要害,所以本該這樣變換但沒走到正確的路。隨便啦,都耗了我三個小時,解出來就好。

2022年5月3日 星期二

一道排列組合的問題(2020,印度理工學院,Main, 5 Sep II)

 ==問題== 

某份試卷分為3節,每節都有5道試題,考生只需在全卷中選擇5道試題作答,但規定每節都至少要選擇1題作答,請問作答方法數為

(a) 3000        (b) 1500        (c) 2255        (d) 2250

(原文)

There are 3 sections in a question paper and each section contains 5 questions. A candidate has to answer a total of 5 questions, choosing at least one question from each section. Then the number of ways, in which the candidate can choose the questions, is

(a) 3000        (b) 1500        (c) 2255        (d) 2250

==解答==

首先假定考卷所分的3節分別為A, B, C,而其中各節所含的題目又可以編號為A1, A2, ..., A5, B1, B2, ..., B5, C1, C2, ..., C5。

我們可以寫出幾種可能的作答方式,如:

A1, A3, B4, B5, C1,

A1, B5, C1, C3, C5.

等等。在處理排列組合問題時,我們必須搞題目問的是什麼,要能夠在心中明確寫出幾個題目所求的具體結果。

暫且先忽略具體選的是哪幾題,只注意所選的題目是從哪裡選出、以及選了多少題,那麼所有作答方式可分為兩大類:

(甲) 某節選3題、某節選1題、某節選1題;

(乙) 某節選2題、某節選2題、某節選1題。

對於(甲)而言,顯然又可再細分3類:

A選3題,B選1題,C選1題

B選3題,A選1題,C選1題

C選3題,A選1題,B選1題

(先處理哪部分選3題比較好做)因此能計算出(甲)的方法數一共有C53×C51×C51×3=750種方法。

對於(乙)而言,顯然又可再細分3類:

A選1題,B選2題,C選2題

B選1題,A選2題,C選2題

C選1題,A選2題,B選2題

(先處理哪部分選1題比較好做)因此能計算出(甲)的方法數一共有C51×C52×C52×3=1500種方法。

綜上所述,所求方法數為750+1500=2250種,選(d)。

(解答終了)