1

そこで、下の図に示すシステムをモデル化しようとしています。以前、Matlab を使用して ODE45 でモデル化したことがありますが、何らかの理由で、次のコードで奇妙な出力が得られます。逆バイアスされたときにダイオードの間違った方向に電流が流れるのを防ぐのに役立つ if( (Vin - Vd - Vc)/Ll >= 0 ) の条件を入れたためかどうかはわかりません。

写真を添付し​​ようとしましたが、私の担当者は十分に高くないので、説明する必要があります.

source -> Inductor(Ll) -> diode -> Capacitor -> ground
                               |->Inductor(Lline)->resistor(Rf) -> ground

どんな助けでも大歓迎です！

function dy = ThreeLamp(t,y)

I1 = y(1);
I2 = y(2);
Vc = y(3);

Is = 1E-15;                                     % these are parameters for the diode
n = 1.05;
Vt = .025;

D = 24;                                        % this number is from 0 to 50
f = 100E3;                                     % this is the frequency of the input waveform
Ll = 5E-7;
Lline = 5.3E-6;
C = 10E-6;
Rf = 10;
Vin = (465/2)*(6/266)*square(2*pi*f*t,D);       % Input voltage waveform

if( Vin < 0 )                                   
   Vin = 0;                                     % get rid of negative half cycle
end

%%

Vd = n*Vt*log(I1/Is + 1)/log(10);               % This is the voltage across the diode

if( (Vin - Vd - Vc)/Ll >= 0 )                   % Is the diode conducting when forward biased?
    dy(1) = (Vin - Vd - Vc)/Ll;                 % Sounds good, what is the current?
else
    dy(1) = 0;                                  % force current to 0 when reversed biased
end
dy(2) = (Vc - I2*Rf)/Lline;                     % KCL for second loop
dy(3) = (I1 - I2)/C;                            % Voltage on cap

dy = dy';                                       % ODE45 likes output to be in columns

end

% [t, y] = ode45('ThreeLamp',[0 .001],[0 0 0 ]);
4

1 に答える 1

1

最近、DCDC ブースト コンバーターをモデル化しようとして、まったく同じ問題が発生しました。am304 が言うように、スティッフ ソルバーがこの問題の解決策になる可能性があります。ダイオードがオン/オフするときの不連続性は非常に硬く、可変ステップ ソルバーで抽出するのは非常に困難です。私のアプリケーションでは、固定ステップ ソルバーがうまく機能することがわかりました。

編集

ただし、horchler が指摘しているように、許容誤差の収束やその他の数値的なセーフティ ネットは存在しないため、出力を非常に慎重に解釈する必要があります。非常に小さなステップ サイズを使用し、出力を非常に慎重に扱うことで、出力が予想どおりの形状になると、多くの場合、この問題が軽減されます。しかし、スティッフなシステムを解こうとするときは、horchler の警告に注意する必要があります。固定ステップ法がうまくいかなかったとしても、ソリューションに与えられた洞察は努力する価値があることがわかりました.

編集終了

また、ダイオードを流れる電流の方向も確認する必要があります。ダイオードに大きな逆バイアス電流が流れることがわかりましたが、不可能です!

ここに回路モデルを設定します。これは写真です。これは正しいトポロジですか? チッ！写真

編集

さて、回路が理解できたので、微分関数ThreeLamp. ダイオードが逆バイアスされている場合、インダクタを流れる電流はダイオードを流れることができません。したがって、ダイオードの導通/非導通のポイントでの電流の変化は、インダクタを流れる電流を遮断すると同時にダイオードを流れる電流を遮断するために大きくなります。この問題は、インダクタ電流の計算で 25 行目あたりが好きだと思います。

于 2013-06-24T14:36:45.003 に答える