clc;
clear;
close?all;
global?a?b?r?h?n1?n2?x0?y0
u0=4*pi*10^（-7）;
a=input（‘請(qǐng)輸入方形線圈的長:‘）;
b=input（‘請(qǐng)輸入方形線圈的寬:‘）;
r=input（‘請(qǐng)輸入圓形線圈的半徑:‘）;
h=input（‘請(qǐng)輸入兩線圈的豎直高度:‘）;
n1=input（‘請(qǐng)輸入方形線圈的匝數(shù):‘）;
n2=input（‘請(qǐng)輸入圓形線圈的匝數(shù):‘）;
x0=input（‘請(qǐng)輸入圓形線圈中心點(diǎn)的橫坐標(biāo):‘）;
la=input（‘請(qǐng)輸入圓形線圈中心點(diǎn)縱坐標(biāo)的初始值:‘）;
dl=input（‘請(qǐng)輸入縱坐標(biāo)移動(dòng)的步長:‘）;
lt=input（‘請(qǐng)輸入圓形線圈的中心點(diǎn)縱坐標(biāo)的終值:‘）;
a=a/2;
b=b/2;
y0=la:dl:lt;
k=0;
for?y0=la:dl:lt
?????P_1=dblquad（@integral102*pi-aa）;
?????k=k+1;
?????P_2=dblquad（@integral202*pib-b）;
?????P_3=dblquad（@integral302*pia-a）;
?????P_4=dblquad（@integral402*pi-bb）;?
?????mmm=n1*n2*u0/（4*pi）;
?????M=mmm*（P_1+P_2+P_3+P_4）;?
?????disp（[‘原線圈中心點(diǎn)縱坐標(biāo)在‘?num2str（y0）]）
?????disp（[‘處時(shí)兩線圈的互感M=‘?num2str（M）]）;
?????kk（k）=M;??
?????x（k）=la+dl*k;
end
plot（xkk）
xlabel（‘y0/m‘）
ylabel（‘M/H‘）

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件????????145??2012-10-25?21:08??integral2.m

?????文件????????136??2012-10-25?21:11??integral3.m

?????文件????????126??2012-10-25?21:12??integral4.m

?????文件????????941??2012-10-26?09:28??Coil_Utual_Inductance_test.m

?????文件????????128??2012-10-25?21:07??integral1.m

-----------?---------??----------?-----??----

?????????????????1476????????????????????5