x_force_inductance_relation.m

变压器仿真

%计算电感系数、电磁力与位移x之间的关系曲线
clc
clear
% x和L的数据
xdata=[0,0.2,0.4,0.6,0.8,1.0,1.2,1.4,1.6,1.8,2.0];
Ldata=[2.8,2.26,1.78,1.52,1.34,1.26,1.20,1.16,1.13,1.11,1.10];
%换算成国际单位制
x=xdata*1.e-2;
L=Ldata*1.e-3;
len=length(x)
xmax=x(len);
%曲线拟合
a=polyfit(x,L,4);
%拟合检查
for n=1:101
    xfit(n)=xmax*(n-1)/100;
    Lfit(n)=a(1)*xfit(n)^4+a(2)*xfit(n)^3+a(3)*xfit(n)^2+a(4)*xfit(n)+a(5);
end
%绘图
plot(xdata,Ldata,'*')
hold
plot(xfit*100,Lfit*1000)
hold
xlabel('x[cm]')
ylabel('L[mH]')
fprintf('\n Paused.Hit any key to plot the force.\n')
pause;
% 电磁力计算与绘图
%电流数值
I=0.75;
for n=1:101
    xfit(n)=0.002+0.016*(n-1)/100;
    F(n)=4*a(1)*xfit(n)^3+3*a(2)*xfit(n)^2+2*a(3)*xfit(n)+a(4);
    F(n)=(I^2/2)*F(n);
end
plot(xfit*100,F)
xlabel('x[cm]')
ylabel('Force [N]')