clc;
clear;
close all;
global a b r h n1 n2 x0 y0
u0=4*pi*10^（-7）;
a=input（‘请输入方形线圈的长:‘）;
b=input（‘请输入方形线圈的宽:‘）;
r=input（‘请输入圆形线圈的半径:‘）;
h=input（‘请输入两线圈的竖直高度:‘）;
n1=input（‘请输入方形线圈的匝数:‘）;
n2=input（‘请输入圆形线圈的匝数:‘）;
x0=input（‘请输入圆形线圈中心点的横坐标:‘）;
la=input（‘请输入圆形线圈中心点纵坐标的初始值:‘）;
dl=input（‘请输入纵坐标移动的步长:‘）;
lt=input（‘请输入圆形线圈的中心点纵坐标的终值:‘）;
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（[‘原线圈中心点纵坐标在‘ num2str（y0）]）
disp（[‘处时两线圈的互感M=‘ num2str（M）]）;
kk（k）=M;
x（k）=la+dl*k;
end
plot（xkk）
xlabel（‘y0/m‘）
ylabel（‘M/H‘）属性 大小 日期 时间 名称
----------- --------- ---------- ----- ----
文件 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
