1. 程式人生 > >有限差分法的MATLAB程式!!!!

有限差分法的MATLAB程式!!!!

設矩形金屬槽,長為10cm,寬為6cm,其側板和底板接地,蓋板電壓為5V(相對值)。 求槽內場(電位)的分佈? 迭代解法:(取最小步長為0.1cm) 在這裡插入圖片描述執行效果圖如下:在這裡插入圖片描述

原始碼如下: lx=61;ly=101; v1=zeros(ly,lx); for j=2:lx-1 v1(ly,j)=5; end v2=v1;m=1;t=0;k=0; while(m>0.00001) k=k+1 m=0; for i=2:ly-1 for j=2:lx-1; v2(i,j)=(v1(i,j+1)+v1(i+1,j)+v1(i-1,j)+v1(i,j-1))/4; t=(v2(i,j)-v1(i,j)); if(t>m)m=t;end end end v1=v2; end subplot(1,2,1),mesh(v2) axis([0,61,0,101,0,5]) subplot(1,2,2),contour(v2,32)

這個實驗我今天上午電磁場課程的原實驗,我課上沒做出來,所以下午自己找教程安裝了MATLAB 2016A,自己做了一遍各方面還是很有收穫。我是一個硬體技術菜鳥,希望能和大家一起持續學習。