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
function y=lagr(x0,y0,x) %x0,y0為節點 %x是插值點 n=length(x0); m=length(x); for i=1:m z=x(i); s=0.0; for k=1:n p=1.0; for j=1:n if j~=k p=p*(z-x0(j))/(x0(k)-x0(j)); end end s=p*y0(k)+s; end y(i)=s; end
標簽： lagr
上傳時間： 2020-06-09
上傳用戶：shiyc2020
VIP專區-嵌入式/單片機編程源碼精選合集系列（58）資源包含以下內容：1. veriloghdl語言工具書.2. 嵌入式詞典查詢.3. 基于ZigBee的無線網絡技術.4. IO 口模擬232通訊程序使用兩種方式的C程序.5. 用Turbo C 實現動畫的一個簡單的實現方法 這種簡單方法利用cleardevice()和delay()函數相互配合.6. dsp2812系列開發辦的原理圖.7. 此程序為網絡通信的16位校驗累加和的經典程序。雖然簡單但是很經典。代碼簡潔.8. 觸摸式人行燈控制器原理圖 觸摸式人行燈控制器原理圖.9. 3310LCD的C51驅動程序,還不錯哦.10. 一些常用的LCD的C51驅動程序,還不錯哦.11. NXP示例編碼集,Software that is described herein is for illustrative purposes only which provides custom.12. W78E516B在系統編程的上位機程序.13. 使用TI TMS2407 實現SVPWM同步電機控制算法。源代碼在聞廷科技的2407開發板上已經調試成功。.14. Nucleus Kernel Demo Full source code.15. Altera CPLD矩陣鍵盤一個很不錯的Altera芯片程序.16. 基于Nios II的串口通信.17. The combined demo is dedicated for S1C33L05, so DMT33L05 should be used to load and run the demo. F.18. 用CPLD驅動SJA1000 CAN控制器.19. arm uart program lpc2104.20. 嵌入式PLC,基于8051F020,在KEIL上開發,可實現自定義ＰＬＣ功能.21. SMDK2410原理圖protel格式和PCB庫元件.22. 功能：調用IAP服務程序實現Flash ROM的讀/寫操作向Flash ROM地址1000H處寫入10個數據.23. 功能：由KEY1鍵使系統進入空閑狀態.24. WDT功能試驗。使用WDT的定時模塊.25. 功能：使用CCU的4個模塊.26. s3c2440開發板can控制器sja1000驅動.27. 一個3分頻器。可進一步改裝成實際需要的分頻器使用.28. 通過軟件的仿真,硬件的調試,實現串口數據的傳輸,開發環境包括C-C++與匯編語言..29. ICD2原理圖（protel格式、帶usb接口和串口）.30. 免費分享版網路硬碟 01.創意風格首頁 02.申請會員 03.密碼查詢 04.會員容量限制 05.上傳檔案支援 Persits.Upload Dundas.Upload Lyf.31. 是一個ARM214系連的最小系統原理圖,希望能多愛好嵌入是開發的人員有幫助.32. ds1621的c51驅動,帶12864液晶驅動..33. ds1302時鐘代碼.34. LCD顯示屏驅動sed1335驅動.35. 漢字顯示字庫及程序源碼14*16標準字庫.36. 14*16字庫生成器升級版.37. This utility downloads code to Philip LPC 2000 series chips.38. 測試ARM的外圍SDAM的狀態和參數.39. 三星44B0上的MP3源碼.40. CPLD 的程序,分頻器.
上傳時間： 2013-07-23
上傳用戶：eeworm
asp.net數據庫程序范例
上傳時間： 2013-12-23
上傳用戶：wqxstar
動態聚類k-means演算 將輸入在程式中的數據資料 給予適當的分群
上傳時間： 2015-03-16
上傳用戶：離殤
VC odbc 連接Access 數據庫 的例子,比較簡單
上傳時間： 2014-01-27
上傳用戶：1051290259
大型電子製造業后臺數據庫的詳細設計!絕對好東西!(部分刪減)
上傳時間： 2015-05-17
上傳用戶：aa54
數據庫轉資料程式用SQL SERVER DTS 轉表到另一能SQL DB上
上傳時間： 2014-01-26
上傳用戶：caiiicc
一個數據庫頁面,用jsp寫成。希望給點意見
標簽： jsp
上傳時間： 2013-12-14
上傳用戶：shanml
修改dbf資料庫的好工具.有源碼.對那些仍在使用dos或數據庫為dbase的用戶非常實用.
上傳時間： 2017-01-04
上傳用戶：duoshen1989
