橢圓的生成課件

＞橢圓的角度微分法

橢圓的角度微分法是利用橢圓的內(nèi)接多邊形來(lái)逼近橢圓

設(shè)橢圓的參數(shù)方程為

x=acosg

(1-31)

Iy=bsin^p

其中a》0〉0,旋轉(zhuǎn)角中的起始角、終止角分別為*/?,且滿足關(guān)系式0&a〈p&360

當(dāng)把橢圓的圓周角〃等分之后,橢圓的內(nèi)接多邊形所對(duì)應(yīng)的中心角Ag為

△W=2TT/〃(1-32)

此時(shí)第i條直線的兩端點(diǎn)坐標(biāo)分別是

(i=O,l(1-33)

aCOSQ+1=x,cosA^

(1-34)

+J—ycosAw

-

cosA*卡sin△中j.-

L#」

卜n△卓區(qū)q中

劣-1

y」Ml」

sin△中cosAg

a..5

IccsA^p一coM中―齊inA0

/sinAF

iyi+i'■加-1■^sin△中cos△中ucisA?I;;]

oos△中一號(hào)win△平OOsAfTJi-sinA??

cosA勺-?sinMcusA手

0

=2

ir)sAyJ』

1011*「

=2oosA?0JLv.-'

-(2cceA$>),Jr-JT,-I

Kf+1l(f=1,2,…，n-I)(1-35)

ly+i二(2co4中)“-y!

a^b時(shí).取H=0.5a+20*

更為一般的橢圓是橢圓的長(zhǎng)軸與坐標(biāo)系的X軸具有一夾角6,

見(jiàn)圖1-8。

為了討論方便，令其留心坐標(biāo)為（0，。），利用旋轉(zhuǎn)公式

IJC=、/?cos。~-sin?

(1-36)

ly=j：r-SLK0-y'ctxO

可得此時(shí)橢圓的參數(shù)方程

x="?cosgyus。~〃，sing,sin9

(1-37)

y=a,cos*,sinU+b?sin儀,cos5

令△伊=2兀/%其中〃為其內(nèi)接多邊形的邊數(shù),旋轉(zhuǎn)角儀從起始角

a旋轉(zhuǎn)至終止角,8處，這時(shí)其內(nèi)接多邊形的第：點(diǎn)坐標(biāo)是

<f>t=a?,平.i

<Xj~a'cos5,cas(fi-b,sin。?sin?(i=0,l,2,…，4)(1-38)

(y,-a?sing,cos傷+b-cos8,sinR

其第￡+1點(diǎn)坐標(biāo)是

Jx；+1=(a?oas^)*cos(<p?A夕)~(^,sint?)*sin(^>-+4g)

tJ(1-39)

(汕」i=(a(衿+△—+(b?axif?),sin($>,?+Acp)

Ico?(仍+△$?)=cos△g'COS物-sinAg,sin價(jià)

而|sin(<f>i?>甲、=cos△3,sin^pr+InAg'cos弘

若設(shè)A-a?co姐，B=/Lsin6,C=a，sinO,D=?oose

EccsAa.E-sind中、S：=cos%,Rsin仍

則上述遞推公式可簡(jiǎn)化為

5“=E￡-CT.

力?i=EF+F?S；

,_儀口，(f=0,1?…，”-I)(1-40)

工,-i-八A6.[-B?rr(,?

,??i=C?5：ji+D?7/

注意此時(shí)

So=c

T&=sina

(1-41)

xnA-&B,/'1)

No=C'-S(|?D?丁o

可見(jiàn)斷轉(zhuǎn)后的確圓其遞推公式并不簡(jiǎn)單，它僅是把大量的三角函數(shù)運(yùn)算轉(zhuǎn)化成實(shí)數(shù)運(yùn)算

而已,但它卻是生成橢mi的基本公式之一;

設(shè)橢圓的圓心位于坐標(biāo)系原點(diǎn)，我方程為

111T桶園的對(duì)稱性.因此只需時(shí)論在坐標(biāo)系第-象限中橢園?。槙r(shí)針?lè)较蛐D(zhuǎn)）的生成即nJ

達(dá)到畫一個(gè)完整橢圓的目的

對(duì)橢回方程進(jìn)行微分，并令--1.得N點(diǎn)坐標(biāo)

（元仔2*'潦1玄）.見(jiàn)圖J以根據(jù)橢網(wǎng)的性質(zhì).

對(duì)于俞橢㈱弧.顯然八〉0,△.、，＜0,且誓2?I.所以

W,\ANI.這就是說(shuō)為了逼近而弧.各點(diǎn)的坐標(biāo)步

進(jìn)規(guī)律整X坐標(biāo)每次推算加1,而坐標(biāo)利會(huì)推算是圖1-9四分之一橢圓弧勺具分界點(diǎn)N

否減】由偏差判別式來(lái)確定。同理,逼近母弧的冬

點(diǎn)坐標(biāo)步進(jìn)規(guī)律是了坐標(biāo)每次推算減I.而*也標(biāo)是否加1中一偏差判別式來(lái)確定，

現(xiàn)對(duì)橢I劇方程進(jìn)仃適當(dāng)變形井?把它定義成函數(shù)4（1.））.有

4（1，3＞）色62比2+?2/_.“2,/（142）

可以發(fā)現(xiàn)：

當(dāng)點(diǎn)3”，）位于橢圓之外時(shí)，函數(shù)〃（7,、＞＞0；

當(dāng)點(diǎn)（工4?）位于橢圓之上時(shí)，函數(shù)d（*G）=0；

當(dāng)點(diǎn)Q"）位于橢圓之內(nèi)時(shí)，函數(shù)＜/（.r,,y）＜0.

首先解決八N弧封點(diǎn)偏差的遞推計(jì)算問(wèn)題.

設(shè)第i-I點(diǎn)從病弧上方遇近煙潮I,根據(jù)中點(diǎn)偏差判別法則.這意味后第?：-1點(diǎn)的中點(diǎn)

偏差￡1<0.因此可得第？-1點(diǎn)與第r點(diǎn)的中點(diǎn)坐標(biāo),見(jiàn)圖I-10(a)o那么第i-1點(diǎn)與第r

點(diǎn)的中點(diǎn)偏M分別為

di?-<1(M,?尸t<r|MI2

dt一(I(M)=Z>?2；1+1產(chǎn)+a°a'b2

故&=dj-l4?4工，:『(:一？，??”)(I47)

設(shè)第￡1點(diǎn)從寂弧下方遇近橢圓，根據(jù)中點(diǎn)偏差判別法則，這意味有乩I》。，因此可

得第i-1點(diǎn)與第i點(diǎn)的中點(diǎn)坐標(biāo)，見(jiàn)圖l”O(jiān)(b)n那么第點(diǎn)與第i點(diǎn)的中點(diǎn)偏景分別為

44ytrb-

d,—b'(.；r,-i+1尸+a].y；?

2

故di=dt?la''y,?+2b''x,?+6(/=2,(148)

需要說(shuō)明：對(duì)于U-47)與EL4于兩式中的2川.者_(dá)］運(yùn)算,由于它發(fā)生于中點(diǎn)偏差遞推升

算的每步，因此可利用在速報(bào)過(guò)程中底次彳都加1這一機(jī)會(huì),迭加格城以實(shí)現(xiàn)該運(yùn)算;

而對(duì)于2“Z.y1運(yùn)算.由于它僅發(fā)生在4RO.即》,=丫一I期間，因此可以利用當(dāng).y斌1

時(shí),從常址2個(gè)、”(注意常員可能是一個(gè)巨型性數(shù))中遞減常量2“？以實(shí)現(xiàn)該運(yùn)算經(jīng)

過(guò)這樣處理之后.遞推公式(1-47)'j(148)就只剩下簡(jiǎn)單的加減法運(yùn)弟了

AN強(qiáng)的初始條件：

由于梆圓晶從八（0.8）點(diǎn)處開(kāi)蛀州點(diǎn)，所以.O.M,一〃.㈣第一點(diǎn)的中點(diǎn)坐標(biāo)為MU1.

621相應(yīng)俏三小為

出二標(biāo)（1尸+儀4一告）

=Z>：a'h,n7Z4

對(duì)切，弧也挑行類似的處理并有納論：

當(dāng)比iXiat

當(dāng)M,?)n-!

?/)~d(i,2Zi"?.1t2M:?、？?tu(1-51)

NB啊勺初始條竹t

顯然茄，弧的起點(diǎn)應(yīng)從X點(diǎn)處開(kāi)始，僧、?點(diǎn)的坐標(biāo)計(jì)算發(fā)朵目為耶惟數(shù).使川極為小使.

可丐虛

芝(/尸+a2yz-(rfj2)~l>

科瓠第

4y___2bzjmIfRI-A'IMHVJfi

dr2a-y

令&;「2〃'/.d_y=2“'*.用有帖論："fAx—Ay

時(shí)，此時(shí)的.r.,V坐硬應(yīng)于N點(diǎn).；/A.r<Ay時(shí)，此時(shí)

的1?一坐料理應(yīng)上而如；當(dāng)Ar>生時(shí).此時(shí)的」q，

弁一郎，森其\e>9l

生你對(duì)應(yīng)「工6強(qiáng)乂因2方。燈2a〃的什尊在M菽

犯時(shí)已解決，皿里只常簡(jiǎn)單比較一卜兩部的大小就健

達(dá)到糾斷所摭京派是否到達(dá)或雄過(guò)N點(diǎn)：出III倘定YHflt附切蠟條件

NB弧中點(diǎn)偏左初始的處理如F-

卿i-1點(diǎn)仍同余版,技眾弧的遞推公式推一

口』1點(diǎn)已超過(guò)N點(diǎn).因此應(yīng)把第，'點(diǎn)的中心M”的儡冷轉(zhuǎn)換成3%起點(diǎn)的中點(diǎn)以的就

卷.如陽(yáng)LI1所示所以勺

M.融途：d,~b'[^xi?v4)十。'(乂17)'-a.(I52)

.%點(diǎn)偏Zhd'「人彳gi-I)'+u'(y,「J)a'b2(t?53)

故</,(7,*3((.''/r)/4-i：/?1??；?ir??*v,)

二￡廠卜3((？-。2)/4(Ar-Ay)/2(|-54)

當(dāng)把式U-54)中的用作為篇弧巾點(diǎn)偏差的初始化并運(yùn)用(150).(1-51),(145).(l-46)?q

式推笠.就旋到達(dá).0)點(diǎn)處，從而'完成整個(gè)；而弧的繪制

1.4自由曲線的生成

所謂自由曲線，可以把它理解為曲線的形狀可由人們

隨意指定多個(gè)控制點(diǎn)（型值點(diǎn)、插值點(diǎn)）而確定。為了避免

高次曲線計(jì)算復(fù)雜，性能不穩(wěn)定等缺陷，通常的做法是采

用分段的低次曲線來(lái)構(gòu)造整條曲線，這樣的曲線稱為樣條

曲線。

自由曲線分兩類：

一類用于曲線的數(shù)值分析，它要求曲線準(zhǔn)確地通過(guò)插

值點(diǎn)，這類曲線主要有三次樣條曲線與三次參數(shù)樣條曲線

等；

另一類用于曲線的綜合，它不要求曲線一定通過(guò)其控

制點(diǎn)，這類曲線主要有貝齊埃曲線與B樣條曲線等，它們

主要用于曲線、曲面的形狀設(shè)計(jì)中。

-次樣條函數(shù)的基本內(nèi)容可表述如.

已知〃個(gè)型值點(diǎn)P（r,.x）,i=12…，，八且口<4<一<.（…苦>=53）滿足下列

條件：

①型僮點(diǎn)在函數(shù).丫=凱工）上；

②SQ）在整個(gè)伏皿叫以上二次連續(xù)可導(dǎo)；

③在每個(gè)于國(guó)網(wǎng)［小T，/《；-1,27）l：,SCr）都是」的三次多項(xiàng)式"

則稱SJ）是過(guò)型然點(diǎn)的三次樣條函數(shù)，由三次樣條函數(shù)構(gòu)成的曲線稱三次樣條曲線

由上述定義可知，5⑺在每個(gè)F區(qū)間二都站三次多期式，因此.第，段的、s.

（龍）可寫成

S,（x）=a,I/>,（?<?-^7）+C（J'.r,）2+C/,（J:.z；）1（I55）

其中.為特定系數(shù)二分段函數(shù).S（r）與型值點(diǎn)/”,「之間的關(guān)系見(jiàn)圖112.

現(xiàn)用型值點(diǎn)處的二階導(dǎo)數(shù)來(lái)表達(dá)三次樣條函數(shù)二

由（155）式叮求得分段函數(shù)%（丁）的一階一階號(hào)數(shù)如F:

2

S1（M）=b,+2cf（.ra,）+3（/f（x,x；）

,S”,（i）=2*+6d,（f-&）

由于股值點(diǎn)匕（u,，M）在函數(shù)$上，即y=S（r,）.故

圖1-12分段函數(shù)總⑴與取值

.=加

點(diǎn)匕“」1之間的關(guān)系

考慮點(diǎn)處的二階導(dǎo)致去達(dá)式,有年（七）二2仃.今此時(shí)分

段函數(shù)S.{.r）仁型強(qiáng)點(diǎn)P,處的二階汴數(shù)為M,是另?待定系數(shù)，有

由于函數(shù)S"）二階導(dǎo)數(shù)連續(xù)，即第I段分段函數(shù)Sf在其終點(diǎn)處的二階導(dǎo)數(shù)等于第j+I

段分段函數(shù)S.7在其起點(diǎn)處的二階導(dǎo)數(shù),見(jiàn)圖1-12,宥

,竽式即，1*=離,i（曲+r）

得4?=泰（C什I-。）=蔣2（由4+1-M；）

其中"i=H：11-4c

同理，由于函數(shù)5（H）是連續(xù)的，即第i段分段函數(shù)S；在其終點(diǎn)的函數(shù)值等于第￡+1段

分段函數(shù)字+i在其起點(diǎn)的函數(shù)值，有

把已求得的系數(shù)分別代入上式，得

■三土G+.那）丁斯（當(dāng)七埠U）

則用型值點(diǎn)處二階導(dǎo)數(shù)表達(dá)的三次樣條函數(shù)為

Si⑺=%+看-與（必+^^卜％-的）

十學(xué)?（*一力/+自（Mr（1-56）

現(xiàn)可利用函數(shù)S（z）的一階導(dǎo)數(shù)連續(xù)來(lái)求解各型值點(diǎn)處的二階導(dǎo)數(shù)M：的數(shù)值c

乂=5。（工：）

令k=,%T._&-八__6/%+「yw-i）

Di

'11,-t+h,'內(nèi)兒—+兒?~h1.l+hi\h,一）一.1

則有左?Mj-i+2?M+/vM,+j=D,（1-57）

注意,此時(shí)￡=2,3,…，（附—1）出+修=1。

（1-57）式稱為三次樣條函數(shù)S（H）的“M連續(xù)性方程"，它反映廣力學(xué)上的"三彎矩關(guān)系:

不難發(fā)現(xiàn)（1-57）式共有〃-2個(gè)線性方程，還不能惟一地求解出各型值點(diǎn)處的二階導(dǎo)數(shù)M.a

但如果在實(shí)際應(yīng)用中，根據(jù)實(shí)際問(wèn)題要求給出Mi與或-階導(dǎo)數(shù)等已知邊界條件,即在曲

纜兩端給出約束條件，就能惟一地求解出一組M：參數(shù)，從而惟一地確定過(guò)型位點(diǎn)的三次樣條

函數(shù)。

三次樣條曲線常用的邊界條件有夾持端、自由端、拋物端等三種，現(xiàn)分別敘述如下。

①夾持端：根據(jù)實(shí)際問(wèn)題的需要限定曲線兩端的切線,即已知曲線在始端與終端的一階

導(dǎo)致3與京,從而給出約束M]與M”的約束條件,由于

S|（x1）=y'1

21

故2M|+M2=6（^-y；）Aq（1-58）

又S：

故Mn^2Mn=6（y\-）1hn..,（1-59）

加上這兩個(gè)條件與（I-57）式就能確定出惟一的一組H.G=；,2,3,…切）來(lái)二

②自由端：即曲線在始端與終端的二階導(dǎo)數(shù)為零，它說(shuō)明曲線在始端與終端不受外力約

束,其切線方向僅受該曲線其他型值點(diǎn)的影響而變化。由于已知二階導(dǎo)數(shù)與Mn的大小,

則它與（1-57）式聯(lián)立在一起共有n個(gè)條件,能求解出惟一的一組=為了

使該已知條件與（1-57）式的表達(dá)形式相同，令

[l,M0+2?Mi+0?M2=2M]

10,M…42?M,+1-M/i=2M“（1?）

則Ma-i=0。

③拋物端：認(rèn)為曲線在第1段S|（和第HI段S“|（工）（即末段）為拋物端、也就是此

二段曲線的二階導(dǎo)數(shù)為常數(shù)（非零），因此，令Ji?二“，即有

p-Mn+-2-M,-2Mi0

〕2iM”.廣2?M“+3?A4“，i=0(1-61)

為了簡(jiǎn)化表達(dá)，將這:種邊界約束條件寫成統(tǒng)?表達(dá)式上述M連續(xù)性方

[2,M|+i）i程是由型值點(diǎn)處

(1-62)

["?陸-1+2?1討產(chǎn)。一階導(dǎo)數(shù)連續(xù)而

表11三次樣條曲線的邊界約束條件導(dǎo)出的，它反映

邊界條件九了相鄰型值點(diǎn)處

央持用產(chǎn)1=1.七二1二階導(dǎo)數(shù)的關(guān)系。

如果在推導(dǎo)型值

自由端，“廣。區(qū)=0D「2MI,D“=2M“，其中JW：M..-0

他物端2兒=2D|O.D?0點(diǎn)關(guān)系式時(shí)改用

把這三個(gè)邊界條件與（I-57）式結(jié)合在一起即構(gòu)成?個(gè)完整的線性方程組,可川矩陣表示s”M（x尸

如下

作為條件，可以

2”[Mi-Dr

類似地推導(dǎo)出反

入22/0M5

2映相鄰型值點(diǎn)處

432”3M力3

3zsz(163)

??????……一階導(dǎo)數(shù)關(guān)系的

0h-12/G-1M,..?D,m連續(xù)性方程式。

兒

L2,一M“.從

山JYI-63）式是一個(gè)主對(duì)角線元素為2且占優(yōu)的矩陣，因此它存在裱一解

①用n個(gè)型值點(diǎn)a逼近一條已知曲線且滿足其*坐標(biāo)遞增這一條件；

②根據(jù)實(shí)際情況確定曲線的邊界條件；

③求解M.的值，用追趕法求解而不用一般的消元法，這樣能節(jié)省大量的時(shí)間和存儲(chǔ)空

間。

④將求解出的M分段代人s,a）的表達(dá)式中。按照繪圖要求，計(jì)算各段內(nèi)的若干插值點(diǎn),

并依次用直線相連，可畫出所求的三次樣條曲線。

①由于它要求嚴(yán)格保證其…V/，這就決定了該曲線不能繪制具有垂直切線之

類的曲線。

②如果用三次樣條曲線代替人工放樣，當(dāng)曲線的導(dǎo)數(shù)曠區(qū)1時(shí)，這兩者吻合的效果較好;

但當(dāng)其曠|>1時(shí)，雖然三次樣條曲線過(guò)指定的型值點(diǎn)并保證其二階導(dǎo)數(shù)連續(xù)，但它仍與

人工放樣所畫出的曲線有較大的差異。

③當(dāng)曲線中夾有直線時(shí)，如直線與圓弧的接合處，如硬要讓其保證二階導(dǎo)數(shù)連續(xù)，則

將使直線或曲線產(chǎn)生波動(dòng)。解決這個(gè)問(wèn)題的方法是采用分段擬合，找出直線與圓弧之

間的切點(diǎn)和切點(diǎn)處的斜率，把該切點(diǎn)作為型值點(diǎn)，其斜率作為樣條曲線的邊界條件，

分段繪出直線與曲線。顯然此時(shí)切點(diǎn)處只能保證兩者一階導(dǎo)數(shù)連續(xù)，而不能保證二階

導(dǎo)數(shù)連續(xù)。

給定9個(gè)型值點(diǎn)：(0,0),(45,71),(90,

100),(135,71),(180,0),(225,-71),(270,

-100),(315,-71),(360,0),畫出過(guò)這組型值

點(diǎn)的三次樣條曲線。

u?

(

二Q

=_

s二

一U

〕=

m，d

iW5

含3

BEM

AmA

(-)\

二

「二;i

M—u

4上.

「:

二*

”{

3l-l/

乂±B

二

〕《

?

二

5M「

三

￡〕/

〔u

〔-

二/

〕一(G

K>t(3

3〔《

二ql

36x〔

A〔一

E!n

—

￡*十-+

一HU

20〔

q「/

百

二I

二X￡

=￡m《

一1

=(t〔|

K〕-

導(dǎo)

—L

工

?+X

U二w?

=1一26Mm

工

二TI+)

q-XS》I七-

一

二

內(nèi)

—〔

》

=〕「U

LtH*HZl

胃I

y5巳X9/