光電成像器件計(jì)算機(jī)輔助設(shè)計(jì)與總結(jié)

1、光電成像器件CAD報(bào)告 1120101048 陸原介光電成像器件計(jì)算機(jī)輔助設(shè)計(jì)與總結(jié)實(shí)驗(yàn)報(bào)告 84摘要：摘要：在現(xiàn)實(shí)生活工作中，我們經(jīng)常要計(jì)算特定電場(chǎng)中的電位分布。在這一過(guò)程中，我們不斷思考，希望能有一種開(kāi)發(fā)程序，能夠完成其特定功能。該程序應(yīng)運(yùn)而生。本程序的基本算法，基于多種數(shù)學(xué)理論以及C語(yǔ)言基礎(chǔ)的環(huán)境支持，應(yīng)用VS編程環(huán)境，最終完成。在經(jīng)歷這個(gè)過(guò)程后，我們大致了解了程序開(kāi)發(fā)的過(guò)程，從中收獲了經(jīng)驗(yàn)與知識(shí)。在未來(lái)，我們將受益匪淺。關(guān)鍵字：電場(chǎng)分布，程序，計(jì)算機(jī)輔助設(shè)計(jì)Abstract:In our real life,we always need to calculate electric po

2、tential in a special electric field.In this processing,we need think continually and develop a program to finish its function.And this program was born.Base on mathematic and C Language and VS develop environment,I do finished it.In this processing,I know the knowledge of program developing,gaining

3、a lot of experience.Key words: electric potential,program,CAD目錄前言3理論基礎(chǔ)4設(shè)計(jì)思想8主要變量8流程圖9使用說(shuō)明9結(jié)論與討論10致謝11參考文件12源程序13前言：該課程是我光電工程系科學(xué)與技術(shù)專業(yè)主干專業(yè)課程之一的重要組成部分。其主要任務(wù)是結(jié)合比較實(shí)際情況的像管設(shè)計(jì)計(jì)算課題，綜合運(yùn)用光電成像原理課程中的相關(guān)基本原理和其他相關(guān)知識(shí)，將本課程設(shè)計(jì)基本理論所給出的光電成像器件電子光學(xué)系統(tǒng)計(jì)算機(jī)輔助設(shè)計(jì)的物理模型與數(shù)學(xué)模型轉(zhuǎn)化為可以實(shí)際進(jìn)行數(shù)值計(jì)算的程序軟件系統(tǒng)，已達(dá)到培養(yǎng)運(yùn)用綜合知識(shí)的能力，編制，調(diào)試，開(kāi)發(fā)實(shí)際工程軟件的能力，提高運(yùn)

4、用現(xiàn)代設(shè)計(jì)方法與手段及計(jì)算機(jī)應(yīng)用開(kāi)發(fā)能力的目的。計(jì)算機(jī)輔助設(shè)計(jì)(CAD-Computer Aided Design)指利用計(jì)算機(jī)及其圖形設(shè)備幫助設(shè)計(jì)人員進(jìn)行設(shè)計(jì)工作。 在設(shè)計(jì)中通常要用計(jì)算機(jī)對(duì)不同方案進(jìn)行大量的計(jì)算、分析和比較，以決定最優(yōu)方案；各種設(shè)計(jì)信息，不論是數(shù)字的、文字的或圖形的，都能存放在計(jì)算機(jī)的內(nèi)存或外存里，并能快速地檢索；設(shè)計(jì)人員通常用草圖開(kāi)始設(shè)計(jì)，將草圖變?yōu)楣ぷ鲌D的繁重工作可以交給計(jì)算機(jī)完成；由計(jì)算機(jī)自動(dòng)產(chǎn)生的設(shè)計(jì)結(jié)果，可以快速作出圖形，使設(shè)計(jì)人員及時(shí)對(duì)設(shè)計(jì)作出判斷和修改；利用計(jì)算機(jī)可以進(jìn)行與圖形的編輯、放大、縮小、平移和旋轉(zhuǎn)等有關(guān)的圖形數(shù)據(jù)加工工作。常用的CAD軟件，也就是所謂

5、的三維制圖軟件，較二維的圖紙和二維的繪圖軟件（比如AutoCAD）而言，三維CAD軟件能夠更加直觀、準(zhǔn)確地反映實(shí)體和特征。對(duì)于專業(yè)企業(yè)，因?yàn)槔L制目標(biāo)不同，還常存在有多種CAD系統(tǒng)并行的局面，那么就需要配置統(tǒng)一的、具備跨平臺(tái)能力的零部件數(shù)據(jù)資源庫(kù)，將標(biāo)準(zhǔn)件庫(kù)和外購(gòu)件庫(kù)內(nèi)的模型數(shù)據(jù)以中間格式（比如通用的有IGS、STEP等）導(dǎo)出到三維構(gòu)型中去，如主流的Autodesk Inventor，SolidWorks,CATIA,SolidEdge,Pro/E,AutoCAD,UG NX,Onespace等。國(guó)外，這種網(wǎng)絡(luò)服務(wù)被稱為“零部件圖書(shū)館”或“數(shù)據(jù)資源倉(cāng)庫(kù)”。目前，航天航空領(lǐng)域使用較多的為Pro/E

6、，飛機(jī)和汽車等復(fù)雜產(chǎn)品制造領(lǐng)域則使用Catia居多，而在中小企業(yè)使用Solidworks較多。在歐美和日本的PLM用戶中，基于互聯(lián)網(wǎng)的PLM零部件數(shù)據(jù)資源平臺(tái)LinkAble PARTcommunity（簡(jiǎn)稱PCOM）的知名度一點(diǎn)都不亞于今天我們所熟知的BLOG和SNS這樣的網(wǎng)絡(luò)平臺(tái)。理論基礎(chǔ)：靜電場(chǎng)的三個(gè)基本定理唯一性定理若已經(jīng)給定系統(tǒng)中所有電極的形狀和排列，并給定每一電極的電位，那么由這些電極所產(chǎn)生的靜電場(chǎng)將由拉普拉斯方程唯一地確定。該定理在實(shí)際中的用處為：不論用什么方法找到一個(gè)函數(shù)，若它既能滿足拉普拉斯方程，又能在區(qū)域的邊界上符合給定的電位值，那么，它就一定是真正的解，而且也是唯一的解。

7、相似性定理若電極系統(tǒng)中各電極的電位都增大為K倍，當(dāng)電位零點(diǎn)不變時(shí)，則空間各點(diǎn)的電位也都增大為K倍，從而系統(tǒng)中等位面的形狀不變。當(dāng)系統(tǒng)各電極的電位保持不變，而電極尺寸按比例相似增大為K倍，則原系統(tǒng)中任一點(diǎn)上的電位和放大了的系統(tǒng)中對(duì)應(yīng)點(diǎn)（Kz，Kr）上的電位完全相同，只要系統(tǒng)坐標(biāo)零點(diǎn)不變即可。因之，新系統(tǒng)中的等位面可以看作是原系統(tǒng)中的等位面保持幾何相對(duì)形狀和電位數(shù)值不變，只是尺寸放大為K倍似的。相似性定理只對(duì)拉普拉斯方程成立。多電極系統(tǒng)的電位疊加定理當(dāng)各電極的形狀、相對(duì)位置確定后，電場(chǎng)分布滿足下述疊加定理：其中為第i個(gè)電極上所加的電位值，為=1、其它=0（ji）時(shí)的系統(tǒng)內(nèi)的電位分布函數(shù)，稱為相應(yīng)電

8、極的單位電位分布函數(shù)。對(duì)一個(gè)各電極形狀、相對(duì)位置確定的系統(tǒng)，當(dāng)需要不斷通過(guò)改變一個(gè)或部分電極的電位值獲得系統(tǒng)新的電位分布時(shí)，可依據(jù)該定理由已先行計(jì)算、存儲(chǔ)的單位電位分布函數(shù)的簡(jiǎn)單疊加來(lái)獲得，而不必每次都重新進(jìn)行電位分布的迭代計(jì)算。計(jì)算電場(chǎng)的拉普拉斯方程：該方程是橢圓型偏微分方程，求解時(shí)必須用到電位所滿足的不同類型的封閉邊界條件。通常都是給定邊界區(qū)域和電極上的電位。這類問(wèn)題稱為第一類邊界值問(wèn)題，或稱為荻利赫萊問(wèn)題。由微分方程知，這方程有唯一確定解。拉普拉斯方程為：u=d2u/dx2+d2u/dy2=0，其中 為拉普拉斯算子，此處的拉普拉斯方程為二階偏微分方程。三維情況下，拉普拉斯方程可由下面的形

9、式描述，問(wèn)題歸結(jié)為求解對(duì)實(shí)自變量x、y、z二階可微的實(shí)函數(shù) ：其中 稱為拉普拉斯算子.拉普拉斯方程的解稱為調(diào)和函數(shù)。如果等號(hào)右邊是一個(gè)給定的函數(shù)f(x,y,z)，即：則該方程稱為泊松方程。 拉普拉斯方程和泊松方程是最簡(jiǎn)單的橢圓型偏微分方程。偏微分算子或 （可以在任意維空間中定義這樣的算子）稱為拉普拉斯算子，英文是Laplace operator或簡(jiǎn)稱作Laplacian。狄利克雷問(wèn)題拉普拉斯方程的狄利克雷問(wèn)題可歸結(jié)為求解在區(qū)域D內(nèi)定義的函數(shù)，使得在D的邊界上等于某給定的函數(shù)。為方便敘述，以下采用拉普拉斯算子應(yīng)用的其中一個(gè)例子熱傳導(dǎo)問(wèn)題作為背景進(jìn)行介紹：固定區(qū)域邊界上的溫度（是邊界上各點(diǎn)位置坐標(biāo)

10、的函數(shù)），直到區(qū)域內(nèi)部熱傳導(dǎo)使溫度分布達(dá)到穩(wěn)定，這個(gè)溫度分布場(chǎng)就是相應(yīng)的狄利克雷問(wèn)題的解。諾伊曼邊界條件拉普拉斯方程的諾伊曼邊界條件不直接給出區(qū)域D邊界處的溫度函數(shù)本身，而是沿D的邊界法向的導(dǎo)數(shù)。從物理的角度看，這種邊界條件給出的是矢量場(chǎng)的勢(shì)分布在區(qū)域邊界處的已知效果（對(duì)熱傳導(dǎo)問(wèn)題而言，這種效果便是邊界熱流密度）。方程的解稱為調(diào)和函數(shù)，此函數(shù)在方程成立的區(qū)域內(nèi)是解析的。任意兩個(gè)函數(shù)，如果它們都滿足拉普拉斯方程（或任意線性微分方程），這兩個(gè)函數(shù)之和（或任意形式的線性組合）同樣滿足前述方程。這種非常有用的性質(zhì)稱為疊加原理。可以根據(jù)該原理將復(fù)雜問(wèn)題的已知簡(jiǎn)單特解組合起來(lái)，構(gòu)造適用面更廣的通解。接下來(lái)

11、我們?cè)O(shè)任意網(wǎng)格點(diǎn)0的坐標(biāo)為（），其電位為；相鄰的四點(diǎn)1，2，3，4與它的間距分別為（均取絕對(duì)值），各點(diǎn)的電位分別為。利用Taylor公式，并忽略高階小量。代入到Laplace方程中得到旋轉(zhuǎn)對(duì)稱場(chǎng)“十”字形不等距五點(diǎn)差分公式（如下圖）：（不等距“十”字形五點(diǎn)差分格式） （1-3）式中各系數(shù)是相鄰的1，2，3，4點(diǎn)網(wǎng)格間距（步長(zhǎng)）及r的函數(shù)： ， （1-4） ， （1-5） （1-6）對(duì)軸上點(diǎn)=0，其差分公式的形式與（1-3）式相同，且其系數(shù)也與（1-4）式相同，但不同： ， （1-7）當(dāng)然也不同： （1-8）對(duì)于實(shí)際網(wǎng)格點(diǎn)，對(duì)軸外點(diǎn)， ， ， ， ， ， ，其中i，j分別為0點(diǎn)所在處的行號(hào)和列號(hào)

12、。對(duì)于軸上點(diǎn)，i=0，不出現(xiàn)，=0?？梢钥闯觯鲜霾罘止骄褪菍^(qū)域內(nèi)任何點(diǎn)的點(diǎn)位與其周圍相鄰點(diǎn)的電位聯(lián)系起來(lái)的線性代數(shù)方程，且方程的系數(shù)在電子光學(xué)系統(tǒng)結(jié)構(gòu)尺寸確定以及網(wǎng)格一旦劃定后都是已知的。如果在電場(chǎng)區(qū)域內(nèi)，自軸線開(kāi)始，沿z軸方向逐點(diǎn)逐點(diǎn)地、按r方向逐行逐行地（即在z方向從左到右、在r方向從下到上），把區(qū)域網(wǎng)格點(diǎn)的電位依次編為（網(wǎng)格點(diǎn)總數(shù)為N），就可以在所有網(wǎng)格點(diǎn)上都建立起相應(yīng)的差分方程。于是便得到一個(gè)N 元一次線性方程組：= （1-9）其中常數(shù)矩陣不為零矩陣，是因?yàn)閷⑦吔缣幍碾娢恢党?shù)代入的緣故。由于差分方程的系數(shù)矩陣是稀疏矩陣，對(duì)于這類方程組，一般采用迭代法求出近似解。連續(xù)超張弛迭代法

13、：u=u*(1-w)+(c1*u1+c2*u2+c3*u3+c4*u4)/c0;對(duì)于荻利赫萊問(wèn)題，從理論上證明當(dāng)w在0和2之間時(shí)，迭代過(guò)程是收斂的。當(dāng)w大于2時(shí)，迭代將是發(fā)散的。在1和2之間時(shí)，將加速收斂過(guò)程，成為超張弛迭代過(guò)程。拉格朗日插值法：在數(shù)值分析中，拉格朗日插值法是以法國(guó)十八世紀(jì)數(shù)學(xué)家約瑟夫·路易斯·拉格朗日命名的一種多項(xiàng)式插值方法。許多實(shí)際問(wèn)題中都用函數(shù)來(lái)表示某種內(nèi)在聯(lián)系或規(guī)律，而不少函數(shù)都只能通過(guò)實(shí)驗(yàn)和觀測(cè)來(lái)了解。如對(duì)實(shí)踐中的某個(gè)物理量進(jìn)行觀測(cè)，在若干個(gè)不同的地方得到相應(yīng)的觀測(cè)值，拉格朗日插值法可以找到一個(gè)多項(xiàng)式，其恰好在各個(gè)觀測(cè)的點(diǎn)取到觀測(cè)到的值。這樣的多項(xiàng)

14、式稱為拉格朗日（插值）多項(xiàng)式。數(shù)學(xué)上來(lái)說(shuō)，拉格朗日插值法可以給出一個(gè)恰好穿過(guò)二維平面上若干個(gè)已知點(diǎn)的多項(xiàng)式函數(shù)。拉格朗日插值法最早被英國(guó)數(shù)學(xué)家愛(ài)德華·華林于1779年發(fā)現(xiàn)，不久后（1783年）由萊昂哈德·歐拉再次發(fā)現(xiàn)。1795年，拉格朗日在其著作師范學(xué)校數(shù)學(xué)基礎(chǔ)教程中發(fā)表了這個(gè)插值方法，從此他的名字就和這個(gè)方法聯(lián)系在一起。定義：一般地，若已知y=f(x)在互不相同 n+1 個(gè)點(diǎn)x0,x1,.,xn處的函數(shù)值y0,y1,.,yn, 則可以考慮構(gòu)造一個(gè)過(guò)這n+1 個(gè)點(diǎn)的、次數(shù)不超過(guò)n的多項(xiàng)式y(tǒng)=Pn(x),使其滿足：Pn(xk)=yk, k=0,1,2,.,n （*）要估計(jì)任一

15、點(diǎn)，xi,i=0,1,2,.,n,則可以用Pn（）的值作為準(zhǔn)確值f()的近似值，此方法叫做“插值法”。稱式（*）為插值條件（準(zhǔn)則），含xi(i=0,1,.,n)的最小區(qū)間a,b(a=minx0,x1,.,xn,b=maxx0,x1,.,xn)定理：滿足插值條件的、次數(shù)不超過(guò)n的多項(xiàng)式是存在而且是唯一的設(shè)計(jì)思想：上述的理論支持，我們可以這樣設(shè)計(jì)程序:先設(shè)計(jì)數(shù)據(jù)的輸入部分，這里主要依靠C語(yǔ)言的語(yǔ)言環(huán)境進(jìn)行完成。然后完成初始化，將給定的邊界條件和電極電位賦值到初始化空間網(wǎng)格中。之后，我根據(jù)迭代的特點(diǎn)，決定每迭代一個(gè)點(diǎn)就計(jì)算一次步長(zhǎng)。計(jì)算步長(zhǎng)的公式大致為：H1=Zx/NxH2=Zx/NxH4=r1/m

16、1H3=r1/m1在完成這些計(jì)算之后，我們進(jìn)行第一次迭代。完成后在進(jìn)行十二次全局迭代。之后應(yīng)用繪圖軟件，完成圖形輸出，并將處理后的數(shù)據(jù)輸出。程序主要變量與數(shù)組：int x,y,t,g1,g2;doublew,u,u1,u2,u3,u4,c0,c1,c2,c3,c4,r0,H1,H2,H3,H4 ,zong1,zong2,langmu,w1,w2=1.375,miu,o=0,wsx=0,rfv=0,tgb=0,ujm=0; double derta=0;double ypxl;double Z10;int n,zxc,edc=0,zxc1;double qaz,d,e;inti,ik,qwe,q

17、wer,qwe1,qwert1,qwe2,qwe3,qwe4,qwe5,qwe6,qwe7,qwe8,qwe9,qwe10,qwert2,qwert3,qwert4,qwert5,qwert6,qwert7,qwert8,qwert9,qwert10,asd=0;double r1=0;double r2=0;double m1=0;double m2=0;double N10=0;double V10=0;int NST=0;int INS=0;int DERTAV=0;程序流程框圖：使用說(shuō)明：輸入數(shù)據(jù)時(shí)，必須按照給定的數(shù)據(jù)順序，依次為，電極數(shù)（整數(shù)型），電極厚度（小數(shù)型），相鄰電極距離（小

18、數(shù)型），相鄰電極間隔步長(zhǎng)數(shù)（整數(shù)型），電極電位（整數(shù)型），電極半徑（小數(shù)型），迭代控制精度（小數(shù)型），打印間隔（整數(shù)型），軸上電位步長(zhǎng)數(shù)（整數(shù)型），差值間隔（整數(shù)型）。若電場(chǎng)分布，電極數(shù)量有不同，應(yīng)該適當(dāng)修改后，再次輸入，以免計(jì)算出現(xiàn)失誤。結(jié)論與討論：在經(jīng)過(guò)全局迭代之后，可以看出來(lái)有些值并不是十分準(zhǔn)確，說(shuō)明算法在設(shè)計(jì)的時(shí)候存在一些問(wèn)題。而且由于個(gè)人能力問(wèn)題，未能設(shè)計(jì)出有鞍點(diǎn)時(shí)的情況。我想這是我的不足。但是經(jīng)理了這個(gè)學(xué)習(xí)，設(shè)計(jì)的過(guò)程之后我想我對(duì)于程序的設(shè)計(jì)有了一定的認(rèn)識(shí)。并且在時(shí)間的過(guò)程中，我也在不斷努力不斷進(jìn)步，最終也設(shè)計(jì)出了自己的程序我想這是我的收獲。相信如果以后還有類似的程序設(shè)計(jì)項(xiàng)目，我能

19、做得更好。同時(shí)，雖然有一些算法雖然計(jì)算結(jié)果沒(méi)問(wèn)題，但是顯得過(guò)于繁瑣。這也是需要我以后注意和改進(jìn)的。致謝：在這一個(gè)學(xué)期中，我最終能完成程序，離不開(kāi)老師的講解，指導(dǎo)和幫助。同時(shí)還有學(xué)長(zhǎng)學(xué)姐在我遇到困難的時(shí)候能幫助我解決困難。參考文獻(xiàn)：周立偉，倪國(guó)強(qiáng)，方二倫。圖像無(wú)旋轉(zhuǎn)的電磁聚焦移向系統(tǒng)的研究。電子學(xué)報(bào) 1984向世明，倪國(guó)強(qiáng)編著。光電子成像器件原理。北京：國(guó)防工業(yè)出版社 1999附錄：源程序：#include<stdafx.h>#include <conio.h>#include <graphics.h> #include <stdio.h>#inc

20、lude <malloc.h>#include<math.h>void main() double derta=0;double ypxl;double Z10;int n,zxc,edc=0,zxc1;double qaz,d,e;int i,ik,qwe,qwer,qwe1,qwert1,qwe2,qwe3,qwe4,qwe5,qwe6,qwe7,qwe8,qwe9,qwe10,qwert2,qwert3,qwert4,qwert5,qwert6,qwert7,qwert8,qwert9,qwert10,asd=0;double r1=0;double r2=0;d

21、ouble m1=0;double m2=0;double N10=0;double V10=0;int NST=0;int INS=0;int DERTAV=0;FILE *fp;fp=fopen("D:/1.txt","r+"); / 文件輸入fscanf(fp,"%lf",&derta);printf("%lfn",derta);fscanf(fp,"%d",&n);printf("%dn",n);for(i=0;i<n;i+)fscanf(fp,

22、"%lf",&Zi);for(i=0;i<7;i+)printf(" %lf",Zi);for(i=0;i<n;i+)fscanf(fp,"%lf",&Ni);for(i=0;i<7;i+)printf(" %lf",Ni);for(i=0;i<n;i+)fscanf(fp,"%lf",&Vi);for(i=0;i<7;i+)printf(" %lf",Vi);fscanf(fp,"%d",&r

23、1);printf("%dn",r1);fscanf(fp,"%d",&m1);printf("%dn",m1);fscanf(fp,"%d",&r2);printf("%dn",r2);fscanf(fp,"%d",&m2);printf("%dn",m2);fscanf(fp,"%lf",&ypxl);printf("%lfn",ypxl);fscanf(fp,"%d&q

24、uot;,&NST);printf("%dn",NST);fscanf(fp,"%d",&INS);printf("%dn",INS);fscanf(fp,"%d",&DERTAV);printf("%dn",DERTAV);double *p,*qw;p=(double *)malloc(5001*sizeof(double); qw=(double *)malloc(5001*sizeof(double); / 初始化，網(wǎng)格化*p=0;int k=0;int q,j;

25、for(q=0;q<=15+23+30;q+)if(q=0)for(i=0;i<7;i+)for(j=Ni;j>=0;j-)if(i=0)*(p+k)=0+4*(Ni-j);k=k+1;elseif(j=Ni)*(p+k)=Vi-1;k=k+1;else*(p+k)=Vi-1+(Ni-j)*(Vi-Vi-1)/(Ni);k=k+1;if(q>0&&q<=15)*(p+k)=0;k=k+1;for(i=0;i<7;i+)for(j=Ni;j>0;j-)if(i=0)if(j=1)*(p+k)=Vi;k=k+1;else*(p+k)=0;k

26、=k+1;elseif(j=1)*(p+k)=0;k=k+1;*(p+k)=Vi;k=k+1; break;if(j=Ni)*(p+k)=Vi-1;k=k+1; else*(p+k)=0;k=k+1;if(q>15)*(p+k)=0;k=k+1;for(i=0;i<7;i+)for(j=Ni;j>0;j-)if (i=6&&j=1)*(p+k)=0;k=k+1;*(p+k)=0;k=k+1;*(p+k)=0;k=k+1;*(p+k)=0;k=k+1;*(p+k)=0;k=k+1;*(p+k)=0;k=k+1;*(p+k)=100;else*(p+k)=0;k=

27、k+1;int x,y,t,g1,g2;double w,u,u1,u2,u3,u4,c0,c1,c2,c3,c4,r0,H1,H2,H3,H4,housan3,zong1,zong2,langmu,w1,w2=1.375,miu,o=0,wsx=0,rfv=0,tgb=0,ujm=0; / 計(jì)算H1,H2,H3,H4t=0;j=0;for(y=38;y>=0;y-)for(x=0;x<7;x+)for(j=0;j<=Nx-1;j+)if(y<15)if(x=0)H3=r2/m2;H4=r2/m2;if(j=Nx)H1=Zx/Nx;H2=0.5;elseH1=Zx/Nx

28、;H2=Zx/Nx;if(x=6)H3=r2/m2;H4=r2/m2;if(j=0)H2=Zx/Nx;H1=0.5;elseH1=Zx/Nx;H2=Zx/Nx;else H3=r2/m2;H4=r2/m2;if(j=Nx)H1=Zx/Nx;H2=0.5;if(j=0)H2=Zx/Nx;H1=0.5;elseH1=Zx/Nx;H2=Zx/Nx;r0=(r2/m2)*(y+1);if(y=15)if(x=0)H3=r1/m1;H4=r2/m2;if(j=Nx)H1=Zx/Nx;H2=0.5;elseH1=Zx/Nx;H2=Zx/Nx;if(x=6)H3=r1/m1;H4=r2/m2;if(j=0)

29、H2=Zx/Nx;H1=0.5;elseH1=Zx/Nx;H2=Zx/Nx;else H3=r1/m1;H4=r2/m2;if(j=Nx)H1=Zx/Nx;H2=0.5;if(j=0)H2=Zx/Nx;H1=0.5;elseH1=Zx/Nx;H2=Zx/Nx;r0=(r2/m2)*(y+1);if(y>15)if(x=0)H3=r1/m1;H4=r1/m1;if(j=Nx)H1=Zx/Nx;H2=0.5;elseH1=Zx/Nx;H2=Zx/Nx;if(x=6)H3=r1/m1;H4=r1/m1;if(j=0)H2=Zx/Nx;H1=0.5;elseH1=Zx/Nx;H2=Zx/Nx;e

30、lse H3=r1/m1;H4=r1/m1;if(j=Nx)H1=Zx/Nx;H2=0.5;if(j=0)H2=Zx/Nx;H1=0.5;elseH1=Zx/Nx;H2=Zx/Nx;r0=15*(r2/m2)+(y-15)*(r1/m1); if(j=1)continue; if(j=(Nx-1) break;if(y=0)break;if(y=38)c1=2/(H1*(H1+H2); /第一次迭代c2= 2/(H2*(H1+H2);c3=0;c4=4/(H4*H4);c0=c1+c2+c3+c4;w=1;t=0;for(i=0;i<=x;i+)t=t+Ni+1; t=t-Nx;u1=*

31、(p+69*y+t+j-1);u2=*(p+69*y+t+j+1);u3=*(p+69*(y+1)+t+j);u4=*(p+69*(y-1)+t+j);*(p+69*y+t+j)=(c1*u1+c2*u2+c3*u3+c4*u4)/c0;elsec1=2/(H1*(H1+H2);c2= 2/(H2*(H1+H2);c3=(2*r0-H4)/(r0*H3*(H3+H4);c4=(2*r0+H3)/(r0*H4*(H3+H4);c0=c1+c2+c3+c4;w=1;t=0;for(i=0;i<=x;i+)t=t+Ni+1; t=t-Nx;u1=*(p+69*y+t+j);u2=*(p+69*

32、y+t+j+2);u3=*(p+69*(y+1)+t+j+1);u4=*(p+69*(y-1)+t+j+1);*(p+69*y+j+t)=(c1*u1+c2*u2+c3*u3+c4*u4)/c0;for(zxc=0;zxc<1;zxc+) /后迭代for(zxc1=0;zxc1<=11;zxc1+)ik=0;zong1=0;zong2=0;for(y=38;y>=0;y-)for(x=0;x<7;x+)for(j=0;j<=Nx-1;j+)if(y<15)if(x=0)H3=r2/m2;H4=r2/m2;if(j=Nx)H1=Zx/Nx;H2=0.5;els

33、eH1=Zx/Nx;H2=Zx/Nx;if(x=6)H3=r2/m2;H4=r2/m2;if(j=0)H2=Zx/Nx;H1=0.5;elseH1=Zx/Nx;H2=Zx/Nx;else H3=r2/m2;H4=r2/m2;if(j=Nx)H1=Zx/Nx;H2=0.5;if(j=0)H2=Zx/Nx;H1=0.5;elseH1=Zx/Nx;H2=Zx/Nx;r0=(r2/m2)*(y+1);if(y=15)if(x=0)H3=r1/m1;H4=r2/m2;if(j=Nx)H1=Zx/Nx;H2=0.5;elseH1=Zx/Nx;H2=Zx/Nx;if(x=6)H3=r1/m1;H4=r2/m

34、2;if(j=0)H2=Zx/Nx;H1=0.5;elseH1=Zx/Nx;H2=Zx/Nx;else H3=r1/m1;H4=r2/m2;if(j=Nx)H1=Zx/Nx;H2=0.5;if(j=0)H2=Zx/Nx;H1=0.5;elseH1=Zx/Nx;H2=Zx/Nx;r0=(r2/m2)*(y+1);if(y>15)if(x=0)H3=r1/m1;H4=r1/m1;if(j=Nx)H1=Zx/Nx;H2=0.5;elseH1=Zx/Nx;H2=Zx/Nx;if(x=6)H3=r1/m1;H4=r1/m1;if(j=0)H2=Zx/Nx;H1=0.5;elseH1=Zx/Nx;H

35、2=Zx/Nx;else H3=r1/m1;H4=r1/m1;if(j=Nx)H1=Zx/Nx;H2=0.5;if(j=0)H2=Zx/Nx;H1=0.5;elseH1=Zx/Nx;H2=Zx/Nx;r0=15*(r2/m2)+(y-15)*(r1/m1); if(j=1)continue; if(j=(Nx-1) break;if(y=0)break;if (y=38)c1=2/(H1*(H1+H2);c2= 2/(H2*(H1+H2);c3=0;c4=4/(H4*H4);c0=c1+c2+c3+c4; t=0; for(i=0;i<x;i+)t=t+Ni+1; t=t-Nx; u1=

36、*(p+69*y+t+j-1); u2=*(p+69*y+t+j+1); u3=*(p+69*(y+1)+t+j); u4=*(p+69*(y-1)+t+j); if(fabs(wsx/2)-w2)/(2-w2)<0.05) w=w2;d=*(p+69*y+t+j); *(p+69*y+t+j)=(*(p+69*y+t+j)*(1-w)+(c1*u1+c2*u2+c3*u3+c4*u4)/c0; e=*(p+69*y+t+j); if(fabs(*(p+69*y+t+j)-d)>0.0005) ik=ik+1;else w=w2; d=*(p+69*y+t+j);*(p+69*y+

37、t+j)=(*(p+69*y+t+j)*(1-w)+(c1*u1+c2*u2+c3*u3+c4*u4)/c0;zong1=zong1+fabs(*(p+69*y+t+j)-d)/2691; if(zxc1!=11) zong2=zong1; elsec1=2/(H1*(H1+H2);c2= 2/(H2*(H1+H2);c3=(2*r0-H4)/(r0*H3*(H3+H4);c4=(2*r0+H3)/(r0*H4*(H3+H4);c0=c1+c2+c3+c4;t=0;for(i=0;i<x;i+)t=t+Ni+1;t=t-Nx;u1=*(p+69*y+t+j-1);u2=*(p+69*y+

38、t+j+1);u3=*(p+69*(y+1)+t+j);u4=*(p+69*(y-1)+t+j);u=(c1*u1+c2*u2+c3*u3+c4*u4)/c0;if( fabs(wsx/2)-w2)/(2-w2)<0.05) w=w2;d=*(p+69*y+t+j); *(p+69*y+t+j)=(*(p+69*y+t+j)*(1-w)+(c1*u1+c2*u2+c3*u3+c4*u4)/c0; if(fabs(*(p+69*y+t+j)-d)>0.0005) ik=ik+1;else w=w2; d=*(p+69*y+t+j);*(p+69*y+t+j)=(*(p+69*y+t+

39、j)*(1-w)+(c1*u1+c2*u2+c3*u3+c4*u4)/c0;zong1=zong1+fabs(*(p+69*y+t+j)-d)/2691; if(zxc1!=11) zong2=zong1; if(zxc1=11)langmu=zong1/zong2;miu=(langmu+w2-1)/(sqrt(langmu)*w2);w1=2/(1+sqrt(1-miu*miu);wsx=w2;w2=1.25*w1-0.5;*qw=langmu; qw=qw+1;wsx=wsx+w2;if(ik=0)break;fp=fopen("D:/3.txt","w+&

40、quot;);fprintf(fp," %lf",&derta);fprintf(fp," %d",n);fprintf(fp," rn");for(i=0;i<n;i+)fprintf(fp," %lf",Zi);fprintf(fp," rn");for(i=0;i<n;i+)fprintf(fp," %lf",Ni);fprintf(fp," rn");for(i=0;i<n;i+)fprintf(fp," %lf

41、",Vi);fprintf(fp," rn");fprintf(fp," %d",r1);fprintf(fp," rn");fprintf(fp," %d",m1);fprintf(fp," rn");fprintf(fp," %d",r2);fprintf(fp," rn");fprintf(fp," %d",m2);fprintf(fp," rn");fprintf(fp," %lf&quo

42、t;,ypxl);fprintf(fp," rn");fprintf(fp," %d",NST);fprintf(fp," rn");fprintf(fp," %d",INS);fprintf(fp," rn");fprintf(fp," %d",DERTAV);fprintf(fp," rn");for(i=0;i<2700;i+)if(NST=2)fprintf(fp," %8.5lf",*(p+i);if(i%68=0&am

43、p;&i/68=1)fprintf(fp,"n");if(i-68)%69=0&&(i-68)/69>=1) fprintf(fp,"n");elseif(i/69)%2=0) fprintf(fp," %lf",*(p+i);else;for(i=0;i<2700;i+)fprintf(fp," %8.5lf",*(qw+i);for(i=0;i<2700;i+)if(10>*(p+i-1)&&(10<*(p+i)qwer=i;break;qw

44、e1=qwer;for(i=0;i<2700;i+)if(10>*(p+i-1)&&(10<*(p+i)qwer=i;if(i>2691)break;qwert1=qwer;for(i=0;i<2700;i+)if(20>*(p+i-1)&&(20<*(p+i)qwer=i;break;qwe2=qwer;for(i=0;i<2700;i+)if(20>*(p+i-1)&&(20<*(p+i)qwer=i;if(i>2691)break;qwert2=qwer;for(i=0;i<2700;i+)if(30>*(p+i-1)&&(30<*(p+i)qwer=i;break;qwe3=qwer;for(i=0;i<2700;i+)if(30>*(p+i-1)&&(30<*(p+i)qwer=i;if(i>2691)break;qwert3=qwer;for(i=0;i<2700;i+)if(40>*(p+i-1)&&(40<*(p+i)qwer=i;break;qwe4=qwer;for(i=0;