コルピッツ発振回路の発振周波数を求める(2)
2015-05-17
続きです。回路方程式を使って発振周波数を求めます。
wxMaximaを使いますが、まず下準備です。
下図の回路の伝達関数を求めておきます。これぐらいなら手計算で簡単に出来ます。
位相が180度回転する周波数を求めたい回路です。
wxMaximaに以下の式を入力します、コピーペーストでOKです(のはずです)、Enterで処理が始まります。
変数の付け方が回路図では大文字です、長年の習慣です。しかしwxMaximaでは小文字にしています、理由はwxMaximaでは頻繁に打ち込まなくてはいけないので大文字だとシフトキーを押すのが面倒だからです。ただし例外はLです。これを小文字にすると数字の1と区別しにくくなります、例外規定は設けたくなかったのですが、仕方ありません。
-------------この下の行から
s:%i*2*%pi*freq;
r1:100;
r2:10^7;
r3:1;
r4:1/10;
r5:1/10;
c1:10^-7;
c2:10^-7;
L:10^-3;
z1:r1;
z2:r4+1/(s*c1);
z3:r3+s*L;
z4:r2*(1+s*c2*r5)/(1+s*c2*r2+s*c2*r5);
f:z2*z4/(z1*z2+z1*z3+z2*z3+z1*z4+z2*z4);
ratsimp(%);
imagpart(%);
solve(%,freq);
float(%);
-------------この上の行まで
処理結果です(クリックで拡大)
一番最後の行にfreqの計算結果が出ています。
約16MHzと22.5KHzで位相が180度回転することになります。
続く
wxMaximaを使いますが、まず下準備です。
下図の回路の伝達関数を求めておきます。これぐらいなら手計算で簡単に出来ます。
位相が180度回転する周波数を求めたい回路です。
wxMaximaに以下の式を入力します、コピーペーストでOKです(のはずです)、Enterで処理が始まります。
変数の付け方が回路図では大文字です、長年の習慣です。しかしwxMaximaでは小文字にしています、理由はwxMaximaでは頻繁に打ち込まなくてはいけないので大文字だとシフトキーを押すのが面倒だからです。ただし例外はLです。これを小文字にすると数字の1と区別しにくくなります、例外規定は設けたくなかったのですが、仕方ありません。
-------------この下の行から
s:%i*2*%pi*freq;
r1:100;
r2:10^7;
r3:1;
r4:1/10;
r5:1/10;
c1:10^-7;
c2:10^-7;
L:10^-3;
z1:r1;
z2:r4+1/(s*c1);
z3:r3+s*L;
z4:r2*(1+s*c2*r5)/(1+s*c2*r2+s*c2*r5);
f:z2*z4/(z1*z2+z1*z3+z2*z3+z1*z4+z2*z4);
ratsimp(%);
imagpart(%);
solve(%,freq);
float(%);
-------------この上の行まで
処理結果です(クリックで拡大)
一番最後の行にfreqの計算結果が出ています。
約16MHzと22.5KHzで位相が180度回転することになります。
続く
- 関連記事
スポンサーサイト
コメント
edyさん、お早うございます。
コルピッツ発振回路の解析、随分と極まった感じですね。ご苦労様でした。この路線を進めて行かれる中で、発振条件のもう一つの、増幅度(A)の解析は、如何でしょうか? 果たして、Maximaが爆発的な式を吐き出すのでしょうか?(笑)
コルピッツ発振回路の解析、随分と極まった感じですね。ご苦労様でした。この路線を進めて行かれる中で、発振条件のもう一つの、増幅度(A)の解析は、如何でしょうか? 果たして、Maximaが爆発的な式を吐き出すのでしょうか?(笑)
コメントの投稿
トラックバック
この記事へのトラックバックURL
https://edycube.blog.fc2.com/tb.php/780-ef6ed45d
コルピッツ発振回路の解析、随分と極まった感じですね。ご苦労様でした。この路線を進めて行かれる中で、発振条件のもう一つの、増幅度(A)の解析は、如何でしょうか? 果たして、Maximaが爆発的な式を吐き出すのでしょうか?(笑)