SCILAB&SCICOSについて教えてください

>>246
-->A=-1 ; B=1 ; C=1; D=0;
-->sys=syslin("c",A,B,C,D);
-->bode(sys)
-->ss2tf(sys)
ans =

　　　 1
　　-----
　　1 + s

-->[aa,bb,cc,dd]=abcd(sys)
dd =

0.
cc =

1.
bb =

1.
aa =

- 1.

>>247
もう少し言葉も含めて詳しくおねがいします。

>>248
->A=-1 ; B=1 ; C=1; D=0;
-->sys=syslin("c",A,B,C,D);
-->bode(sys)

-->A=-1 ; B=1 ; C=1; D=0;
-->sys=syslin("c",A,B,C,D);
-->ssprint(sys)
.
x = |-1 |x + | 1 |u

y = | 1 |x

-->bode(sys)

状態空間の場合、evansによるグラフ表示はどうやったらできるかわかる人いますか？

>>251
なんで、まずやってみないかなー？

A=-1 ; B=1 ; C=1; D=0;
sys=syslin("c",A,B,C,D);
ssprint(sys)
evans(sys)

エディタにプログラムを書いて保存するとき、.sciで保存するのと.sceで保存するのとで
何か違いはありますか？自分で実際にやってみた限りでは特に違いがわかりません。

むだ時間要素を扱う方法ってある？

>>254
パーデ近似。

s=%s;
Ts=1/100;
G =(12-6*s*Ts+s*s*Ts*Ts)/(12+6*s*Ts+s*s*Ts*Ts)
sys=syslin("c",G)
bode(sys)

サンプリング周期の1/4くらいまでならOK?

Scilabの起動後にrandで乱数を出すと同じ数が出ます
これを避けて違う乱数を発生させることは出来ますか？

>>256
dt=getdate() ; dt(10)
が毎回変わる処理系ならば（dt(10)はOSのミリセコンド値）
dt=getdate() ; rand("seed" , dt(10) ) っていうのはどうだろう

>>257
これ良いですね、一気に解決しました。
ありがとうございました（・∀・)

>>257
rand("seed", sum(getdate()))

の方が周期が長いかもｗ

>>259
ヒントありがとうございます、自分でも考えてみました。
rand('seed')は8桁のようなのでNも8桁になるようにしてみました。
これだと1時間周期�_秒単位で乱数設定ができると思います。

dt=getdate();
N=60*60*1000*dt(7)+60*1000*dt(8)+1000*dt(9)+dt(10);
rand('seed',N)

>>260
24時間周期の間違いでした
N=0〜89,999,999？

scilab のデータを c# の DLL に渡したり、c# DLL のデータを scilab が受け取ったりできる？

Scilabってdec2bin, bin2decないんだね(´；ω；｀)

>>263
片方向ならある。


-->base2dec("1000",2)
ans =

8.

-->base2dec("1000",10)
ans =

1000.

-->base2dec("1000",16)
ans =

4096.

>>264
片方向でもありがたいです、どうも(｀･ω･´)
でもdec2hexがあるのになんでdec2binが無いんだろうね・・・
来月のV5.0に期待しましょうか

179とかで出てきてるint2dを使って
試しに半径1の半球の体積=2.0944を求めてみようということで、

clear
N=700;

for i=1:N

x0(i)=cos((i-1)*2*%pi/N);
y0(i)=sin((i-1)*2*%pi/N);

x1(i)=cos(i*2*%pi/N);
y1(i)=sin(i*2*%pi/N);

x2(i)=0;
y2(i)=0;

end

X=[x0';x1';x2'];
Y=[y0';y1';y2'];

deff('z2=g(x,y)',['z2=sqrt(1-x.^2-y.^2)'])
int2d(X,Y,g)

というようなソースを書いてみましたが、さらにNを
増やしていって(積分領域を真円に近づけていっても)
結果は2.0918程度となって厳密解と差が出てしまいます。

どなたか、原因がお分かりの方がいらっしゃったら御教授願いたいです。

Scicosで変数を使うことって出来ないですか？
結果のグラフを見ながら何度もゲインKを変更するんですけど
そのKが何箇所かにあって、その全部のブロックを変更するのって
やたらと時間がかかってしまって面倒なんです。
もしご存知の方がいたら教えてください、お願いします。

>>267
ttp://atmori.info/engineer/scilab/scicos/tips.htm

Edit / Context だってさ

パーデ近似でexp(-Ls)を表現するとどうなりますか？
むだ時間が必要です。
お願いします。

一次
-(s*L-2)/(s*L+2)

二次
(s^2*L^2-6*s*L+12)/(s^2*L^2+6*s*L+12)

>>269
>>255
過去スレとは言わないが、せめてこのスレくらい検索してみるくらいの余裕はないのかね？

>>271
ちゃんと見てますよ。
だから、ここで再質問したんです。
255で回答されてますが、この式をどう使っていいかわからないから、質問しました。

>>270
ありがとうございます。
一次と二次で形が違うんですね。
>>255だけでは、わかりませんでした。

Lをむだ時間として、exp(-Ls)の替わりに、これを使えばいいのですね？

>>273
横レスでアレだが、

>>255だけでは、わかりませんでした。

が理解できない。　>>255を実行すればこれが無駄時間を表すことは
一目瞭然だろ。

現に、わからなかったから、質問したのですが？
わかっていたら、質問していませんね。
それがわからない方が、どうかしています。

しかも、一次と二次とで形が違ってるし。

>>276
>しかも、一次と二次とで形が違ってるし。

そんなに自分がアホって自慢しなくても....

>>268
ありがとうございます！！
おかげ様で作業がかなりやり易くなりました。