Scilabを触ってみるｽﾚ

ちょっと質問が正しくなかったので訂正します。

変数に計算結果の行列を入れたいのですが、それに条件を入れたいのです。

たとえば、ｘの範囲を-πからπまで101ポイント取るとして、
y=sin(x)として
ynの結果が0以下であったならynの値を強制的に0とする．．

といった感じにしたいのです。

よろしくお願いします。

y .* (y<0) みたいなので

できますた。

ありがとうございます。

ちなみにWindowsのユーザー名が２バイト文字だと起動しないようですね。

修正ができず、新たに1カウント作っちゃいましたよ。

度々すみません・・・

演習として、ガウシアン分布の３Ｄプロット図を作ろうとしましたが、
ｘとｙをマトリクスで演算させる方法がわかりません．．．。



d=0.5
x=[-2:0.01:2]';
y=[-2:0.01:2]';

z=exp(-2*((x/d)^2+(y/d)^2));

clf()
plot3d(x,y,z,)

よろしくご指南ください m(__)m

追記

要するに、ｘとｙの行列同士のマトリクスの計算結果をｚに入れたいのです。

http://kiwi.emse.fr/SCILAB/sci-bot/x330.htm
のExample 2-1

あとユーザー名２バイト文字にすんなよ

うーん、せっかくですがわかりません．．．

d=0.5
x=[-2:0.1:2];
y=[-2:0.1:2];


z=exp(-2*((x/d)^2+(y/d)^2));

clf()
plot3d(x,y,z,)

for iy=1:length(y)
for ix=1:length(x)
z(iy,ix)=exp(-2*((x(ix)/d)^2+(y(iy)/d)^2));
end
end

plot3d(x,y,z)

>161
ありがとうございます。
参考にさせて頂きます。

一応努力してみました。

clear

d=0.5;

x = linspace(-2,2,101);
y = linspace(-2,2,101);

[X,Y] = ndgrid(x,y);

deff("z=f(x,y)","z=exp(-2*(x.^2/d+y.^2/d ))");

Z=f(X,Y);

plot3d(x,y,Z)

>161

要約すると

clear

d=0.5
x=[-2:0.1:2];
y=[-2:0.1:2];

clf()

for iy=1:length(y)
for ix=1:length(x)
z(iy,ix)=exp(-2*((x(ix)/d)^2+(y(iy)/d)^2));
end
end

plot3d(x,y,z)


こうですね。

申し訳ないのですが
Ｎクイーン問題をＣ言語ではなくScilabを使って解くやり方を教えてください。。。

この３月に発刊されたこの本は、有用でしょうか？
Scilab入門　フリーソフトで始める数値シミュレーション　計測・制御シリーズ　大野修一　CQ出版　2009/03/15　\2,310
http://www.cqpub.co.jp/hanbai/books/40/40991.htm

おお。まさか機械工学板に立ってるとは。

>>166
昔は、電電板とかシミュレーション板がなかったんでなー

>>167
そうだったんすか。古いスレですね
今期から授業で使う事になったんで9月頃から導入しました
先日何となくググったら2000年くらいのスレが見つかって、
そこからこのスレをようやく見つけたところです

>>168
SCILAB&SCICOSについて教えてください
http://science6.2ch.net/test/read.cgi/sim/980480122/

こっちのほうがにぎやかかも

>>169
ありがとうございます。色々見て回ってみます。
まだ数理統計の簡単な計算をすこし始めたところで、
質問や手助けできる事もあまりないですが、大変参考になります。

>>170
統計処理なら

統計解析フリーソフト　Ｒ　【第３章】
http://science6.2ch.net/test/read.cgi/math/1224142396/

こっちもよろしくｗｗ

>>171
ありがとうございます！
さうは一人で持ちきれぬ
ヨイヨイ。

SIVPでAVIファイル作成するときavifileを使用すると
エラーが出てしまうのですが対応方法ないでしょうか。

MatlabからScilabに書き換えてるのですが，以下のエラーの解決策が分かる方がいましたら，教えてください．

<以下エラーの内容>
ii =

2.
!--error 144
与えられたオペランドのための不確定の操作です。

負荷をかけすぎる関数 %s_3_s 定義、または、チェックです。

at line 34 of function bilinear called by :
-d2(ii),d1(ii-1)-d2(ii-1),F1(ii-1),in1,st1,va1)
at line 413 of exec file called by :
exec(%fileToExec); disp(msprintf(gettext("Execution done.\n")));
while executing a callback

どうやら，whileループでfunctionを実行しようとしたら1回目は回るのですが，2回目のループで引っかかるようで，Matlabでは上手く回るので，functionやmainに間違いは無いと思います．

どなたか分かる方がいましたらよろしくお願いします．

>>174
元のプログラムを載せないで、わかるやつがいたとしたら、そいつは悪魔と思わないのか？

過疎ってますね...。質問なんですが、小数点で数字を0.1ではなく0.10のように表示するにはどうすればいいんでしょうか？

>>176
-->printf("%7.2f",0.1)
0.10

>>175
すみません，プログラムの権利の問題上全て載せるわけにはいかないですが，
引っかかってるところを載せます．

abs_dd = abs(dd)

is = 2;

if (ds <= 0) then
is = 1;
end;

sn = (is+is-3);

select (st.ll)

Whileループの中の関数のプログラムの一部で，どうやらifの条件式のところで引っかかってるようです．

>>177
遅れましたが、ありがとうございました。無事解決しました。

>>178
解決したかい？

-->a = sqrt(-1)
a =

i

-->a <= 0
!--error 144
与えられたオペランドのための不確定の操作です。

負荷をかけすぎる関数 %s_3_s 定義、または、チェックです。



-->A = cell(2,2)
A =


!{} {} !
! !
!{} {} !

-->A <= 0
!--error 144
与えられたオペランドのための不確定の操作です。

負荷をかけすぎる関数 %ce_3_s 定義、または、チェックです。


とかあやしいいかも！

A=[1,2,3,4,5,6,7,8,9]
の行列があるとき、三つごとの平均を出すにはどうすればいい？
B1=2
B2=5
B3=8
みたいに出したい。

>>181
エレガントにやってもいいが、
てっとりばやくやるのなら、 for ループで回すほうが速いかな

>>180
すみません遅くなりましたが，無事解決いたしました．
原因は，dsが複素数になっていたらしく，if文の条件式で判定できなかったようです．
real(ds)にすることで無事に解決できました．
ありがとうございました．

Scilab6はいつころ出るかな。。。 予想では1.5年後くらいか。。

http://www.scilab.org/products/scilab/roadmap
Major Scilab 6 Family Objectives

より、期待できるのは、次の改良点か。

・New kernel (interpretor, syntax tree, memory management ...)

Fortran風の言語仕様を排除するらしい。

・Generator of C code
Enable Scilab to generate C code from Scilab script for embedded systems

ScilabからC言語コードを吐き出してくれると、確かにうれしい。

と思ったら、ScilabからC言語の生成は既にあるのか。

http://atoms.scilab.org/toolboxes/scilab2c/
Scilab 2 C
Translate Scilab code into C code

使ったことないから、知らねえ。

>>185
Scilab5.2.0winでScilab2Cを試してみたんだがエラーでCコードは生成できなかった
おまいらもやってみてよ

高階関数って使えます？

f(x) = ax+b
g(x) = gx+b

とあるときに
f(g(x))を関数として出力したいんですけど

できました有難うございました

>>187
horner関数で代入する。

5.2.2アゲ
http://www.scilab.org/products/scilab/download

複数の変数に一度に代入する方法はないだろうか？


関数だったら
-------------------------
function [x,y] = myfunc( a )
x = a*2;
y = a*3;
endfunction

[m, n ]= myfunc(5)
-------------------------
でＯＫだけど、

-------------------------
[x, y] = [3, 8];
-------------------------
はエラーなる。
どうしてだろう？


ベクトルとして、一度に複数の変数に値をセットできると、すっきりするのだが。

えーっと、ス、ス、スシラブ

Ｓｃｉｌａｂは、サイラブ と呼びます。 もしくはサイラボ。

scilabの正当後継は　scilabなのかscicoslabなのか、
どっちなんだ？

誤差の２乗面積が最小になる関数をscilabで数値探索で求めたいのですが
やり方が分かりません。

>>195
> 誤差の２乗面積が最小になる関数をscilabで数値探索で求めたいのですが
> やり方が分かりません。

具体的に、どんな問題？

式を立てられないのはScilabの問題ではないよね。matlabでも同じだろう。

5.3.0のβ1が公開
http://www.scilab.org/products/scilab/download/5.3.0-beta-1

>>198

初期のVer5,0のころに比べ
ver5.3 は動作が速くなった気がする。

Scilab 5.3.0とか、Windows7　Professionl 64bitにいれているんだけど、
64bit版も、32bit版もキーボードの入力を受け付けなくて、
なんにもできない状態なんですが、何が間違っていますか？

>>200

とりあえず、当方 WinXP(32bit)で Scilab5.3 beta2 は、それなりに動いています。

まだ、動作不安定なところがあるけど、Ver5.0のころと比べると、確実に安定かつ速くなっている。

>>200
scicoslab　にしろ。
これで全部解決だ