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
VIP專區-嵌入式/單片機編程源碼精選合集系列（132）資源包含以下內容：1. siemens PLC Prodave 通訊測試.2. 嵌入式系統下IC卡讀寫的源碼.3. 數字音樂盒的完整畢業論文（無意中發現的）
課程設計題目：數字音樂盒.4. 單片機程序 51頻率計 電子琴 電子鐘程序
按鍵顯示程序最終成功.5. C和C++嵌入式系統編程,一本很嵌入式編程很不錯的書..6. 一本很好的介紹yaffs文件系統的書籍 一本很好的介紹yaffs文件系統的書籍.7. 一份不錯的介紹linux 字符設備驅動的文檔.8. 工程源文件.9. 用holtek3000開發環境編寫的單片機C語言鍵盤掃描程序．.10. 8255A芯片的中文資料（8255中文資料）.11. mmc_sd_memory_card_driver_sample_project_files spi接口的sd card driver 示例.12. 單片機驅動PS2鍵盤的c程序
希望對大家有用.13. a voice guide client ,it is the second part of voice guide center.14. msc map program to control client of tools,use for map guide.15. C#開發win ce代碼程序.16. 2個Ds18b20溫度顯示電路,12Mhz晶振.17. 5分頻cpld原理圖實現.18. 可實現對ST公司的SPI flash的控制.19. VB實現的文字圖片原理,VB實現的文字圖片原理.20. 嵌入式C精華, 嵌入式C精華,嵌入式C精華.21. 一個簡單的電梯控制系統,主要實現一個八層樓的電梯控制..22. 這是自己做的一個gui
及供初學者參考
如有問題可以一起討論
可發到我的郵箱里
完成的功能就是
對uicontrol的應用.23. 觸摸屏校正算法,主流嵌入式設計都使用此算法,個人覺得經典..24. 基于51單片機的源碼,內有液晶顯示(IIC)接口,還實現串口和定時操作等功能..25. 鏈接器和加載器,這是一本全面介紹我們編寫的程序是如何被鏈接的,個人覺得做嵌入式開發的軟件工程師必備的一本書..26. 使用Embest IDE for ARM 開發環境及ARM 軟件模擬器；
通過實驗掌握簡單ARM 匯編指令的使用方法。.27. AD7714是一款24位高速AD轉換芯片.28. Cortex-M3是首款基于ARMv7-M架構的處理器.29. 剛剛學習ne2000驅動開發.30. ATMEL的芯片資料pdf文檔.31. GPS衛星定位模塊的驅動程序.32. 還可以.33. 基于CAN總線的分布式三維步進電機控制系統.34. 自己做的交通燈~~~用c編寫的.35. 基于PCI總線的高速數據傳輸.36. [轉載]嵌入到系統進程中,檢查并報告主板型號.37. Freescale_071112_1公司內部的器件資料.38. 關于2410開發板的u-boot的說明.39. CX1-22R 可編程控制器(簡稱 CX1)是 22 路 I/O 一體化小型可編程控制器（PLC）.40. PIC匯編各個時間延時小程序.
標簽：
SD卡
上傳時間：
2013-04-15
上傳用戶：eeworm
