第一類邊界問(wèn)題的有限差分法探討
《第一類邊界問(wèn)題的有限差分法探討》由會(huì)員分享,可在線閱讀,更多相關(guān)《第一類邊界問(wèn)題的有限差分法探討(25頁(yè)珍藏版)》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。
1、.第一類邊界問(wèn)題的有限差分法探討摘要:本次重點(diǎn)是對(duì)于第一類邊界問(wèn)題的兩種不同方法的對(duì)比研討,通過(guò)計(jì)算機(jī)仿真有限差分法和計(jì)算分離變量法對(duì)同一問(wèn)題的求解,對(duì)結(jié)果進(jìn)行對(duì)比,能夠發(fā)現(xiàn)有限差分法更加快捷簡(jiǎn)便,只要迭代次數(shù)足夠多就能使誤差趨于零。而分離變量法則是準(zhǔn)確的計(jì)算出結(jié)果,只是運(yùn)算相對(duì)復(fù)雜。關(guān)鍵字:有限差分法,分離變量法,加速收斂因子,迭代次數(shù),邊界條件。引言:在給定的三類邊界條件下求解標(biāo)量位或矢量位的泊松方程或拉普拉斯方程的解一般的理論依據(jù)是唯一性定理和得加原理,由此而得出的解題方法有很多。主要分為兩大類:一是解析法(如分離變量法,鏡像法等),二是數(shù)值法(如有限差分法,有限元法等)。這兩種方法各有
2、優(yōu)點(diǎn)和不足,相比較而言在許多實(shí)際問(wèn)題中由于邊界條件過(guò)于復(fù)雜而無(wú)法求得解析解。這就需要借助于數(shù)值法來(lái)求電磁場(chǎng)的數(shù)值解。有限差分法便是一種比較容易的數(shù)值解法。本次研討就以第一類邊界問(wèn)題進(jìn)行為例來(lái)分析研究有限差分法。一、 有限差分法的定義:微分方程和積分微分方程數(shù)值解的方法為有限差分法?;舅枷胧前堰B續(xù)的定解區(qū)域用有限個(gè)離散點(diǎn)構(gòu)成的網(wǎng)格來(lái)代替, 這些離散點(diǎn)稱作網(wǎng)格的節(jié)點(diǎn);把連續(xù)定解區(qū)域上的連續(xù)變量的函數(shù)用在網(wǎng)格上定義的離散變量函數(shù)來(lái)近似;把原方程和定解條件中的微商用差商來(lái)近似, 積分用積分和來(lái)近似,于是原微分方程和定解條件就近似地代之以代數(shù)方程組,即有限差分方程組 , 解此方程組就可以得到原問(wèn)題在離
3、散點(diǎn)上的近似解。然后再利用插值方法便可以從離散解得到定解問(wèn)題在整個(gè)區(qū)域上的近似解。精品.二、 有限差分法解題的基本步驟:(1)、區(qū)域離散化,即把所給偏微分方程的求解區(qū)域細(xì)分成由有限個(gè)格點(diǎn)組成的網(wǎng)格; (2)、近似替代,即采用有限差分公式替代每一個(gè)格點(diǎn)的導(dǎo)數(shù); (3)、逼近求解。換而言之,這一過(guò)程可以看作是用一個(gè)插值多項(xiàng)式及其微分來(lái)代替偏微分方程的解的過(guò)程。三、 有限差分法公式的推導(dǎo):把求解的區(qū)域劃分成網(wǎng)格,把求解區(qū)域內(nèi)連續(xù)的場(chǎng)分布用網(wǎng)絡(luò)節(jié)點(diǎn)上的離散的數(shù)值解來(lái)代替。網(wǎng)格劃分的充分細(xì),才能夠達(dá)到足夠的精度。應(yīng)用有限差分法計(jì)算靜態(tài)場(chǎng)邊值問(wèn)題時(shí),需要把微分方程用差分方程替代。用圖形法解釋如下:精品.有
4、限差分網(wǎng)格劃分在由邊界L界定的二維區(qū)域D內(nèi),電位函數(shù)滿足拉普拉斯方程且給定第一邊界條件,則:如圖將區(qū)域D劃分為正方形網(wǎng)格,網(wǎng)格線的交點(diǎn)稱為節(jié)點(diǎn),兩相鄰平行網(wǎng)格線間的距離稱為步距 h。然后,拉普拉斯方程離散化,對(duì)于任一點(diǎn)0,有一階偏導(dǎo)數(shù):而后,對(duì)于二階偏導(dǎo)數(shù):精品.對(duì)于Y軸同理:因此拉普拉斯方程的差分格式為:緊鄰邊界節(jié)點(diǎn)的拉普拉斯方程的差分格式為:其中p、q為小于1的正數(shù);1、2為邊界上的節(jié)點(diǎn),其值為對(duì)應(yīng)邊界點(diǎn)處的值,是已知的。具體如圖:緊鄰邊界節(jié)點(diǎn)的網(wǎng)格劃分精品.應(yīng)用數(shù)值計(jì)算解釋(泰勒公式展開(kāi)法):1點(diǎn)電位的泰勒公式展開(kāi)為3點(diǎn)電位的泰勒公式展開(kāi)為,當(dāng)h很小時(shí),忽略4階以上的高次項(xiàng),得 同樣可得
5、將上面兩式相加得在上式中代入,得精品.對(duì)于,即F=0的區(qū)域,得到二維拉普拉斯方程的有限差分形式通過(guò)以上兩種方法的推導(dǎo),可得任意點(diǎn)的電位等于圍繞它的四個(gè)點(diǎn)的電位的平均值。當(dāng)用網(wǎng)格將區(qū)域劃分后,對(duì)每一個(gè)網(wǎng)絡(luò)點(diǎn)寫出類似的式子,就得到方程數(shù)與未知電位的網(wǎng)絡(luò)點(diǎn)數(shù)相等的線性方程組。已知的邊界條件在離散化后成為邊界上節(jié)點(diǎn)的已知電位值。四、 差分方程組的解法方法一:高斯賽德?tīng)柕ǎê?jiǎn)單迭代法)其步驟是先對(duì)每一網(wǎng)格點(diǎn)設(shè)一初值。然后按一個(gè)固定順序(一般點(diǎn)的順序按“自然順序”,即:從左到右,從下到上)如圖:之后,利用二維拉普拉斯方程的有限差分形式用圍繞它的四個(gè)點(diǎn)的電位的平均值作為它的新值,當(dāng)所有的點(diǎn)計(jì)算完后,用它
6、們的新值代替舊值,即完成了一次迭代。然后再進(jìn)行下一次迭代,如此循環(huán)。如下式:精品.(迭代公式1)其式中的上角標(biāo)(k)表示k次近似值,下腳標(biāo)i,j表示節(jié)點(diǎn)所在位置,即第i行第j列的交點(diǎn)。其中要特別注意:在迭代過(guò)程中遇到邊界點(diǎn)式,需將邊界條件 帶入。循環(huán)迭代時(shí)當(dāng)所有內(nèi)節(jié)點(diǎn)滿足以下條件時(shí)停止迭代:其中,W是預(yù)定的最大允許誤差。方法二:逐次超松弛法 簡(jiǎn)單迭代法在解決問(wèn)題時(shí)收斂速度比較慢,實(shí)用價(jià)值不大。實(shí)際中常采用逐次超松弛法(又稱高斯賽德?tīng)柕ㄗ冃危?,相比之下它有兩點(diǎn)重大的改進(jìn) ,第一是計(jì)算每一網(wǎng)格點(diǎn)時(shí),把剛才計(jì)算得到的臨近點(diǎn)的新值代入,即在計(jì)算(i,j)點(diǎn)的電位時(shí),把它左邊的點(diǎn)(i-1,j)和下面
7、的點(diǎn)(i,j-1)的電位用剛才算過(guò)的新值代入,即:(迭代公式2)第二,是引入“加速收斂因子”。上式中的即為“加速收斂因子”,且1 2。特別關(guān)注的是逐次超松弛法收斂的快慢與有明顯關(guān)系。并且最佳的取值隨著條件的不同而不同,如何選擇最佳,是個(gè)復(fù)雜問(wèn)題。精品.在計(jì)算時(shí)可以嘗試求取最佳值,以使計(jì)算快速。五、 應(yīng)用計(jì)算機(jī)仿真有限差分法解決具體問(wèn)題本次討論我選擇第四題為具體實(shí)例進(jìn)行研究。題目:如圖所示,有一長(zhǎng)方形的導(dǎo)體槽,a = 20,b = 15,設(shè)槽的長(zhǎng)度為無(wú)限長(zhǎng),槽上有一塊與槽絕緣的蓋板,電位為100V,其他板電位為零,求槽內(nèi)的電位分布。 通過(guò)MATLAB進(jìn)行仿真,運(yùn)用有限差分法,源代碼如下:u=ze
8、ros(15,20);i=2:14;u(i,20)=100;for j=1:20 for i=2:14 u(i,j)=100/19*(j-1); endenda=input(please input a(1a2);a=);m=input(please input m(1m);m=);精品.for n=1:m; for i=2:14; for j=2:19; b=u(i,j); c=u(i,j+1); d=u(i+1,j); e=u(i-1,j); f=u(i,j-1); g=0.25.*(c+d+e+f); u(i,j)=b+a.*(g-b); end endendmesh(u);通過(guò)對(duì)相關(guān)資
9、料的查詢及學(xué)習(xí),我認(rèn)為運(yùn)用MATLAB編程就是要將有限差分法的基本公式實(shí)現(xiàn)。以此為出發(fā)點(diǎn)可以讓解析和編寫過(guò)程相對(duì)簡(jiǎn)單化。首先,先要對(duì)電場(chǎng)的邊界條件進(jìn)行設(shè)定:u=zeros(15,20);i=2:14;u(i,20)=100;使其滿足題目的要求右側(cè)邊界電位為100,其余三邊為0,由此也可以判斷出此題目為第一類邊界條件的問(wèn)題。精品.同時(shí),此處,也要特別注意,因?yàn)樵贛ATLAB程序中,編寫一個(gè)矩陣的默認(rèn)順序是從左到右,從上到下。另外,命名矩陣時(shí)是先編寫行,后編寫列,因此,對(duì)于u(i,j)中,i代表的是第i行即對(duì)應(yīng)縱坐標(biāo),j代表的是第j列即對(duì)應(yīng)橫坐標(biāo),與坐標(biāo)表示有一定的差別,需要區(qū)分和注意。之后是相對(duì)
10、關(guān)鍵的編寫部分,就是對(duì)電場(chǎng)內(nèi)部的電位設(shè)定,此處,我通過(guò)分離變量法將u分解為u(x)和u(y),分別求解。為讓計(jì)算簡(jiǎn)便,我設(shè)定其內(nèi)電場(chǎng)為線性變化的,即u(x)=kx+b,u(y)=0作為初值進(jìn)行計(jì)算(應(yīng)用分離變量法得到的結(jié)果是u(x)=Bsinh m/15u(y)=Asinm/15。可以通過(guò)此處在得出結(jié)果后進(jìn)行對(duì)比和誤差分析。)又因?yàn)橐贛ATLAB中進(jìn)行運(yùn)算,所以首先要轉(zhuǎn)化為i和j的形式,因此可得:u(j)=kj+b和u(i)=0。從而得到方程組:u(1)=k+b=0;u(20)=20k+b=100。解得:u(j)=100/19*(j-1),即:u(i,j)=100/19*(j-1)。由此得到
11、了內(nèi)電場(chǎng)除去上下邊界(衡為零)的電位方程,在MATLAB中我運(yùn)用了for循環(huán)語(yǔ)句來(lái)實(shí)現(xiàn)其內(nèi)電場(chǎng)的賦值,源代碼如下:for j=1:20 for i=2:14 u(i,j)=100/19*(j-1); end精品.end將電場(chǎng)各點(diǎn)的電位賦初值完成后,就需要帶入公式進(jìn)行計(jì)算了。在之前先要對(duì)“加速收斂因子” 的最佳值和預(yù)定的最大允許誤差進(jìn)行確定。關(guān)于加速因子的最佳值確定問(wèn)題,我通過(guò)學(xué)習(xí)和查閱相關(guān)資料得到了兩種方法,如下:一是逐步搜索法;將的取值區(qū)間(1,2)進(jìn)行M等分,分別取1+1/M,1+2/M1+(M-1)/M,通過(guò)上面提到的迭代公式2依次對(duì)同一個(gè)精度要求求出迭代次數(shù)k的值,并從中選出“加速收斂
12、因子” 的最佳值,具體步驟如下:步驟1 給定等分?jǐn)?shù)M和精度要求的值,令的初始值為1;步驟2 令p=1,2,3,M-1,重復(fù)步驟3-5;步驟3 p=1+p/M;步驟4 按照如下公式迭代找出符合精度要求的迭代次數(shù)kp;步驟5 比較找出kp值最小的p為最佳的加速收斂因子的值。二是黃金分割法:依據(jù)黃金分割法的思想,通過(guò)計(jì)算機(jī)自動(dòng)選取最佳加速收斂因子的近似值,具體步驟如下:步驟1 對(duì)(1,2)區(qū)間第一次0.618的分割,區(qū)間分界a精品.1=1,b1=2,在(a1,b1)區(qū)間分割出黃金點(diǎn)point1= a1+0.618(b1- a1),進(jìn)行迭代公式2的迭代,求出迭代次數(shù)K值,如果迭代次數(shù)沒(méi)有超出規(guī)定的發(fā)散
13、常數(shù),迭代結(jié)束,否則轉(zhuǎn)為步驟2。步驟2 在(1,1.618)和(1.618,2)之間進(jìn)行第二次的黃金分割,找出分割點(diǎn)point2=a2+0.618(b2- a2),其中a2和b2是新分割區(qū)間的左右邊界。找出迭代次數(shù)最少的。以上兩種方法是比較具有實(shí)踐和研究?jī)r(jià)值進(jìn)行的,在這里介紹和分享給大家探討。在我分析過(guò)后認(rèn)為第二種較好,因?yàn)槠淇梢酝ㄟ^(guò)計(jì)算機(jī)自動(dòng)選取到最佳的加速收斂因子,但是實(shí)施起來(lái)比較復(fù)雜,要通過(guò)MATLAB或其它語(yǔ)言進(jìn)行編碼實(shí)驗(yàn)得到,我初步試驗(yàn)過(guò),沒(méi)能成功,就運(yùn)用了第一種方法進(jìn)行了計(jì)算,但是,我發(fā)現(xiàn)其運(yùn)算過(guò)于復(fù)雜,筆算正確的成功率比較低,而運(yùn)用電腦編程要重復(fù)多次運(yùn)算,雖然,計(jì)算機(jī)運(yùn)算快速且方
14、便,但是,這樣的源代碼同樣過(guò)于復(fù)雜,可以進(jìn)一步進(jìn)行深入的研究。本次的題目解答中,由于是應(yīng)用計(jì)算機(jī)計(jì)算,所以計(jì)算速度和計(jì)算效率在不同的取值中沒(méi)有明顯區(qū)別。因此,我就只是設(shè)定了一個(gè)12之間的一個(gè)數(shù)字:1.8為“加速收斂因子” 的值了,進(jìn)行了近似計(jì)算,實(shí)驗(yàn)誤差在可接受的范圍。然后,是對(duì)于迭代次數(shù)的確定。當(dāng)?shù)鷿M足停止條件:時(shí),對(duì)應(yīng)的迭代次數(shù)即為最終的迭代次數(shù),這里的W為設(shè)定好的允許誤差,為從簡(jiǎn)方式,我選擇運(yùn)用輸入迭代次數(shù),并通過(guò)多次輸入對(duì)比輸出結(jié)果選取最佳迭代次數(shù)的方法來(lái)選取,存在了一定的誤差。但是,相比于我運(yùn)用上面的公式求解計(jì)算迭代次數(shù),進(jìn)而得出的圖形的誤差比要小許多,我一直嘗試進(jìn)行改進(jìn),但一直沒(méi)
15、能很好的解決誤差過(guò)大的問(wèn)題:最主要的問(wèn)題是邊界值一直不滿足我的設(shè)定值。希望通過(guò)進(jìn)一步的研討和大家的交流能過(guò)得到最優(yōu)方式。精品.通過(guò)以上的分析,我得出了進(jìn)一步的編程代碼:a=input(please input a(1a2);a=);m=input(please input m(1m);m=);其中a和m分別為最佳“加速收斂因子”和迭代次數(shù),均為我們輸入數(shù)值計(jì)算的方式。先要預(yù)算出一些取值范圍,然后通過(guò)分別輸入,對(duì)比輸出結(jié)果進(jìn)而來(lái)選取最優(yōu)解。之后的源代碼的編寫思想就是通過(guò)迭代計(jì)算公式來(lái)實(shí)現(xiàn)對(duì)電場(chǎng)內(nèi)部各點(diǎn)電位值的計(jì)算,因此,在輸入“加速收斂因子”和“迭代次數(shù)”的數(shù)值后,我通過(guò)for的循環(huán)語(yǔ)句實(shí)現(xiàn)了對(duì)
16、電場(chǎng)內(nèi)部的電位迭代求解,并用圖形形式表示出來(lái)。其源代碼如下:for n=1:m; for i=2:14; for j=2:19; b=u(i,j); c=u(i,j+1); d=u(i+1,j);精品. e=u(i-1,j); f=u(i,j-1); g=0.25.*(c+d+e+f); u(i,j)=b+a.*(g-b); end endendmesh(u);上述代碼中,b,c,d,e,f分別場(chǎng)內(nèi)一點(diǎn)及其左右上下各點(diǎn)的電位值,之后的g和u(i,j)則實(shí)現(xiàn)了逐次超松弛法的迭代,最后用mesh實(shí)現(xiàn)圖形的輸出。六、 結(jié)果和誤差的對(duì)比分析由分離變量法得到的初步結(jié)果是u(x)=Bsinh m/15和u
17、(y)=Asinm/15,從而應(yīng)有u=Csinh(m/15)sinm(/15)的結(jié)果,從其表達(dá)式的形式,我們可以簡(jiǎn)單推斷出電場(chǎng)的電位在x軸上的分布是成虛指數(shù)函數(shù)的形式呈現(xiàn),而在y軸上是成正弦函數(shù)呈現(xiàn)的,這和通過(guò)有限差分法的計(jì)算機(jī)仿真的圖形基本符合,如下圖:精品.因此,可以基本確定通過(guò)計(jì)算機(jī)仿真得到的結(jié)果的正確性。另外,可能存在較大誤差的方面是跌代次數(shù),因?yàn)榈螖?shù)是我們?nèi)斯ぽ斎氲?,所以判斷其正確與否要通過(guò)輸出的結(jié)果來(lái)判定、來(lái)選擇,下面幾幅圖是當(dāng)“加速收斂因子” =1.8時(shí),不同的迭代次數(shù)m的對(duì)應(yīng)圖形:m=1時(shí):m=10時(shí):精品.m=20m=30m=32精品.m=35m=40通過(guò)上面迭代次數(shù)m取
18、不同的值對(duì)應(yīng)的不同圖形分析可以看出,在m=1,10,20時(shí)期計(jì)算結(jié)果都存在較大誤差,圖形與分離變量法得到的結(jié)果相差很大,有明顯的不規(guī)則部分,不成現(xiàn)光滑的曲線。而當(dāng)精品.m=30后圖形與分離變量法計(jì)算到的圖形近似,而且圖形規(guī)則,曲線光滑,而且隨著m值的增大圖形的變化非常微小了,因此,我選擇m=32為迭代次數(shù),所得到的結(jié)果誤差比較小,結(jié)果相對(duì)準(zhǔn)確。七、 研究結(jié)論學(xué)習(xí)了運(yùn)用有限差分法計(jì)算第一類邊界條件的電位值,通過(guò)計(jì)算機(jī)仿真進(jìn)行了實(shí)際例題的求解,通過(guò)與分離變量法求解的對(duì)比,證明可行性和有效性,適用于解決第一類邊界問(wèn)題,而且更加方便和快捷。同時(shí),“加速收斂因子”最佳解問(wèn)題比較復(fù)雜,不過(guò)通過(guò)計(jì)算機(jī)仿真的
19、過(guò)程中,不同數(shù)值“加速收斂因子”的加速效果不易體現(xiàn)出來(lái)。跌代次數(shù)對(duì)計(jì)算結(jié)果又較大影響,迭代次數(shù)越多,結(jié)果越準(zhǔn)確,誤差越小,但相應(yīng)的計(jì)算也會(huì)更加復(fù)雜,因此只要迭代次數(shù)得到得結(jié)果滿足一定的允許誤差范圍即可。附錄: 三類邊界條件:第一類:已知電位在邊界上的數(shù)值。第二類:已知電位的導(dǎo)數(shù)在邊界上的數(shù)值。第三類:已知電位在一部分邊界上的數(shù)值和在其余邊界上的導(dǎo)數(shù)的數(shù)值。 解析法包括:分離變量法,鏡像法,復(fù)變函數(shù)法,保角變換法,格林函數(shù)法。 數(shù)值法包括:數(shù)值積分法,有限差分法,有限元法,矩量法。 解析法的優(yōu)點(diǎn):由有限個(gè)項(xiàng)構(gòu)成閉合解或無(wú)窮級(jí)數(shù),通常具有鮮明的物理意義。精品.缺點(diǎn):解題范圍窄小,對(duì)邊界形狀十分挑剔。數(shù)值法的優(yōu)點(diǎn):解題范圍寬廣,對(duì)邊界的形狀沒(méi)有限制。缺點(diǎn):數(shù)據(jù)離散,物理意義深藏其中。心得體會(huì):有限差分法的公式套用很簡(jiǎn)單,記住結(jié)果很簡(jiǎn)單,但是要想深刻理解其內(nèi)容和靈活應(yīng)用則比較困難,需要一定時(shí)間的深刻學(xué)習(xí)和做題應(yīng)用,同時(shí),這個(gè)方法確實(shí)是一個(gè)非常方便和有效地解決三類邊界條件問(wèn)題的工具。MATLAB在電磁場(chǎng)方面的應(yīng)用相當(dāng)廣泛,其內(nèi)擁有的很多函數(shù)都可以模擬靜電場(chǎng)中的實(shí)例,就比如這次的研討就應(yīng)用了矩陣,通過(guò)有限差分法來(lái)模擬靜電場(chǎng)的電位分布,相信將來(lái)還會(huì)有更多的有意思,有研究?jī)r(jià)值計(jì)算和模擬電磁場(chǎng)的相關(guān)知識(shí)。是學(xué)好電磁場(chǎng)必須要牢固掌握的工具。如有侵權(quán)請(qǐng)聯(lián)系告知?jiǎng)h除,感謝你們的配合!精品
- 溫馨提示:
1: 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
2: 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
3.本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
5. 裝配圖網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 對(duì)標(biāo)管理概述和實(shí)施
- 市場(chǎng)營(yíng)銷信息系統(tǒng)與營(yíng)銷調(diào)研
- 市場(chǎng)經(jīng)濟(jì)的一般特征(優(yōu)質(zhì)課)
- 粒子物理學(xué)的發(fā)展PPT
- 四年級(jí)下冊(cè)語(yǔ)文園地三口語(yǔ)交際(檀麗萍)
- 四年級(jí)下冊(cè)品德與社會(huì)《創(chuàng)造文明的環(huán)境》
- 創(chuàng)業(yè)設(shè)立公司流程介紹
- XX綜合經(jīng)濟(jì)開(kāi)發(fā)區(qū)戰(zhàn)略咨詢
- 零基礎(chǔ)對(duì)外漢語(yǔ)第一課
- 寫作+寫出人物精神+(共23張PPT)(教育精品)
- 螺桿擠出機(jī)
- 考點(diǎn)4-價(jià)格變動(dòng)對(duì)經(jīng)濟(jì)生活的影響
- 美國(guó)專利查詢介紹
- GP綠色產(chǎn)品環(huán)境管理物質(zhì)知識(shí)培訓(xùn)
- 簡(jiǎn)潔明快的優(yōu)秀模板之商業(yè)管理技巧suervisory