基于MATLAB计算两平行线圈的互感

本代码计算了平行的方形线圈和圆形线圈在任意位置的互感，用户可以设定线圈参数。代码开发环境为MATLAB软件

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