function [alpha,N,U]=youxianchafen2(r1,r2,up,under,num,deta)
%[alpha,N,U]=youxianchafen2(a,r1,r2,up,under,num,deta)
%該函數用有限差分法求解有兩種介質的正方形區域的二維拉普拉斯方程的數值解
%函數返回迭代因子、迭代次數以及迭代完成后所求區域內網格節點處的值
%a為正方形求解區域的邊長
%r1,r2分別表示兩種介質的電導率
%up,under分別為上下邊界值
%num表示將區域每邊的網格剖分個數
%deta為迭代過程中所允許的相對誤差限
n=num+1; %每邊節點數
U(n,n)=0; %節點處數值矩陣
N=0; %迭代次數初值
alpha=2/(1+sin(pi/num));%超松弛迭代因子
k=r1/r2; %兩介質電導率之比
U(1,1:n)=up; %求解區域上邊界第一類邊界條件
U(n,1:n)=under; %求解區域下邊界第一類邊界條件
U(2:num,1)=0;U(2:num,n)=0;
for i=2:num
U(i,2:num)=up-(up-under)/num*(i-1);%采用線性賦值對上下邊界之間的節點賦迭代初值
end
G=1;
while G>0 %迭代條件:不滿足相對誤差限要求的節點數目G不為零
Un=U; %完成第n次迭代后所有節點處的值
G=0; %每完成一次迭代將不滿足相對誤差限要求的節點數目歸零
for j=1:n
for i=2:num
U1=U(i,j); %第n次迭代時網格節點處的值
if j==1 %第n+1次迭代左邊界第二類邊界條件
U(i,j)=1/4*(2*U(i,j+1)+U(i-1,j)+U(i+1,j));
end
if (j>1)&&(j U2=1/4*(U(i,j+1)+ U(i-1,j)+ U(i,j-1)+ U(i+1,j));
U(i,j)=U1+alpha*(U2-U1); %引入超松弛迭代因子后的網格節點處的值
end
if i==n+1-j %第n+1次迭代兩介質分界面(與網格對角線重合)第二類邊界條件
U(i,j)=1/4*(2/(1+k)*(U(i,j+1)+U(i+1,j))+2*k/(1+k)*(U(i-1,j)+U(i,j-1)));
end
if j==n %第n+1次迭代右邊界第二類邊界條件
U(i,n)=1/4*(2*U(i,j-1)+U(i-1,j)+U(i+1,j));
end
end
end
N=N+1 %顯示迭代次數
Un1=U; %完成第n+1次迭代后所有節點處的值
err=abs((Un1-Un)./Un1);%第n+1次迭代與第n次迭代所有節點值的相對誤差
err(1,1:n)=0; %上邊界節點相對誤差置零
err(n,1:n)=0; %下邊界節點相對誤差置零
G=sum(sum(err>deta))%顯示每次迭代后不滿足相對誤差限要求的節點數目G
end
標簽:
有限差分
上傳時間:
2018-07-13
上傳用戶:Kemin
VIP專區-嵌入式/單片機編程源碼精選合集系列(94)資源包含以下內容:1. 本文介紹了嵌入式系統的概念, 分析了μC/OS 的內核結構, 并詳細介紹了在具有ARM 體系結構的S3C44B0 微處理器
上進行μC/OS 操作系統的移植和應用程序及驅動程序的開發。.2. 本文介紹了嵌入式系統的概念, 分析了μC/OS 的內核結構, 并詳細介紹了在具有ARM 體系結構的S3C44B0 微處理器
上進行μC/OS 操作系統的移植和應用程序及驅動程序的開發。.3. 介紹嵌入式開發適合初學者的學習及有一定工作經驗的人作為參考手冊.4. 通過C++和GLUT.5. ht芯片通過IO口讀寫I2C芯片的匯編源代碼.6. 該文檔為dsp c6000系列的原理圖資料,對于嵌入式開發者很有使用價值..7. 12864液晶的驅動程序.8. 一個完整的ASM程序.9. 單相電子式液晶電能表源程序已經是成熟產品的程序。.10. 一個用C語言控制的讀寫7022.11. 51開發板的源程序.12. sofia-sip-1.12.4.13. CPLD EPM7256原理圖PCB圖.14. 嵌入式LINUX 的驅動程序。采用2410的開發板全部可以通用(如使用引腳不同只要重新改腳定義).15. 著名EDA工具軟件VCS得技術資料。pdf格式。.16. 非常好完的游戲.17. keil c中io的編程.18. keil c中BUZZ的運用.19. Lcd的編成.20. 射頻卡讀卡電路和程序,以及網絡芯片8019的電路和程序,功能是實現一個射頻卡讀卡,讀出數據傳輸到上位機.通過網絡..21. 這是一個串口通信程序.22. 該源碼實現了為現場人員創建擁有GUI的嵌入式數據庫,現場人員通過獲得數據可以進行薪水調查,數據包括職位,企業類型,年收入,閱歷.23. * 一、功能: Timestamp驅動演示代碼.
* 二、該源碼需要硬件開發板的支持,因為ISS對Timestamp定時器的模擬還不夠精確
* 如果將該源碼運行于ISS模式下,將得不到精確的結.24. rtl8019驅動程序及其main函數的源碼.25. C++嵌入系統實例不是很全,總共7個分別是2,3,5,6,7,8,9.26. C++嵌入系統實例不是很全,總共7個分別是2,3,5,6,7,8,9.27. C++嵌入系統實例不是很全,總共7個分別是2,3,5,6,7,8,9.28. C++嵌入系統實例不是很全,總共7個分別是2,3,5,6,7,8,9.29. WINCE MFC COM sample. (from EVC高級編程及其應用開發).30. 三星ARM9的LCD驅動板原理圖和PCB圖.31. 我設計的CAN總線模塊.32. 三星ARM9 S3C2410 核心板原理圖與PCB圖.33. 我用三星ARM9 S3C2410 做的掌上電腦 原理圖與pcb圖.34. 一個大公司的ARM9開發板原理圖.35. cypress fx2 firmware代碼示例.36. cypress ezusb driver 代碼模板.37. 不錯的X86匯編代碼示例.38. 匯編代碼示例.39. 非常好的匯編代碼示例.40. 本人以前做的些東西.
標簽:
金屬材料
性能
對照
上傳時間:
2013-06-09
上傳用戶:eeworm