tas-net-2の日記

自動車整備、修理対応

ExcelデータをScilabで読みこむ:例G-scan>CSV:変換:>Excel>Scilab:wav file read:FFT:

Scilab

ExcelデータをScilabで読みこむ:例G-scan>CSV:変換:>Excel>Scilab:wav file read:FFT: - CQF - Yahoo!ブログ

 

イメージ 1

ExcelデータをScilabで読みこむ例G-scan>CSV:変換:>Excel>Scilab:wav file read:FFT:

スキャンデータの多くが基本的にはCSVなどのテキストにbinデータから変換されて、Excelなどの
処理ソフトに読み込まれる。
例:
G-scanもCSVから読み込む。
Excelで読み込んだデータは以下の手順で Scilabに移してみる。そこでの、メリット（自分のスキルのっ内で使えるのか、多少スキルアップも期待して）があるのかも
その過程で見極めたい。


今確認できているのは、*.xls形式で保存したもの。　これ以降の拡張子ではまだ対応してないようだが
未確認。

テキストで読み込む方法も試しておきたい。


-->[fd,SST,Sheetnames,Sheetpos] = xls_open('wfy11-047-04.xls')
警告: ファイル 'C:\DOCUME~1\user\LOCALS~1\Temp\SCI_TMP_3832_\wfy11-047-04.xls\Workbook' は既に Scilabで開かれています。
Sheetpos =

1652. 98379. 98642. 
Sheetnames =

!Sheet1 Sheet2 Sheet3 !
SST =


[]
fd =

2. 




-->[Value,TextInd] = xls_read(fd,Sheetpos(1));

-->Value(:,1)
ans =

0.09 
0.78 
0.73 
0.39 
0.05 
0.14
.
.
.
-->data1=Value(:,1)
data1 =

0.09 
0.78 
0.73 
0.39 
0.05 

.
.
.

-->plot(data1)

で何とか数値データが取り込めた。

ここから、何らかの分析とか評価が出来るものを摸索する。




http://tokida.sakura.ne.jp/Hobbies/computer/Scilab.htm

ベクトル・行列の一部を呼び出す 
関数を定義する 
関数を呼び出す 
グローバル変数とローカル変数 
関数への引数の与え方 
プログラミング 
スクリプトの実行 
変数をファイルに保存する 
参考図書 
役に立つページ


http://www.bekkoame.ne.jp/~ponpoko/index.html

Scilabを中心としたMATLABクローン入門 
SCICOSで遊ぼう


http://www.ecl.hiroshima-u.ac.jp/~ohno/scilab/introscilab/node7.html

6.1 二次元プロット 
6.1.1 LineSpec 
6.1.2 タイトルと凡例 
6.1.3 GUIグラフィックエディタ 
6.1.4 図の出力




:wav file read:FFT:[20160609]

scilab でオーディオ信号を扱ってみる

http://blogs.yahoo.co.jp/susanoo2001_hero/6652650.html


wav sample files
http://www.wavsource.com/sfx/sfx.htm

http://www.wavsource.com/snds_2016-06-05_6174362464058618/sfx/chime_big_ben.wav

chime_big_ben.wav > sound.wav

現在の作業ディレクトリを示すコマンドは、present working directoryの略で
pwd

-->pwd

 

 コメント（91）

  

顔アイコン

scilab：scilab の素晴らしさを実感
http://blogs.yahoo.co.jp/puccinisan/16494435.html

モータートルクと負荷トルク
http://blogs.yahoo.co.jp/puccinisan/folder/598353.html?m=lc&p=2削除

2010/11/26(金) 午後 10:41[ ogw*og*3 ]返信する

  

顔アイコン

http://user.numazu-ct.ac.jp/~mochizuki-k/scilab/index.html

http://atmori.info/engineer/scilab/scicos/more01.htm
ちょっとだけ応用
マルチプレクサ
セレクター
SciFunc
ReadFromInputFile削除

2010/11/26(金) 午後 10:49[ ogw*og*3 ]返信する

  

顔アイコン

６４点高速フーリエ変換回路 設計仕様書 (Ver 1.0)
http://www.ie.u-ryukyu.ac.jp/~wada/design07/spec_j.html

http://dsp.cse.kyutech.ac.jp/lec/h22_dsp/No11.pdf

Scilabの魅力を伝えるべく、「Scilab超入門」と題した連載を続けようと思います。工学系の学生や研究者の方に興味を持っていただけるきっかけにでもなれば幸いです。
Scilab超入門 第1回 Scilabのインストールと簡単な使い方（この記事）
Scilab超入門 第2回 行列とベクトルの基本
Scilab超入門 第3回 行列式/逆行列/固有値/固有ベクトルを求める関数
Scilab超入門 第4回 ベクトルの大きさ/正規化/連結/等間隔のベクトル
Scilab超入門 第5回 plot関数でグラフを描く
http://networkprogramming.blog18.fc2.com/blog-entry-45.html削除

2010/11/27(土) 午前 0:40[ ogw*og*3 ]返信する

  

顔アイコン

[scilab] fft()関数とフーリエ級数
2008年7月13日 18:38 滑稽者 | 個別ページ | コメント(0) | トラックバック(0) 
scilabのfft()関数による、フーリエ変換と、そのフーリエ級数が収められる配列についての素朴な疑問です。

試験にもでたので自分なりにまとめて見ました。 
タグ: プログラミング 
[scilab] fft()関数とフーリエ級数の続きを読む
http://networkprogramming.blog18.fc2.com/blog-entry-45.html
http://blog.delta50.net/2008/07/scilab-fft.html削除

2010/11/27(土) 午前 0:42[ ogw*og*3 ]返信する

  

顔アイコン

3章 どのようにして成分を知るのか？ 2章にて総ての波形は位相や振幅の異なる各周期成分の寄せ集めである事を述べたが、フーリエ変換の目的は総ての波形から各周期成分の振幅を知る事にある。それには、関数の直交性という関係を利用する。
例えば、COS波とSIN波を互いに掛けて合成した波形の面積は0になる。下図ではSIN波に各周期のSIN波やCOS波を掛けているが、同じ周期のSIN波同士の積以外は面積が0になる事が判る。
つまり、合成波のan*cos nxの成分を知りたければ、合成波にcos nxを掛けて足せばan*cos nxの成分だけ残る。
https://sites.google.com/site/yoshinoritakeuchi/Home/dejitaru-shingou-shori/furie-henkan/3shou--dono-you-ni-shi-te-seibun-wo-shiru-no-ka削除

2010/11/27(土) 午前 0:57[ ogw*og*3 ]返信する

  

顔アイコン

(1)フーリエ変換
1章 フーリエ変換とは
2章 総ての波形はSin波とCos波の組み合わせで再現できる
3章 どのようにして成分を知るのか？
4章 フーリエ級数から成分表示へ
5章 エクセルのシート上でフーリエ変換を手動で計算してみる
6章 フーリエ変換の問題点・良い点・仕事での使い方など
7章 リーケージ誤差
(2)デジタル信号処理
1章 簡単なフィルター(移動平均・移動差分)
2章 フィルターの分類
(3)Scilabによる信号解析
1章 Scilabの用法
(1)初歩の使い方
(2)変数と配列
(3)配列の計算
(4)配列どうしの計算
(5)配列の連結
(6)配列の一部利用
(7)配列情報
2章 Scilabによるフーリエ変換
3章Scilabによるフィルター処理
(1)ScilabによるIIRフィルター
(2)ScilabによるFIRフィルター
(3)ScilabによるFFTフィルタ
4章伝達関数
https://sites.google.com/site/yoshinoritakeuchi/Home/dejitaru-sh削除

2010/11/27(土) 午前 1:05[ ogw*og*3 ]返信する

  

顔アイコン

A semi-colon at the end of a line will stop Scilab output for that line. This is
useful if you are working with large vectors or matrices. Try
y = A*x ;
Note that the answer is not printed out. To see y just type
y

http://www.blogdopapeleiro.com.br/biblioteca/AplicativosLivres/Scilab/tutorial-all.pdf

p8削除

2013/11/2(土) 午後 6:31ogw*og*2返信する

  

顔アイコン

x = 0 : 0.1 : 10;
y = sin(x);
z = cos(x);
plot(x, y, x, z);

http://spiff.rit.edu/classes/phys317/lectures/plot/plot.html削除

2013/11/2(土) 午後 11:56ogw*og*2返信する

  

顔アイコン

Simple Plotting in Scilab - Michael Richmond's computer
spiff.rit.edu/classes/phys317/lectures/plot/plot.html

Mar 6, 2007 - You can use Scilab not only to calculate quantites, but to plot them as well. In this course, we'll do only simple x-y plots. It's easiest to explain ...削除

2013/11/2(土) 午後 11:59ogw*og*2返信する

  

顔アイコン

Scilab - Community - Google+
https://plus.google.com/communities/111361848707994482299

The training course of the Openeering Scilab Professional experts is now .... A new Openeering tutorial for beginners is available: learn how to use Scilab as a ...削除

2013/11/3(日) 午前 0:09ogw*og*2返信する

  

顔アイコン

The Scilab MySQL toolbox. - Google Code
code.google.com/p/scilab-mip/wiki/MySqlToolbox
scilab // In Scilab, in the mysql directory, build the toolbox (once). exec builder. sce; // Load the toolbox (each time you need to use the mysql toolbox). exec ...削除

2013/11/3(日) 午前 0:12ogw*og*2返信する

  

顔アイコン

[PDF]
Scilab - GPIB toolbox installation - Reveyrand
www.reveyrand.com/products/GPIB/GPIB_visa.pdf

statement 'Scilab'.??? VISA_read – Read text from selected instrument trough VISA. Character 0 ends reading of the text. Note: It is not possible to use this ...

http://www.reveyrand.com/products/GPIB/GPIB_visa.pdf削除

2013/11/3(日) 午前 0:22ogw*og*2返信する

  

顔アイコン

Index of /RTAI

https://www.rtai.org/RTAI/削除

2013/11/3(日) 午前 5:17ogw*og*2返信する

  

顔アイコン

// Define your x-values
x = 0 : .1 : 10*%pi;
// Define your function
y = x .* cos(x);
// Plot as used in Matlab
plot(x, y)
// Add title, labels and legend
title('x * cos(x)')
xlabel('x'); ylabel('y');
legend('y = x * cos(x)');


http://www.matrixlab-examples.com/scilab-examples.html削除

2013/11/3(日) 午前 5:26ogw*og*2返信する

  

顔アイコン

Polar Plots

polarplot - Plot polar coordinates
polarplot(theta, rho, <opt_args>) 

// Define your variable
t = linspace(0, 2*%pi, 100);
// Define your polar function
r = sqrt(abs(2 * cos(5*t)));
// Plot in polar coordinates
polarplot(t, r);削除

2013/11/3(日) 午前 5:31ogw*og*2返信する

  

顔アイコン

Histograms in Scilab 

histplot - plot a histogram
histplot(bins, data, <opt_args>)


// Define your data or get it somehow
y = 10 * rand(500, 1);
// Create the histogram. The 5 at the beginning means 5 bins.
// The 2 at the end means lines in blue.
histplot(5, y, 2)削除

2013/11/3(日) 午前 5:40ogw*og*2返信する

  

顔アイコン

Pie Plots


pie - draw a pie
pie(x [, sp [,txt]]) 

x is a scalar or a vector of real, positive numbers.
sp is a real scalar or vector.
txt means a cell or a vector of strings.





pie([1 3 4 5], [1 0 1 0],...
["Part1 Exp.","Part2","Part3 Exp.", "P4"]);
// [1 3 4 5] is the number and size of sectors
// [1 0 1 0], 1 for exploded sector, 0 for nonexploded
// ["Part1 Exp.","Part2","Part3 Exp.", "P4"] is the associated
// text削除

2013/11/3(日) 午前 5:46ogw*og*2返信する

  

顔アイコン

Signal processing for sound & music
www.mnt-aq.it/english/dsp_tut.htm

Here you can find the site from where you can download Scilab. These Scilab and TclTk scripts are widely commented to show how things work, and the use of ...削除

2013/11/3(日) 午前 5:52ogw*og*2返信する

  

顔アイコン

A HID USB Host Driver - AHID and Scilab
ahidlib.com/index.php/ahid-and-scilab

Make AHidDemoScilab to your current folder in Scilab and type: ... In the following sections, we describe the usage of the wrapper functions from wrapper. sci.削除

2013/11/3(日) 午前 6:12ogw*og*2返信する

  

顔アイコン

Sound Processing with Short Time Fourier Transform - Ceremade
www.ceremade.dauphine.fr/~peyre/numerical-tour/.../audio_1_processin...

For Scilab user: you must replace the Matlab comment '%' by its Scilab ... around a point t0 of x, which is the FFT (use the function fft) of the windowed signal x.削除

2013/11/3(日) 午前 6:22ogw*og*2返信する

  

顔アイコン

SIP(Scilab Image Processing Toolbox) + Scilab
www.google-melange.com/gsoc/proposal/review/google/.../8001

Firstly, It involves usage of SIP (Scilab Image Processing Toolbox),which leverages the extremely simple Scilab programming environment for prototyping ...削除

2013/11/3(日) 午前 6:28ogw*og*2返信する

  

顔アイコン

The working directory assignment
On windows system, type pwd at the command prompt to determine current working directory. Usually it is c:\documents and settings\user. It would be a good idea to change it to other directory. This is accomplished by following the steps given below.削除

2013/11/3(日) 午前 6:35ogw*og*2返信する

  

顔アイコン

1. Create a new folder in the directory C:\Program Files\scilab-4.1\ and name it work. This will create a new working folder with directory address C:\Program Files\scilab-4.1\work
2. The working directory can be changed to this new director by typing the following command at the command prompt
chdir('C:\Program Files\scilab-4.1\work');削除

2013/11/3(日) 午前 6:35ogw*og*2返信する

  

顔アイコン

Youngstown State University
College of Engineering & Technology
Civil & Environmental/Chemical Engineering Program
ENGR6924: Computer Based tools for Engineers

http://www.eng.ysu.edu/~jalam/engr6924s07/sessions/session23/session23.pdf削除

2013/11/3(日) 午前 6:36ogw*og*2返信する

  

顔アイコン

[PDF]
download - Engineering Technology - Youngstown State University
www.eng.ysu.edu/~jalam/engr6924s07/sessions/session23/session23.pdf

Introduction to SCILAB interface and usage. The first chapter of the book “ Introduction to SCILAB” that is discussed in the session 23 can be can be downloaded ...削除

2013/11/3(日) 午前 6:37ogw*og*2返信する

  

顔アイコン

Think@SysFera: Distributing Scilab computations with DIET
blog.sysfera.com/2012/10/distributing-scilab-computations-with.html
Scilab is a freely distributed open source scientific software package, ... a high level of transparency for the user and limits the wrong usage of the software.

http://blog.sysfera.com/2012/10/distributing-scilab-computations-with.html削除

2013/11/3(日) 午前 6:40ogw*og*2返信する

  

顔アイコン

RTC-scilab

http://ysuga.net/rtc/122?lang=en

http://ysuga.net/rtc/122
RTC-scilabは，scilabで開発したコントローラと，RTミドルウエアによって通信が標準化されたロボットシステムとの接続を容易にするツールボックスです．削除

2013/11/3(日) 午前 6:52ogw*og*2返信する

  

顔アイコン

Introduction to lp_solve 5.5.2.0

http://lpsolve.sourceforge.net/5.5/削除

2013/11/4(月) 午前 0:37ogw*og*2返信する

  

顔アイコン

lp_solveは、混合整数線形計画法 (Mixed Integer Linear Programming, MILP)のソルバーであり、純粋な線形、(混合)整数／バイナリ、半連続、特別な指示によるセット(SOS)モデルなどへ解を与えます。lp_solveは、ANSI Cで記述され、LinuxやWindowsのような多くの異なるプラットフォーム上でコンパイルすることができます。
http://sourceforge.jp/projects/sfnet_lpsolve/削除

2013/11/4(月) 午前 0:44ogw*og*2返信する

  

顔アイコン

Scilabを用いたcubeSPの周波数特性補正処理
http://blogs.yahoo.co.jp/cubesp_speaker/66662133.html
scilab: 伝達関数を使ってＲ－Ｌ直列回路の過度電流を描く。
http://blogs.yahoo.co.jp/puccinisan/16754336.html
数値処理＞SciLab>scicos()では:EXCEL file読込
http://blogs.yahoo.co.jp/ogw2ogw2/20879561.html削除

2013/11/4(月) 午前 11:25[ ogw*og*3 ]返信する

  

顔アイコン

フリーアーム式プリズムホルダーUD
http://store.shopping.yahoo.co.jp/c-task/04-0123.html削除

2013/11/5(火) 午前 6:53[ ogw*og*3 ]返信する

  

顔アイコン

2D plot
Calling Sequence
plot(y,<LineSpec>,<GlobalProperty>)
plot(x,y,<LineSpec>,<GlobalProperty>)
plot(x1,y1,<LineSpec1>,x2,y2,<LineSpec2>,...xN,yN,<LineSpecN>,<GlobalProperty1>,<GlobalProperty2>,..<GlobalPropertyM>)
plot(<axes_handle>,...)

http://help.scilab.org/docs/5.4.1/en_US/plot.html削除

2013/11/10(日) 午後 1:35ogw*og*2返信する

  

顔アイコン

x = -1:0.1:1; y = -1:0.1:1; [X,Y] = meshgrid(x,y); for i=1:size(X,1) for j=1:size(X,2) Z(i,j) = sinc(2*%pi*X(i,j)*Y(i,j)); end end surf(X,Y,Z)

http://www.physicsforums.com/showthread.php?t=577993

Read more: http://www.physicsforums.com削除

2013/11/10(日) 午後 6:05ogw*og*2返信する

  

顔アイコン

x = -1:0.1:1;
y = -1:0.1:1;
[X,Y] = meshgrid(x,y);
for i=1:size(X,1) 
for j=1:size(X,2)
Z(i,j) = sinc(2*%pi*X(i,j)*Y(i,j));
end
end
surf(X,Y,Z)



http://www.physicsforums.com/showthread.php?t=577993

Read more: http://www.physicsforums.com削除

2013/11/10(日) 午後 6:16ogw*og*2返信する

  

顔アイコン

x = -1:0.1:1;
y = -1:0.1:1;
[X,Y] = meshgrid(x,y);
for i=1:size(X,1) 
for j=1:size(X,2)
Z(i,j) = atan*1+atan*2/X(i,j));
end
end
surf(X,Y,Z)削除

2013/11/10(日) 午後 6:21ogw*og*2返信する

  

顔アイコン

[PDF] 
Scilab古典制御プログラム
home.hiroshima-u.ac.jp/saeki/Appendix/ClassicalControl_Scilab.pdf
2014年4月7日 ... 制御性能の解析や制御系設計では，時間応答のシミュレーションや周波数応答などを 計 ... [4] 佐伯正美、Scilab ボード線図関数プログラム bode2.sci、[2]のホームページに 掲載 .... P(s)のボード線図を [0.01,10][rad/s]の周波数区間で描け.削除

2016/6/8(水) 午前 2:13[ ogw*og*3 ]返信する

  

顔アイコン

PID制御EXCEL版
www.geocities.jp/mytecana/pidcon3.html
PID制御をEXCELでシミュレーションする技法を解説し,I-PD制御,PIDセルフ チューニング,システム同定などの技術についても紹介している。

http://www.geocities.jp/mytecana/pidcon3.html削除

2016/6/8(水) 午前 6:34[ ogw*og*3 ]返信する

  

顔アイコン

PIDとI-PD制御シミュレーション
www.geocities.jp/mytecana/pidsimulation.html
Excel版. 自動制御の数値実験技法. ＰＩＤとＩ-ＰＤ制御のシミュレーション. ◇◇◇ ＰＩＤ（図 １）とＩ-ＰＤ（図２）制御の違いを体験しよう！ ... 手順１： PID制御 or I-PD制御のどちらか を選択します。 ＰＩＤ制御 Ｉ-ＰＤ制御 手順2： Ｐ、Ｉ、Ｄ定数を入力し、スタートボタンを押し ...

http://www.geocities.jp/mytecana/pidsimulation.html削除

2016/6/8(水) 午前 6:36[ ogw*og*3 ]返信する

  

顔アイコン

Engineers-Excel.com: PID Loop Simulator
www.engineers-excel.com/Apps/PID_Simulator/Description.htm
PID Loop Simulator. ... This app simulates a control loop consisting of a first order process with deadtime and a PID controller. An advanced version of this app suitable ... Integrating Process · Simulation of a Model Predictive Controller (MPC ) ...

http://www.engineers-excel.com/Apps/PID_Simulator/Description.htm削除

2016/6/8(水) 午前 6:40[ ogw*og*3 ]返信する

  

顔アイコン

エクセルのちょっとフェチな使い方。

http://www.irs-japan.com/?p=541削除

2016/6/8(水) 午前 6:45[ ogw*og*3 ]返信する

  

顔アイコン

エクセルのちょっとフェチな使い方。 - 株式会社 アイ・アール・エス ...
www.irs-japan.com/?p=541
2013年3月27日 ... インジェクター洗浄とは直接は関係のない、ある２次遅れ系の制御に関する簡単な シミュレーションだったので、エクセル上で２次遅れ系の関数とＰＩＤ制御器を作って全部 やってしまったんですが、使い方によってはエクセルって実はかなりの ...削除

2016/6/8(水) 午前 6:50[ ogw*og*3 ]返信する

  

顔アイコン

PID制御と二自由度制御のシ ミュレーション

http://www.02.246.ne.jp/~m-adachi/sim/sim3.htm

PID制御のシミュレーションへ
www.02.246.ne.jp/~m-adachi/sim/sim3.htm
PID制御と二自由度制御のシミュレーション. 一次遅れ ... 注）この場合の目標r(t)は常に 零なので、二自由度制御とPID制御は、同じ応答を示す。 ..... 注） m-adachi@02.246. ne.jp へメールで依頼があれば、上記エクセルのソフトをメールの添付ファイルで送り ます。削除

2016/6/8(水) 午前 6:56[ ogw*og*3 ]返信する

  

顔アイコン

エクセルでシミュレーション Vol.10 二次元非定常熱伝導問題シミュレーション＋P(ID)制御エクセルシートを作ってみよう


http://www.hirax.net/keywords/log/エクセルでシミュレーション/latest削除

2016/6/8(水) 午前 7:01[ ogw*og*3 ]返信する

  

顔アイコン

Public Sub Rotation()

Dim t, theta As Double
Const Pi = 3.1415926535

Worksheets("Sheet1").Activate
With Worksheets("Sheet1")

.Range("C3") = 0
.Range("C4") = 0
.Range("C5") = 0

For i = 1 To 2
For t = 0 To 2.001 Step (1 / 32)

theta = Pi * t
.Range("C3") = theta


.ChartObjects(1).BringToFront

Next t

.Range("C4") = -0.5
.Range("C5") = -0.5

Next i

.Range("C3") = 0

End With
End Sub削除

2016/6/8(水) 午前 8:00[ ogw*og*3 ]返信する

  

顔アイコン

7.5 例題のVBAスクリプトの解説

以下に，例題で用いたVBAのステートメントや関数の簡単な説明を載せる．詳しくはVBAのヘルプを参考にしなさい．
Const Pi = 3.14 ...

Pi が3.14...という値を持つ定数(constant)であることを示している．削除

2016/6/8(水) 午前 8:01[ ogw*og*3 ]返信する

  

顔アイコン

Worksheets("Sheet1").Activate

Sheet1を前面に表示し，アクティブにする．
With Worksheets("Sheet1") ... End with

ワークシートSheet1にあるセルを参照するには，通常，Worksheets("Sheet1").Range("C3")としなければならないが，毎回これを書いていては表記が長くなり，読みづらい上に面倒なので，Work...の部分を省略して，.（ドット）から始めて良いようにする．この様に．（ドット）から始まる場合は，それ以前にWith...以下で設定したものが適用される．削除

2016/6/8(水) 午前 8:03[ ogw*og*3 ]返信する

  

顔アイコン

.Range("C3")

Withで設定されたワークシート上のセルC3を参照している．この他に，.Cells(3,3) とする方法がある．
For t = 0 To 2.001 Step (1 / 8) ... Next t

For ... Next で反復（ループ）を表す．ここではループカウンタに浮動小数(Double)の t を用いており，0 から 2 まで，1/8 きざみで（すなわち16回）ループが廻る（反復する）ことを示している．
.Range("C3") = theta

セルC3にthetaの値を代入している．これにより，ワークシート上に作った表のθの値が書き換えられ，結果，回転行列と，それを用いた行列・ベクトルの演算が行われる．
ChartObjects(1).BringToFront

グラフを最前面に最前面に移動し，再描画を行う．

また，例題プログラムの中では使われていないが，
Application.Wait (Now + TimeValue("0:00:01"))

とすると，実行を1秒間停止させて，計算過程をゆっくり表示させ削除

2016/6/8(水) 午前 8:04[ ogw*og*3 ]返信する

  

顔アイコン

半日で作る倒立振子

http://www.instructables.com/id/半日で作る倒立振子/削除

2016/6/8(水) 午後 1:18[ ogw*og*3 ]返信する

  

顔アイコン

Scilabで物理シミュレーション（手習い）放物線ボール投げ - Qiita
qiita.com/HirofumiYashima/items/e5f8660281cbf7728fa7
2014年7月1日 ... 教材ウェブサイト （リンク）亀山 建太郎（福井大学）「Scilabで楽々数値計算」 （実行 スクリプト） Scilab console h = 1; // ボール初期高度v = 20; // 投げ上げ初速度g = 9.8; // 重力加速度dt = 0.1; // 時間刻みt = 0; // 時刻STEP = 1; ...

http://qiita.com/HirofumiYashima/items/e5f8660281cbf7728fa7


ball-test.sce削除

2016/6/8(水) 午後 2:44[ ogw*og*3 ]返信する

  

顔アイコン

2. スクリプトファイルの作成


http://www.phys.u-ryukyu.ac.jp/~myogi/member/yogi/scilab/calc2.htm

scipadと打ち込んで実行（Enterキーを押す）か，窓の上の方にあるEditorをクリックして下さい。そうすると下記のような窓が立ち上がると思います。これがエディタです。削除

2016/6/8(水) 午後 2:57[ ogw*og*3 ]返信する

  

顔アイコン

ScilabによるIIRフィルタのシミュレーション

http://bfin.sakura.ne.jp/?p=319削除

2016/6/8(水) 午後 3:47[ ogw*og*3 ]返信する

  

顔アイコン

FFTしてみる
Xcos上ではFFTを直接実行するモジュールはなさそう．
なので，実行結果をScilabに持って行ってFFTする．

FFTしたい信号に "ToWorkspace"モジュールを接続する．
Size of bufferに保存する時刻数，variable nameに変数名を指定する．

Xcosでシミュレーションを実行する．削除

2016/6/8(水) 午後 5:30[ ogw*og*3 ]返信する

  

顔アイコン

次に，Scilabで以下の関数を定義する．

function MyFFT(Variable)

val=Variable.values;
time=Variable.time;
dt=time(2)-time(1);
val_t=val';
N=size(val,'r');
freqs=1/dt*(0:N/2)/N;
n=size(freqs,"*");

scf(1);
clf(1);
plot(time, val);

scf(2);
clf(2);
fft_res=fft(val_t);
plot2d(freqs, abs(fft_res(1:n)));

endfunction

Scilabのコマンドで MyFFT( 変数名 )とすると，FFTの結果がWindowに表示される．
表示されるWindow番号は，Functionのscf(1), clf(1)あたりを変更することで変えられる．削除

2016/6/8(水) 午後 5:30[ ogw*og*3 ]返信する

  

顔アイコン

FFTしてみる

http://freeplanets.ship.jp/Scilab/Scilab-TIPS.html削除

2016/6/8(水) 午後 5:31[ ogw*og*3 ]返信する

  

顔アイコン

休止期間の無いTSPを用いたインパルス・レスポンス測定法のScilabシミュレーション・プログラム

http://www.cepstrum.co.jp/download/cont_tsp/cont_tsp.html削除

2016/6/8(水) 午後 5:42[ ogw*og*3 ]返信する

  

顔アイコン

【プログラムの主なな概要】
2行目ではサンプリング周波数1000Hzとします。
3行目ではサンプリング時間を0.5(s)で時間ステップは0.001(s)です。
4行目ではサンプリング数で501個となります。
6行目でノイズがない場合の50,90Hzの合成波を求めています。
7行目でサブウインドウを指定しています。
8行目では合成波を表示しています。
9行目ではFFTの計算を行っています。
10行目で周波数を求めています。
11行ではサブウインドウを指定しています。
12行目ではfを横軸にFFTの値を表示します。

http://scilab.kani33.com/2015/05/post-523/削除

2016/6/8(水) 午後 9:23[ ogw*og*3 ]返信する

  

顔アイコン

syslin

A=[0,1;0,0];B=[1;1];C=[1,1];
S1=syslin('c',A,B,C) //線形システムを定義
S1("A") //A行列を表示
S1("X0"), S1("dt") // X0と時間領域を表示
s=poly(0,'s');
D=s;
S2=syslin('c',A,B,C,D)
H1=(1+2*s)/s^2, S1bis=syslin('c',H1)
H2=(1+2*s+s^3)/s^2, S2bis=syslin('c',H2)
S1+S2
[S1,S2]
ss2tf(S1)-S1bis
S1bis+S2bis
S1*S2bis
size(S1)

https://help.scilab.org/docs/5.4.1/ja_JP/syslin.html削除

2016/6/9(木) 午前 5:36[ ogw*og*3 ]返信する

  

顔アイコン

// s131 - FFT

Fs=44100
t=0:1/Fs:5;
y=0.05*sin(2*%pi*440*t)./(0.1+(t-0.3).^2) ...
+ 0.05*sin(2*%pi*770*t)./(0.1+(t-1.1).^2) ;

//[y,Fs,bits]=wavread('sound.wav');

playsnd(y,Fs,bits);

m=8192;
d=0;
clf();
for n=1:7削除

2016/6/9(木) 午前 6:03[ ogw*og*3 ]返信する

  

顔アイコン

yy=y(1,1+(n-1)*m:n*m);
//wavsave('sound2.wav'，yy);
playsnd(yy,Fs,bits);

scf(d);
// d=d+1;

subplot(2,1,1);
plot2d*3-200*n,style=n);

end

//end

http://user.numazu-ct.ac.jp/~mochizuki-k/scilab/削除

2016/6/9(木) 午前 6:03[ ogw*og*3 ]返信する

  

顔アイコン

scilab のワークフォルダ
の確認コマンド


-->pwd
ans =

C:\Documents and Settings\user\My Documents 



http://www.wavsource.com/sfx/sfx.htm削除

2016/6/9(木) 午前 7:07[ ogw*og*3 ]返信する

  

顔アイコン

// sound.wavというファイルを読み込んで、その波形とＦＦＴ結果を表示する。

[y,Fs,bits]=wavread('sound.wav');
y=y(1,:)
//片方のチャネルのみ使う

disp(Fs)
disp(bits)
// Fs=サンプリング周波数、bits=量子化ビット数を表示

m=32768;
// 入力したデータサイズの内輪でもっとも大きい２のべき乗の値をいれる。

clf();

yy=y(1,1:m);

scf(0);

subplot(2,1,1);
plot2d*4-200,style=1);
//修正しました。

http://blogs.yahoo.co.jp/susanoo2001_hero/6652650.html削除

2016/6/9(木) 午前 7:08[ ogw*og*3 ]返信する

  

顔アイコン

[6]ログファイルの利用

どのような計算をしてたかを確かめておくことは重要です。そのために、重要な計算の前では必ずログを取っておくべきでしょう。
ログをとり始めるコマンドは、
diary('ファイル名')
ログを終了させるコマンドは、
diary(0)
です。


このようにしておくと、(今の場合)z:\wd_scilabフォルダに、test.logというファイルが保存され、中身はこんな感じでテキストで保存されます。後で見返すときに便利です。

http://www.ic.daito.ac.jp/~tkadoda/initScilab/initScilab.htm削除

2016/6/9(木) 午前 7:37[ ogw*og*3 ]返信する

  

顔アイコン

フリーソフト Scilab 活用講座資料(2014)

https://scilabinfo.wordpress.com/scilab-入門/フリーソフト-scilab-活用講座資料2014/削除

2016/7/12(火) 午前 7:08[ ogw*og*3 ]返信する

  

顔アイコン

scilab つかいませんか - 制御情報工学科 - 沼津高専
user.numazu-ct.ac.jp/~mochizuki-k/scilab/
Last Update : '16.4.24 「応用1a」内のミスプリ１個所の修正と，「外部リンク 応用1a- 2016」を追加; '07.5.8言い回しの修正; '07.2.15 例題の番号を振りなおした; '07.1.31 Filter を Scilab 化できた; '07.1.24 = ミスプリ訂正（Change Directry のメニューが「edit 」 ...削除

2016/7/12(火) 午前 7:11[ ogw*og*3 ]返信する

  

顔アイコン

[PDF] 
はじめての Xcos - Scilab
https://www.scilab.org/content/.../1141/.../Xcos_beginners_jp.pdf
実験を行わなくても、複雑な現象のシミュレーション(物理、機械、電子機器など)は、 それらの振. る舞いや結果などの研究を可能にしています。また産業界では幅広く活用 されていますし、エン. ジニアや科学者を志望する若者にはモデリングやシミュレーション の ...削除

2016/7/12(火) 午前 7:12[ ogw*og*3 ]返信する

  

顔アイコン

Scilab と PHP を使った効率的な数学処理 - IBM
www.ibm.com/developerworks/jp/opensource/.../os-php-scilab/
2010年3月16日 ... ... できます。コミュニティー主導の科学を実現するための 1 つの方法として、PHP を Scilab に結びつける方法について調べましょう。 ... 最近のサーバー・サイド Web 開発 には、Ruby、Python、PHP などのスクリプト言語が活用されています。


http://www.ibm.com/developerworks/jp/opensource/library/os-php-scilab/削除

2016/7/12(火) 午前 7:23[ ogw*og*3 ]返信する

  

顔アイコン

Scilab で楽々数値計算
亀山 建太郎（機械工学科）

http://crf.flib.u-fukui.ac.jp/dspace/bitstream/10461/4628/1/vol.50.p31-36.pdf削除

2016/7/12(火) 午前 7:51[ ogw*og*3 ]返信する

  

顔アイコン

// Throw ball
clear;

h = 1; // ボール初期高度
v = 20; // 投げ上げ初速度
g = 9.8; // 重力加速度

dt = 0.1; // 時間刻み
t = 0; // 時刻

STEP = 1; // カウンタ
STORE_H = h; // 高さh を保存するベクトル
STORE_T = t; // 時間t を保存するベクトル

while(h>0)
dh = v*dt - g*t*dt; // ボール高度の更新
h = h + dh;

STEP = STEP + 1; // カウンタを進める
t = t + dt; // 時刻を進める

STORE_H(STEP+1) = h; // データの保存
STORE_T(STEP+1) = t;
end

plot(STORE_T, STORE_H); // グラフ表示削除

2016/7/12(火) 午前 7:52[ ogw*og*3 ]返信する

  

顔アイコン

scilabの使用法6-シミュレーション例（比熱） - 琉球大学
www.phys.u-ryukyu.ac.jp/~myogi/member/yogi/scilab/calc6.htm
シミュレーション例（比熱）. 電気抵抗と同じように，比熱のシミュレーションを行って みたいと思います。格子振動による比熱はアインシュタインの式とデバイの式があります 。 アインシュタインの式. デバイの式. 特性温度をそれぞれ300Kとして，比熱の温度変化 を ...

http://www.phys.u-ryukyu.ac.jp/~myogi/member/yogi/scilab/calc6.htm削除

2016/7/12(火) 午前 9:17[ ogw*og*3 ]返信する

  

顔アイコン

2012-05-04
Scilab を使ってLPC1343とシリアル通信 & データ処理 ADCデータのリアルタイムデモ

scilab

http://d.hatena.ne.jp/hira-rin/20120504/1336112893削除

2016/7/13(水) 午前 1:13[ ogw*og*3 ]返信する

  

顔アイコン

ベイズの定理（入門編） - Pashango's Blog
d.hatena.ne.jp/pashango_p/20090809/1249805193
2009年8月9日 ... ベイズの定理」とは、まさしく今まで説明してきた「条件付き確率」を求めるための計算式 なのです。 まずは「ベイズの定理」を見て .... この「例外」（←数学ではあまり見かけない 言葉ですが）の具体例を教えてください。 私は「P(A∧B)=P(A)」の方 ...削除

2016/7/13(水) 午後 2:22[ ogw*og*3 ]返信する

  

顔アイコン

[PDF] 
ベイズの定理
antlers.cis.ibaraki.ac.jp/KAKURITU/331.pdf
ベイズの定理. 問題１. 次の確率を考察せよ。 (a)袋1に赤玉が3個、白玉が2個入って いる。 (b)袋2に赤玉が8個、白玉が4個入っている。 (c)袋1と袋2が選ばれる割合は、1:2 とする。 （１）赤玉が選ばれる確率. （２）赤玉が選ばれたときに袋1から選ばれている確率 .

http://antlers.cis.ibaraki.ac.jp/KAKURITU/331.pdf削除

2016/7/13(水) 午後 2:33[ ogw*og*3 ]返信する

  

顔アイコン

[PDF] 
ベイズの定理と確率
school.gifu-net.ed.jp/ena-hs/ssh/H23ssh/sc3/31111.pdf
いうものがあると知ったので、ベイズ確率やベイズの定理について研究を進めていった。 ベイズの定理. は事前確率や事後確率などと深く関係があることが分かり、「モンティ ホール問題」や、身近なところ. ではスパムメールの検出に使われていることも分かった 。削除

2016/7/13(水) 午後 2:38[ ogw*og*3 ]返信する

  

顔アイコン

待ち合わせに遅れる彼女、ベイズの定理、そして例題 - MyEnigma
myenigma.hatenablog.com/entry/20110128/1296220831
2011年1月28日 ... はじめに. 自律ロボットを研究する上で、. ベイズの法則は欠かすことのできない理論だと 思います。 ロボットの位置推定や物体認識の一般的な技術である. カルマンフィルタや パーティクルフィルタは. このベイスの法則を元に定式化されています ...削除

2016/7/13(水) 午後 2:40[ ogw*og*3 ]返信する

  

顔アイコン

ベイズの定理 - So-net
www004.upp.so-net.ne.jp/s_honma/probability/bayes.htm
いま、ベイズの定理が熱い。古典的確率論の一つの定理であるが、現在、検索エンジン でよく利用される Ｇｏｏｇｌｅ の高いヒット率を支えたり、Ｉｎｔｅｌ や Ｍｉｃｒｏｓｏｆｔ における アプリケーション開発の数学的基礎として注目を集めているらしい。 どのような形で ベイズの ...削除

2016/7/13(水) 午後 2:41[ ogw*og*3 ]返信する

  

顔アイコン

ベイズ推定 集合論 例題 解法比較 市民のための環境学ガイド
www.yasuienv.net/BayesExamples.htm
2013年2月10日 ... 前回から始めたベイズ推計ですが、そこで提示した例題を、（ａ）ベイズの定理を用いた 解法、（ｂ）古典的・頻度主義的な解法の二種類の解を示 ... 解説：この例題から分かる ように、この問題をベイズの定理を使って解く理由は特にない。普通に ...削除

2016/7/13(水) 午後 2:44[ ogw*og*3 ]返信する

  

顔アイコン

[PDF] 
ベイジアンになろう
www.iba.t.u-tokyo.ac.jp/iba/SE/shiryou/Baysian-print.pdf
例題. •あなたは乳ガンの検査（マンモグラフィ）で陽性になった。 •この検査では、乳ガン になっている人を見つける正確さ. は９０％、 ... ベイズの定理. ▫ 用途. 結果が起きたとき 、その原因を調べる. ▫ 結果の確率（事前確率）. ▫ 原因の確率（事後確率）. ▫ 原因.

http://www.iba.t.u-tokyo.ac.jp/iba/SE/shiryou/Baysian-print.pdf削除

2016/7/13(水) 午後 3:00[ ogw*og*3 ]返信する

  

顔アイコン

[PDF] 
確率論（第4講）
www.k.hosei.ac.jp/~yana/pr09/probability04-09.pdf
例題. 喫煙者. 非喫煙者. 肺癌になる. 1200. 300. 肺癌になる. 1200. 300. 肺癌に ならない 1900. 6600. 上記統計をもとに次の問いにこたえなさい. 喫煙者の割合 ... ベイズ則とよばれ統計的決定理論、推定理論の基礎となる。 イズ則とよ ... ベイズの 定理例題. ▫.


http://www.k.hosei.ac.jp/~yana/pr09/probability04-09.pdf削除

2016/7/13(水) 午後 3:11[ ogw*og*3 ]返信する

  

顔アイコン

C:\Users\aka\Documents\SCiLab\SAMPL削除

2016/10/1(土) 午後 11:56[ ogw*og*3 ]返信する内緒

  

顔アイコン

clear

// 入出力ファイル名
infile='C:\Users\aka\Documents\SCiLab\SAMPL\waipa.au';
outfile='C:\Users\aka\Documents\SCiLab\SAMPL\waipa2.au';

// 音声読み込み
[dat, fs] = auread(infile);

// グラフ表示
clf();
len=length(dat);
time=1000*[0:1/fs:(len-1)/fs];
plot(time, dat);
title(infile);
xlabel('time[ms]');

// 音声書き出し
auwrite(dat, 2*fs, outfile);削除

2016/10/2(日) 午前 0:50[ ogw*og*3 ]返信する内緒

  

顔アイコン

// 入出力ファイル名
infile='oden.au';
outfile='oden2.au';

// 音声読み込み
[dat, fs] = auread(infile);

// グラフ表示
clf();
len=length(dat);
time=1000*[0:1/fs:(len-1)/fs];
plot(time, dat);
title(infile);
xlabel('time[ms]');

// 音声書き出し
auwrite(dat, 2*fs, outfile);

http://vision.kuee.kyoto-u.ac.jp/lecture/dsp/?menu=mat_io削除

2016/10/2(日) 午前 0:51[ ogw*og*3 ]返信する

  

顔アイコン

-->im = imread('C:\Users\aka\Documents\lena.png');

-->imshow();削除

2016/10/2(日) 午前 1:45[ ogw*og*3 ]返信する内緒

 

全2ページ

[1] [2]

 

*1:Y(i,j)-2)/X(i,j

*2:2-Y(i,j

*3:1+(n-1)*m:n*m)/Fs, yy,style=n);

zz=fft(yy);

subplot(2,1,2);
plot2d((1:m/16)*Fs/m+60*n, ...
abs(zz(1:m/16

*4:1:m)/Fs, yy,style=1);

zz=fft(yy);
// yyをfftします。

subplot(2,1,2);
plot2d((1:m/16)*Fs/m, abs(zz(1:m/16