測(cè)控技術(shù)與儀器課程設(shè)計(jì)電流環(huán)磁場(chǎng)的分布

1、測(cè)試與光電工程學(xué)院課程設(shè)計(jì)任務(wù)書 測(cè)控技術(shù)與儀器 系 090812 班 學(xué)號(hào) 09081209 姓名 王若玨 課題名稱： 電流環(huán)磁場(chǎng)的分布 課題要求：1. 建立電流環(huán)磁場(chǎng)模型2. 掌握?qǐng)D形的繪制3. 編寫圖形顯示界面4. 顯示電流環(huán)磁場(chǎng)分布圖課題進(jìn)程：第一周：確定使用哪種計(jì)算機(jī)語言編程； 編寫主要功能。第二周：程序的完善、運(yùn)行與調(diào)試；撰寫課程設(shè)計(jì)報(bào)告。 指導(dǎo)教師： 袁麗華 【摘要】 matlab是由美國(guó)mathworks公司發(fā)布的主要面對(duì)科學(xué)計(jì)算、可視化以及交互式程序設(shè)計(jì)的高科技計(jì)算環(huán)境。它將數(shù)值分析、矩陣計(jì)算、科學(xué)數(shù)據(jù)可視化以及非線性動(dòng)態(tài)系統(tǒng)的建模和仿真等諸多強(qiáng)大功能集成在一個(gè)易于使用的視窗

2、環(huán)境中，為科學(xué)研究、工程設(shè)計(jì)以及必須進(jìn)行有效數(shù)值計(jì)算的眾多科學(xué)領(lǐng)域提供了一種全面的解決方案，并在很大程度上擺脫了傳統(tǒng)非交互式程序設(shè)計(jì)語言（如c、fortran）的編輯模式，代表了當(dāng)今國(guó)際科學(xué)計(jì)算軟件的先進(jìn)水平。 本課程設(shè)計(jì)根據(jù)畢奧一薩伐爾定律推導(dǎo)出圓環(huán)電流磁場(chǎng)分布的積分表示，利用matlab的符號(hào)積分得出計(jì)算結(jié)果，并繪制磁場(chǎng)分布曲線。并通過gui制作可交互式界面進(jìn)行仿真與演示。本課程設(shè)計(jì)選用matlab r2008a軟件。 關(guān)鍵詞： 圓環(huán)電流；磁場(chǎng)；matlab；符號(hào)積分；gui 目錄一 引言11.1 課題簡(jiǎn)介11.2 課題制作要點(diǎn)11.3 matlab軟件1二 理論分析22.1 畢奧一薩伐爾

3、定律22.2 電流元在圓環(huán)外的磁場(chǎng)22.3 電流元在圓環(huán)內(nèi)的磁場(chǎng)32.4 距離圓心z處的磁感應(yīng)強(qiáng)度3三 程序設(shè)計(jì)43.1 gui界面設(shè)計(jì)43.1.1 創(chuàng)建gui文件43.1.2 gui布局與控件設(shè)計(jì)43.1.3 gui菜單欄設(shè)計(jì)63.2 功能實(shí)現(xiàn)63.2.1 菜單功能63.2.2 繪圖設(shè)置相關(guān)控件功能實(shí)現(xiàn)7四 gui運(yùn)行效果114.1 gui初始化效果114.2 菜單欄效果114.3 繪圖效果12五 心得體會(huì) 13六 參考文獻(xiàn)131 引言1.1 課題簡(jiǎn)介圓電流是一個(gè)經(jīng)典的物理模型。對(duì)圓電流內(nèi)磁場(chǎng)分布特性的研究，在宏觀上可以改善線圈繞組在電機(jī)中的應(yīng)用。微觀上?？梢越o分子電流提供理論基礎(chǔ)，對(duì)進(jìn)一步

4、研究磁介質(zhì)的磁性特征也有重要參考價(jià)值。本課題主要著重于對(duì)電流環(huán)磁場(chǎng)進(jìn)行建模與仿真，根據(jù)畢奧一薩伐爾定律通過設(shè)定電流環(huán)半徑、電流強(qiáng)度，模擬不同情況下磁場(chǎng)與磁感線的分布情況，并繪制成圖像。 1.2 課題制作要點(diǎn) 電流環(huán)磁場(chǎng)參數(shù)的設(shè)定 電流環(huán)磁場(chǎng)的建模與仿真 可交互式gui界面的設(shè)計(jì)1.3 matlab軟件簡(jiǎn)介matlab是由美國(guó)mathworks公司發(fā)布的主要面對(duì)科學(xué)計(jì)算、可視化以及交互式程序設(shè)計(jì)的高科技計(jì)算環(huán)境。它將數(shù)值分析、矩陣計(jì)算、科學(xué)數(shù)據(jù)可視化以及非線性動(dòng)態(tài)系統(tǒng)的建模和仿真等諸多強(qiáng)大功能集成在一個(gè)易于使用的視窗環(huán)境中，為科學(xué)研究、工程設(shè)計(jì)以及必須進(jìn)行有效數(shù)值計(jì)算的眾多科學(xué)領(lǐng)域提供了一種全面

5、的解決方案，并在很大程度上擺脫了傳統(tǒng)非交互式程序設(shè)計(jì)語言（如c、fortran）的編輯模式，代表了當(dāng)今國(guó)際科學(xué)計(jì)算軟件的先進(jìn)水平。matlab是一個(gè)高級(jí)的矩陣/陣列語言，它包含控制語句、函數(shù)、數(shù)據(jù)結(jié)構(gòu)、輸入和輸出和面向?qū)ο缶幊烫攸c(diǎn)。用戶可以在命令窗口中將輸入語句與執(zhí)行命令同步，也可以先編寫好一個(gè)較大的復(fù)雜的應(yīng)用程序（m文件）后再一起運(yùn)行。而且這種語言可移植性好、可拓展性極強(qiáng)，這也是matlab能夠深入到科學(xué)研究及工程計(jì)算各個(gè)領(lǐng)域的重要原因。matlab是一個(gè)包含大量計(jì)算算法的集合。其擁有600多個(gè)工程中要用到的數(shù)學(xué)運(yùn)算函數(shù)，可以方便的實(shí)現(xiàn)用戶所需的各種計(jì)算功能。函數(shù)中所使用的算法都是科研和工程

6、計(jì)算中的最新研究成果，而前經(jīng)過了各種優(yōu)化和容錯(cuò)處理。在通常情況下，可以用它來代替底層編程語言，如c和c+ 。在計(jì)算要求相同的情況下，使用matlab的編程工作量會(huì)大大減少。matlab自產(chǎn)生之日起就具有方便的數(shù)據(jù)可視化功能，以將向量和矩陣用圖形表現(xiàn)出來，并且可以對(duì)圖形進(jìn)行標(biāo)注和打印。新版本的matlab對(duì)整個(gè)圖形處理功能作了很大的改進(jìn)和完善，使它不僅在一般數(shù)據(jù)可視化軟件都具有的功能而且對(duì)于一些其他軟件所沒有的功能，matlab同樣表現(xiàn)了出色的處理能力。同時(shí)對(duì)一些特殊的可視化要求，另外新版本的matlab還著重在圖形用戶界面（gui）的制作上作了很大的改善，對(duì)這方面有特殊要求的用戶也可以得到滿足

7、。 2 理論分析2.1 畢奧一薩伐爾定律畢奧一薩伐爾定律是以實(shí)驗(yàn)為基礎(chǔ)經(jīng)過科學(xué)抽象而得到的，描述的是電流元在空間任一點(diǎn)產(chǎn)生的磁感應(yīng)強(qiáng)度。原則上利用畢奧一薩伐爾定律并結(jié)合磁感應(yīng)強(qiáng)度疊加原理，可以計(jì)算任意形狀的電流所產(chǎn)生的磁場(chǎng)。 畢奧薩伐爾定律指出：載流導(dǎo)線上的電流元在真空中某點(diǎn)p的磁感度db的大小與電流元的大小成正比，與電流元和從電流元到p點(diǎn)的位矢r之間的夾角 的正弦成正比，與位矢r的大小的平方成反比，即： 上式中，為位矢的單位矢量。此即畢奧薩伐爾定律的公式表述。2.2 電流元在圓環(huán)外的磁場(chǎng) 圖1 電流元在電流環(huán)外的磁場(chǎng)如圖1所示，根據(jù)畢奧一薩伐爾定律，任一電流元在p點(diǎn)產(chǎn)生的磁感應(yīng)強(qiáng)度 ，考慮到

8、，, 有（1）。2.3 電流元在圓環(huán)內(nèi)的磁場(chǎng) 圖2 電流元在電流環(huán)內(nèi)的磁場(chǎng) 如圖2所示，電流元在p點(diǎn)所產(chǎn)生的磁感應(yīng)強(qiáng)度，即（2）。2.4 距離圓心z處的磁感應(yīng)強(qiáng)度比較電流元在圓環(huán)內(nèi)、外的磁場(chǎng)表達(dá)式，兩式僅相差一個(gè)負(fù)號(hào)。該負(fù)號(hào)體現(xiàn)了兩種情況下磁場(chǎng)方向正好相反。若規(guī)定垂直紙面向外為磁感應(yīng)強(qiáng)度的正方向，則電流元在圓環(huán)內(nèi)、外產(chǎn)生的磁感應(yīng)強(qiáng)度可統(tǒng)一表示為( 2 )式的形式。結(jié)合磁感應(yīng)強(qiáng)度疊加原理，整個(gè)圓環(huán)電流在距圓心處產(chǎn)生的磁感應(yīng)強(qiáng)度為： 3 程序設(shè)計(jì)3.1 gui界面設(shè)計(jì)3.1.1 創(chuàng)建gui文件打開matlab軟件后，點(diǎn)擊菜單欄中的“guide”按鈕，進(jìn)入“guide quick start”，選擇

9、“create new gui”選項(xiàng)卡，選擇“blank gui（default）”創(chuàng)建一空白gui文檔。并命名為“cichang1.fig”如圖3所示。圖3 創(chuàng)建一空白gui3.1.2 gui布局與控件設(shè)計(jì) gui布局如圖4，所用控件羅列于表1。圖4 gui布局控件名稱控件屬性控件用途fontsizestring（title）tagpanel15繪圖設(shè)置uipanel4放置繪圖設(shè)置的其他控件panel15視角切換uipanel3放置用于視角切換的控件axesaxes1顯示磁場(chǎng)及磁感線分布圖像static text25帶電電流環(huán)周圍空間的磁場(chǎng)分布text1標(biāo)題static text15電流值t

10、ext3顯示static text15atext4電流值單位static text15環(huán)半徑text5顯示static text15mtext6環(huán)半徑單位listbox15電流環(huán)的磁場(chǎng)分布圖;電流環(huán)的磁感線分布圖listbox1選擇繪圖類型push button15繪制圖像pushbutton2繪制圖像push button15重置參數(shù)pushbutton7重置參數(shù)edit text15edit3設(shè)置電流參數(shù)值edit text15edit4設(shè)置環(huán)半徑數(shù)值edit text15測(cè)試圖edit1顯示圖像名稱radio button12自由視角radiobutton1開啟三維圖像自由視角radio

11、 button12鎖定視角radiobutton1關(guān)閉三維圖像自由視角表1 所用控件3.1.3 gui菜單欄設(shè)計(jì)為使用戶獲得更方便的使用體驗(yàn)，故增加菜單欄。點(diǎn)擊“menu editor”進(jìn)行制作。如圖5。圖5 菜單欄編輯3.2 功能實(shí)現(xiàn)gui界面設(shè)計(jì)過后，對(duì)相關(guān)控件的“view callbacks”中的函數(shù)進(jìn)行編輯，以實(shí)現(xiàn)既定功能。3.2.1 菜單功能（1） 菜單欄中的退出功能 點(diǎn)擊“menu editor”中的“退出”選項(xiàng)，設(shè)置tag屬性為“untitled_1”，并對(duì)回調(diào)函數(shù)進(jìn)行編輯。如圖6.圖6 退出功能編輯程序代碼如下：function untitled_1_callback(hobj

12、ect, eventdata, handles)selection=questdlg(是否關(guān)閉,get(gcf,name),窗口？, . close ,get(gcf,name),.,是,否,是);if strcmp(selection,否) return; else clc;%當(dāng)選擇關(guān)閉時(shí)，清空所有matla輸入面上的所有錯(cuò)誤信息，同時(shí)關(guān)閉圖像窗口 clear all;delete(gcf); end（2） 菜單欄中的說明功能 點(diǎn)擊“menu editor”中的“說明”選項(xiàng)，設(shè)置tag屬性為“untitled_2”，并對(duì)回調(diào)函數(shù)進(jìn)行編輯。程序代碼如下：function untitled_2_

13、callback(hobject, eventdata, handles)msgbox(利用畢奧薩伐爾定律作圖,說明)3.2.2 繪圖設(shè)置相關(guān)控件功能實(shí)現(xiàn) （1） 選擇繪制圖像類型 當(dāng)選擇“l(fā)istbox1”中的某一選項(xiàng)時(shí)，坐標(biāo)軸下方的“edittext”也顯示相應(yīng)的文字說明。鼠標(biāo)右鍵點(diǎn)擊“l(fā)istbox1”選擇“view callbacks”中的“callback”進(jìn)行編輯。如圖7。圖7 程序代碼如下：function listbox1_callback(hobject, eventdata, handles)selected.index=get(handles.listbox1,value

14、);str=get(handles.listbox1,string);set(handles.edit1,string,strselected.index);（2） 繪制圖像按鈕功能實(shí)現(xiàn) 當(dāng)使用者選擇好所要繪制的圖像類型后，點(diǎn)擊“繪制圖像”按鈕，則在坐標(biāo)軸中出現(xiàn)相應(yīng)的圖像。圖像下方的文本框中顯示圖像名稱。對(duì)“pushbutton 2”的回調(diào)函數(shù)進(jìn)行編輯。程序代碼如下：function pushbutton2_callback(hobject, eventdata, handles)axes(handles.axes1);cla;popup_sel_index = get(handles.lis

15、tbox1, value);rh=str2double(get(handles.edit4, string); %通過edit4輸入環(huán)半徑值i0=str2double(get(handles.edit3, string); %通過edit3輸入電流值switch popup_sel_index case 1mu0=4*pi*1e-7; c0=mu0/(4*pi)*i0; nh=40; %將電流分段為40段x=linspace(-3,3,nh);y=x; %確定觀察點(diǎn)的范圍 thta0=linspace(0,2*pi,nh+1); %環(huán)周的分段數(shù) thta1=thta0(1:nh);y1=rh*

16、cos(thta1); %環(huán)各段矢量的起點(diǎn)坐標(biāo)y1,z1;z1=rh*sin(thta1);thta2=thta0(2:nh+1); y2=rh*cos(thta2); %環(huán)各段矢量的起點(diǎn)坐標(biāo)y2,z2; z2=rh*sin(thta2); xc=0;yc=(y2+y1)/2;zc=(z2+z1)/2; %計(jì)算環(huán)各段矢量中點(diǎn)的3個(gè)坐標(biāo)分量 dlx=0;dly=y2-y1;dlz=z2-z1; %計(jì)算環(huán)各段矢量dl的3個(gè)長(zhǎng)度分量 ngx=nh;ngy=ngx; %確定網(wǎng)格數(shù)for i=1:ngy %循環(huán)計(jì)算各網(wǎng)格點(diǎn)上的b(x,y)值 for j=1:ngx rx=x(j)-xc;ry=y(i)-

17、yc;rz=0-zc; r3=sqrt(rx.2+ry.2+rz.2).3; dlxr_x=dly.*rz-dlz.*ry; dlxr_y=dly.*rx-dlx.*rz; bx(i,j)=sum(c0*dlxr_x./r3); %把環(huán)各段產(chǎn)生的磁場(chǎng)分量累加 by(i,j)=sum(c0*dlxr_y./r3); b=(bx.2+by.2).0.5; %計(jì)算b的大小 endendsurf(x,y,b); %繪制磁場(chǎng)分布圖像shading interpgrid on case 2mu0=4*pi*1e-7; c0=mu0/(4*pi)*i0;nh=40; x=linspace(-3,3,nh);

18、y=x; thta0=linspace(0,2*pi,nh+1); thta1=thta0(1:nh); y1=rh*cos(thta1); z1=rh*sin(thta1); thta2=thta0(2:nh+1);y2=rh*cos(thta2); z2=rh*sin(thta2); xc=0;yc=(y2+y1)/2;zc=(z2+z1)/2; dlx=0;dly=y2-y1;dlz=z2-z1; ngx=nh;ngy=ngx; for i=1:ngy for j=1:ngx rx=x(j)-xc;ry=y(i)-yc;rz=0-zc; r3=sqrt(rx.2+ry.2+rz.2).3

19、; dlxr_x=dly.*rz-dlz.*ry; dlxr_y=dly.*rx-dlx.*rz; bx(i,j)=sum(c0*dlxr_x./r3); by(i,j)=sum(c0*dlxr_y./r3); b=(bx.2+by.2).0.5; end end quiver(x,y,bx,by,0.5,linewidth,1); %繪制磁感線圖線 grid on end（3） 重置參數(shù)按鈕功能實(shí)現(xiàn) 初始化gui圖形界面：坐標(biāo)軸中的圖像還原為初始時(shí)的“測(cè)試圖”，坐標(biāo)軸下方的文本框顯示文字“測(cè)試圖”，edittext3、edittext4無顯示，listbox1中選項(xiàng)顯示為“電流環(huán)的磁場(chǎng)分布圖

20、”，視角切換中的單選框選中“視角鎖定”。首先對(duì)cichang1.fig中的“openfunction”進(jìn)行編輯。程序代碼如下：function cichang1_openingfcn(hobject, eventdata, handles, varargin)guidata(hobject, handles); if strcmp(get(hobject,visible),off) surf(peaks); rotate3d off;end然后對(duì)pushbutton7的回調(diào)函數(shù)進(jìn)行編輯。程序代碼如下：function pushbutton7_callback(hobject, eventdat

21、a, handles)cla;surf(peaks);rotate3d off;set(handles.radiobutton2,value,1);set(handles.radiobutton1,value,0);set(handles.edit3,string, );set(handles.listbox1,value,1);set(handles.edit4,string, );set(handles.edit1,string,測(cè)試圖);（4） 視角切換功能實(shí)現(xiàn) 當(dāng)選擇“自由視角”時(shí)，用戶可通過鼠標(biāo)自由旋轉(zhuǎn)坐標(biāo)軸中的圖像；當(dāng)選擇“鎖定視角”時(shí)，用戶可以鎖定任意視角時(shí)的圖像。 首先對(duì)rad

22、iobutton2的回調(diào)函數(shù)進(jìn)行編輯。程序代碼如下：function radiobutton2_callback(hobject, eventdata, handles)rotate3d off; 然后對(duì)radiobutton1的回調(diào)函數(shù)進(jìn)行編輯。程序代碼如下：function radiobutton1_callback(hobject, eventdata, handles)rotate3d on; 最后再保存一次“cichang1.fig”。4 gui運(yùn)行效果4.1 gui初始化效果打開“cichang1.fig ”，顯示如圖8效果。圖84.2 菜單欄效果點(diǎn)擊菜單欄“退出”選項(xiàng)，顯示一對(duì)話框，如點(diǎn)擊“是