!SET bitmap SIZE 600,600 !画面を大きくする
SET WINDOW -0.5,6.5, -50,5 !表示領域
DRAW grid(1,5) !左端の目盛り
FOR f=1 TO 6 !x軸が対数
PLOT TEXT ,AT f-0.3,-0.15: mid$("10 100 1k 10k 100k1M ",4*(f-1)+1,4)
NEXT f
FOR f=10 TO 100000 STEP 100 !周波数[Hz]
CALL routine
LET vv=vo(1)/vi(1)
PLOT LINES: LOG10(f),20*LOG10(ABS(vv)); !利得[dB]
NEXT f
PLOT LINES
SET TEXT COLOR 2
FOR i=5 TO -45 STEP -5 !右端の縦軸目盛り
PLOT TEXT ,AT 6,i: STR$(i*2)&"°"
NEXT i
SET LINE COLOR 2
FOR f=10 TO 100000 STEP 100 !周波数[Hz]
CALL routine
LET vv=vo(1)/vi(1)
PLOT LINES: LOG10(f),DEG(ATN(Im(vv)/Re(vv)))/2; !位相θ ※1/2倍
NEXT f
PLOT LINES
END
周波数プログラム
投稿日:2008年11月27日(木)15時48分15秒周波数のプログラムを有難うございました。
サブルーチンにつなげたら直ぐ軽快に動きました。
(1)最後のDB表示の#3と位相表示の#2のグラフ
を一緒に表示したいのですが、
今のままで、「!」マークを外し、そのまま
つなげても動きませんでした。
DBの縦目盛りは左側に、位相の0-90度目盛りは
右側にといった具合です。
こうすると、本にある「ボード線図」ずばりに
なります。
(2)位相表示は、なにか周波数の文字の上が表示されて
ません。
私の バージョンは、VER 5.02 ですが、
最新はVER 7.2.8 のようです。
これが原因でしょうか。
私は、ハードデスクの「D」に本「ブルーバックス」
に添付のCD ROM を指示に従い展開したのですが、
インターネット上の「(仮称)10進BASIC」にある
VER 7.2.8 のを、そのまま同じ「フォルダー」に
コピー「展開?}したら駄目でしょうか。
前のVER 5.02 は、この本のサンプルも一緒に
入ってるのでどこまでがVER 5.02 本体なのか
分らずVER 5.02 だけ消すのは困難なのです。
VER 5.02をVER 7.2.8 で上書きするのはどう
するのですか。
(3)突然ですが、同じ行数と列数のMATには、
固有値成るものが、あるそうですが、
この