版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)
文檔簡介
1、精選優(yōu)質(zhì)文檔-傾情為你奉上橢圓形方程的有限元 上機實習(xí)報告專心-專注-專業(yè)兩點邊值問題有限元法(必做)從Galerkin原理出發(fā)用線性元解兩點邊值問題:精確解:。1.1變分形式從Galerkin原理出發(fā)推導(dǎo)出兩點邊值問題的變分形式,將積分區(qū)間等分為N份,則步長,記為。寫出有限元方程及系數(shù)矩陣元素。解: 由題可知p(x)=1,q(x)=1,f(x)= 所以有限元方程為,j=1,2,.,n 其中, , , , 計算有: 1.2利用MATLAB求解問題的過程依次取用MATLAB求解并圖形比較數(shù)值解與精確解,用表格列出不同剖分時的誤差。程序如下:function u= bianzhi(p,q,N)h=
2、1/N;x=0:h:1;A=zeros(N-1);for i=2:N-1 a3=(t)-p./h+h.*q.*t.*(1-t); a2=(t)p./h+h.*q.*(t.2)+p./h+h.*q.*(1-t).2); a1=(t)-p./h+h.*q.*t.*(1-t); A(i,i-1)=quad(a1,0,1); A(i,i)=quad(a2,0,1); A(i-1,i)=quad(a3,0,1);endA(1,1)=quad(a2,0,1); f=zeros(N-1,1);for i=2:N f1=(t)(x(i-1)+h.*t).2.*t+(x(i)+h.*t).2.*(1-t); f
3、(i-1)=h.*quad(f1,0,1);endu=inv(A)*f;precise_value=(exp(2)-1)(-1).*(2-3*exp(1)*exp(x)-(2*exp(1)-3)*exp(1-x)+x.2+2;plot(x,0;u;0,b-,x,precise_value,r-);legend(數(shù)值解,精確解);err=norm(0;u;0-precise_value)end N=4 N=8 N=16 N=32 N=64 N=128 N=256不同N對應(yīng)的誤差表格如下:N48163264128256err2.46e-048.56e-053.01e-051.06e-053.76e-061.33e-064.70e-07 1.3 方法總結(jié)及分析由上圖和表格可以看出,從Galerkin原理出發(fā)推導(dǎo)的兩點邊值問題的解和真實解的誤差
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 如何在小學(xué)生中樹立正確的學(xué)習(xí)價值觀
- 青少年對職業(yè)的選擇與發(fā)展
- 小學(xué)數(shù)學(xué)模擬練習(xí)試卷模擬考試
- 體育訓(xùn)練:練習(xí)簡單的球類運動
- 如何正確面對失敗
- 高中學(xué)業(yè)水平考核
- 小學(xué)中國生命之源測試
- 古老的建筑材料和青磚瓦石
- 九年級《道德與法治》(上冊)教學(xué)計劃及教學(xué)進度
- 百師聯(lián)盟2025屆高三10月一輪復(fù)習(xí)聯(lián)考語文試卷(含答案)
- GB/T 617-2006化學(xué)試劑熔點范圍測定通用方法
- GB/T 3215-2019石油、石化和天然氣工業(yè)用離心泵
- GB/T 25113-2010移動消防指揮中心通用技術(shù)要求
- GA 450-2013臺式居民身份證閱讀器通用技術(shù)要求
- 世界十大航空公司課件
- 松江新城-總體城市設(shè)計2021-2035年
- 最新國家開放大學(xué)電大《課程與教學(xué)論》形考任務(wù)4試題及答案
- lu《雨巷》 (共45張)課件
- 中石油公司資質(zhì)管理相關(guān)規(guī)定及工作要求課件
- 體檢服務(wù)禮貌用語課件
- 幼兒園大班社會:《四大發(fā)明》 課件
評論
0/150
提交評論