7、2,z2,選項(xiàng)2,…,xn,yn,zn,選項(xiàng)n)
其中每一組x,y,z組成一組曲線的坐標(biāo)參數(shù)，選項(xiàng)的定義和plot函數(shù)一樣。當(dāng)x,y,z是同維向量時(shí)，那么x,y,z 對(duì)應(yīng)元素構(gòu)成一條三維曲線。當(dāng)x,y,z是同維矩陣時(shí)，那么以x,y,z對(duì)應(yīng)列元素繪制三維曲線，曲線條數(shù)等于矩陣列數(shù)。
5.legend命令：為繪制的圖形加上圖例
調(diào)用格式：legend('string1','string2',...)
例如：legend('電信161班','學(xué)號(hào)：05401111','三','Location','best');
6.xlabel命令：給X軸加標(biāo)題
調(diào)用格式：xlabel('string
8、')
例如：xlabel('x');
三、實(shí)驗(yàn)容
1. 在命令窗口中運(yùn)行一個(gè)加法程序；
2. 在命令窗口中練習(xí)幫助命令〔doc命令〕的使用。
3. 建立第一個(gè)M文件，并運(yùn)行，觀察并保存運(yùn)行結(jié)果。
四、實(shí)驗(yàn)步驟
1. 在命令窗口中運(yùn)行一個(gè)加法程序
〔1〕點(diǎn)擊桌面上matlab7.0快捷方式圖標(biāo)，如圖1.1所示，啟動(dòng)該軟件。
圖1.1 matlab7.0快捷方式圖標(biāo)
〔2〕在翻開(kāi)的界面右方，是命令窗口〔mand Windows〕，如圖1.2所示，在閃動(dòng)光標(biāo)處可以寫入命令；
圖1.2 Matlab的命令窗口
〔3〕在光標(biāo)處寫入如圖1.3所示的命令〔注意：前
9、兩個(gè)語(yǔ)句后面有分號(hào)，最后一個(gè)語(yǔ)句沒(méi)有分號(hào)〕；按回車鍵，那么得到運(yùn)行結(jié)果為50，如圖1.4所示。
圖1.3 在命令窗口輸入命令
圖1.4 按回車鍵執(zhí)行命令得到正確運(yùn)行結(jié)果
2. 在命令窗口中練習(xí)幫助命令〔doc命令〕的使用。
在命令窗口光標(biāo)處輸入命令：doc plot；回車，那么進(jìn)入在線幫助文件，顯示plot命令的使用方法頁(yè)面，如圖1.5所示。
圖1.5 plot命令的在線幫助頁(yè)面
3. 建立第一個(gè)M文件，并運(yùn)行，觀察并保存運(yùn)行結(jié)果。
〔1〕點(diǎn)擊圖標(biāo)，如圖1.6中紅色圓圈所示，即創(chuàng)立了一個(gè)新的M文件，如圖1.7所示。
圖1.6 紅色圓圈的圖標(biāo)用于創(chuàng)立新的M文件
10、
圖1.7創(chuàng)立的空白M文件
〔2〕在空白M文件中輸入“二、實(shí)驗(yàn)原理〞例子的程序，保存，運(yùn)行，得到運(yùn)行結(jié)果如圖1.8所示。
要求：在E盤建立新文件夾，命名為Fiele_Wave_simulation_2012_10_27；將M文件保存在Fiele_Wave_simulation_2012_10_27目錄下，命名為Exp_1.m；
特別說(shuō)明兩點(diǎn)：a. M文件名及保存的路徑名均應(yīng)為英文，否那么運(yùn)行出錯(cuò)；
b.程序中的所有字符均應(yīng)為英文狀態(tài)下輸入，特別注意單引號(hào)，逗號(hào)，空格，這些細(xì)節(jié)會(huì)導(dǎo)致運(yùn)行報(bào)錯(cuò)，又極難發(fā)現(xiàn)。
圖1.8 M文件的保存、運(yùn)行按鍵及運(yùn)行結(jié)果
五、實(shí)驗(yàn)工具
1．計(jì)
11、算機(jī) 1臺(tái)
2．MATLAB 仿真軟件 1套
六、實(shí)驗(yàn)報(bào)告要求
1．寫出仿真程序源代碼。
2．在同一窗口用不同的線性繪制y=sinx，y=cosx在[0,2*pi]上的圖像，并加標(biāo)注。
3. 在同一窗口用不同的線性繪制y=sin2x，y=cos2x在[-2*pi,2*pi]上的圖像，并加標(biāo)注?！惨趫D中繪制出與學(xué)號(hào)〕
實(shí)驗(yàn)二 單電荷的場(chǎng)分布
一、實(shí)驗(yàn)?zāi)康?
1. 掌握 MATLAB 仿真的根本流程與步驟；
2. 學(xué)會(huì)繪制單電荷的等位線和電力線分布圖。
二、實(shí)驗(yàn)原理
1．根本原理
單電荷的外部電位
12、計(jì)算公式：
等位線就是連接距離電荷等距離的點(diǎn)，在圖上表示就是一圈一圈的圓，而電
力線就是由點(diǎn)向外輻射的線，比擬簡(jiǎn)單，這里就不再贅述。
2. 參考程序
theta=[0:0.01:2*pi]';%創(chuàng)立向量theta，確定theta的圍為0到2*pi，步距為0.01
r=0:10; %半徑分別為r=0,1,2,...,10
x=sin(theta)*r;
y=cos(theta)*r;
plot(x,y,'b') %繪制10個(gè)圓
x=linspace(-5,5,100);%創(chuàng)立線性空間向量x,從-5到5,等間距分為100個(gè)點(diǎn)
13、
for theta=[-pi/4 0 pi/4]
y=x*tan(theta); %分別繪制y=x*tan(theta)的三條直線，其中theta分別取-pi/4，0，pi/4
hold on;%Hold住繪制的圖形
plot(x,y); %繪制y=x*tan(theta)的三條直線
end
grid on
% axis tight
% legend('電信本162班','學(xué)號(hào)：16401111','三','Location','best'); legend('boxoff');%加上圖例
% xlabel('x');%加上橫
14、坐標(biāo)標(biāo)題
% ylabel('y');;%加上縱坐標(biāo)標(biāo)題
% hold on;
3.程序參考運(yùn)行結(jié)果
運(yùn)行程序，獲得圖像大致如圖2.1所示。
圖2.1 單電荷的等位線和電力線分布圖
三、實(shí)驗(yàn)容
繪制單電荷的等位線和電力線分布圖。
四、實(shí)驗(yàn)步驟
1．在E盤建立新文件夾，命名為Fiele_Wave_simulation_2012_10_27；
2.翻開(kāi)Matlab軟件，新建一個(gè)空白的M文件，保存在Fiele_Wave_simulation_2012_10_27目錄下，命名為Exp_2.m；
3.將源程序拷貝到M文件中，保存；
4.點(diǎn)擊運(yùn)行按鈕，觀察程序運(yùn)行結(jié)果
15、。
五、實(shí)驗(yàn)工具
1．計(jì)算機(jī) 1臺(tái)
2．MATLAB 仿真軟件 1套
六、實(shí)驗(yàn)報(bào)告要求
1．寫出仿真程序源代碼。
2．繪制單電荷的等位線和電力線分布圖。〔要在圖中繪制出與學(xué)號(hào)〕
實(shí)驗(yàn)三 點(diǎn)電荷電場(chǎng)線的圖像
一、實(shí)驗(yàn)?zāi)康?
學(xué)會(huì)由解析表達(dá)式進(jìn)展數(shù)值求解的方法。
二、實(shí)驗(yàn)原理
1.根本原理
考慮一個(gè)三點(diǎn)電荷系所構(gòu)成的系統(tǒng)。如下圖，其中一個(gè)點(diǎn)電荷-q位于坐標(biāo)原點(diǎn)，另一個(gè)-q位于y軸上的點(diǎn)，最后一個(gè)+2q位于y軸的-點(diǎn)，那么在xoy平面，電場(chǎng)強(qiáng)度應(yīng)滿足
任意
16、條電場(chǎng)線應(yīng)該滿足方程：
〔1〕
求解〔1〕式可得：
〔2〕
這就是電場(chǎng)線滿足的方程，常數(shù)C取不同值將得到不同的電場(chǎng)線。
2.參考程序
解出y=f(x) 的表達(dá)式再作圖是不可能的。用Matlab語(yǔ)言即能輕松的做到這
一點(diǎn)。其語(yǔ)句是：
syms x y % 設(shè)置x，y變量；
for C=0:0.1:3.0
ezplot(2*(y+1)/sqrt((y+1)^2+x^2)-y/sqrt(y^2+x^2)-
(y-1)/sqrt((y-1)^2+x^2)-C, [- 5,5,0.1]); %其中取了a=1，C=
17、0，0.1，0.2，……，3.0
hold on;
end
3.程序參考運(yùn)行結(jié)果
運(yùn)行程序，獲得圖像大致如圖3.1所示。
圖3.1 點(diǎn)電荷電場(chǎng)線的圖像
三、實(shí)驗(yàn)容
根據(jù)給出的三點(diǎn)電荷系所構(gòu)成的系統(tǒng)電場(chǎng)線滿足的方程，繪制其圖像。
四、實(shí)驗(yàn)步驟
1．在E盤建立新文件夾，命名為Fiele_Wave_simulation_2012_10_27；
2.翻開(kāi)Matlab軟件，新建一個(gè)空白的M文件，保存在Fiele_Wave_simulation_2012_10_27目錄下，命名為Exp_3.m；
3.將源程序拷貝到M文件中，保存；
4.點(diǎn)擊運(yùn)行按鈕，觀察程序運(yùn)行結(jié)果
18、。
五、實(shí)驗(yàn)工具
1．計(jì)算機(jī) 1臺(tái)
2．MATLAB 仿真軟件 1套
六、實(shí)驗(yàn)報(bào)告要求
1．寫出仿真程序源代碼。
2．繪制三點(diǎn)電荷系所構(gòu)成的系統(tǒng)電場(chǎng)線的圖像。〔要在圖中繪制出與學(xué)號(hào)〕
實(shí)驗(yàn)四 線電荷產(chǎn)生的電位
一、實(shí)驗(yàn)?zāi)康?
理解交互式程序運(yùn)行的過(guò)程。
二、實(shí)驗(yàn)原理
1.根本原理
設(shè)電荷均勻分布在從z=-L到z=L,通過(guò)原點(diǎn)的線段上，其密度為q(單位C/m),
求在xy平面上的電位分布。
點(diǎn)電荷產(chǎn)生的電位可表示為
是一個(gè)標(biāo)量。其中r為
19、電荷到測(cè)量點(diǎn)的距離。線電荷所產(chǎn)生的電位可用積分或疊
加的方法來(lái)求。為此把線電荷分為N段，每段長(zhǎng)為dL。每段上電荷為q*dL,看作
集中在中點(diǎn)的點(diǎn)電荷，它產(chǎn)生的電位為
然后對(duì)全部電荷求和即可。
2.參考程序
把xy平面分成網(wǎng)格，因?yàn)閤y平面上的電位僅取決于離原點(diǎn)的垂直距離R，所
以可以省略一維，只取R為自變量。把R從0到10米分成Nr+1點(diǎn)，對(duì)每一點(diǎn)計(jì)算其電位。
clear all;
L=input(‘線電荷長(zhǎng)度L＝：’);
N=input(‘分段數(shù)N＝：’);
Nr=input(‘分段數(shù)Nr＝：’);
q=input(‘電荷密度q=:’);
E0=8
20、.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);
for k=1:Nr+1
Rk=sqrt(Lm.^2+R(k)^2);
Vk=C0*dL*q./Rk;
V(k)=sum(Vk);
end
[max(V),min(V)]
plot(R,V),grad
3.程序參考運(yùn)行結(jié)果
輸入：
線電荷長(zhǎng)度L＝：5
分段數(shù)N＝：50
分段數(shù)Nr＝：50
電荷密度q=:1
21、可得最大值和最小值為:
ans =
1.0e+010 *[9.3199 0.8654]
圖像大致如圖4.1所示。
圖4.1 點(diǎn)電荷電場(chǎng)線的圖像
三、實(shí)驗(yàn)容
根據(jù)給出的線電荷產(chǎn)生的電位方程，繪制電位圖像。
四、實(shí)驗(yàn)步驟
1．在E盤建立新文件夾，命名為Fiele_Wave_simulation_2012_10_27；
2.翻開(kāi)Matlab軟件，新建一個(gè)空白的M文件，保存在Fiele_Wave_simulation_2012_10_27目錄下，命名為Exp_4.m；
3.將源程序拷貝到M文件中，保存；
4.點(diǎn)擊運(yùn)行按鈕，觀察程序運(yùn)行結(jié)果。
五、實(shí)驗(yàn)工具
1．
22、計(jì)算機(jī) 1臺(tái)
2．MATLAB 仿真軟件 1套
六、實(shí)驗(yàn)報(bào)告要求
1．寫出仿真程序源代碼。
2．繪制線電荷產(chǎn)生的電位的圖像?！惨趫D中繪制出與學(xué)號(hào)〕
實(shí)驗(yàn)五 有限差分法處理電磁場(chǎng)問(wèn)題
一、實(shí)驗(yàn)?zāi)康?
理解有限差分法處理電磁場(chǎng)問(wèn)題。
二、實(shí)驗(yàn)原理
1.根本原理
在很多實(shí)際情況下，往往我們不知道電荷的分布情況，而只知道邊界的電位。
例如在靜電場(chǎng)中的導(dǎo)體，邊界是等位面，但是外表上的電荷分布往往不一樣，而且很難求得。
如果我們想求導(dǎo)體附近的電場(chǎng)，這里就介紹用差分
23、法解電場(chǎng)的方法。
如果我們想求導(dǎo)體附近的電場(chǎng)，這里就介紹用差分法解電場(chǎng)的方法。
解決這個(gè)問(wèn)題的關(guān)鍵是對(duì)電位使用以下結(jié)論：在一個(gè)沒(méi)有電荷的區(qū)域，給出一個(gè)
點(diǎn)的電位等于周圍點(diǎn)的電位數(shù)值的平均值。我們使用高斯定律以及以下公式來(lái)證
明這個(gè)結(jié)論：
〔根據(jù)電壓V的電場(chǎng)分量〕〔5.1〕
我們將集中討論該情況，其中電位只取決于兩個(gè)坐標(biāo)，x和y。一個(gè)例子是一個(gè)帶電的長(zhǎng)圓柱體。在一個(gè)點(diǎn)的電勢(shì)只依賴于這一點(diǎn)在平面垂直于圓柱體的軸線，而不是z坐標(biāo)。對(duì)于這樣一個(gè)二維的情況下，考慮一個(gè)點(diǎn)P的坐標(biāo)〔x，Y，Z〕，并在一個(gè)由高斯外表封閉的立方體的一面長(zhǎng)度是2Δl，中心在P〔圖T8.1〕。如果立方體的部沒(méi)有電荷，
24、通過(guò)立方體的電通量ΦE等于零。由方程 (T8.1) 可知Z軸的電場(chǎng)分量為零，因?yàn)殡妱?shì)V并不是z的函數(shù)。因此，并沒(méi)有通過(guò)高斯外表的平行于xy -平面的電通量。由于是一個(gè)小的立方體，通過(guò)其他四個(gè)面每通量有一個(gè)良好的逼近，等于在每面的中心和每一面的〔2Δl〕2的E的垂直分量的乘積 ?？偭髁俊驳扔?〕可以表示為；
〔5.2〕
在一個(gè)沒(méi)有電荷的區(qū)域，點(diǎn)P的電位數(shù)值等于P點(diǎn)周圍電位值的平均值。
圖5.1 沒(méi)有電荷的區(qū)域，點(diǎn)P的電位的求解示意圖
使用式5.1，我們可以近似的寫出電場(chǎng)的各分量：
〔5.3〕
我們可以得到P點(diǎn)的電位是：
〔5.4〕
總之，P點(diǎn)的電位值等于P點(diǎn)周圍點(diǎn)的電位值的
25、平均值，前提是非常小。
2.計(jì)算機(jī)繪圖算法
3.參考程序
計(jì)算場(chǎng)點(diǎn)電壓值程序如下：
m=12
for k=1:m
for j=1:m
if k==1
V(j,k)=1;
elseif((j==1)|(j==m)|(k==m))
V(j,k)=0;
else
V(j,k)=0.5;
end
end
end
cha=0.01;
delta=0;
n=0;
wh
26、ile(1)
n=n+1;
for k=2:m-1
for j=2:m-1
Vnew(j,k)=1/4*(V(j+1,k)+V(j-1,k)+V(j,k+1)+V(j,k-1));
d=abs((Vnew(j,k)-V(j,k))/V(j,k));
if d>delta
delta=d;
end
V(j,k)=Vnew(j,k);
end
end
if delta100)
27、
break;
end
delta=0.;
end
代入m＝22
繪圖程序：
k=1:m;
j=1:m;
[DX,DY] = gradient(V,.4,.4);
hold on
quiver(k,j,DX,DY,2)
hold off
k=1:m;
j=1:m;
[DX,DY]=gradient(V,.4,.4)
A=(DX.^2+DY.^2).^0.5;
[DA,DB]=gradient(A,.4,.4);
hold on
quiver(k,j,DA,DB,2)
hold off
3.程序參考運(yùn)行
28、結(jié)果
電場(chǎng)圖和電力線圖分別如圖5.2、圖5.3所示。
圖5.2 電場(chǎng)線的圖像
圖5.3 電力線的圖像
三、實(shí)驗(yàn)容
有限差分法繪制帶電長(zhǎng)圓柱體的電位和電力線圖像。
四、實(shí)驗(yàn)步驟
1．在E盤建立新文件夾，命名為Fiele_Wave_simulation_2012_10_27；
2.翻開(kāi)Matlab軟件，新建一個(gè)空白的M文件，保存在Fiele_Wave_simulation_2012_10_27目錄下，命名為Exp_5.m、Exp_5_1.m；
3.將源程序拷貝到M文件中，保存；
4.點(diǎn)擊運(yùn)行按鈕，觀察程序運(yùn)行結(jié)果。
五、實(shí)驗(yàn)工具
1．計(jì)算機(jī) 1臺(tái)
2．MATLAB 仿真軟件 1套
六、實(shí)驗(yàn)報(bào)告要求
1．寫出仿真程序源代碼。
2．繪制帶電長(zhǎng)圓柱體的電位和電力線圖像?！惨趫D中繪制出與學(xué)號(hào)〕
27 / 27
2016【電磁場(chǎng)和電磁波】仿真實(shí)驗(yàn)
