電磁場的Matlab仿真_第1頁
電磁場的Matlab仿真_第2頁
電磁場的Matlab仿真_第3頁
電磁場的Matlab仿真_第4頁
電磁場的Matlab仿真_第5頁
已閱讀5頁,還剩3頁未讀, 繼續(xù)免費閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)

文檔簡介

……………最新資料推薦…………………Matlab與電磁場模擬一單電荷的場分布:單電荷的外部電位計算公式: q40r外輻射的線。MATLAB程序:theta=[0:.01:2*pi]';r=0:10;x=sin(theta)*r;y=cos(theta)*r;plot(x,y,'b')x=linspace(-5,5,100);fortheta=[-pi/40pi/4]y=x*tan(theta);holdon;plot(x,y);endgridon單電荷的等位線和電力線分布圖:……………最新資料推薦…………………二多個點電荷的電場情況:模擬一對同號點電荷的靜電場,其帶電量分別為+Q1和Q2>0)距離為2aP(x,y)處產(chǎn)生的電勢為:由電場強度可得E=-U,在xOy平面上,電場強度的公式為:為了簡單起見,對電勢U做如下變換:。Matlab程序:q=1;xm=2.5;ym=2;……………最新資料推薦…………………x=linspace(-xm,xm);y=linspace(-ym,ym);[X,Y]=meshgrid(x,y);R1=sqrt((X+1).^2+Y.^2);R2=sqrt((X-1).^2+Y.^2);U=1./R1+q./R2;u=1:0.5:4;figurecontour(X,Y,U,u)gridonlegend(num2str(u'))holdonplot([-xm;xm],[0;0])plot([0;0],[-ym;ym])plot(-1,0,'o','MarkerSize',12)plot(1,0,'o','MarkerSize',12)[DX,DY]=gradient(U);quiver(X,Y,-DX,-DY);surf(X,Y,U);同號電荷的靜電場圖像為:11.5250 2.533.540 43020100-2 -2-10 012 2……………最新資料推薦…………………同理,將程序稍作修改,便可以得到異號電荷的靜電場圖像:4040302011.522.533.54100-10-20-30-40-2.52-2-1.5-1-0.500.511.522.5-20三、線電荷產(chǎn)生的電位:z=-Lz=L,q(C/m),xy平面上的電位分布。dVqdLr00/dVqdLr00然后對全部電荷求和即可。xyxy平面上的電位僅取決于離原點的垂直距離R維,只取R為自變量。把R010Nr+1點,對每一點計算其電位。Matlab程序:clearall;L=input(L=’);N=input(N=’);……………最新資料推薦…………………Nr=input(‘分段數(shù)Nr=’);q=input(‘電荷密度q= E0=8.85e-12;C0=1/4/pi/E0;L0=linspace(-L,L,N+1);L1=L0(1:N);L2=L0(2:N+1);Lm=(L1+L2)/2;dL=2*L/N;R=linspace(0,10,Nr+1);fork=1:Nr+1Rk=sqrt(Lm.^2+R(k)^2);Vk=C0*dL*q./Rk;V(k)=sum(Vk);end[max(V),min(V)]plot(R,V),grid線電荷產(chǎn)生的靜電位分布圖:……………最新資料推薦…………………四計算平面上N個電荷之間的庫倫引力1建模:FFqq /r2Fqq(xx)/r3F qq(yy)/r3120x12 2 10yr12 210(xx)2(yy)22 1 2 1先輸入電荷的數(shù)目,各電荷的坐標及電荷量,再選一個電荷,求其它電荷對它的作用力,疊加求合力。再選下一個電荷,依次類推。Matlab程序:clearall;N=input('ê?è?μ?oéêy??N=:');foric=1:N/n??μ?oé#%g\n',ic);rc=input('ê?è?μ?oé????[x,y]£¨?×£?:');x(ic)=rc(1);y(ic)=rc(2);q(ic)=endE0=8.85e-12;C0=1/(4*pi*E0);foric=1:NFx=0.0;Fy=0.0;forjc=1:Nif(ic~=jc)xij=x(ic)-x(jc);yij=y(ic)-y(jc);Rij=sqrt(xij^2+yij^2);Fx=Fx+C0*q(ic)*q(jc)*xij/Rij^3;Fy=Fy+C0*q(ic)*q(jc)*yij/Rij^3;endendend……………最新資料推薦…………………五有限差分法處理電磁場問題Matlab程序:m=40fork=1:mforj=1:mifk==1V(j,k)=1;elseif((j==1)|(j==m)|(k==m))V(j,k)=0;elseendendendcha=0.01;delta=0;

V(j,k)=0.5;……………最新資料推薦…………………n=0;while(1)n=n+1;fork=2:m-1forj=2:m-1Vnew(j,k)=1/4*(V(j+1,k)+V(j-1,k)+V(j,k+1)+V(j,k-1

溫馨提示

  • 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. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論