Scilabを触ってみるｽﾚ

108 ：bo!：2006/05/26(金) 13:33:40 ID:1QGwhKJy

scilabのプログラムを考えてほしいんです。お願いします！！
「100000までの素数をすべて求めよ。」
誰かほんとにお願いします＞＜。

109 ：名無しさん＠３周年：2006/05/28(日) 22:57:11 ID:gNjOxOR3

>>108 効率悪いけど。まぁループ終端をsqrt(N)にするとか、適当に直して。

lines(0)
N=100;

printf("N=%d\n\n",N)
p(1:N)=1;
p(1)=0;
for i=2:N
for j=i+i:i:N
p(j)=0;
end
end
for i=1:N
if p(i)<>0 then
printf("%d\n",i)

end
end

disp "finish"

110 ：名無しさん＠３周年：2006/05/28(日) 22:57:39 ID:gNjOxOR3

-->N=100

2
3
5
7
11
13
17
19
23
29
31
37
41
43
47
53
59
61
67
71
73
79
83
89
97

finish

111 ：名無しさん＠３周年：2006/05/29(月) 13:07:31 ID:7GVgeN7V

すみませんがこの問題教えてください。。。。

x,yを1000までの自然数として
y^3=x^2+4 をみたすx,yの組をもとめよ

今日の１７時が提出期限なんです↓

112 ：プジョル：2006/05/29(月) 13:14:37 ID:7GVgeN7V

申し訳ないのですが、この問題をなるべく早く解いてもらえませんか？今日の5時までが期限なんです。お願いします。
自然数nに対し、
n=a^2+b^2+c^2＋d^2
を満たす整数a,b,c,dを求めるプログラムを作れ。
です。お願いします

113 ：名無しさん＠３周年：2006/05/29(月) 13:52:20 ID:g7vqhy2y

ほいっ

114 ：名無しさん＠３周年：2006/09/04(月) 14:27:53 ID:9TlKPD9D

MacOS10.4 でScilabやってマス
playsnd()等音声出力ができないんですけど、何かほかに設定いりますか？
エラーメッセージは
>Audio capabilities not available
です。

115 ：114：2006/09/04(月) 14:32:57 ID:9TlKPD9D

↑
Scilab 3.11 と 4.0でやってみますた
上のエラーメッセージは 3.11のもので 4.0では
Writing Wave file: Microsoft PCM format, 1 channel, 22050 samp/sec
44100 byte/sec, 2 block align, 16 bits/samp
Finished writing Wave file, 132302 data bytes
と正常らしいですが、音はでませんです。

116 ：名無しさん＠３周年：2006/09/04(月) 23:30:28 ID:fBT/fxRq

playsndのsciファイルを見るとPlaySoundという関数?で実際の再生を行なっているので
こいつをscilabのunix関数とMacOSXのopen関数の組み合わせにすれば
音を鳴らすことはできると思う。

117 ：114：2006/09/05(火) 22:02:34 ID:JROw1EaZ

即レス感謝です。
if MSDOS then
savewave(TMPDIR+'/_playsnd_.wav',y,rate);
PlaySound(TMPDIR+'/_playsnd_.wav');
return
else
if aplay<>'/dev/audio' then
// should send aplay to PlaySound
savewave(TMPDIR+'/_playsnd_.wav',y,rate);
PlaySound(TMPDIR+'/_playsnd_.wav',aplay);
これですね。
とりあえずTMPファイルをsavewave で作成してからそれをPlaySound関数で読み出しているっていう手順みたいです。
エラーメッセージからはとりあえずwavefileは作成できてるみたいですねー。
>こいつをscilabのunix関数とMacOSXのopen関数の組み合わせにすれば
実際どうすればいいか分かりません、スンマセン。octaveでは試してないですけど音出るみたいなんで、インスコしてoctaveのも見てみます〜。

118 ：名無しさん＠３周年：2006/09/05(火) 23:28:17 ID:n7WoAD8u

PlaySoundの行をを以下に書き換えたらできないかな?
（playsnd.sciからplaysnd.binを作り直す必要があるかも。
.sciファイルと.binファイルの関係は忘れました。）

unix('open '+TMPDIR+'/_playsnd_.wav');

ただし、デフォルトではiTunesが立ち上がると思いますので
ちょっと不便です。

Macのコマンドラインから音を再生するプログラム
（例えばplaysoundってのがあります）がインストール
してあれば

unix('playsound '+TMPDIR+'/_playsnd_.wav');

のように書き換えれば（playsoundというプログラムが
PATHの通っているディレクトリに入ってる必要が
あります）できると思います。

119 ：名無し：2006/11/11(土) 04:11:29 ID:dpcozmDI

tan(x)=3をニュートン法で解くって、どうやったらいいですか？？

120 ：名無しさん＠３周年：2006/11/11(土) 09:56:25 ID:v4O4lKQR

>>119
1)
atan(3)


2)
clear fx
function y=fx(x)
y=tan(x)-3
endfunction

fsolve(0,fx)

121 ：名無し：2006/11/11(土) 11:54:09 ID:3TMpmDdE

早速のお返事、ありがとうございますっ

122 ：名無しさん＠３周年：2006/11/23(木) 12:08:33 ID:/va2RJ8W

scilab 4.0でint16の挙動がWindows版とLinux版(どちらもバイナリ版をインス
トールしました)で異なっているんですが、どちらが正しいんでしょうか。

Linux版
------------------------------------------------------------------
-->int16([32767, 32768, 32769])
ans =

! 32767 -32768 -32768 !
------------------------------------------------------------------

Windows版
------------------------------------------------------------------
-->int16([32767, 32768, 32769])
ans =

! 32767 -32768 -32767 !
------------------------------------------------------------------

123 ：名無しさん＠３周年：2006/11/23(木) 12:33:44 ID:HXGqiJ6s

>>122
32769っていう時点で想定外なので、どーなっても構わないということになっているので
気にするな、っていうかそんな数入れるな。

124 ：名無しさん＠３周年：2006/11/23(木) 18:56:27 ID:/va2RJ8W

>>123
範囲外の値に対する動作は未定義ってことですね。

「Scilabで学ぶディジタル信号処理」という書籍で、sign bitが立つ値
は負値に変換されるという解説があり、Linux版のバグなのかと思ったもので。

125 ：GCofSS：2006/12/07(木) 22:57:11 ID:kmVMOtUV

Scilabを画像処理に使っているのですが、以前Matlabで画像処理をやっていて急遽freeの
Scilabに書き換えなきゃならなくなって書き換えを行っているのですが、Matlabで使用している
imagescに替わるScilabの関数を教えてください！

126 ：名無しさん＠３周年：2006/12/08(金) 12:31:19 ID:MlXhpj8S

>>125
imagescって、ナニやるものですか？
SIP (Scilab Image Processing)とか？

127 ：名無しさん＠３周年：2006/12/08(金) 13:50:08 ID:6AyRveCp

>>126
そうですね。imagescとはデータをスケーリングし Image を表示するものです。
SIPにあるいろいろな関数を試してみたのですが、なかなかうまくいかなくて・・・。

128 ：名無しさん＠３周年：2006/12/23(土) 11:26:59 ID:9J/FO075

4.1age

129 ：名無しさん＠３周年：2007/01/01(月) 20:27:10 ID:LWLGrPE6

SCILAB-4.1をダウンロードしたのですが英語だらけで分かりません(>_<)すいませんが使い方を教えていただけないでしょうか？

130 ：名無しさん＠３周年：2007/01/02(火) 18:13:41 ID:G1wRsfB6

>>129
scilab の検索結果のうち 日本語のページ 約 81,200 件中 1 - 50 件目 (0.24 秒)
http://www.ecl.sys.hiroshima-u.ac.jp/scilab/introscilab/introscilab.html

131 ：名無しさん＠３周年：2007/01/06(土) 00:24:26 ID:4J7LQeFf

根号の部分が小数表示されるのですが、根号のまま表示させることは可能でしょうか？

132 ：名無しさん＠３周年：2007/01/06(土) 00:50:27 ID:714NIKRq

かたつむり

133 ：名無しさん＠３周年：2007/01/06(土) 10:33:38 ID:tAuRVBI1

>>131
MAXIMAじゃないんだから、計算しちゃうんじゃない？

134 ：名無しさん＠３周年：2007/01/06(土) 11:44:46 ID:M/ef4eID

-->a=0　と入力すると
a =

0.　と自動的に表示されますが

これを表示させない方法はありますか？つまり
-->a　と入力しない限り
a =

0.　が表示しないようにしたいのですが。

よろしくお願いしま。

135 ：名無しさん＠３周年：2007/01/06(土) 11:45:44 ID:M/ef4eID

す

136 ：名無しさん＠３周年：2007/01/06(土) 13:12:59 ID:lormVx9K

>>134
-->a=0;
ではダメなの?

137 ：名無しさん＠３周年：2007/01/06(土) 13:45:11 ID:M/ef4eID

>>136
ありがとうございました、それで出来ました（・∀・）

138 ：名無しさん＠３周年：2007/01/11(木) 19:27:30 ID:YWhcPYHW

scilabでクラスは作ることはできますか？
ヘルプで検索してもJavaｲﾝﾀﾌｪｰｽとmarkovなんちゃらのことしかでてきません・・・

139 ：名無しさん＠３周年：2007/01/11(木) 19:31:15 ID:YWhcPYHW

>>138
tlistというのがあることがわかりました。
もうちょっと自分で調べてみます。

140 ：名無しさん＠３周年：2007/03/11(日) 17:59:19 ID:l8+Q/qiG

ScilabはCやFortranとの連携は出来るようですがC++はどうでしょうか？

141 ：名無しさん＠３周年：2007/12/20(木) 01:23:30 ID:tn69q6wz

ちょいメモ：
グラフのＸＹ軸のラベル名を変更する場合はtlistを使用する
tlistでないと既定されているmatrixサイズが変更できない

a=gca();
a.x_ticks=tlist(['ticks' 'locations' 'labels'],[1:5],['o' 'b' 'a' 'k' 'a']);

142 ：名無しさん＠３周年：2008/01/15(火) 01:40:12 ID:/68kNl3s

これを解いていただけませんか？

(1)漸近安定性、可制御性、可観測性を調べる。
(2)極を(0.31±0.85*i)に配置するための、状態フィードバックゲインを計算。
(3)状態フィードバックの有無による状態変数の時間変化をシミュレーション。

◆状態方程式
x[i+1] = A*x[i] + b*u[i]
y[i] = c*x[i]
x[0]=[1;2]
A = [3, 8; 1, 5]
b = [1;2]
c = [1,0]

※不可制御となってしまう場合、A行列の数値の並びを変えてよい。

143 ：名無しさん＠３周年：2008/03/20(木) 20:04:30 ID:8ykkfMqd

保守

144 ：名無しさん＠３周年：2008/04/03(木) 23:54:52 ID:FSOW6Acv

保守

145 ：名無しさん＠３周年：2008/04/17(木) 00:23:13 ID:cm9SwHDQ

すいません。。自力で調べたのですがいまだに解決できないので、教えてください。

SIVPを用いてグレイスケール化された画像をg_imとしたとき、
位置P(i,j)とg_im(i,j)の値の関係を３次元プロットしたいと考えています。

plot3d関数を使おうと思ったのですが、できませんでした。

146 ：名無しさん＠３周年：2008/05/02(金) 21:52:54 ID:9/SQ+8wz

どうしてもjavaから呼べないんですけど、、教えてクンですみません。教えてください。

147 ：名無しさん＠３周年：2008/05/29(木) 23:46:44 ID:29CjSwco

>>142
1.rankを使って可制御行列、可観測行列のランクを調べる
2.ppolを使う
3.scicos使えば簡単

148 ：名無しさん＠３周年：2008/06/18(水) 03:39:36 ID:yC9nfMi+

寂れてますね・・・
みんなoctaveに引っ越したのかな

149 ：名無しさん＠３周年：2008/06/20(金) 00:09:55 ID:+OyVGt3V

はじめまして
突然質問でごめんなさい

うちの環境だと
csimで1000ポイント以上計算させると誤差が大きくなるんですけど
ポイント数が多くても精度を保って計算させる方法とかないですか？

150 ：名無しさん＠３周年：2008/06/20(金) 10:42:47 ID:kgx2Gyj2

>>149
csimのパラメータに、誤差を指定するパラメータ : tol あるんだけど、

y=csim(u,t,sl,x0,tol)

試してみた？

151 ：名無しさん＠３周年：2008/06/21(土) 21:58:27 ID:OBAHoZyz

>>150
どうもありがとうございます。

いろいろいじってみましたけど変わりません。


精度が悪くなるというかポイント数を増やすと
最終的に計算結果がめちゃくちゃになるんですよね。

152 ：昨日からはじめた人：2009/01/11(日) 22:32:52 ID:eCuJyVOV

エディターで、

x=-π〜π
y=sin(x)

で、もしｘ>0 ならyを1にしたい場合、どういった記述に
するのが妥当なのでしょうか？

よろしくお願いします。

153 ：昨日からはじめた人：2009/01/12(月) 00:14:39 ID:+ATel+I1

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

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

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

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

よろしくお願いします。

154 ：名無しさん＠３周年：2009/01/12(月) 00:17:14 ID:YmhqvdsU

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

155 ：昨日からはじめた人：2009/01/12(月) 00:47:17 ID:+ATel+I1

できますた。

ありがとうございます。

156 ：155：2009/01/12(月) 17:34:35 ID:+ATel+I1

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

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

157 ：昨日からはじめた人：2009/01/12(月) 19:53:58 ID:+ATel+I1

度々すみません・・・

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



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

158 ：昨日からはじめた人：2009/01/12(月) 20:03:44 ID:+ATel+I1

追記

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

159 ：名無しさん＠３周年：2009/01/12(月) 23:26:20 ID:YmhqvdsU

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

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

160 ：昨日からはじめた人：2009/01/13(火) 00:37:49 ID:Gh1DJUVm

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

161 ：名無しさん＠３周年：2009/01/13(火) 20:03:49 ID:p85x2yXD

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)

162 ：名無しさん＠３周年：2009/01/13(火) 20:32:39 ID:Gh1DJUVm

>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)

163 ：名無しさん＠３周年：2009/01/13(火) 20:35:16 ID:Gh1DJUVm

>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)


こうですね。

164 ：名無しさん＠３周年：2009/06/07(日) 18:15:58 ID:DOrg7rdg

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

165 ：名無しさん＠３周年：2009/08/04(火) 22:26:22 ID:ulGsLK4m

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