用matlab仿了一下前半截的电流, 感觉略吓人, 当线圈的电流为sin形,电流峰值1000A,125us上升到顶时,弹丸电流有2*10^10安........ 是不是我传递函数搞错了。。。
代码如下:
R2 = 87*10^(-6);%弹丸电阻
L2 = 7.2*10^(-9);%弹丸电感
M = 0.3;%耦合系数
trans_I2 = tf([ M 0],[L2 R2]);%线圈电流到弹丸电流传递函数
T = 500*10^(-6);%线圈电流周期
f = 1/T;%线圈电流频率
T_seq_div = 5000;%时间周期细分数
T_seq = 0:T/T_seq_div:T/4;%时间序列
I1_seq = 1000*sin(2*pi*f.*T_seq);%线圈电流序列
I2_seq = lsim(trans_I2, I1_seq,T_seq);
plot(T_seq, I1_seq*10000000,'.');%线圈电流, 乘1000000放大好观察
hold on;
plot(T_seq, I2_seq);%弹丸电流