MATLABプログラミング 質問箱 その1
レス数が1000を超えています。これ以上書き込みはできません。
0001デフォルトの名無しさん
NGNGhttp://www.cybernet.co.jp/matlab/
MATLAB 本家
http://www.mathworks.com/
Yahoo MATLAB
http://dir.yahoo.co.jp/Science/Mathematics/Software/MATLAB/
MATLAB入門
http://feedback.mech.fukui-u.ac.jp/matlab/matlab0.htm
MATLAB チュートリアル
http://www-isys.eei.metro-u.ac.jp/~sho/pub/ttmatlab/Index.html
0002デフォルトの名無しさん
NGNG0003デフォルトの名無しさん
NGNGこのスレにお世話になることあると思うんで、盛り上げていきましょう。
0004デフォルトの名無しさん
NGNGって、俺はほとんど使ったこと無いけど。
普段は octave 使ってる。
0005デフォルトの名無しさん
NGNGもしかしたら >3 と知り合いかも。
・・・って、音響やっててMatlab使っているなんてたくさんいるか。
0007デフォルトの名無しさん
NGNGかくいう俺も音響関係の研究室に所属。信号処理関連のライブラリは魅力的だなあ。
0008デフォルトの名無しさん
NGNGその研究室の先生が教授を辞めてしまった。
0009デフォルトの名無しさん
NGNGFORTRANと比べたらそりゃそうだろう。
けどC/C++/Javaなどと比べたら(ry
0010デフォルトの名無しさん
NGNG東京23区内か？
0011デフォルトの名無しさん
NGNG0012デフォルトの名無しさん
NGNG思ったらPerl程度のものしかできなくて萎えた
0013デフォルトの名無しさん
NGNGソースをぱっと見ただけで
それが関数であるのか、配列であるのか、ベクトルであるのか
解らないところだね。
C/C++/Javaならすぐにわかるのに
0014デフォルトの名無しさん
NGNGいくら音響やっててMATLAB使ってるのがたくさんいるからといって
MATLABだけで音響ソフトを作るのはﾊｽﾞｶｼｲ
0015デフォルトの名無しさん
NGNG思いついたアイディアをすぐに形にできて、実験できるのがいいところ。
その後に自分で書くなり外部のプログラマを雇うなりしてC/C++で最適化します。
0016デフォルトの名無しさん
NGNG今まで音響の研究室見てきたけれど
雇ってもろくなもんができない罠。
本格的にやろうとすればするほど金がかかるからね。
100～500万じゃ雇っても無駄。
プログラマという職についてみればあのソフトには何100万じゃ足りないということがわかてしまう
0017デフォルトの名無しさん
NGNG自分の研究成果を組み込み機器に取り入れたがるうえに
かつ考え方が古いので
いつまでもC言語にこだわる。
0018デフォルトの名無しさん
NGNG1人月30万で3人でチーム組んだとすると
一ヶ月でその100万が消えてしまう。
その三人が無能だったら使えないソフトしか作れない。
0019デフォルトの名無しさん
NGNG対話形式やウィザード、スクリプトが大好きで
設計が苦手な傾向にあるようだ。
彼らには型の違いに厳しいC/C++/Javaという言語が苦手だという傾向もある。
0020デフォルトの名無しさん
NGNGMATLAB付属以外のものを
使ってる香具師いる？
0021デフォルトの名無しさん
NGNG個人としてはMatlabよりはCを使うことが多いんだけど、
Matlabの優秀さは認めています。
基本的な数学がわかっていればプログラミングの訓練を
しなくても短期間でそれなりのことができちゃうからね。
Matlabの初(?)スレとしてはお勧め書籍とか値段とかを書くのがよいのかも。
0022デフォルトの名無しさん
NGNGこの本、ほとんどがDSPのことばかり
でMATLABのことやソースコードは
ほとんど載っていなかった。
0023デフォルトの名無しさん
NGNG熱心に使おうという気が起きない。
0024デフォルトの名無しさん
NGNGプログラムはCとかで書いた方が速いし
0025デフォルトの名無しさん
NGNG> matlabは主にデータのプロット用だね
あはは。ちょっと贅沢なgnuplot、ってところか。
0026デフォルトの名無しさん
NGNG学会発表とか授業のプレゼン準備とかにMatlabでプロット作成してます。
EPSに吐けるし、いろいろと細かいところまでいじれるのがいいですね。
0027デフォルトの名無しさん
NGNG同じことをやってるひとがいるとは
0028デフォルトの名無しさん
NGNGデータ形式はどのようなものを用意して、どのようなコマンドでプロットすればよいのでしょうか？
0029デフォルトの名無しさん
NGNGMATLABの「LMI Control Toolbox」のLMIエディッタというのはどう言う風に使えばいいのか・・・
初めて触ったけど全然わからない；；
これ使ったら行列不等式とけるってわけではないのかな？Orz
0030デフォルトの名無しさん
NGNGを表示したい場合はどのようにすればよいのでしょうか？
0031デフォルトの名無しさん
NGNGx=0:0.1:pi;
y=0:0.1:2*pi;
[X,Y]=meshgrid(x,y);
surf(X,Y,sin(X).*cos(Y));
0032デフォルトの名無しさん
NGNGWarning: Unable to load Java Runtime Environment: /usr/local/matlab/sys/java/jre/glnx86/jre/lib/i386/native_threads/libjava.so: symbol __libc_wait, version GLIBC_2.0 not defined in file libc.so.6 with link time reference.
Warning: Disabling Java support.
ってエラーが出て GUI が起動しません
だれか fedora 2 で GUI を使えている人いますか？
0033デフォルトの名無しさん
NGNGたとえば
A=[1 0;0 1]
B=[0 0;1 1]
C=[0 1;1 0]
A*B*X+C>0
になるようなX
とか><
ぜんぜんわかんない～；；
0034デフォルトの名無しさん
NGNG0035デフォルトの名無しさん
NGNG0036デフォルトの名無しさん
NGNG左辺が対称行列じゃなきゃダメ。
0037大学生くん
NGNGどなたか教えてください。
DCTを使ってどうにかできないものなのでしょうか？
0038デフォルトの名無しさん
NGNGCやC++などでやるのがよろしいんじゃなくて？
0039デフォルトの名無しさん
NGNG0040デフォルトの名無しさん
NGNGもちろん使っているコマンドだかオプションは
deps2cとかいうやつだよな？
カラーでな。
みんなeps使わず出力された
ポップアップウィンドウから
マウスでコピペしてる奴ばかり。
一つのグラフ書くたびに
毎回Xtick, Ytickを指定しset（’gca' とか
をコピペして
効率悪すぎ。
クラス作れなくても関数くらい作れや。
攻めてforとかwhileループとか
使って効率化くらいしたらどうだ
っていいたくなる奴ばかり。
0041デフォルトの名無しさん
NGNGシミュレーションに使うには
重たくて使いにくそうだな
0042デフォルトの名無しさん
NGNGつーか数値計算や信号処理の本でも嫁
0043デフォルトの名無しさん
NGNGシェルで
demoと入力して
エンターを押せ
そこに答えが見つかる
0044デフォルトの名無しさん
NGNG他の作業ができないのですが、matlabでディスプレイ上に表示させないで、aviムービーを作成することは可能でしょうか？
0045Mach Go! Go! ◆6opQoo0UZg
NGNG表示内容はなんですか？
画像データとかであれば表示させずにimage変数内で用意して書き出せるんですが。
004644
NGNG例えばこのようなものです
Z = peaks; surf(Z)
axis manual % Axes の範囲を固定します。
set(gca,'nextplot','replacechildren');
M = moviein(20);
for j = 1:20
surf(sin(2*pi*j/20)*Z,Z)
M(:,j) = getframe;
end
0047Mach Go! Go! Go! ◆6opQoo0UZg
NGNGttp://www.mathworks.com/support/tech-notes/1200/1204.html#Section%2010
二次元ではできてももしかしたら三次元だとできない、とかそういう制約があるかもしれませんね。
0048Mach Go! Go! Go! ◆6opQoo0UZg
NGNG>調べてみたら、画面表示なしにMatlab Tech Noteに載っていました。
調べてみたら、画面表示なしにムービーを作る方法がMatlab Tech Noteに載っていました。
0049デフォルトの名無しさん
NGNGところで、皆は行列計算時に Out of memory が出たら
どう対応してる？よく repmat でメモリ不足して困る。
今のところ行列を適当に小行列に分解し、for ループで
まわしてるんだが、やはり遅さは否めない。
なんとかならんものかのぅ...
0050デフォルトの名無しさん
NGNG＞どう対応してる？
そっち面がシビアなものにはMatlabは使わない。
0051デフォルトの名無しさん
NGNG以下を参照。
http://www.linuxquestions.org/questions/showthread.php?s=&forumid=7&threadid=183462
www.blackdown.org からj2re1.3.1をとってきて入れて、適当にパスを張って
ln -s /where-you-install-blackdownjava/j2se/1.3.1/jre /local-matlab-dir/sys/java/jre
環境変数を
export LD_ASSUME_KERNEL=2.4.19
とすれば、fedora core 2でmatlab12が動くようになる。
多分、上のリンクの
export LD_PRELOAD=~/libcwait.so
の部分はやらなくてもうまくいった。
0052デフォルトの名無しさん
NGNGどうもありがとうございます。無事できました。
が、サイズがとても小さいのですが、大きくする方法はないのでしょうか。
GUIの場合は、windowを大きくすれば良かったのですが...
0053デフォルトの名無しさん
NGNGそうか...
そう大規模でもないんだがなぁ。たかだか400サンプルを
もとに計算しているだけだが、如何せん途中で大規模な
行列が生じるためにOut of memory...
ネックになっている行列操作だけ C で書いて走らせるか
なぁ。Optimization Toolbox が使いたいから、MATLABを
諦めきれん。
0054Mach Go! Go! Go! ◆6opQoo0UZg
NGNGメモリ管理に関してのTech Noteなんてのを見つけましたが、参考になりますでしょうか？
ttp://www.mathworks.com/support/tech-notes/1100/1106.html
「変数は連続したメモリ空間にしか確保できないので、メモリを大量に使用することが分かっている場合は、あらかじめhoge = zeros(n,m)で確保しておくべし」ということらしいですが。
0055デフォルトの名無しさん
NGNGそんなフィルタ、何使えばいい？
0056Mach Go! Go! Go! ◆6opQoo0UZg
NGNG[b, a] = butter(3, [0.4 0.6]);
y = filter(b, a, x);
0057デフォルトの名無しさん
NGNGわざわざ探してくれてありがとう．
扱う行列がでかすぎるのが問題なんだよね．
WS.調べたけど，1Gメモリ積んでても連続領域が460M
ちょいしかなかったのが原因だろうな．Windows 起動
時でそれなので，460M を越えないサイズで計算させる
しかないね．
せっかく計算をベクトル化したのになぁ(´･ω･｀)ｼｮﾎﾞｰﾝ
0058デフォルトの名無しさん
NGNG0059デフォルトの名無しさん
NGNGMATLABでJavaをよみこんでいるのか？
だったらJava単体コードを実行したときには
java.lang.OutOfMemoryErrorはでないのか？
0060デフォルトの名無しさん
NGNG＞MATLABでJavaをよみこんでいるのか？
GUI周りはJava使ってるよ。
0061Mach Go! Go! Go! ◆6opQoo0UZg
NGNG>49さんは「ところで」の部分から違う話題に入ったっぽいので、
Javaとメモリ不足とは関係ないみたいです。
R13くらいからのMatlabではJavaのGUIコマンドを実行することもできるようになってるみたいです。
Swingを使ってウィンドウを出したりして遊んだことがあります。
0062デフォルトの名無しさん
NGNG0063matlab help!!
NGNGhighpass butterworth filterをデザインしようとしてます。
同じ設定のlowpass filterをデザインすることはできたんですが、
そこから、highpass filterに変換することができません。誰かご存知だったら、
お願いします。lowpassのコードはこれです。
>> [B1,A1]=butter(2,1000,'s');
>> f=logspace(2,5);
>> h1=freqs(B1,A1,f);
>> mag1=20*log10(abs(h1));
>> semilogx(f,mag1),grid;
0064デフォルトの名無しさん
NGNG0065デフォルトの名無しさん
NGNG0066matlab help!!
NGNG0067choix ◆pnCm4bpG3Y
NGNG>> [B1, A1] = butter(2, 1000, 'high', 's');
じゃダメですか？
0068デフォルトの名無しさん
NGNGdx(i)/dt = P(i) - k(i)x + a(i)x(i) - b(SUMx(j))
ここでxは2つあり、x＝１の時は以下の式で解を求めます。
When D <> 0
x１(t) = (Ps/2S + Pd/2D) + (x１(0) - Ps/2S)e^(-St) - (Pd/2D)e^(e-Dt)
When D = 0
x１(t) = xd(0) + Pdt/2 + Ps/2S + (xs(0) - Ps/2S)e^(-St),
where
S = k + B;
D = k - B;
xs=x１+x2;
xd=x１- x2;
Ps = P１+ P2;
Pd = P１- P2;
x=2の場合、Pd -> -Pd, x1 -> x2, xd -> -xdとし、同じ様に解を求めます。
0069デフォルトの名無しさん
NGNG平均については1を足せばいいとわかりました。
0071デフォルトの名無しさん
NGNG正規化できる。 すなわち、Zは正規分布N(0, 1)に従う。
であるから、 Z～N(0, 1)なる正規乱数に対して、 X = σZ+μ と置けば平均μ、分散
σ×σの正規乱数が生成できる。
要するに、出来た乱数に２をかけて1を足せばいいのだ。
基本中の基本と思うけど。
0072デフォルトの名無しさん
NGNG0073デフォルトの名無しさん
NGNG実装されているACF(AutoCorrelationFunction)の使いかたがわかりません。
[ACF]=acf(Z,N)
とあり引数はZ=Signal,N=# of coefficientsと書いてあります。
Nには何を入れれば良いのでしょう…
0074デフォルトの名無しさん
NGNG約-21億から約21億までの値を
一個づつNに代入しながら
確かめてみよう。
0075デフォルトの名無しさん
NGNG時間軸にある相関関数の数
0076デフォルトの名無しさん
NGNGMatlab一般って本はあまりないと思うぞ。
坂和先生の本とかは一般論に近いけど、僕はそんなに良いとは思わない。
分野を限定したら色々ある。
大学とかだと先輩が書いたマニュアルが便利なことが多いと思う。
007773
NGNGレスさんくすこ
相関関数の数ってデータ系列から判断できるものなの？
勉強不足か…逝ってきます orz
0078デフォルトの名無しさん
NGNG付属のマニュアルとオンラインヘルプとデモ
まじめに書いてあるし、下手な入門書よりわかりやすいと思う
0079デフォルトの名無しさん
NGNGこれはhist(x1),hist(x2)の後や中に、何のコマンドを書けばいいんでしょうか？
0080デフォルトの名無しさん
NGNGどうやるんでしょうか？
0081デフォルトの名無しさん
NGNG質問の意味がいまいちはっきりしないんだけど
x = rand(100,1) * 4 - 2;
y = x .^ 3;
っていうことなのかな。
0082デフォルトの名無しさん
NGNG同時に二つ書くというのはどういう状態がほしいんですか？
色違いの棒グラフが並んでいるようなイメージですか？
それだったらhistじゃなくて以下のような感じでもいいかも。
x1 = randn(1000, 1);
x2 = randn(1000, 1) * 2 + 1.5;
edges = -4:0.5:4;
h1 = histc(x1, edges);
h2 = histc(x2, edges);
bar(edges, [h1 h2]);
0083デフォルトの名無しさん
NGNG同一のグラフに描くをいうのがありますね。
それをx^3でやれというんです。(ちょっと違いますが)
乱数なら簡単なんですが、任意の関数からサンプルを取るとなるとわかんなくて。
0084デフォルトの名無しさん
NGNG>>80
の方です。
0085デフォルトの名無しさん
NGNG参考になりました。ありがとうございました。
重ねることは出来ないんですね。
0086デフォルトの名無しさん
NGNG重ねようと思えば
hist(x1)
hold on
hist(x2)
でできますが、意図した結果にならないかもしれません。
0087デフォルトの名無しさん
NGNGやっぱりいまいち分からんのですが、こういうことですか？
x1 = linspace(-2, +2, 100);
y1 = x1 .^ 3;
[n, center] = hist(y1);
plot(center, n, 'k--');
hold on;
x2 = rand(100, 1) * 4 - 2;
y2 = x2 .^ 3;
hist(y2);
0088デフォルトの名無しさん
NGNG乱数を使えばいいんですね～。
TAに質問したら、ｙ＝e(-ｘ)の場合を考えてみろとかかれました。
e(-x)の関数の分布の形に乱数を発生させるのってどうすればいいんですか？
どう乱数を変換すればいいんですかね？
0089デフォルトの名無しさん
NGNGxに均一な乱数を生成して、それをexp(-x)に入れれば、
y=exp(-x)からの出力yはその分布に従うと思うのですが。
x = rand(10000, 1) * 3; % [0, 3]の一様乱数生成
y = exp(-x);
hist(y, 100); % 分布を表示
>>86の後半でやっていることも同じですね。
0090迷える大学生
NGNG大学でscilabを使う問題が出されたのですがまったく無知ゆえさっぱりです
どうか住人の方のお力をお借りしたいです
aシューティングメソッドを使って解け
b有限差分法を使って解け
y''=-2y'-4y
y(0)=2
y(1)=2
0091デフォルトの名無しさん
NGNG授業で出された問題だったら「まったくの無知」というのもおかしな話なので、
大学の先生に質問してみましょう。
0092デフォルトの名無しさん
NGNG宿題スレとちゃうって。
0093デフォルトの名無しさん
NGNG0095デフォルトの名無しさん
NGNGしかもスレ違いなんだし。
0096デフォルトの名無しさん
NGNG下のMATLABの簡単なプログラムをOCTAVEで走らすのに苦闘中です。
function bdot = heatup(t,b)
global a h q w;
bdot = (q-h*(b-a))/w;
endfunction
global a h q w;
a = 18; h = 50; q =1000; w=10;
[time,bt] = ode23('heatup',0,1,30)
plot(tme,bt)
これを下の様にしたんですがグラフが見事にちがうんです。どこがいけないですか?
function bdot = heatup(t,b)
global a h q w;
bdot = (q-h*(b-a))/w;
endfunction
global a h q w;
a = 18; h = 50; q =1000; w=10;
t= linspace(0,1,100);
tb = lsode('heatup',[0],t)
plot(tb)
0097デフォルトの名無しさん
NGNGなんでだめなんすかね?
0098matlab help!!
NGNG可能な限りのsquare waveを維持しながら、
可能な限りのノイズを吸収するするフィルターを作ろうとしてます。
バターワースもしくはChebyshev1 のフィルターのハイパス、ローパス
バンドパスかnotchを使ってフィルターを作りたいんですが、下のコードの
A,Bに入れるcodeがわかりません。 どなたかわかりますか？
0099matlab help!!
NGNGt = 0:0.00001:0.005; % Generate time index.
x = 100*square(2*pi*1000*t); % Generate square wave.
var(x) % Determine the variance of x.
pause
plot(t, x), axis([0 0.005 -150 150]); % Axis command means to plot the
% signal from 0 to 0.005 for t axis
% and from -150 to 150 for x axis.
noise = 100.0*randn(1, length(t)); % Generate gaussian (normal) noise
% signal.
var(noise) % Determine the variance of the noise.
pause
plot(t, noise), axis([0 0.005 -300 300]); % Plotting noise
pause
xn = x + noise; % Add square wave and noise.
var(xn) % Determine the variance of xn.
pause
plot( t, xn); % Plot square wave plus noise signal.
y1 = filter(B, A, xn); % Filter the noisy signal. You get B & A
% from your filter design using butter
% or cheby1, or?
plot(t, y1); % Plot filtered noisy signal.
pause
plot(t, x, t, y1); % Plot noiseless signal and filtered
% noisy signal together.
----------------------------------------------------------------------------------
0100matlab help!!
NGNG---------------------------------------------------------------------------------
t = 0:0.00001:0.005;
x = 100*square(2*pi*1000*t);
var(x)
pause
plot(t, x), axis([0 0.005 -150 150])
noise = 100.0*randn(1, length(t));
var(noise)
pause
plot(t, noise), axis([0 0.005 -300 300]);
pause
xn = x + noise;
var(xn)
pause
plot( t, xn);
y1 = filter(B, A, xn);
plot(t, y1);
plot(t, x, t, y1);
0101choix
NGNG奇数倍の周波数に倍音成分が入っているので、
単純なbutterやchebyフィルタじゃ取り出せないと思います。
各倍音成分を取り出せるようなフィルタを考えてみましょう。
0102matlab help!!
NGNGbutterworh/chebyshefフィルタしか使えるのが、ないんですよね。
(まだその２つしか習ってない・・。）
例えば、square wave+noise signalっていうのは、butterやchebyでは
ちょっとのノイズでも取り出せないんですかね？
0103choix
NGNG>>101 で書いた「各倍音成分を取り出す」ということを
Butterworthを使って実装してみました。
もっといい方法はあるんでしょうが、
自分は信号処理の専門家じゃないので
この程度しか思いつきませんでした・・・。
fs = 100000;
t = 0 : 1/fs : 0.005;
x = 100 * square(2*pi * 1000 * t);
noise = 100.0 * randn(1, length(t));
xn = x + noise;
width = 20;
y = zeros(size(xn));
for n=1:20
fc = (2*n-1)*1000;
[B, A] = butter(1, [fc-width fc+width] / (fs/2));
y = y + filter(B, A, xn)/(2*n-1);
end
y = y/max(abs(y)) * 100;
plot(t, xn, 'r', t, y, 'k');
0104matlab help!!
NGNGちょっと見たことがないコマンドがいくつか
ありますけど、打ち込んだらしっかりとしたフィルターグラフが出てきました。
すいません、ここでまたまた質問なんですが、butterworthを使った場合
1st order butter worthでないと、いけないんでしょうか？
0105matlab help!!
NGNGみてください。
-------------------------------------------------------------------------------------
t = 0:0.00001:0.005;
x = 100*square(2*pi*1000*t);
var(x) % Determine the variance of x.
plot(t, x), axis([0 0.005 -150 150]);
xlabel('Period(Second)'), ylabel('Amplitude')
title('Noiseless Signal');
noise = 100.0*randn(1, length(t));
var(noise)
plot(t, noise), axis([0 0.001*5 -300 300]);
xlabel('Period(Second)'), ylabel('Amplitude')
title('Noise');
pause
xn = x + noise;
var(xn)
plot( t, xn);
xlabel('Period(Second)'), ylabel('Amplitude')
title('Noised Signal');
pause
Wp = (5/500); Ws = (250*0.707/500); Gp=-2; Gs=20;
[n,Wn] = buttord(Wp,Ws,-Gp,-Gs);
[B, A]=butter(n,Wn);
y1 = filter(B, A,xn);
plot(t, y1)
xlabel('Period(Second)'), ylabel('Amplitude')
title('Filtered Noisy Signal');
plot(t, x, t, y1);
xlabel('Period(Second)'), ylabel('Amplitude')
title('Noiseless Signal and Filtered Noisy Signal');
0106デフォルトの名無しさん
NGNG次数を求める故の阻止域端減衰量、通過域端減衰量、阻止域エッジ周波数の求め方がわかりません。
これは、ノイズのグラフから求められるものなのですか？それとも、何か
公式があるのでしょうか？
ちなみに、フィルター使用としてるノイズは１kHzのsquare waveにSNR値を0dBに
したノイズです。
0107choix
NGNGいくつか変えて実行してみて、
それの出力が良さそうだったので選んだまでです。
buttordで次数を決めるのはいいかもしれませんね。
buttordに指定する係数は自分がどういうフィルタを作りたいかに関係するので、
この場合はこうしろ、という公式があるわけではなさそうです。
いくつか試して、最適なものを選ぶと良いでしょう。
使ったことなければMatlab Signal Processing Toolboxの FDATool というのも
グラフィカルにフィルタが作れるので、
フィルタデザインを理解するのに良いかもしれません。
0108デフォルトの名無しさん
NGNG0109choix
NGNGノイズの種類などが分かっていれば除去できるアルゴリズムを見たことがあります。
もちろん、完全に除去できるわけではありません。
たぶん音声信号でも同じでしょう。
多くの場合は「耳で聞いた限りはノイズがない状態」に持っていくことはできても、
「ノイズが加えられる前の状態」にするのは相当難しい（あるいは不可能）だと思います。
0110デフォルトの名無しさん
NGNG0111デフォルトの名無しさん
NGNGfpritfで保存すると参考書に書いてあるんですが，
どのフォルダのどこに保存されているかも解りません．
cやFortranのように，データファイルの保存は出来ないのでしょうか！？
0112111
NGNG保存されている場所が解りました．workフォルダだった．
すみませんでした．
0113デフォルトの名無しさん
NGNGfopenのときフルパス出開くこともできるんでないか？
0114111
NGNGレス，ありがとうございました．
あと，自分のMatlabでは，
randint()やmodmap()などの関数が認識していません．
どうやったら，使えるようになりますか！？
0115111
NGNGcommunications Toolboxが入っていないからかな？？
誰か教えてください．
0116choix ◆pnCm4bpG3Y
NGNG一番簡単な方法は「save -ascii」だと思われます。
0117choix ◆pnCm4bpG3Y
NGNGfid = fopen([filename, 'w');
format = '%9.4f';
for row=1:size(data, 1)
for column=1:size(data, 2)
fprintf(fid, format, data(row, column));
end
fprintf(fid, '¥n');
end
fclose(fid);
0118デフォルトの名無しさん
NGNG選ぶしかないんですか?
0119choix ◆pnCm4bpG3Y
NGNG特にこれといった要求がなければ1でいいんじゃないでしょうか。
もしそれが学校の宿題の問題なんだったら、
先生に聞くのが一番手っ取り早いですよ。
0120デフォルトの名無しさん
NGNGraindintもmodmapもCommunications Toolboxだね。
入ってないなら使えない。
どうしても使いたいなら買ってください
0121デフォルトの名無しさん
NGNG0122choix ◆pnCm4bpG3Y
NGNG完全互換にするのは面倒ですが、
似たものを自作という手もありますよ。
たとえばrandint()なら
round(rand*m)
なんていうぐあいに。
0123111
NGNG認識出来ませんでした。
研究室から全ツールが入ってものを入れて、解決しました。
みなさん、色々なアドバイスありがとうございました。
そこで、もうひとつ質問があります。
ある参考書にQPSK信号を作るのに、
% 信号作成
Pd = 2000; % Pd : 送信されたシンボル数
Fd =1;% Fd：シンボルレート；１Hzに正規化
Fs = 4*Fd; % Fs : シンボルあたりの標本数
R = 0.5;% R : roll-off factor
Delay = 5;% Delay : フィルタの遅延時間 [symbol]
No = 1;% No : 雑音の電力密度スペクトル [W/Hz]
M=4;% M：多値数
% 信号系列発生(M値（０，１，２，３）の乱数発生)
x1 = randint(Pd,1,M);
% 信号配置変調のためのマッピング
y1 = modmap(x1,Fd,Fs,'qask',M);
%フィルタリング(raised cosine filter)
[rcv_a1,ti] = rcosflt(y1, Fd, Fs, 'fir/sqrt/Fs', R, Delay);
% 変調波信号の発生複素振幅変調
s1 = amodce(rcv_a1,Fs,'qam');
のように作ってありました。
でも一番下の行のamodce(rcv_a1,Fs,'qam')はQAM信号を作るソースですよね！？
上のプログラミングソースはQAM信号、QPSK信号どちらを作るものですか？
分かる方がいらしたら、是非教えてください！
お願いします
0124デフォルトの名無しさん
NGNGいいなぁ。ウラヤマスィ。
0125デフォルトの名無しさん
NGNG0126デフォルトの名無しさん
NGNGどうせ研究室外に出たら使えないし、ライセンス数も同時使用数の形でサーバー側で管理してるんだから。
それか、自分のっていっても占有できるだけで、研究室の備品だとか。
ﾏｲｸﾛｿﾌﾄのキャンバスアグリーメントは契約次第で
学生や教職員の個人のＰＣにインストールしてもＯＫってのがあるらしい
0127デフォルトの名無しさん
NGNG0128デフォルトの名無しさん
NGNGネットワークライセンスでも…技…にはでき…なんだな。
実際の使用条件ってどうだっけ？オレ、鯖のセットアップは自分でやった
けど研究室でしか使わないからその辺りを意識したことなかったな。
＞学生や教職員の個人のＰＣにインストールしてもＯＫってのがあるらしい
うちの学科はそれやってる。Officeはついてこないけどね。
0129デフォルトの名無しさん
NGNGＣで計算し、出力したプログラム（txt)をmatlabを用いて
グラフ化するやり方がわかりません。
出力したプログラムは
時間Ｔに関して、それぞれx y zの三つの値がでて、それが５千行ある結果なのですが。。
matlabにうつコマンドは
plot3(ファイル名(:,2),ファイル名(:,3),ファイル名(:,4))
だとおもうのですが、、
ちなみに、matlab5.3を使っています。これは９８年くらいで
結構ふるいのですが、
よろしければおしえてください。
0130デフォルトの名無しさん
NGNGＣで計算し、出力した結果（txt)をmatlabを用いて
です。
0131デフォルトの名無しさん
NGNGp = load('filename.txt');
plot3(p(:,2),p(:,3),p(:,4));
でいいんじゃないでしょうか。(時間軸を無視しちゃいますけどね)
0132デフォルトの名無しさん
NGNGれすありがとうございます。
まだうまくできないのですが。。。
ほんとうに初心者ですいません。
》 p = load(''filename.txt');
??? エラー: ==> load
result5.txt: Can't open file.
0133デフォルトの名無しさん
NGNG>> help fopen
>> help fread
0134デフォルトの名無しさん
NGNG0135デフォルトの名無しさん
NGNG0136choix ◆pnCm4bpG3Y
NGNG使ってもいいんじゃないでしょうか？
0137デフォルトの名無しさん
NGNG0138choix ◆pnCm4bpG3Y
NGNG了解しました。
0139デフォルトの名無しさん
NGNGそれを教授に言ったら、「複素数なんて、実際の通信ではありえない」と
怒られました。
複素数の信号から実際の使われている値に変換することって出来ませんか？
論文の提出期限まで時間がないので、かなりビンチですぅ
0140デフォルトの名無しさん
NGNGｗｗﾍ√ﾚｖｖ～（ﾟ∀ﾟ）─ｗｗﾍ√ﾚｖｖ～─ !!!
0141デフォルトの名無しさん
NGNGbitmap 2階調でデータを保存したいのですが、
imwrite( data, gray, 'savename', 'BMP' ); %%（データは0~64の整数値）
とすれば、一応RGBで保存してくれるのですが・・・・。
白黒で良いので0or1の2階調で保存したいのです。
いいアイデアあれば教えてください。
>>139
時間領域の信号に複素数はありえないよ。
周波数領域なんでしょ、FM変調なんだし。
教授が怒るのも分かる気がする。
0142choix ◆pnCm4bpG3Y
NGNG複素数から実数にするには一般的には
「実数成分だけ取ってくる real(x)」とか
「振幅だけ取ってくる abs(x)」とかを使いますが、
それは「実際に使われている値」に変換する方法ではありません。
そもそもFM変調の仕方がおかしいような気がしますので、
その部分を見直した方が良いでしょう。
0143choix ◆pnCm4bpG3Y
NGNGbw = data > 128;
などとやってムリヤリ０と１に変換してしまうのが手っ取り早いかと思います。
0144choix ◆pnCm4bpG3Y
NGNG>>143 は >>141 への答えでした。
>>139 なんですが、Signal Processing Toolboxには以下の
FM変調コマンドがあります。
y = modulate(x, Fc, Fs, 'fm');
入力であるxが実数であれば出力yも実数値が出てくるはずです。
139さんの用途にあうかどうか分かりませんが、試してみてください。
0145choix ◆pnCm4bpG3Y
NGNGたびたびすみません。
X = dither(I);
というコマンドがImage Processing Toolboxにありました。
もし白黒2値化にディザが必要であればチェックしてみてください。
0146ken
NGNG全く二値化のあたりからわからないです。だれかアドバイスください。
0147139
NGNGレスありがとうございます。
今から、試してみます。
0148デフォルトの名無しさん
NGNG灰色にして、どこかで切ればいいんじゃないの？
0149139
NGNG今、教えられたようにy = modulate(x, Fc, Fs, 'fm',kf);
をhelpで調べ、使うと思っているんですが、
xには一般的ににどのような実数を入れれば良いのでしょうか？
試しに、ｘにsin関数の値を入れて実行しましたが、、、。
質問ばかりでごめんなさい。
でも、本当に感謝しています。
0150139
NGNG今、FM信号とAM信号を生成し、シミュレーションしました。
良いと思われる結果が出ました。
これからFM,AM信号の他に、PSK信号やFSK信号も作りたい思います。
また、何かありましたらよろしくお願いします！
0152choix ◆pnCm4bpG3Y
NGNG画像の二値化の簡単な方法は >>143 >>145 に。
>>148 の「どこかで切る」の切る場所が知りたければ
大津先生の方法とかが参考になるかも。
理解するには統計知識もちょっと必要になったりします。
0153デフォルトの名無しさん
NGNG2値化は既に出来ているんです。（ディザにする必要はないです）。
０と１の値からなる行列データを白黒(2bit)ビットマップとして保存したいのです。
iwriteでは無理なのでしょうか。HELPみてもそこいらへんについて書いていなかったので・・・。
なにか良い方法があるのでしょうか
0154デフォルトの名無しさん
NGNGaとnは任意変数で、
aのｎ 乗根を計算するコマンドがないのですか？
sqrt(n)しかしらないので。
0155デフォルトの名無しさん
NGNGa^n
0156154
NGNG0157デフォルトの名無しさん
NGNG0158154
NGNGどうもです、吊ってきます。
0159choix ◆pnCm4bpG3Y
NGNGimwriteに'BitDepth'パラメータを指定することで
1ビットの白黒データにできますよ。
もしくはTRUE/FALSEデータをimwriteに
突っ込めば勝手に白黒だと認識してくれます。
ただ、PNG、BMP、TIFFくらいしか1ビットには対応してない気がします。
0160choix ◆pnCm4bpG3Y
NGNGimwrite(IMG, 'hoge', 'png', 'BitDepth', 1);
です。
0161デフォルトの名無しさん
NGNGご丁寧に有難うございます。
試してみます。
0162デフォルトの名無しさん
NGNGシミュレーションのサンプルってありませんかね？
たとえば、１ｋHz正弦波をだした時の音圧の濃密がカラーで
わかるような。
0163デフォルトの名無しさん
NGNGやってみました。
imwrite( kami, savename, 'bmp', 'BitDepth', 1 );
これだと、
?? エラー: ==> c:/matlab6p5/toolbox/matlab/iofun/private/writebmp
入力引数が多すぎます
エラー: ==> C:\MATLAB6p5\toolbox\matlab\iofun\imwrite.m
行番号: 310 ==> feval(fmt_s.write, data, map, filename, paramPairs{:});
エラー: ==> C:\MATLAB6p5\work\func2bmp3.m
行番号: 67 ==>imwrite( kami, savename, 'bmp', 'BitDepth', 1 );
と怒られます。bmpをpngにするとうまくいきます。
うーん、bmpはだめなのか・・・。
ちなみにバージョンは6.5(R13)です。
ImageProcessingToolboxは持っていません。
SignalProcessingToolboxはあります。
0164デフォルトの名無しさん
NGNG256点の領域に60点放して等強度（強度１）のスペクトルを仮定し、（その他のデータ点は０）これを①、
これとは別に129点の真ん中にスプレッドの中心がくるような全面積１のラインスプレッド関数を
つくります。これを②
この二つをコンボリューションした後を③とします。
これをそれぞれフーリエ変換すると③＝②＊①が成り立つので、①＝③÷②にして逆フーリエ変換
したのですが、グラフが正しく出てこず、ためしに絶対値表示にしたらうまく元の①の図が出てきました。
絶対値表示にしなければいけなかった理由を教えてください。スペクトルは波の強さを表す必要が
あるからですか？？
よろしくお願いします。
0165デフォルトの名無しさん
NGNGfftをN点で行ったら値はN倍されて算出されるので、N割りしないといけないから。
0166デフォルトの名無しさん
NGNGありがとうございます。
でも、勉強不測のためわからんかったっす。
高速フーリエ変換と離散的フーリエ変換を勉強すれば上記のことはわかりますか？
0167デフォルトの名無しさん
NGNG分かると思います。
2次元信号ではなく1次元信号で確かめてみてはいかが？
簡単にN倍されていることが分かるはず。
0168choix ◆pnCm4bpG3Y
NGNGPNGでは試したんですが、BMPでは試してなかったのがバレバレですね。
すみませんでした。
>>159 で書いた別の方法（TRUE/FALSEを使う）だとうまくいきました。
以下の文ではlogical()を使って0と1だけが含まれているAをT/Fに変換しています。
imwrite(logical(A), 'test.bmp', 'bmp');
0169139
NGNGを使い、4-PSK信号を作るにはどうした良いですか！？
0170choix ◆pnCm4bpG3Y
NGNGdmodceを使うのがラクでしょう。
0171choix ◆pnCm4bpG3Y
NGNG用途に合わせてどうぞ。
0172ken
NGNG0173デフォルトの名無しさん
NGNGカネはいくら出すんだ？
0174139
NGNGどうもです！
さっそく，自宅に帰り次第試してみます．
0175ビギナー
NGNG例えば十進数を二進数に直すときは
clear;close all
input1=[
1
2
3
4
5
6
7
];
target1=[
0 0 1
0 1 0
0 1 1
1 0 0
1 0 1
1 1 0
1 1 1
];
p=input1';
t=target1';
net=newff(minmax(p),[5,3],{'logsig','purelin'});
net.trainParam.epochs=500;
net.trainParam.lr=0.1;
net=train(net,p,t);
a=sim(net,p)
といった感じなんですけれど、こんな感じのをベースにして、足算をさせるにはどうするのでしょうか？
どうしても分からないので、教えていただけたら幸いです。
0176デフォルトの名無しさん
NGNG3=2+1=(1+1)+1
4=3+1=((1+1)+1)+1
0177デフォルトの名無しさん
NGNGっていうんだったらこのスレでいいけど、理論のことはスレ違い。
卒論や宿題のお手伝いスレじゃないぞ。
0178デフォルトの名無しさん
NGNG返事が遅くなりすみません。
早速、行列をlogicalに変換してやってみます。
0179デフォルトの名無しさん
NGNG出来ました。
有難うございました。
0180デフォルトの名無しさん
NGNG具体的に言うと3次元でのグラフがz軸で対象なためとりあえず2次元で解いて
それをz軸において回転させて三次元のプロットにできたらいいなと思ってます。
世の中そんなに甘くないのでしょうか?；；
0181180
NGNG0182デフォルトの名無しさん
NGNG0183デフォルトの名無しさん
NGNGhelp mesh
0184180
NGNG>183 meshじゃちょっと違うんです。X(i,j) とZ(i,j)があってそれをstreamsliceで
vector fieldでplotしてるんですよ。今は二次元なんですが、matrixを三次元にする
というよりは二次元(i,j)を回転させて三次元のplotにしたいだけで。。
いろいろやってみましたが唯一できるかなって思ったのはX(i,j)とZ(i,j)をcylindrical
coordinateと考えてthetaをいれてそれをpol2cartでかーてじあんに変換してplotかなって。。
どうでしょうか？わかりずらい文章ですいません。。
0185デフォルトの名無しさん
NGNG古いバージョンのマニュアルに出ていた例ですが、z軸からの距離をrで表し、
二次元の関数をZ=f(r)=sin(r)/rとします。
x, yともに-8から+8までの範囲を描くとすれば
[X,Y] = meshgrid(-8:.5:8);
R = sqrt(X.^2 + Y.^2) + eps;
Z = sin(R)./R;
mesh(Z);
これでf(r)=sin(r)/rをz軸中心に回転させた図形が描けます。
0186180
NGNGなるほど、その方法参考にしてどうにか当てはめてみようと考えましたが、
私の場合はFx=f(x,y,z)、Fy=f(x,y,z)、Fz=f(x,y,z)のvector fieldをplotしたいんです；；
それでもってFx=FyなのでFx=Fr=f(x,z)、Fz=f(x,z)で解いてz軸に回転させたvectorfieldを
plotできるかなーって思ってました。
185さんの意見を参考に考えぬいてFrにcos(theta)をかけてFxとし、sin(theta)をかけて
Fyとしておいてみました。これでFx= f(x,z,theta)、Fｙ= f(x,z,theta)、Fz= f(x,z,theta)と
一応三次元にしてみたものの、、自分でよくわからなくなりました。(汗
もうちょっと考えて見ます。
アドバイスありがとうございます。
0187デフォルトの名無しさん
NGNGどなたかお願いします！
例えば1～5の数があるとして、
1のでる確率が0.3、5のでる確率が0.1
といった感じにしたいのですが…。
詳しい方、アルゴリズムだけでもいいのでお願いします。
0188choix ◆pnCm4bpG3Y
NGNGその簡単な例を実現するだけであれば
a = rand(1,1);
switch(a)
case a<0.3
%確率0.3 [0.0, 0.3)
b = 1;
case a<0.4 [0.3, 0.4)
%確率0.1
b = 2;
case a<1.0 [0.4, 1.0)
%確率0.6
b = 3;
end
みたいにできますが、もっと複雑な乱数生成がしたいのであればStatistics Toolboxに含まれているrandom関数を調べてみましょう。
0189choix ◆pnCm4bpG3Y
NGNG0190デフォルトの名無しさん
NGNGもう少し複雑に乱数を生成したいのでrandom関数をいじってみます。
0191デフォルトの名無しさん
NGNGmatlab用の実験toolboxを使って心理実験等を作って、それをWeb上で動かして
より多くの人のデータを取りたいと思っているのですが、可能でしょうか？
基本的な作業の流れ等、教えていただければうれしいです。
0192デフォルトの名無しさん
NGNGCGIは書けますか？できればC言語で。
それができるならばUNIX用 or Linux用のMatlabでは可能だと思います。
インターフェイス部分はWebでは利用できないでしょうが、エンジン部分なら
使えそう。
＞matlab用の実験toolboxを使って心理実験等を作って
ここを詳しく書かないと有用な情報を得るのは難しいでしょ…
0193choix ◆pnCm4bpG3Y
NGNGでも、WebであればCGIをPHPやPerlで作ってしまって、集めたデータ処理だけをMatlabで行う方が簡単だと思いますよ。
どうしてもWeb上でGUIを使いたいのであれば、JavaアプレットかFlashを使わないといけませんね。
0194デフォルトの名無しさん
NGNG> MATLABでFM信号を作ると、複素数の値で出てきます。
> それを教授に言ったら、「複素数なんて、実際の通信ではありえない」と
> 怒られました。
> 複素数の信号から実際の使われている値に変換することって出来ませんか？
> 論文の提出期限まで時間がないので、かなりビンチですぅ
いいよ、そのまま論文出しちまえ。
そんなことを言う教授は戒めのつもりで言っているつもりなんだろうが
今の時代にその教授の学生に対するやり方、態度は時代遅れで通用しない。
「複素数なんて、実際の通信ではありえない」
なんて言ってきたら
「だから何だってんだ？ あんたは何が言いたいんだ？全然意味がわからん。
先生が何を言っているのか先生が何を考えているのか全くわからん。
先生が何のためにそんなことを言っているのかわからん。」
と言ってやれ。教官によってはかなりへこむぞ。
わからねえならとりあえずabs()で絶対値をとっちまえ。
それでも教授が文句をいうなら
「先生の指導能力が悪く、先生の意図していること、
目的がわからず時間がないのでこうさせてもらいました。」
とでも言っておけ。
0195デフォルトの名無しさん
NGNG> 自分は音響心理実験をMatlabのGUIでやってますが、実際にPCの前でヘッドホンをしてもらう、というのが中心で、Web版の経験はないです。
MatLabでいちいちGUIとはおめでたい。
> でも、WebであればCGIをPHPやPerlで作ってしまって、集めたデータ処理だけをMatlabで行う方が簡単だと思いますよ。
> どうしてもWeb上でGUIを使いたいのであれば、JavaアプレットかFlashを使わないといけませんね。
いまどきJava Appletなんかつかわねえよ。
これからはリッチクライアントたる Java Web Startの時代だ。
そもそも、CGIやPHPといったサーバサイドを意識しているなら
そこでクライアントサイド専用のJava AppletやFlashを持ち出すのは
何か不釣り合いだ。Servlet/JSPを忘れてないか？
0196デフォルトの名無しさん
NGNG> >>191
> CGIは書けますか？できればC言語で。
なぜのろまなCGIや時代遅れなC言語に拘るのか理解できんな。
CではなくPerlやPHPやServlet使えば何だって可能だぞ。
MatLabの出力結果をサーバサイド言語に入力して
Webに出力すればいいだけなんだからよ。
そもそも>>191が目指しているWeb上で動かすものとやらが
どの程度もものなんだかな。
サーバ側にあるプログラムをクライアント側に
ダウンロードしてから実行する程度のものならCGIもServletも不要なんだがな。
HTMLとあまりかわらんな。そこでAppletやらFlashやらJavaWebStartやら
ﾄﾞﾄﾈﾄのノータッチデプロイメントやらが出てくるわけだが。
サーバ側で処理してから結果をクライアント側に返す程度だけの場合も
いまいちだな。
サーバ側で計算しクライアント側に応答を送るだけでなく
クライアントからサーバにリクエストを送信する場合のときのみ
CGIやPHP, JavaServletやJava RMI
といったサーバサイドテクノロジが必要だ。
サーバ側にServlet、クライアント側にAppletという組み合わせで動く
システムもあるがな。
Flashをサーバサイドと連携する場合はFlash RemotingMXなんてものを
使う手もある。ColdFusionというものも使う手もありだ。
0197choix ◆pnCm4bpG3Y
NGNG質問者のことを理解していない低俗な煽りですが、とりあえずマジレスしておきます。
>MatLabでいちいちGUIとはおめでたい。
Matlabを使えば複数言語を覚える必要なく、実験からデータ処理まで一貫して行えるので便利です。
特に心理実験では一回実験をやったら用無しになるものが多く、販売用アプリケーションを作るわけでもないので、思ったものが素早く作れるMatlabが自分にとっては最適です。
>いまどきJava Appletなんかつかわねえよ。
>これからはリッチクライアントたる Java Web Startの時代だ。
JWSも十分理解していますよ。でもアプレットの方が現在でも情報が入手しやすいのでおすすめしただけです。
>そもそも、CGIやPHPといったサーバサイドを意識しているなら
>そこでクライアントサイド専用のJava AppletやFlashを持ち出すのは
>何か不釣り合いだ。Servlet/JSPを忘れてないか？
サーバーだのクライアントだの技術面で切り分けるのではなく、実験者と被験者にとって最適なものを提案したいだけです。
Java系のSEをやっていたこともありますので、そういった技術も十分に承知していますが、
ServletやJSPだとリッチクライアントがまだ作れない上に、プログラムもCGIに比べて大規模になりがちです。
また、理解すべきことが多すぎますので、プログラミング初心者だという>191さんには向かないと思いました。
0198デフォルトの名無しさん
NGNGCGIが言語だと勘違いしていないか？
いまどきServletやJSPが難しいって、いつの時代の人間なんだが。
本当にリッチクライアントの意味わかって言ってるのか？
StrutsもJSFもTapestryもVelocityもXDocletも知らないのか。
これらを使いこなせば大規模に見えるJava開発もCGIよりも簡単に
なるわけだが。
0199デフォルトの名無しさん
NGNG質問者のことを理解していないという以前に
質問者の質問がアホ過ぎるほうが低俗。
0200choix ◆pnCm4bpG3Y
NGNG自分はJava言語が一番のお気に入りですし、Webサービス系の仕事をしていました。
StrutsなどのMVCフレームワークやO/RマッピングのためのフレームワークもSEをやっていたときに使いましたし、様々な技術に関してフィージビリティ・スタディのための調査もしました。
その経験からもServletやJSPが初心者にとって難しいのは事実だと思っています。
少なくとも質問者は、Webプログラムの書き方をマスターしたいのではなく、効率のよい実験をしたいだけのように思えますし、
Javaのサーバーサイド技術に関してはMatlabスレで議論する内容じゃないでしょう？ 他のところでやりませんか？
0201デフォルトの名無しさん
NGNGやりたい実験の内容を書いていないとしょうがなかったですね。
とりあえずMatlabに慣れようと思って、簡単な刺激提示（文字／絵）→キーを押すまでの
反応時間の計測という程度の実験（ストループ実験のような）を作ってみようかなと思った次第です。
そこで、それがなるべくウェブを通じて多くの人にやってもらえる方法はないかなと、
いう質問でした。
皆さんにとっては大分しようもない質問だと思いますが、なにぶん私が全くの門外漢なもので、宜しくお願いします。
0202デフォルトの名無しさん
NGNGとにかく藻前はCGIが言語だと勘違いしている時点で(ry
0203デフォルトの名無しさん
NGNGただの体感ソフトを作りたいだけなら
CGIもServletもいらんな。
Webを通して他人に使って貰いさらに
そこからデータを取得したいなら
多少のサーバサイド技術を使うが。
0204choix ◆pnCm4bpG3Y
NGNGMatlabを使うのであれば、オフラインで（被験者にPCの前に座ってもらって）行うようなプログラムがいいと思います。
また、オンラインで行いたいのであれば、Java Applet、Java Web Start、Flashなど、リアルタイム性のあるクライアントで動作するものがいいでしょう。計測結果はプログラムがメールで発送するようにすれば、サーバーを用意する必要もありません。
どうしても何らかのかたちでMatlabを使い、かつオンラインで行いたい場合、Matlab Web Server （http://www.mathworks.com/products/webserver/
） を導入し、リアルタイム性のあるクライアント・ソフトと組み合わせるという手もありますが、むりやり難しい解決法を選ぶ必要もないでしょう。
0205choix ◆pnCm4bpG3Y
NGNGCGIをプログラミング言語だなんて言った覚えはありません。
そう読み取れる文章を書いたこともないはずです。
ツッコミを入れる時はちゃんと文章を読んでからにしてもらいたいです。
0206デフォルトの名無しさん
NGNGところでプラットフォームは何？
0207デフォルトの名無しさん
NGNG0208デフォルトの名無しさん
NGNGだったら誤解されないよう文章を書くこった。
誰もが自分と同じ考え方を持っていると勘違いしないようになw
0209デフォルトの名無しさん
NGNG> 計測結果はプログラムがメールで発送するようにすれば、サーバーを用意する必要もありません。
sendmailはどうやって動かすんだい？
0210デフォルトの名無しさん
NGNGもう一度自分の文を読み返してみ。
> でも、WebであればCGIをPHPやPerlで作ってしまって
これはチミがPHPのことをろくにわかってない証拠になるものですな。
これで「自分の思い通り読みとれ」などとよく言えますな。
こういうのを厚顔無恥というのですよ。
0211デフォルトの名無しさん
NGNGSEらしくない発言が目立つのだが。
顧客と交渉する能力もたいしたことなさそうだし
0212choix ◆pnCm4bpG3Y
NGNG仕事で使ったのはJavaのサーバーおよびクライアントが中心で、PHPでは小規模システムしか作ったことはないですし、何十年もSEをやったわけでもありませんから。
でも、この場を荒らすようなことはしたくないというのを理解してもらえませんか？ さっきも書きましたが、Javaの話はスレ違いですし。
0213デフォルトの名無しさん
NGNG多分、自分が偉いと勘違いしてるおぼっちゃま君なんだろう。
0214デフォルトの名無しさん
NGNG俺はえらぞ
0215デフォルトの名無しさん
NGNGじゃあPHPの話もスレ違いだな坊主
0216デフォルトの名無しさん
NGNGなるほど、よくわかりました。
そうですね、Matlab練習用にと考えていたので、とりあえずオフラインの実験で
やってみようかと思います。Javaの方も少し考えてみます。
アドバイスありがとうございました。
0217デフォルトの名無しさん
NGNGMATLABは金がかかるからなあ。
気をつけることだよ。
いろいろやりたいことがあるなら
自前でサーバ構築くらいできるようになっとけ。
0218デフォルトの名無しさん
NGNG0219デフォルトの名無しさん
NGNGMが作ったあの穴だらけの古いJavaでさ。
0220choix ◆pnCm4bpG3Y
NGNGそれとは違うJava VM使ってますね。
0221助けてください
NGNGしているのですが何かお勧めのコマンドとかありますか？教えてください(＞_＜)
それとMATLABコマンドのODEでのメインファイルとサブファイル（関数ファイル）
のやり取りは具体的にどのように計算されているのかについても教えていただきたいです。
お願いします。
0222choix ◆pnCm4bpG3Y
NGNG申し訳ないのですが、こっちから質問させてください。
微分方程式の数値解析的な計算方法が知りたいのでしょうか？
それとも対象関数との間の値の受け渡し方法が知りたいのでしょうか？
＃スライディングモードに関しては無知なので
＃力になれそうにありません。
0223助けてください
NGNGode45を使っているのですが、これを使うときにメインファイルとサブ関数ファイル
を使う必要があり、実際にシミュレーションできたのですが、この2つのファイルが
どのようにして受け渡しされているか教えてください。お願いします。
0224choix ◆pnCm4bpG3Y
NGNG関数名の先頭に@マークをつけて引き渡されるので、ハンドラとして扱われます。
ode45内ではfevalでその関数が処理されます。
fevalはハンドラとして与えられた関数を実行するための命令です。
一言で言えば、ハンドラを使って関数そのものを受け渡し、fevalを使って呼び出す、ということをしているわけです。
0225助けて
NGNGがエラーが出ます。使っているMATLABがVer5.1なので…。
私が勘違いしているのかもしれません。。
ode45というコマンドがこれらの行為を自動で行っているのですか？
0226choix ◆pnCm4bpG3Y
NGNG自動で行なわれる、という言い方でいいと思います。
ただ、ODEの計算可能性は定義域や初期値などに左右されますし、
ode45（ルンゲ・クッタ法）で解けない微分方程式の種類もあります。
0227助けて
NGNG微分方程式で初期値や許容範囲も大丈夫です。
ode45の内部で関数がfevalによって処理されていることを確認するためには
どうすればいいでしょうか？ ちなみにMATLABはVer6.5にしました。
0228choix ◆pnCm4bpG3Y
NGNG自分で定義した関数の中にdisp命令を入れておく
というのがあります。
function dydt = myfunc(t,y)
dydt = hogehogehoge.....;
disp("myfunc executed");
disp(["y = " num2str(y)]);
disp(["t = " num2str(t)]);
disp(["dydt = " num2str(dydt)]);
のようなかんじにすると引数や戻り値なども表示できます。
0229助けて
NGNG時間を0秒から3秒で0.1秒刻みでシミュレーションしたのですが、関数ファイル
がode45に読み込まれた回数が半端じゃなく多いのはどうしてなんでしょうか？
いろいろと聞いてしまってすみません。。
0230choix ◆pnCm4bpG3Y
NGNG何度もfeval関数を使って違ったtとyの組み合わせを計算させているのです。
ここらへんはode45のもとになっているルンゲ・クッタ法の
教科書を読み返していただけると理解が早いと思います。
学部生のときに課題でやらされましたが、
手で計算してみるとかなり大変なのが身にしみて分かります ;-)
0231助けて
NGNG[t,y]=ode45('関数名',tspan,初期X) でode45が何度もfeval関数を使い異なる組み合わせ値を[t,y]に返しているということですね。ode45で実際は凄い小さい時間刻みで計算
されているけど、例えばtの刻み幅を0.1に指定してt=[0,0.1,0.2,0.3,…]とするときにはそれぞれの時間tのときの状態yを返しているのですよね？
t=0.1ではt=[0.0098,0.101,0.110,…]のような値(状態値)を近似しているのですか？
0232デフォルトの名無しさん
NGNGSimulink を使用しているのですが、途中でフレームサイズを変更することは
出来ますでしょうか？
モデルは以下のようになっています ( 簡略化してます )
Ramdom Signal Generator => Block1 => Scope
ここで、Random Signal はフレームサイズ 2048 で出力するのですが
Block1 はフレームサイズ 1024 しか受付ません。
そこで Block1 の前に何かブロックを挿入して
1. 2048 のフレームを 1024 にする
2. 1024 ずつ Block1 に渡す
3. Block1 の後で 1024 のフレームを結合して 2048 にする
ということをしたいと思っております。
こういう事はできますでしょうか？
よろしくお願いいたします。
0233デフォルトの名無しさん
NGNG自己レスです。DSP Blockset の buffer を使うことでフレームサイズの
分割と結合は出来ましたが、今度は必ず最初に値が 0 のフレームが負荷される
ようになってしまいました。
この 0 フレームが付かないようにするにはどうするのでしょうか？
0234choix ◆pnCm4bpG3Y
NGNG>されているけど、例えばtの刻み幅を0.1に指定してt=[0,0.1,0.2,0.3,…]とするときにはそれぞれの時間tのときの状態yを返しているのですよね？
そうですね。
>t=0.1ではt=[0.0098,0.101,0.110,…]のような値(状態値)を近似しているのですか？
数値解析的な解き方をしているので、ある程度の近似はしているでしょうが、ode45内部でどの状態値で近似しているのかはわかりません。
詳しく知るためにはソースコードを読むしかないでしょう。
0235デフォルトの名無しさん
NGNGどうにか実数で出力することは出来ませんか？
このままだと卒業できない・・・
0236助けて
05/01/17 17:40:53わかりました。本当にいろいろありがとうございました。
また何かあればよろしくお願いします。
0237デフォルトの名無しさん
05/01/18 17:49:21MATLABでUMLを使えるプラグインないかな。
0238デフォルトの名無しさん
05/01/20 11:31:290239助けてください・・・・
05/01/21 12:03:52とある大学生です。
自分は今卒論でＭＡＴＬＡＢ使ってるんですが、
いかんせんプログラムが苦手で困ってます。
どれくらい困ってるかと言うと、
就職決まってるのに卒研のせいでダブる危機があるくらい、
リアルにピンチです。
そこの教授も放置プレーなので、何も教えてくれません。
本当に聞いても自分で調べてくださいっと言われます。
締め切りまで2週間もないのですが、もし力を貸してくれるなら、
七世代まで神として祭り上げます。
と他力本願な俺ですね・・・叩いてもらってもいいです。すいません。
しかし、就職したいのです。だれか力を貸していただけませんか？
2-Ｄ画像に離散ウェーブレット変換をかけて例えばsをa1、d1に分解し、
さらにa1をa2,d2に分解、さらにa2をa3,d3・・・
といった具合に分解していって、sとa1～anをそれぞれ2乗平均して、
その値にLogをかけて縦軸をlog aとこ軸を2の階上として、グラフを作り、
その傾きを求めるというプログラムを作成しろとの事なのですが、
自分には何がなにやらさっぱりわかりません。
できればプログラム作っていただく、もしくは作り方を教えていただけませんでしょうか？？おねがいします。
0240デフォルトの名無しさん
05/01/21 12:11:310241デフォルトの名無しさん
05/01/21 12:14:45おい、何も教えてくれない先生なんておかしいだろ。
授業料返してもらえ。
0242デフォルトの名無しさん
05/01/21 14:58:20MATLAB使ってる院生に聞け。
0243デフォルトの名無しさん
05/01/21 15:15:46具体的にお題を与えられないとできないんじゃ、
いいとこ実験レベル。とても卒業研究とは言えまいに。
しかもそのお題の持つ意味すらわからないのでは……。
素直に留年しなさい。
0244デフォルトの名無しさん
05/01/21 15:40:250245choix ◆pnCm4bpG3Y
05/01/21 19:28:55卒業研究なんだから、それなりに基礎知識をつけているはずです。もしwaveletの「s」「a?」「d?」ですら意味が分からないのであれば、他の方々がおっしゃるように、もう一年しっかり勉強した方がいいでしょう。
後半部分に書かれているプログラムは、画像処理系の研究で10年ほど前から行なわれているのを見たことがありますし、IEEEをあたれば類似の手法はたくさん出てきますので、それらの論文も読んでみましょう。
0246デフォルトの名無しさん
05/01/21 20:38:530247デフォルトの名無しさん
05/01/21 23:28:53＞2-Ｄ画像に離散ウェーブレット変換をかけて（以下略）
すげえ。ここまで親切に教えてくれる教授なんて、信じられん。
ちょっと甘やかしすぎではないか？
0248デフォルトの名無しさん
05/01/22 03:20:20最近は卒研内容を「習っていません(だからできなくても問題ない)」なんて後輩もいるんだな．
なんで研究室配属を許可しちゃっただろ(というか演習や実験の単位をあげたんだろう)…
0249まと
05/01/23 13:30:38質問です。
ファイル名を絶対パスで取ってきているのですが、
￥が入ってなかなか使いづらいのですが、最後のファイル名だけ
取るのにはどうしたらよいでしょうか？？？
絶対パスでとってくるGUIを手に入れたので便利だなと
思っていたら落とし穴です。
それを使ってファイル名を作れなくなった次第です。
どなたか。お助けを。
0250デフォルトの名無しさん
05/01/23 23:00:28help fileparts fullfile
0251デフォルトの名無しさん
05/01/23 23:17:50やりたいことは，多数のファイルに対して一括処理を施したいだけなのですが．．．作った関数は下記の通りです．
1. DataLoad_Callback()：データ読み込み関数
2. DataProcessing_Callback()：データ処理関数
3. BatchProcess_Callback()：バッチ処理関数
3の関数をメインとして，1,2の順に呼び出して処理を行うものです．
1で読み込まれたデータは，handles構造体に格納して，2の処理関数へと渡されます．
しかし，1で格納したデータに2の関数からアクセスできません（2の関数に入った時点で，handles構造体に格納したデータは消えている）．
ソースは単純化していますが，次のようなものです．
＃改行が多いので，分割します．
なぜこのようなことが起こるのか教えて頂けないでしょうか？
0252デフォルトの名無しさん
05/01/23 23:19:14% データ読み込み
function DataLoad_Callback(hObject, eventdata, handles)
data = load('****.***');
handles.data = data;
guidata(hObject, handles);
% データの処理
function DataProcessing_Callback(hObject, eventdata, handles)
% この部分でhandles.dataは空となっておりエラーが発生する
data = handles.data;
処理部分（省略），結果をファイルに保存
% バッチ処理
function BatchProcess_Callback(hObject, eventdata, handles)
% こんな感じで複数のファイルを読み込み，処理をしたい
filename(1,:) = 'D:\work\dev\matlab\test1.mat';
:
filename(n,:) = 'D:\work\dev\matlab\testn.mat';
% 各ファイルに対して処理を行う
for iRow = 1:length(filename)
handles.filename = filename(iRow,:);
guidata(hObject, handles);
DataLoad_Callback(hObject, eventdata, handles)
DataProcessing_Callback(hObject, eventdata, handles)
end
0253デフォルトの名無しさん
05/01/24 16:48:210254デフォルトの名無しさん
05/01/27 01:05:410255253
05/01/27 03:59:12このスレにゃ答え出せる人居ないよ。自分で考えてください＞251
以下何事もなかった如く進行。
担当者と話をすると、MATLABの発音がマツゥラヴに聞こえるマットラボ？マトラボ？
0256デフォルトの名無しさん
05/01/27 07:41:27Waveletでずいぶん釣れたね！
0257choix ◆pnCm4bpG3Y
05/01/27 09:42:42GUIDEで生成されたひな形を使ってますか？
もしそうだとしたら、コールバック関数を呼ぶときには
hogehoge('DataLoad_Callback', hObject, eventData, handles, ...);
ってやることになっていると思うんですが、
BatchProcess_Callbackの中ではコールバック関数を直接呼んでますよね。
それが原因になっていたりしませんか？
0258251
05/01/28 22:19:53できました！
ご指摘の通り，関数の呼び出しを間違えていたようです．
ありがとうございました．
0259デフォルトの名無しさん
05/01/30 10:24:33一体どんな実装になったの?
0260choix ◆pnCm4bpG3Y
05/01/30 23:31:43ファイルの圧縮も行なわれるようになりました。
0261デフォルトの名無しさん
05/02/01 16:13:36出来ないのでしょうか？教えてください、お願いします。
0262デフォルトの名無しさん
05/02/04 22:19:421
--------
1 + Z^-1
を連続の
1
--------
1 + S
に変換したいのですが、どうしたらよいでしょうか。
関数などあれば、出来れば詳しく教えてください。
お願いします。
0263デフォルトの名無しさん
05/02/05 13:12:22クライアントサイド厨はネットワーキングに疎い知ったかぶりばかりで
このスレはもとからいいスレじゃあないな。
0264デフォルトの名無しさん
05/02/05 13:16:54教えてくれないといより
プログラミングスキルが低い教授なんだろ。
たいていの場合、教授よりもIT業界で
プログラミングだけで喰っていける実践経験のある人間のほうが
そういうことには強い。
教授は新しい技術にはついていけず数学しか知らないというのが実体だよ。
教授の中でも特に50代という団塊の人間にはコンピュータが苦手な連中が多い。
古い技術にとりつかれて新しいものになかなか鞍替えできない変化に弱い人間
なだけに彼ら教授陣の能力も衰えてきている、
0265デフォルトの名無しさん
05/02/05 13:21:08> data = load('****.***');
> filename(1,:) = 'D:\work\dev\matlab\test1.mat';
> :
> filename(n,:) = 'D:\work\dev\matlab\testn.mat';
ちゃんとファイル分割しているかしらんが
もししていないならこんなハードコーディングはおすすめできんな。
ソースコードを読む側の立場になったら読むだけでも腹が立つようになってくるぞ。
0266デフォルトの名無しさん
05/02/06 01:08:30上の方にいるからまた変なのが来た。
0267デフォルトの名無しさん
05/02/07 11:34:01処理が長くて困っています。
で、↓のfor文の処理が長いようなのですが
高速化する方法ないでしょうか？
for i = 1:(X)
tmp(i,1) = moto_file(i+(10*A+3),1);
end
読み込んだデータ(mot_file)の10*A+3番目以降のデータをtmpの中に
そのまま移行させたいのです...
よい方法をご存知の方がいらっしゃいましたら
なにとぞご教授お願いします。m(_ _;)m
0268デフォルトの名無しさん
05/02/07 11:39:34上のmoto_file=mot_fileです
0269デフォルトの名無しさん
05/02/07 13:21:54tmp=moto_file((10*A+3):end);
でどうだ？
tmpは消え去るけどな。
それか、moto_fileをいじってもよければ
moto_file(1:(10*A+2))=[];
で、１０＊Ａ＋２までがデリートされる
0270デフォルトの名無しさん
05/02/07 13:43:57早速試してみます
0271269
05/02/07 14:24:46めちゃめちゃ早くなりましたよ!
数十分かかってたのが数秒ですよ！
ありがとうございました!!!!
0272choix ◆pnCm4bpG3Y
05/02/08 08:40:00英語文章なんですが、Mathworksのサイトに書き換えの指針がまとめてありました。
http://www.mathworks.com/support/tech-notes/1100/1109.html
0273デフォルトの名無しさん
05/02/08 21:07:540274デフォルトの名無しさん
05/02/08 21:11:45おまいはサイバーネット社のまわし者か？
0275デフォルトの名無しさん
05/02/08 21:22:350277デフォルトの名無しさん
05/02/09 10:10:33--------------------------------------------------------------------------------
The frequency response of a linear time invariant system has been
measured and can be seen in the attachment file below.
A feedback element with a transfer function, H(s), is to be
chosen so that the sysytem steady-state error to a unit step input
is equal to 0.07. What is a suitable transfer function to achieve this???
これで、frequencyとmagnitudeとphase angleの数値をざっと各30個ずつ
ぐらいもらってるんだけど、まったく分かりません。どういう過程で
transfer functionを求めていけばいいのでしょうか？ matlabを使ってtransfer
functionを求められると聞いたんですけど、どういう風に数字を入れていけば
いいのでしょうか？
0278269
05/02/09 12:50:33う、英語文章苦手なんですけど
がんばって読んでみます。
情報ありがとうございますm(_ _)m
0279choix ◆pnCm4bpG3Y
05/02/09 14:30:19M系列信号を計算することはMatlabでできますが、
それを再生してマイクで拾って録音して、という部分は別のシステムを使った方がラクでしょう。
そのあとで録音データにFFTをかけて周波数特性を確認するのはMatlabでできます。
こんな流れなんですが、周波数特性の分かっているマイクとアンプを使わないと、スピーカーだけの周波数特性は計測できません。
完全なスピーカー特性をはかるには専用のデバイス類が必要になるのでかなり大変です。
0280２７３
05/02/10 05:09:23騒音計使うので平気なのです
アンプもフラットなので・・・・
その得られた周波数特性から逆フィルタを求めて
ソースなりに逆フィルかませるということは可能ですか？
0281choix ◆pnCm4bpG3Y
05/02/10 05:34:37別ソースにフィルタをかけることも可能です。
0282デフォルトの名無しさん
05/02/10 10:09:51contec製のADボードを使うと、入力が恐らく0.25倍、出力が2倍になってしまいます。
バージョンは2.7ですがMATLAB6.5付属の物でも同じ症状がでます。
インターフェース製のＡＤＤＡボードだと問題無いのですが、これはバグか何かなのでしょうか。
0.25倍と2倍というのもテスターによって測ったあやふやなものなので詳細がわかりません。
知っている方いらっしゃいますか？
0283282
05/02/10 10:11:11使用しているボードはcontecのAD12-16Eです。
0284デフォルトの名無しさん
05/02/12 02:37:18>>232
すごい遅レスだけど，Unbufferで一旦フレームベースのデータを変換して，
その後でBufferで1024のフレームサイズに直せばいけるかな？
漏れはフレームベース使ってないから間違ってるかもしれないけど・・・
最初に0がつくのは初期値が0になってるからだと思われ。
>>235
Complex-to-Real Imagブロックで実数虚数部に分けるか
absブロックで絶対値取るかだね。
MATLABはともかくとして，Simulinker少なっ
0285デフォルトの名無しさん
05/02/21 13:43:18MetlabはVer6.1です。
0286273
05/02/22 17:51:220287デフォルトの名無しさん
05/02/24 01:59:01具体的にと聞かれても・・・（ ´Д｀；)
適当にプログラム組んでみて，上手くいかないところを聞くのが筋では？
全部プログラムを作ってくれ，というように聞こえるけど・・・
0288choix ◆pnCm4bpG3Y
05/02/24 09:40:18やらなきゃいけないことが重なっていっぱいいっぱいで（ｒ
0289デフォルトの名無しさん
05/02/26 09:10:36オフスクリーンでやれば出来ないことはないが、サイズが指定できないので注意
0290デフォルトの名無しさん
05/02/26 18:12:52>> a = [ 1 2 ];
>> plot(a)
??? Undefined function or variable 'lineseriesmex'.
Error in ==> lineseries at 19
lineseriesmex;
とつれないお返事。このエラーにどう対処したらいいのでしょうか？
0291デフォルトの名無しさん
05/03/07 15:10:54遅レスだけど，Ver6.5.1だと普通にプロットできますた。
Ver7だとできないのかな・・・・
0292choix ◆pnCm4bpG3Y
05/03/09 02:03:33遅レス申し訳ない。ようやくちょっと時間ができました。
十分な長さのノイズ信号を再生した録音データには全周波数帯での応答が含まれています。
その録音データにFFTをかけてマグニチュード応答を割り出し、
それをもとにフィルタを作ります。
フィルタを作るにはYule-Walker法（yulewalk.m）なんかがラクチンです。
0293デフォルトの名無しさん
05/03/09 17:18:18Numerical Recipe 印C
をみろ
0294choix ◆pnCm4bpG3Y
05/03/09 22:16:59FFTのところか乱数生成のところ？
0295デフォルトの名無しさん
05/03/14 15:17:530296デフォルトの名無しさん
05/03/14 15:41:38自己解決
getframe
ね。
0297デフォルトの名無しさん
05/03/16 11:16:410298デフォルトの名無しさん
05/03/16 11:42:13難しい？
0299choix ◆pnCm4bpG3Y
05/03/16 13:28:45このスレの過去に似た話が出てきたのでチェックすると良いかもしれません。
>>44-48 あたりです。
>>298
意外に簡単です。
ただ、システム依存する部分はちょっと難しいです。
自分は、多チャンネル音声再生をマルチプラットフォームで
実行できるようなmexファイルを作ろうとして挫折しました。
0300デフォルトの名無しさん
05/03/17 12:55:09Ｍａｔｌａｂじゃなくて普通のＣのプログラムから呼び出すには
どうしたらよろしいのでしょうか
0301choix ◆pnCm4bpG3Y
05/03/17 21:38:46英語版のホームページで言うと
http://www.mathworks.com/access/helpdesk/help/techdoc/matlab_external/ch_java.html
です。これはMatlabそのものを他のプログラムから利用するというものです。
参考になるかもしれません。
0302デフォルトの名無しさん
05/03/19 12:29:360303デフォルトの名無しさん
05/03/19 22:50:510304デフォルトの名無しさん
2005/03/21(月) 08:27:10f = @(x) x
みたいなことは、octave では出来ないんでしょうか?
0305デフォルトの名無しさん
2005/03/30(水) 09:52:44std=sqrt(varience);
x=std*randn(1,10)+mean*ones(1,10);
この公式分かる人いますか？
0306デフォルトの名無しさん
2005/03/30(水) 11:37:25もうちょっと情報ないの？
0307デフォルトの名無しさん
2005/03/30(水) 12:49:26何ですか？
0308デフォルトの名無しさん
2005/03/30(水) 13:08:51んですけど、出ません。どなたか何故だか分かりますか？
>> m=10;
varience=9;
stdd=sqrt(varience);
%length_10
x_10=stdd*randn(1,10)+m*ones(1,10);
subplot(4,1,1);plot(x_10)
%length_100
x_100=stdd*randn(1,100)+m*ones(1,100);
subplot(4,1,2);plot(x_100)
%length_1000
x_1000=stdd*randn(1,1000)+m*ones(1,1000);
subplot(4,1,3);plot(x_1000)
%length_10000
x_10000=stdd*randn(1,10000)+m*ones(1,10000);
subplot(4,1,4);plot(x_10000)
>>
0309choix ◆pnCm4bpG3Y
2005/03/30(水) 19:17:33とりあえず「平均mu分散sigma^2の正規分布に従う乱数列が作りたい」という質問だと想定して答えますと、以下のようなプログラムができます。
mu = 10; % 平均
variance = 9; % 分散
sigma = sqrt(variance); % 標準偏差
x = randn(1, 10000); % 平均0.0分散1.0の正規分布乱数を生成
x = (x * sigma) + mu; % 平均mu分散sigma^2に移動
hist(x, 50); % ヒストグラム表示
この後、std(x)やvar(x)で指定した標準偏差や分散値に近い値になっていることが確認できます。
0310デフォルトの名無しさん
2005/03/31(木) 04:01:250311デフォルトの名無しさん
皇紀2665/04/01(金) 19:32:510312デフォルトの名無しさん
2005/04/12(火) 15:27:190313デフォルトの名無しさん
2005/04/24(日) 16:00:41面白そうじゃん。
0314デフォルトの名無しさん
2005/04/27(水) 09:00:110315科学者の端くれ
2005/04/27(水) 10:37:580316デフォルトの名無しさん
2005/04/29(金) 18:00:17大きな行列を使った関数をたくさん呼んでるんですが、
値渡しのコストが大きそうな感じ。
値渡しを避けるには変数をglobalにするしか方法はないんでしょうか?
0317デフォルトの名無しさん
2005/05/01(日) 14:35:43演算がコピーしてるじゃないの？
コンパイラの最適化の話じゃない？
0318デフォルトの名無しさん
2005/05/08(日) 19:18:37購入を検討しているのですが。
0319デフォルトの名無しさん
2005/05/08(日) 21:08:48サイバネットに聞け。マジで。
toolboxの組み合わせが膨大で、ここで答えられる範疇を超えてる。
0320デフォルトの名無しさん
2005/05/08(日) 23:37:13本国だと学生向けアカデミックパッケージがだいたい100ドルくらいだけど、
グループライセンスはウェブには価格が出てないね。たいていは大学全体での
サイトライセンスなんじゃないかな。
0321デフォルトの名無しさん
2005/05/09(月) 00:26:26MATLABのグループライセンスの場合、
LAN内にライセンスサーバを一台立てて、同時使用数とかのチェックが入るから
大学全体だと数が多すぎて微妙だったりするので
研究室単位でライセンスを買ってるところも結構ある。
0322デフォルトの名無しさん
2005/05/09(月) 02:07:04もれの研究室がライセンス持っているみたいだけど、他の研究室は持ってないみたい。
ちなみに、研究目的で自分のノートPCにインストールして、自宅でするのは合法？
ってライセンスの組み合わせによって変わるか…
0323デフォルトの名無しさん
2005/05/09(月) 06:13:54起動できないんじゃないか。
＃漏れのノートにも入ってるが研究室でしか使えん。
スタンドアロン（認証サーバにつながない）ライセンスのがあるなら
インストールしてあるPCの数がそのライセンスの数までなら合法
それより多くのPCに入れるとたぶんアウトかな？
一度ライセンスをよく読んでみるのお勧め。
0324デフォルトの名無しさん
2005/05/09(月) 08:14:45サンクス
0325デフォルトの名無しさん
2005/05/09(月) 22:56:41どうもサンクス！
教授がどうもそういうのに疎くて…
一度よく読んでみるよ。
0326デフォルトの名無しさん
2005/05/10(火) 07:48:520327デフォルトの名無しさん
2005/05/10(火) 12:16:09ttp://www.cybernet.co.jp/matlab/support.shtml
0328デフォルトの名無しさん
2005/05/10(火) 14:03:26The MathWorks product documentation you have requested is no longer available.
なのだが、ヘルプはどこへ行ったのでしょうか？
0329デフォルトの名無しさん
2005/05/10(火) 17:59:50最近、研究室でMATLABを使っているのですが、エラーが出てきた時に音が出るのを
消したいのですが、どうしたらよいですか？コントロールパネルからではパソコンの本体から出る「プッ」
なんていう音が消せないのです。エラーが出るたびにほかの人の注目を浴びてしまいます。。。
よろしくお願いします。
0330デフォルトの名無しさん
2005/05/10(火) 19:24:43単位円をx軸に投影するプログラムを作るのですが。
ラジアンをｔで設定し、x=cos(t),y=sin(t)をplot(x,y)で円をプロットするところまではできたのですが。
この後、radonというコマンドを使って投影をするそうなのですが、そっから先に進みません。
radon(???)の？の部分に何を入力したらいいのか分からないのです.
すいません、よろしくお願いします・・・。
0331デフォルトの名無しさん
2005/05/10(火) 19:53:02うぇｗｗｗっおｋｗｗｗｗｗｗうぇｗｗｗっうぇっｗｗｗｗｗｗｗｗｗｗｗｗ
っうぇおｋｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗ
うぇｗｗｗ
っうぇｗｗｗおｋｗｗｗ
0332デフォルトの名無しさん
2005/05/12(木) 17:51:39できます？
0333デフォルトの名無しさん
2005/05/13(金) 16:34:04外側を０とするにはどのようなプログラムでできるのでしょうか？
0334choix ◆pnCm4bpG3Y
2005/05/13(金) 23:49:54これは >>330 のradonの使用法が分かった上で、単位円をradonが使用する形式に書きたい、ということでいいですか？
単に「強度」といわれても、一般的にはさっぱりわからないのです。
もっと分かりやすい文章にしてくださいね。
以下のプログラムで半径256点の円が描けます。
N = 256;
I = zeros(N);
for row = 1:N
for col = 1:N
if (row/N)^2 + (col/N)^2 < 1
I(row, col) = 1;
else
I(row, col) = 0;
end
end
end
imshow(I);
0335choix ◆pnCm4bpG3Y
2005/05/13(金) 23:58:21title命令はtext命令に対するラッパーのような動きをするので、
textのプロパティをそのまま使用できます。
h = title('hogehoge');
set(h, 'Position', [x y]);
で位置を変更できます。
0336デフォルトの名無しさん
2005/05/15(日) 14:44:190337choix ◆pnCm4bpG3Y
2005/05/15(日) 17:29:460338デフォルトの名無しさん
2005/05/16(月) 01:12:48x= dir()
みたいに、dir の返り値をとって変数に代入できますが、
octave で同じことをする方法を教えて下さい。
（そのままだと右辺が未定義として怒られる）
0339デフォルトの名無しさん
2005/05/16(月) 16:25:50LMIを扱った問題について質問させてください。
あるスカラ変数τ*（タウスター）があるとします。
この時 (0<=τ<=τ*) の範囲で、システムを安定させる行列変数Pを求めたいのです。
この場合どのようにLMI式を立てればよいのでしょうか。
whileループでひたすらτ*の値を増加させていくしかないのでしょうか。
できればきっちりと限界値がわかるようにMATLABで計算してもらいたいのです。
しかしlmitermでは変数が一つしか扱えません。
何かうまい使い方があるのでしょうか
0340迷える子羊
2005/05/19(木) 21:32:55行列の積分の事について質問させて下さい。
y=4*x+yなどの積分、たとえば ∬(4*x+y)dxdy （1<x<3,3<y<7)
などはq = dblquad(fun,xmin,xmax,ymin,ymax)でもとまるのですが
行列の積分、たとえば ｜∬(4*x+y)dxdy、∬(*x+y)dxdy ｜
｜∬(2*x+y)dxdy、∬(*3x+y)dxdy ｜
全て範囲は（1<x<3,3<y<7)
の2行2列の積分を一気に求める方法（関数）はないのでしょうか。
誰か知っている方ご教授願います。よろしくお願いします。
0341デフォルトの名無しさん
2005/05/21(土) 09:09:47教授の知り合いで信号処理に詳しい他の教授はいないのか？
他の研究室、大学なら何かしらMATLABについて詳しい香具師がいると
思うのだが。
音響工学、画像工学の研究室なら大抵つかっているし
0342デフォルトの名無しさん
2005/05/23(月) 14:21:14(s^2+4*s)/s/(s-th1)*th2
という解が求まりました。
まずはこの式を
(s^2+4*s)/(th2*s^2-th1*th2*s)
というsに対する項べきの順に整理して、その係数を
分母の係数 [1 4 0]
分子の係数 [th2 -th1*th2 0]
のような行列に抽出したいと思っているのですが、うまく行きません。
どなたかご存知の方がいらしたら方法を教えてもらえないでしょうか？
0343Zeno ◆5nZQbNmQPs
2005/05/23(月) 14:24:27吾は Matlab フルバージョンを使用している。
せっかくだから matlab で *.wav の編集や、画像の編集技術を
学びたいと思っているのだが、それ系でなにかいい website はないか？
____________
"Why the future doesn't need us." by BILL JOY, the cofounder and Chief Scientist of Sun Microsystems
http://www.wired.com/wired/archive/8.04/joy.html
0344デフォルトの名無しさん
2005/05/23(月) 21:07:55クラスに静的メンバ関数を実装することはできるでしょうか？
クラスで実装しているメンバ関数の間で、
簡単なユーティリティ関数を共有したいのですが、
どのようにしたらよいでしょうか？
よろしくお願い致します。
0345デフォルトの名無しさん
2005/05/23(月) 21:08:44「MATLAB によるｘｘｘｘ」という本がたくさん出てるよ。
マルチメディア処理の本とかもあったな。
0346デフォルトの名無しさん
2005/05/24(火) 07:11:53OpenGLを使わないときれいに出ます
0347デフォルトの名無しさん
2005/05/24(火) 14:49:19警告: OpenGLのバージョンが認識できません。デフォルト1.0を利用します.
という警告がグラフィックを操作するたびに出ます。この警告を出ないようにする方法はないでしょうか？
0348Zeno ◆5nZQbNmQPs
2005/05/24(火) 16:06:15本を買う金がない。 Website があると助かるのだが。
____________
"Why the future doesn't need us." by BILL JOY, the cofounder and Chief Scientist of Sun Microsystems
http://www.wired.com/wired/archive/8.04/joy.html
0349デフォルトの名無しさん
2005/05/24(火) 16:35:29オブジェクト指向のサポート度数はあがったの？
PHPがver4からver5にアップグレードしたときに
オブジェクト指向のサポートが非常に強力になったように。
0350choix ◆pnCm4bpG3Y
2005/05/24(火) 22:01:22まずは本家のドキュメントが良いかと思います。
http://www.cybernet.co.jp/matlab/support/helpdesk/r14/index.shtml
Signal Processing ToolboxのUser's Guideで
信号処理の基本あたりを学んで、
Image Processing Toolboxで画像処理を勉強する、
なんていうのはどうでしょう？
0351Zeno ◆5nZQbNmQPs
2005/05/25(水) 16:40:36やはり User's Guide が一番能率的であるかもしれないな。
Matlab 自体は好きなのだが、現在魅惑があらゆる角度から迫ってくる故、
なかなか本気になれないのが現実だ。 好きなものがもっと好きなものに押し
潰されるとは皮肉であるだと思わないか？
しかしながら、今度読む予定に入れておこう。
_________________
セオドア・カジンスキー 「産業社会とその未来」 (aka. マニフェスト) の抜粋
http://www.cnn.com/US/9509/unabomber/09-19/am/excerpt.html
0352デフォルトの名無しさん
2005/05/25(水) 16:58:06素直に「ありがとう」と言いなさい。
0353デフォルトの名無しさん
2005/05/25(水) 17:09:080354デフォルトの名無しさん
2005/05/26(木) 10:42:06このメッセージを消すにはどうしたらよいでしょうか？
0355デフォルトの名無しさん
2005/05/26(木) 15:49:160356デフォルトの名無しさん
2005/05/26(木) 16:32:120357デフォルトの名無しさん
2005/05/26(木) 16:36:140358デフォルトの名無しさん
2005/05/26(木) 16:43:090359デフォルトの名無しさん
2005/05/26(木) 16:48:01どうしようもないアフォだな
0360デフォルトの名無しさん
2005/05/26(木) 19:18:550361デフォルトの名無しさん
2005/05/27(金) 08:47:41貧弱なのか・・・・・
0362Zeno ◆5nZQbNmQPs
2005/06/01(水) 22:16:29貧弱である。
_________________
If T is consistent, T !|- G_T.
If T is ω-consistent, T !|- ￢G_T. -- Kurt Gödel
0363デフォルトの名無しさん
2005/06/02(木) 16:05:46ttp://www.petitcolas.net/fabien/software/mpeg/
ここにあるDownloadから落としたプログラムで音声を合成したいのですがなかなかうまくいきません。だれか分かる方がいたら教えてくれませんか？
ちなみにFFTがオーバーラップしていて合成はadd法を使えばいいということまで分かっているつもりなのですが・・・
よろしくお願いいたします。
0364デフォルトの名無しさん
2005/06/03(金) 00:23:23どこかに解説したページはありませんか？
0365デフォルトの名無しさん
2005/06/03(金) 14:39:450366デフォルトの名無しさん
2005/06/04(土) 14:34:40フィルタのコードを書くしかないのでしょうか？
0367デフォルトの名無しさん
2005/06/05(日) 05:38:47conv2とかじゃだめなん？
0368デフォルトの名無しさん
2005/06/05(日) 14:27:49[C,h] = contour(x,y,z,0);
とやるとうまくいきません
[C,h] = contour(x,y,z,[0 3]);
のようにしてやらないとダメなのでしょうか？3のところにも余分な線が出来てしまっていやなのですが
0369デフォルトの名無しさん
2005/06/05(日) 23:55:17マニュアル読めば分かるよ。
loadlibrary だよ。
loadlibrary で、外部ライブラリと、その関数を宣言したヘッダーを読み込んで、
callib するだけだよ。
0370デフォルトの名無しさん
2005/06/05(日) 23:55:59どういうこと？
func(arg1, arg2)
でいいんだろ？
0371デフォルトの名無しさん
2005/06/07(火) 15:09:330372デフォルトの名無しさん
2005/06/14(火) 18:33:04これは誤差としてみなしてよいのでしょうか？
0373デフォルトの名無しさん
2005/06/15(水) 00:30:56その包茎波の大きさは何よ。
そもそも、MATLABとはいえ所詮有限桁数での演算シミュレーションなんだから
こんな得体のしれないところで訊くまでもなく、ちょっと考えりゃ答えは分かるだろ。
0374デフォルトの名無しさん
2005/06/17(金) 23:27:20良いというか，以下の条件を満たすような説明．
・豊富なサンプルの提示
・落し穴や注意点の説明
・Win/LinuxやC/C++/Fortranといった違いに触れてある解説
・できれば日本語
サイバネットのサーバって遅いし，力が入っていない気がするし…
お薦めページを御存じの方，よろしくお願いします．
0375デフォルトの名無しさん
2005/06/21(火) 13:26:370376デフォルトの名無しさん
2005/06/23(木) 18:12:26simulinkでのやり方がわかりません。
わかる方がいらっしゃたら、ご教授お願いします。
0377デフォルトの名無しさん
2005/06/23(木) 23:17:22Matlabが手元にないんで不明瞭な回答になるけど、
「0次ホールド」「zero order hold」「hold」「A/D」などの
文字がSimulinkの要素選択画面にない？
あったならば、それが正解。
0378デフォルトの名無しさん
2005/06/24(金) 02:52:55こちらの手元にもないんですが、
確か「zero order hold」っていうのはありました。
これは、A/D変換ではなく、D/A変換だと聞いたのですが
間違っていますか？
0379デフォルトの名無しさん
2005/06/24(金) 12:02:34漏れSimulinkerだから適当にモデル説明してくれれば分かると思う。
連続系から離散的にっていうのは，例えばサイン波形をサンプリングとかそういう感じ？
0380デフォルトの名無しさん
2005/06/24(金) 14:44:21おぉ、ありがたいです。
僕もよくわからないのですが、Monodモデルです。
連続系から離散的にっていうのはおっしゃる通りです。
0381デフォルトの名無しさん
2005/06/24(金) 16:34:28>>377氏が言うように，zero order holdか，別にGainでもいいんだけど，
ブロックを入れてサンプル時間を設定すればいい。
サンプル時間についてはそっちの設定が分からないからなんともいえないけど。
試しに1くらいでやってみて，それで駄目なら変えてみるとかね。
0382376
2005/06/24(金) 22:38:24適当に入力にランプ入力入れて zero order hold はさんで
Scopeで見てみましたところ、それらしきものがでてきました。
どうもありがとさんです！
で、疑問に思ったのですが、
どっかから適当に入力とってきてそれをA/D変換して
その値に適当に処理（フィルタなど）を施してD/A変換して
入力を取って来たループに戻したいのですが、
こういう事はできますでしょうか？
zero order holdを通したら、
一気にA/D変換→D/A変換をセットでしているような気がして
A/D変換とD/A変換の間に処理を施す事ができないのでは？
と思ったのですが、どうなんでしょう?(´･ω･`)ﾖｸﾜｶﾗﾝ
0383デフォルトの名無しさん
2005/06/25(土) 15:33:00381だが，Ａ/Ｄ，Ｄ/Ａ変換は正直よく分からない部分もあって･･･('A`)
A/D変換はできるけど，D/A変換は完全に出来るかどうかは疑問。
っていうか，Simulink上では一見アナログに見える波形も，実はデジタル波形な希ガス。
サンプル時間を非常に短くとっているので，アナログっぽく見えるだけのような･･･
まぁ，一応モデルを作ってみてそれっぽい動作は確認できた。
適当なメアドかメッセの番号，あとはうｐろだ探してくれれば送るよ。
0384デフォルトの名無しさん
2005/06/25(土) 21:13:35そうなんです。D/A変換がよく分からないんですよね(^ω^；)
モデルわざわざ作っていただきありがとうごじまます。
メール欄のとこのアドにお願いします。
手元にMATLABがないので確認は月曜日になりますが…
0385デフォルトの名無しさん
2005/06/25(土) 21:19:04×ありがとうごじまます
↓
○ありがとうございます
でした。失礼。
0386デフォルトの名無しさん
2005/06/26(日) 02:21:03OK，明日か明後日には送っておく。
ところでふと疑問に思ったんだけど，A/D・D/A変換する必要があるの？
どんな処理か分からないけどさ。
0387デフォルトの名無しさん
2005/06/27(月) 22:07:040388デフォルトの名無しさん
2005/06/27(月) 23:59:32ちと意味がわからないんだが・・・
ネットワークライセンスのことかな？
0389デフォルトの名無しさん
2005/06/28(火) 00:09:14都市伝説かなんかでPhotoshopには個人データを送信しているとかいないとか
噂されていますよね？そういう小細工みたいなことしこんであるのかなと思いまして。
まあ、ライセンス違反をした時ばれるかってことです。
0390デフォルトの名無しさん
2005/06/28(火) 00:16:10ライセンス違反してるの？（ ´Д｀；)
0391デフォルトの名無しさん
2005/06/28(火) 00:20:47いえ、こういうツールって高いじゃないですか？
だから、高価なツールにはみんな変な細工されてるのかなと思いまして(ﾞ ｀-´)/
0392デフォルトの名無しさん
2005/06/28(火) 00:23:02大丈夫ｼﾞｬﾈｰﾉ?
0393デフォルトの名無しさん
2005/06/28(火) 00:28:47了解っす(ﾞ ｀-´)/
では開発者が現れるまで待ち続けるぞー！
0394デフォルトの名無しさん
2005/06/28(火) 12:08:350395デフォルトの名無しさん
2005/06/28(火) 17:57:58There may be a singularity in the solution. If the model is correct,
try reducing the step size (either by reducing the fixed step size or by tightening the error tolerances)
と表示されて翻訳してみたら、「時0に計算されたを含む代数的な輪における代数的な状態は、
InfかNaNです。 解決策には特異性があるかもしれません。
モデルが正しいなら、ステップサイズを減少させてみてください。
(固定ステップサイズを減少させるか、または誤り寛容をきびしくするのによる)
」と表示されたのですがどうすればいいNO＿？
0396デフォルトの名無しさん
2005/06/29(水) 00:20:35Infは無限大なんだけど，値が発散してるってことかな？
0397デフォルトの名無しさん
2005/06/29(水) 00:23:38時間が余ってからCで作り直すか
0398デフォルトの名無しさん
2005/06/29(水) 00:35:10その名の通り行列計算はかなり楽だよ。
たとえば，ある行列Aの逆行列を求める時は，Cは複雑なアルゴリズムが必要だけど
MATLABはinv(A)だけで出る。
グラフ描画も楽だし，複素数も普通に扱える。
欠点はアプリが重いってことかねぇ･･･
0399あああ
2005/06/29(水) 00:35:17それぞれcmdとdataのメモリの内容、メモリから読み出した直後
のレジスタの内容、実際に出力される内容はそれぞれどのような
値になるでしょう？
char cmdreg = CMD_PORT, datreg=DATA_PORT;
unsigned int16 cmd;
unsigned int32 data;
cmd=0x87cd , data=0x20030719;
outw(cmdreg,(cmd & 0xff00)>>8 | (cmd & 0xff)<<8);
outl(datreg, data);
0400デフォルトの名無しさん
2005/06/29(水) 02:19:35そうみたいです。ｔｈｘ。
それで固定ステップサイズを減少させるというのは
どういった意味があるんでしょうか？
>>397
なんで教授にバカにされるの？
0401デフォルトの名無しさん
2005/06/29(水) 02:28:33発散してるってことはどこか係数か何かを小さくしてみるとか。
それくらいしかアドバイスﾃﾞｷﾈ
0402デフォルトの名無しさん
2005/06/29(水) 04:01:25いろいろ試してみます。ドモー(*＾▽＾*)
0403デフォルトの名無しさん
2005/06/29(水) 12:31:33まず、ほかの人の論文をたくさん読んでみ。
結構MATLABで処理した結果を論文に載せてる人は多いし
そこまで気にすることじゃないだろう
0404デフォルトの名無しさん
2005/06/29(水) 15:18:030405デフォルトの名無しさん
2005/06/29(水) 16:26:48レンダラーをOpenGL以外にすれば良いです。
0406デフォルトの名無しさん
2005/06/30(木) 14:24:38ワロス
0407デフォルトの名無しさん
2005/07/04(月) 11:44:580408デフォルトの名無しさん
2005/07/14(木) 03:58:59うちはファイアウォール内でライセンスサーバ動かしていますけど、
何の問題もなく使えていますよ。
もちろんヤバいことはしていない。
0409デフォルトの名無しさん
2005/07/14(木) 04:01:49「数値シミュレーションをやりたい！」という目的を忘れてCの
Segmentation Faultをデバッグするのが好きな人には馬鹿にされるかも。
0410デフォルトの名無しさん
2005/07/18(月) 10:56:57'LineWidth',5などを追加しても、認識してくれません。
0411デフォルトの名無しさん
2005/07/20(水) 16:33:59質量と剛性なら[d,v]=eig(M,K)でできたんですけど、粘性まで
入ってくるとどうもエラーが．．．
詳しい方お願いします。
0412デフォルトの名無しさん
2005/07/21(木) 09:41:28また緯度経度データを読み込んで地図上にプロットしたいのですが。
よろしくお願いいたします。
0413デフォルトの名無しさん
2005/07/27(水) 04:18:090414デフォルトの名無しさん
2005/07/28(木) 15:52:220415デフォルトの名無しさん
2005/08/01(月) 11:49:04関数f(x,y)をfsolve(@f,x0)のようなことをしたいんですが、
yはgivenで、xが変数なんです。
yはどのように指定したらいいでしょうか。
0416415
2005/08/02(火) 01:45:48x = fsolve(fun,x0,options,P1,P2,...)
ですね。
0417デフォルトの名無しさん
2005/08/02(火) 08:40:50自分も分からなくて困っていたことがありました。
ありがとうございました。
0418デフォルトの名無しさん
2005/08/09(火) 09:37:56の特性F(w)に対してその逆関数1/F(w)を求めることとかって可能でしょうか？
0419デフォルトの名無しさん
2005/08/13(土) 11:59:21逆関数だったら複素数(ゲインと位相)の関数になるでしょ。
周波数応答って全単射じゃないから逆関数は一般には存在しない。
0420418
2005/08/13(土) 16:55:36そうするとスピーカの振幅特性の逆関数が判らないとできなくないですか？
0421デフォルトの名無しさん
2005/08/13(土) 18:08:05板違いじゃない？(Matlabじゃなくって数学の質問のような気がする)
例えば、
G(s)=e^{sL}, sは虚軸の上半分, Lは十分大きな正の数
の逆関数なんて存在しないでしょ。周波数を100Hz～16Hzに限定しても。
0422デフォルトの名無しさん
2005/08/14(日) 00:51:030423デフォルトの名無しさん
2005/08/14(日) 07:33:420424デフォルトの名無しさん
2005/08/14(日) 10:49:21どうするつもりなんだろうか。
0425デフォルトの名無しさん
2005/08/14(日) 14:55:14それがダメならPython。最終的にはC言語。
0426デフォルトの名無しさん
2005/08/15(月) 13:39:32を使う
0427デフォルトの名無しさん
2005/08/22(月) 05:26:110428デフォルトの名無しさん
2005/08/22(月) 07:13:00韓国の大手テレビ局が、８月１５日のニュース番組で報道した「旧日本軍による生体実験」と称する映像が、実は真っ赤なウソ、中国映画の盗用でっち上げだった。
韓国のマスコミは反日キャンペーンで、この種のでっち上げ「歴史歪曲（わいきょく）」をよく行うが、今回の事件で 「日本糾弾ならウソも平気」 といった腐敗体質が、改めて確認された。
番組では戦前、旧日本軍の細菌戦研究７３１部隊の生体実験の様子を「白黒の記録フィルム風」に紹介したが、視聴者から「映画の場面と同じだ」と抗議の声が上がり、悪質なでっち上げが露呈した。
この中国映画はカラー作品だが、テレビでは本物に見せかけるために「白黒」に変えられており、意図的なでっち上げであることは明白。
日本の民主党（岡田党首）は、日本国内に数十万人いる反日的な在日韓国・朝鮮人に日本の選挙権を与え、民主党に投票させようと画策している。
在日韓国・朝鮮人の大部分は、戦後のドサクサで密入国した不法滞在者とその子孫だが、強制連行されたとウソをついて被害者のフリをし、日本国内に居座って、毎月二十数万円の生活保護をもらい、遊んで暮らしている。
0429デフォルトの名無しさん
2005/08/22(月) 08:44:310430デフォルトの名無しさん
2005/08/22(月) 11:45:37Statistics Toolbox内にnlinfitっていうのがあるから、それを使うべし。
0431デフォルトの名無しさん
2005/08/24(水) 12:02:45例えば
100から10000までの値を10000分割して格納したいのですが・・
よろしくお願いします。
0432デフォルトの名無しさん
2005/08/25(木) 02:04:57100～10000を500分割。
ちなみにlogでいきたければlogspaceという関数もあり。
0433431
2005/08/25(木) 21:52:41うわ、そういう風な関数あったんですね。
こんな風に書いてました。
nは分けたい数です。
ありがとうございました。
A = min;
max - min
for i = min + (max-min)/n : (max-min)/n : max;
A = [A i];
end
0434デフォルトの名無しさん
2005/08/26(金) 00:41:51今まで通りやっても大丈夫だよ。
0435デフォルトの名無しさん
2005/08/30(火) 10:18:07データの最初の２行の読み込みを飛ばしたい時はどういう風にすればよいですか？
0436デフォルトの名無しさん
2005/08/30(火) 12:44:08「最初の二行を飛ばしたい」と言ってるので、おそらく自分で用意したテキストファイルやCSVファイルを読み込もうとしているんじゃないかな？
そういうときはdlmreadやcsvreadという専用コマンドを使った方がいいよ。
A = dlmread('hoge.txt', ':', 2);
とやれば、hoge.txtの中にあるコロン区切りのデータを最初の二行を飛ばして読み込んでくれる。
0437デフォルトの名無しさん
2005/09/07(水) 16:38:49f(x,y)=ax+byのような平面スライスの方法を教えていただけないでしょうか
0438437
2005/09/08(木) 11:11:570439名無しさん＠そうだ選挙に行こう
2005/09/11(日) 12:00:52でも、サイバネットのサイトって鯖が遅過ぎなんで使えん。
そう思っているのってボクだけですか？
0440名無しさん＠そうだ選挙に行こう
2005/09/11(日) 17:44:21質問してもすぐ返事が返ってくる。
漏れの使ってるハードウェアのサポートが他の会社に移った途端，
返事が返ってこなくなって実感した('A`)
>>440
関係者乙
0441デフォルトの名無しさん
2005/09/12(月) 09:09:21( ´д)ﾋｿ(´д｀)ﾋｿ(д｀ )
0442デフォルトの名無しさん
2005/09/12(月) 09:53:050443デフォルトの名無しさん
2005/09/12(月) 21:59:390444デフォルトの名無しさん
2005/09/17(土) 01:55:40知らんで０から作ってしまった
0445デフォルトの名無しさん
2005/09/24(土) 04:34:35FFTをかけることで周期を抽出しようとしています．
また，途中でテンポが変化するため，窓のサイズを小さく取り，
時間方向に少しずつずらしていくことで，周期の変化を見ようとしています．
FFTの理解は十分でないのですが，こうした解析を行うために，
スペクトログラムが適しているのではないかと考えました．
つまり，パワースペクトルの時間変化をみたいということになります．
こういった経緯で，下記のようなプログラムを書いたのですが，
なぜかゼロ周波数付近に大きなピークが観測されてしまいます．
データは明確な周期性があり，他人のCプログラムで実行したところ，
ゼロ周波数付近にはピークが観測されず，想定された位置にピークが観測されています．．．
問題点がわからず悩んでいます．
なにか心当たりがある方がいらっしゃれば，ぜひ教えてください．
よろしくお願いします．
------[以下，ソース]
signal = load('hoge.txt');% データ数は約200点
FFT_domain = 32;
Hz = 20;
datNum = length( signal );
% specgramの使い方
% specgram(データ，窓長，サンプリング周波数，窓関数，オーバーラップの長さ)
specgram( signal, FFT_domain, Hz, [], FFT_domain-1 );
0446デフォルトの名無しさん
2005/09/24(土) 04:45:54コードで例えば、polyvalやdeconvはどうやって書かれているか見る
方法はありますか？
0447デフォルトの名無しさん
2005/09/24(土) 06:33:56type deconv
0448デフォルトの名無しさん
2005/09/24(土) 06:36:26DC成分って、つまりはx時系列データの平均値みたいなもんでしょ。
x = x - mean(x);
ってやってからspecgramに渡したらどんな感じ？
0449445
2005/09/24(土) 13:36:24アドバイス通りにやった後，absをかけるとなかなかうまく抽出できました．
ありがとうございました！
0450デフォルトの名無しさん
2005/09/24(土) 15:02:03>> help thanks
0451デフォルトの名無しさん
2005/10/04(火) 14:43:18ゼロ割という警告はいったい何を意味してるのですか？
0452デフォルトの名無しさん
2005/10/04(火) 15:19:16(たとえば 5.5/0 は不定なのでゼロ除算エラーになります)
0453デフォルトの名無しさん
2005/10/04(火) 16:56:47これがゆとり教育の（ry
0454デフォルトの名無しさん
2005/10/04(火) 20:54:24不定でもエラーでもないぞ。
不定とは、Cで未初期化メモリにアクセスしたときとかで、
MATLABでは起こらない。
MATLABではIEEE演算を行い5.5/0はInfになるので、
エラーは起こらない。
やってみればわかるが、エラーではなく警告が出てるはず。
0455デフォルトの名無しさん
2005/10/05(水) 15:40:46どうやって表示させればいいでしょうか？
0456デフォルトの名無しさん
2005/10/12(水) 01:29:39どうすればいいのでしょうか…。
具体的にはカルマンフィルタで求まった状態推定を抽出して、それを用いた計算を
行いたいのですが…。
0457デフォルトの名無しさん
2005/10/12(水) 09:54:48バックアップとりたいのだが。
0458デフォルトの名無しさん
2005/10/13(木) 12:52:14ここは工学とか物理とか、そういう専門のことを聞く場所じゃなくて、
Matlabプログラミングの質問をするところです。
もっと具体的に書いてもらわないと何が何だかわかりません。
>>457
バージョンやプラットフォームは？
コピーガードがついてるのはStudent Versionだけだったような気がする。
0459デフォルトの名無しさん
2005/10/13(木) 19:14:420460デフォルトの名無しさん
2005/10/13(木) 21:04:51４年しかつかえないのでしょうか？
0461デフォルトの名無しさん
2005/10/13(木) 21:53:25各ToolBoxもちゃんと使えるし。
起動時にプログラムCD-ROMを入れておかないといけない、
という起動制限があるだけ。
これも回避する方法はあるみたいだしね。
0462デフォルトの名無しさん
2005/10/14(金) 01:32:35どうやんの？
0463456
2005/10/14(金) 02:29:47sys=ss(A,B,C,D)
で表されているシステムの中の、
A行列にかかるx行列を抽出して変数にしたり(「set('statename')」ではなくて)、
要素のみを抜き出すことって出来るのでしょうか。
0464461
2005/10/14(金) 21:25:23自分はMac OS X版のR13 Student Version。
Roxio Toastでディスクイメージを作ってしまえばいい。
OS付属のDisk Utilityだとイメージの作成に失敗するけど、Toastなら大丈夫。
ただ、中途半端なイメージを強制作成することになるので、
一部ファイルが含まれなくなってしまう（PDFドキュメントとか）。
できあがるイメージはDisk Utilityで作成されるものと互換性があるから
Disk Utilityで圧縮イメージに変換してしまえばサイズも小さくなるよ。
0465デフォルトの名無しさん
2005/10/14(金) 21:29:04その「ss」というルーチンを自分は持ってないよ。
なんかのツールボックスのルーチンなんだな、きっと。
自分はお手上げ。
0466デフォルトの名無しさん
2005/10/14(金) 22:54:08状態変数ってのはベクトルの時系列でしょ。
つまり、
\dot{x}=Ax+Bu
y=Cx+Du
の解なわけだ(初期値は適当に設定するとして)。
だから初期値と入力(u)を決めないとどうにもならない。
初期値と入力が与えられたなら、適当なwrapperを関数として作って特定の状態を抽出すればよいだけ。
別に関数を作んなくてもCを決めるだけでOK。
一出力のシステムなら、例えばC=[1 0 0 0 0 .....]とすればx(1)を抽出できる。
0467463
2005/10/15(土) 04:45:01まさしくコントロールツールボックスというモノの中のやつです…。
>>466
ということとは初期値とＣをいじって、システムの出力を扱えばいいんですね…。
ありがとうございました、やってみます！
0468デフォルトの名無しさん
2005/10/15(土) 10:30:020469デフォルトの名無しさん
2005/10/16(日) 14:52:31thx!
0470デフォルトの名無しさん
2005/10/19(水) 00:15:21どうやればよいのでしょうか？
0471デフォルトの名無しさん
2005/10/19(水) 07:33:58どうすればいいでしょう。
audiorecorderオブジェクトは、1秒録音するのに2秒近くかかるので
困っています。
オブジェクトを複数立ち上げて交代で録音しようとしたら、フリーズ
しました。
プラットフォームはWindowsXPです。
0472デフォルトの名無しさん
2005/10/20(木) 00:58:16その動作がわからないとなんともいけないわけで。
0473デフォルトの名無しさん
2005/10/21(金) 08:32:14慣れてしまってMatlabなしには生きられなくなる→
そういう学生が理系就職する→
高額な一般ライセンスを購入→
Mathworksｳﾏｰ→
プログラミングができない学生を再教育しなければならない企業はﾏｽﾞｰ
という構図なんだが、判っててもはまってしまう。
Matlabは便利だからとりあえずの結果を出すには最適だものな。
スレ違いｽﾏｿ。
0474デフォルトの名無しさん
2005/10/25(火) 16:11:110475デフォルトの名無しさん
2005/10/25(火) 18:50:30ユザー名bayes05
パスワードthomas
↑この利子率(TB3M)の時系列データについて
①データをtextreadで読み込む
②横軸に時間、縦軸；利子率 をとったグラフを書く
③corrで自己相関係数を求める
この３つをやるプログラムがわかりません。初歩的な質問ですみませんが、教えてください。
0476デフォルトの名無しさん
2005/10/25(火) 19:14:030477デフォルトの名無しさん
2005/10/25(火) 23:45:24しかし、単純にfft(data)とすると、虚数のデータが吐き出されるようで、うまくいきません。
どうすればパワースペクトルを求められるでしょうか。
初心者にアドバイスお願いします。
0478デフォルトの名無しさん
2005/10/26(水) 02:18:09アルゴリズムの中身を知りたいのでしたら，信号処理の本を見たほうがいいかと思います。
matlabでの離散フーリエ変換は組み込み関数ifftでできます
>>477
>>特定の周波数帯域について積分しようとしています。
この意味がよくわかりませんが･･･
混同されているようですが，まず，時間ｔの関数をフーリエ変換すると，周波数ｆの関数になります。
これは振幅スペクトルです。パワースペクトルとは，振幅スペクトルの2乗になります。
データを data 、データ数を n とすると，
Y = fft(y,ｎ); %振幅スペクトル
Pyy = Y.* conj(Y) / n; %パワースペクトル
>>474,477
参考これのデータ解析とフーリエ変換より
http://www.cybernet.co.jp/matlab/support/helpdesk/r13/toolbox/matlab/ref/func_by__01.shtml
0479デフォルトの名無しさん
2005/10/26(水) 11:49:352次元ヒートマップ、あるいは3次元2DZ (ZはFrequency)を作りたいのですがどうしてもバイナリーの2次元マップしか作れません。
アドバイスお願いします。
（簡単な例）
1)行列を、A=(1,2,3; 1,4,5; 4,3,1), B=(1,4,2;1,5,3;3,4,1)で定義する。
各々の等しいマトリックスエレメントを(x,y)にアサインしてやりそれらを2次元展開する。
つまり、(X,Y) = (Aij,Bij)でプロットする。
この場合、(1,1),(2,4),(3,2),(1,1),(4,5),(5,3),(4,3),(3,4),(1,1)の9つの点ができる。
2)私はこれらをプロットするのに、マトリックスA,Bを一次元展開してやり、Plot(A(:),B(:))で展開してやりました。
この場合グラフはバイナリー表示です。
3)更にここから同座標の数をヒートマップか3次元で表したい。
(1)の例ならば、(1,1)の座標が3個あります。
これをグラフに反映させたいのです。
よろしくお願いします
0480デフォルトの名無しさん
2005/10/26(水) 14:29:04ありがとうございます。
振幅スペクトルとパワースペクトル、そのあたりのもやが晴れました。
ただ、
>>Pyy = Y.* conj(Y) / n; %パワースペクトル
ここまではたどり着いていたんですが、これをプロットすると、
両端にピークを持つグラフが出てきて、
横軸が周波数ではなくデータ数のようなんです。
単純にプロットの仕方がおかしいだけなのでしょうか？
それともデータの時間に関する情報が欠けてるのかなと思ったり。
データは、サーミスタで測定した温度変化です。
アドバイスお願いします。
0482デフォルトの名無しさん
2005/10/26(水) 22:48:49ttp://www.cybernet.co.jp/matlab/support/helpdesk/r13/toolbox/matlab/ref/func_by__01.shtml
で，グラフィックス－特殊なプロットでグラフ表示の関数を見てみましたが
ヒートマップというのは見受けられませんでした('A`)
0483477=480
2005/10/27(木) 14:18:11どなたかアドバイスお願いします。
0484デフォルトの名無しさん
2005/10/27(木) 16:49:01少しは自分で調べようYO
0486デフォルトの名無しさん
2005/10/27(木) 17:43:09式はだいたいわかるのですが・・・。
始めたばかりで、わからないことだらけですので誰かご教授お願いします。
0487デフォルトの名無しさん
2005/10/27(木) 18:27:410488デフォルトの名無しさん
2005/10/29(土) 19:28:290489デフォルトの名無しさん
2005/10/29(土) 19:40:59それはjavaのスレで聞くべきなんじゃ･･･
0490デフォルトの名無しさん
2005/10/29(土) 22:30:0649.5[Hz]-50.5[Hz]の成分を除去する、
といったことがMatlabではできるのでしょうか。
具体的な関数などを教えていただけませんか。
0491デフォルトの名無しさん
2005/10/29(土) 23:29:59バンドエリミネーションフィルタをかけてやればいいと思う。
filter関数かconv関数を使うことになるんだけど，あまり詳しくないです('A`)
参考
ttp://www.cybernet.co.jp/matlab/support/helpdesk/r14/toolbox/signal/basics1.shtml
0492デフォルトの名無しさん
2005/10/29(土) 23:33:110493デフォルトの名無しさん
2005/11/01(火) 14:59:23コマンドライの引数にmファイルを指定して実行することはできるでしょうか。
$ matlab hoge.m
こんな感じです。(このままだと、matlabが立ち上がるだけです)
octaveだとできるのですが。
やりたいことは、シェルスクリプトの中に組み込むことです。
0494493
2005/11/01(火) 18:25:43正:コマンドライン
matlab --helpとやってみても、matlabが立ち上がるだけだし。
コマンドラインオプションはないのかなぁ
0495470
2005/11/02(水) 16:48:40具体的にはある入力を1秒に1回のサンプリングを1000回分して
その入力の最小値と最大値を求めたいのです。
自分が考えたsimulinkのやり方は、
1秒ずつﾃﾞｨﾚｲさせてそれをMinとMaxのブロックにぶちこむというもの
なのですが、それでは、MinMaxの入力の枝が1000個になってしまい
間違っているような気がしまして質問させて頂きました。
環境：MATLAB6.5、xp
0496デフォルトの名無しさん
2005/11/02(水) 17:29:28ただし手元にLinux版がないので間違ってるかも
Windows版は今試してみたがＧＵＩ専門ぽい。
0497デフォルトの名無しさん
2005/11/02(水) 20:17:18つまり、1000サンプル分のデータの内の最大値と最小値を求めたいってことですか？
それでしたらバッファを使えば何とかなるかもしれません。
ＤＳＰブロックセットはインストール済みですか？
ちなみに、ディレイを1000個も使ったらメモリの関係で動作がすごいことになりそうです・・・
0498デフォルトの名無しさん
2005/11/03(木) 19:29:07うちのはWindows版みたいです。
>>497
おっしゃるとおりです！
実際には1000秒のシミュレーション中の変化する値の最小値と最大値を
随時更新していくという事をしたいのです。拙い説明ですいません。
今手元にMATLABがないので、DSPブロックセットの有無が分からないので
確認してからまた連絡させていただきます。
0499デフォルトの名無しさん
2005/11/03(木) 20:18:48シュミレーションは少なくとも1000秒以上行い，
1000秒ごとに最小値と最大値を更新する，ということでいいのかな？
それともシュミレーションは1000秒間行うだけ？
あと，最小値と最大値の更新はリアルタイムで行う必要がありますか？
ある程度の遅延があったら無理？
質問ばかりですが，一応確認させてください。
0500デフォルトの名無しさん
2005/11/03(木) 21:09:25シミュレーションが全体で1000秒ありまして、
最小値と最高値の更新は1秒ごとにしたいと思っています。
最小値と最高値の更新はリアルタイムであればベストですが、
ある程度の遅延があっても構いません。
＞質問ばかりですが，一応確認させてください。
いえいえ、すごくありがたい気持ちです。
よろしくお願いします。
0501デフォルトの名無しさん
2005/11/03(木) 21:33:26↓じゃまずいんですか？
シミュレーションパラメータで
シミュレーション時間1000秒
固定ステップ，ステップ幅1にして
.┌────┐. ┌────┐. ┌────┐
. | データ |＿＿|.zero-order|＿＿|. MIN |＿＿
. | 系列 | | hold | | |. |
. |＿＿＿＿__| |＿＿＿＿__|. | |＿＿＿＿__|
↑ |
sampling time:1. | ┌────┐
|＿|. MAX |＿＿
| |
|＿＿＿＿__|
0502デフォルトの名無しさん
2005/11/03(木) 22:00:58それと全く同じ事を最近試しましたがダメでした(；´Д`)
考え的にはあってるかなと思っていたのですが
MinとMaxのあとをスコープで覗いてみたところ全く同じ結果になってました。
なぜだぁ_|￣|○
0503デフォルトの名無しさん
2005/11/03(木) 22:10:45なら細工しないと無理そう
0504デフォルトの名無しさん
2005/11/03(木) 22:50:46minからの出力を、またminに
フィードバックさせたりできないんですかね？
zero-order holdからの出力とマージさせれば
それまでの最小値がだせますよね？
0505デフォルトの名無しさん
2005/11/03(木) 23:13:04ただ，出力を直につなぐとMATLABがフリーズするので
Delayブロックを噛ます必要があります。
0506デフォルトの名無しさん
2005/11/04(金) 17:35:30MATLABで株価チャートを表示させたいんだが
とりあえずチャートのX軸に日付を出力させたいのよ。
けどdateを使うと始まりの日付だけしか指定できないし、
ネット上から落としてきた株価データファイルの日付と一致させられない。
どっかにサンプルは落ちてないものかね？
救っておくれ。
0507デフォルトの名無しさん
2005/11/04(金) 18:35:25つーか喪前ら、サイバネットのドキュメントくらい見ようYO
次スレのテンプレに追加した方がいいかな・・・
R13
ttp://www.cybernet.co.jp/matlab/support/helpdesk/r13/toolbox/matlab/learn_matlab/gs_colle.shtml
R14
ttp://www.cybernet.co.jp/matlab/support/helpdesk/r14/toolbox/matlab/learn_matlab/ch1intro.shtml
問題は次スレまで相当かかりそうだってことだけど('A`)
0508デフォルトの名無しさん
2005/11/05(土) 10:44:37PDFドキュメントの存在しか知らんかったよ
どうもありがとう
未だにMATLABのトップページからどの道筋を行けば辿り着けるのか分からないが
ドキュメントの内容は充実しているくせに、サイト自体が使いにくい罠
0509デフォルトの名無しさん
2005/11/05(土) 10:54:02技術サポートのページばかり通ってたよ
0510デフォルトの名無しさん
2005/11/08(火) 21:49:08guideを使ったGUIプログラムをスタンドアローンアプリ
にしたいのですが、
（C/C++のグラフィックライブラリを使うのか、
あるいはそうでないのかがよく分かっていませんが）
mcc -m FILENAME
mcc -B splcpp FILENAME
だと生成されませんでした。正しい操作はどうでしょうか？
0511デフォルトの名無しさん
2005/11/08(火) 21:55:25単純なコマンドウィンドウに表示させるソースを上記で
処理すると、きちんとDOS窓表示のexeができます。
ちなみに、上記GUIソースを-mオプションで処理すると
膨大なCソースが生成され、最後はエラーで止まります。
（エラー内容は必要であれば後日貼ります）
これらCソースを無理やり、VCでリビルドをかけても
当然ながらexe生成はできず・・という状況です。
初心者です。ネットでそれなりに調べ、マニュアルも目を通しています。
ぜひ皆さんのご助言をたまわりたく。
0512デフォルトの名無しさん
2005/11/09(水) 20:43:06f0 = 440; %[Hz]
f1 = 880;
Nsec=2; %[sec]
Sf = 44100%[Hz]
Np= Sf*Nsec;
y =zeros(1,Np);
for i =1:Np
y1(i) = 100*sin(2*pi*f0*i/Sf+(1/3*pi));
y2(i) = 50*sin(2*pi*f1*i/Sf+(4/3*pi));
end;
y3 = y1 + y2;
x(1:Np)=(1:Np)/Sf;
axis([0 200 -150 150 ])
hold on;
box on;
plot(y1(1:200),':');
plot(y2(1:200),':');
という、波形のグラフを作ったのですが、
このときｘ軸（時間軸）には数値が現れるわけですが、
これを時間に置きなおすにはどうしたらよいのでしょうか？
plot(y1(1:200),':');を
plot(x,y1(1:200),':');にするとx=0のところに線が書かれてしまうのですが・・・
0513デフォルトの名無しさん
2005/11/09(水) 20:44:41plot(x（1:200）,y1(1:200),':');でした。
0514デフォルトの名無しさん
2005/11/09(水) 21:14:20自己解決
0515デフォルトの名無しさん
2005/11/09(水) 23:24:0005/11/09といったようにスラッシュで区切ってあるデータを
テキストファイルからloadするにはどうしたらよいのでしょうか？
一列目には日付、二列目・三列目には数字が入力されているテキストです。
05/11/09 100 200（←こんなかんじ）
load test.txt; だけでは05/11/09のうちの最初の5しか読み取って
いないようでした。
0516デフォルトの名無しさん
2005/11/10(木) 00:22:00デリミタ指定したデータの読み込みができたはず。
0517デフォルトの名無しさん
2005/11/13(日) 15:48:38スタートアップ時にデフォルト値を設定できる項目はないのでしょうか？
set(0,'defaultAxesFontSize',14);
set(0,'defaultAxesFontName','Helvetica');
set(0,'defaultTextFontSize',14);
set(0,'defaultTextFontName','Helvetica');
0518デフォルトの名無しさん
2005/11/14(月) 23:58:07パワースペクトルを出した後,ある周波数帯域例えば
10～20HZ間の合計パワーを求めるにはどうしたらよいのでしょうか？
0519デフォルトの名無しさん
2005/11/15(火) 00:28:55コマンドラインでhelp sumって打ち込めば説明が出てくるYO
0520デフォルトの名無しさん
2005/11/15(火) 00:37:14他のスレに貼り付けると□□□□■□□□□■□□□
□■■隠された文字が■浮き出るよ□□■■■□□□□
http://pc8.2ch.net/test/read.cgi/tech/1131282679/l50
□□□■■□■■□□□□□■■□□■■□■□□
237 ：デフォルトの名無しさん ：2005/11/14(月) 23:38:29
さあ話をそらしにかかりました□□□□■□□□■■□
□■■□□■□□□□■■□□■□□□□□■■□
0521デフォルトの名無しさん
2005/11/15(火) 00:38:340522デフォルトの名無しさん
2005/11/15(火) 15:01:58俺も知りたい
0523デフォルトの名無しさん
2005/11/17(木) 02:51:370524デフォルトの名無しさん
2005/11/17(木) 03:05:17複雑なプログラムだと結果出るまでに数時間，数日間ってのもざらだし。
その分書くのは簡単なんだが。
0525デフォルトの名無しさん
2005/11/17(木) 16:19:050526デフォルトの名無しさん
2005/11/20(日) 14:26:42ttp://www.cybernet.co.jp/matlab/support/helpdesk/r13/toolbox/matlab/ref/ode15s.shtml
stiif=堅いらしいけど
ステッフの方が概して精度が低いみたい。
ノンスティッフシステム=例:外力のない剛体の挙動なんて記述もあるけど
ググッても良く分らなかった・・・
0527デフォルトの名無しさん
2005/11/21(月) 00:44:11漏れは全く知らんけどとりあえず
stiffｿﾙﾊﾞｰ(陰的方法)、nonstiffｿﾙﾊﾞｰ(陽的方法)とかあった
0528526
2005/11/21(月) 20:13:21有難う。
硬度比とか良く分らなかったけど、
当然だけどソルバによって計算結果が違う割には、
どのソルバが適当なのか理解出来てないって気持ち悪いな・・・
0529デフォルトの名無しさん
2005/11/21(月) 23:19:280530デフォルトの名無しさん
2005/11/24(木) 12:08:03MATLABでエクセルからデータ(たとえば、時間、変位X、変位Y、角変位φのような)を読み込んで
時間が変化していく際の、状態を(X,Y)とφでグラフにプロットして再現したいと考えています。
最初for文でまわしてplotでhold onとoffをつかって表現してたのですが、
かなり時間がかかってしまいアニメーションとしてみたかったのでうまくいきませんでした。
それでmovieinやgetframeというコマンドでアニメーション化をしてしまおうと試みました。
普通のplotで一つのグラフだけだとうまくやることはできました。
ただ今回は片方が極座標系のグラフを使いたいのでpolarのコマンドを使っていて
subplotで２つのグラフ表示をさせています。
両方同時に再生したいと考えているのですが、subplotを使った場合のアニメーション化が
どうしてもうまくいきません。
これは可能なんですかね？もしご存知の方いらしたら教えてください。お願いします。
0531デフォルトの名無しさん
2005/11/24(木) 13:55:47よくわからんが、XDataとYDataをsetで変更すれば、毎回書き直すより
少しは速くなる。
0532デフォルトの名無しさん
2005/11/28(月) 16:10:000533あ
2005/11/28(月) 16:23:34教えてください。
今サンプリング周波数が16000ｋＨｚなのですが
これを80000ポイントで5秒にしてます。
これを時間軸で表したいのです。
お願いします。
0534デフォルトの名無しさん
2005/11/28(月) 17:11:27figure
0535デフォルトの名無しさん
2005/11/28(月) 22:13:11なっが～いシリアルつかえば
コピーガードも何も関係ないよ
0536デフォルトの名無しさん
2005/11/28(月) 22:15:27> いまのStudent Versionには機能制限はないよ。
> 各ToolBoxもちゃんと使えるし。
> 起動時にプログラムCD-ROMを入れておかないといけない、
ver7はそうなってるんか？
ver6までは長い汁アル使えばそんなことしなくてもよかったし
0537デフォルトの名無しさん
2005/11/28(月) 22:28:32> 学生用のライセンスは機能制限なし（つまり何でも使える）なのに安価→
> 慣れてしまってMatlabなしには生きられなくなる→
> そういう学生が理系就職する→
> 高額な一般ライセンスを購入→
> Mathworksｳﾏｰ→
> プログラミングができない学生を再教育しなければならない企業はﾏｽﾞｰ
> という構図なんだが、判っててもはまってしまう。
> Matlabは便利だからとりあえずの結果を出すには最適だものな。
今業界ではオープン化、オープンソース化の動きが出ている。
その動きを予想して漏れはオープンソース版MATLABライクな言語
R言語に興味を示している。グラフ表示にはGNU plot.
ほかの言語としてはC++/Javaで応戦なり。
今まで研究室にいたメンバーは授業でいくら教わっても皆C/C++を
理解できなかった。そこで彼らは研究室に入ってからMATLABという言語に
出会って大半の者が「これが俺の求めていたものだ」とか言いだして感動しちゃってる。
漏れはその流れに逆らって苦労して時間をかけてC++/Javaに手を出した。
自作しなければならないものが多く、手探りがしばらく続き大変だったが、
今では自動化ツールにより開発効率がMATLABを使うよりも確実に向上している。
なんとかMathworksの魔の手から逃れることができた、漏れもすれ違いｽﾏｿ
0538デフォルトの名無しさん
2005/11/28(月) 22:28:37Javaスレでも質問したな？
だが誰も答えられなかったな？
漏れがちょっと答えたがRuntimeクラスを
使う以外わからなかったのだ。
ただしMATLABのマニュアルには
ＭＡＴＬＡＢからJavaオブジェクトを
呼び出せると書いてあるのをみたことがある。
やりたいことはそれか？
Javaコードがすでに存在しているからそれを
どうにかしたいのか、それともMATLABコードが
すでに存在するからそれをどうにかしたいのか、
どっち？
それとも、webで何かをしたいのか？
一体何をしたいのかかかないと
質問に答えにくい。
JavaからMATLABを呼び出す方法については
なにやら>>496がヒントを出しているみたいだぞ。
こいつをJavaのRuntimeクラスをつかって。
0539デフォルトの名無しさん
2005/11/28(月) 23:08:02フーリエ変換が趣味です。
http://hobby8.2ch.net/test/read.cgi/hobby/1133005088/
0540デフォルトの名無しさん
2005/11/29(火) 02:50:25もうフーリエ変換は(ry
0541デフォルトの名無しさん
2005/11/29(火) 17:41:33計算量では、WTはFFTに太刀打ちできんよ。
0542デフォルトの名無しさん
2005/12/01(木) 23:53:550543デフォルトの名無しさん
2005/12/07(水) 16:09:14Y = A;
のかわりに
[X, Y] = ***(A);
という構文があったと思うのですが、***に入るコマンドを教えてください。
0544デフォルトの名無しさん
2005/12/07(水) 16:55:58どのように入力したら良いですか？
0545デフォルトの名無しさん
2005/12/09(金) 00:17:18*** = 自分で作った関数名
0546デフォルトの名無しさん
2005/12/11(日) 14:34:14まじ、悩み中です
0547デフォルトの名無しさん
2005/12/11(日) 15:04:08どういった条件の下でやるのか書いてもらわんことにはね・・・
0548デフォルトの名無しさん
2005/12/11(日) 21:34:25使いたいんですが。。
0549デフォルトの名無しさん
2005/12/11(日) 21:49:31ttp://www.cybernet.co.jp/matlab/support/helpdesk/r14/toolbox/filterdesign/adaptive.shtml
アダプティブフィルタはそれほど大した式で記述されてるわけじゃないし，
5行程度のコードでできるＹＯ
でもディジタル音声ってことは，1，0のバイナリデータだよねぇ･･･
アナログ音声ならアダプティブフィルタ使うのはわかるけど・・・
ディジタルにのったノイズはどうしようもない気がする。
まぁ，間違ってるかもしれないけど一度Ｄ／Ａ変換する必要があると思われ。
0550デフォルトの名無しさん
2005/12/12(月) 21:40:12定在波観測実験にMATLABを用いようとしたのですが
信号データを作成して再生と録音コマンドを実行しても
再生が終わった後録音が始まってしまい、目的の結果が得られませんでした。
スピーカーから音を出すと同時にマイクで録音するためにはどうすればよいでしょうか？
0551デフォルトの名無しさん
2005/12/13(火) 00:22:41「その文字はMATLABで使えません」とか言われました。
0552デフォルトの名無しさん
2005/12/14(水) 08:25:36>オープンソース版MATLABライクな言語R言語
MatlabにはOctaveのほうが近くね？
R言語はS言語だし。
そういえばMaximaというMathematicaっぽいやつもあった。
0553デフォルトの名無しさん
2005/12/19(月) 17:00:02[ly lx] = size(S1)
m =2; i=0; j=0
for i=m+1:ly-m
for j=m+1:lx-m
a=S1(i-m:i+m,j-m:j+m)
b1=mean2(a)
b2=std2(a)
varx=((b2^2)+(b1^2))/(((b2/b1)^2)+1)-(b1^2)
b3=varx/(((b1^2)*(b2/b1)^2)+varx)
FS(i-m,j-m)=b1+b3*(S1(i,j)-b1)
end
end
S1が元の画像、FS1が処理後の画像です。
なんかやけに計算に時間がかかるんですが、プログラムに無駄があるんでしょうか？
0554デフォルトの名無しさん
2005/12/19(月) 17:16:59以下のような行列で、a,b,cそれぞれの系列について総当たりで単回帰分析を行いたいのです。
a→b,a→c,b→c,b→a,c→b,c→aというふうに。
--------------
data =
a b c
0 1 2
1 5 4
2 3 7
2 4 5
--------------
(ちなみに実際のデータはdouble配列で、a,b,cは便宜上付加しているだけです)
現状では、
for i = 1:3
for i = 1:3
k(i,j) = regress(data(:,i),data(:,j));
end
end
というようにforを駆使して(笑)求めているんですが、ループを使わずに
記述する方法は無いでしょうか?
・単回帰分析に特化されたregress関数
・それぞれの列に対して処理を行うことを明示する関数
などがあれば出来そうだと思い、探しているんですが見つけられません。
素人質問で申し訳ありませんが、よろしくお願い致します。
0555554
2005/12/19(月) 17:18:53for j=1:3
の誤りでした。
0556必死
2005/12/20(火) 12:38:020557デフォルトの名無しさん
2005/12/20(火) 14:11:20ある程度できたんですけど、因果律を満たしていないのか、うまく出来ません。
おそらく入力信号に遅延を付加していない体と思うのですが、その遅延の方法
自体分からず困ってます。どなたか入力信号に遅延を付加する方法を教えてください。
ちなみに、入力信号は、3,000点のランダム信号。タップ数は70でやっています。
0558デフォルトの名無しさん
2005/12/20(火) 20:39:48fft
0559デフォルトの名無しさん
2005/12/22(木) 22:36:390560デフォルトの名無しさん
2005/12/23(金) 02:25:13ものすごい初歩的な質問なのかもしれませんが、仮想メモリについてです。
画像処理を行っていまして、480×640×600のマトリックスを作りたいのですが、「out of memory」のメッセージが出力されてしまいます。
マトリックスサイズは480×640ピクセルの画像が時間軸方向に600枚ある画像を解析していこうと思うのですが、できれば、1000枚くらいまで解析したいのです。
ところが、480×640×200くらいまでのマトリックスしか作ってくれません。
画像はグレースケール画像なのですが、画像マトリックスをAとすると、A=zeros(640,480,600);と宣言して、そこに実際に撮像した画像のグレースケール値を突っ込んでやろうと思うのですが、必ず「out of memory」が出力されてしまうのです。
PCには物理メモリが１G載っていて、仮想メモリを５G当てているのですが、タスクマネージャで見てみると、コミットチャージで１Gちょっとしか使ってくれないのです。
MATLABの環境はR12とR14のどちらかを使用しています。
MATLABのほうで使ってくれるメモリの量は決まっているのでしょうか？もし、それを解除できるような方法がありましたら、どなたかご教授願います。
0561デフォルトの名無しさん
2005/12/23(金) 10:01:32メモリ上限を増やす方法は知らないが、
グレイスケールなら
A = repmat(int8(0), [640 480 600])
にしてみては?
0562デフォルトの名無しさん
2005/12/23(金) 14:16:54http://www.mathworks.com/support/tech-notes/1100/1106.html
0563デフォルトの名無しさん
2005/12/26(月) 00:57:42ありがとうございます。参考にさせていただきたいと思います。
0564デフォルトの名無しさん
2005/12/29(木) 02:09:39Octave
http://www.octave.org/
基本文法は同じですが関数名が違うものも多いので、
Matlab用のM-ファイルがそのまま動かないことも多いです。
0565デフォルトの名無しさん
2005/12/29(木) 02:12:50Signal Processing Toolboxでよければ、
periodgramをはじめ、pwelch、pcov、pyulear、pmusic、peig...
多数のPower Spectrum Density Estimationができます。
0566565
2005/12/29(木) 02:14:53私はこれを参考にして音響信号をプロットしています。
http://www.mathworks.com/support/tech-notes/1700/1702.shtml
0567デフォルトの名無しさん
2006/01/01(日) 00:15:520568デフォルトの名無しさん
2006/01/02(月) 12:14:43できるけど、起動時にCD-ROMが必要になるので、いつも持ち歩いてないといけないぜ。
0569デフォルトの名無しさん
2006/01/02(月) 14:03:33Matlab Student Version R14は起動時のCD-ROMが不要になっているらしいです。
あらためて、R14では、ﾃﾞｽｸﾄｯﾌﾟと2台目のﾉｰﾄﾌﾞｯｸとの2台にｲﾝｽﾄｰﾙして使うことは出来ますか？
0570デフォルトの名無しさん
2006/01/03(火) 01:29:48自分はR13だけど、R14を使ってる友人は起動時にいつもCD-ROMを入れてる。
0571デフォルトの名無しさん
2006/01/03(火) 16:15:15すなわち、FFTや窓関数について、実際のスペクトラムスコープ
との動作と照らし合わせて説明している書籍（和書、洋書)とか、
測定器（Agilentなど)のマニュアルとか何処かにありませんか？
デジタル信号処理の書籍を読めば良いことかもしれませんが、
Simlinkを使う上で、アナログスペアナの概念で捉えると勘違い
してしまう部分など・・・
0572デフォルトの名無しさん
2006/01/03(火) 16:25:43だとすると、Matlab Student Version R14 SP3 以降からなのか、
つい最近から、ｲﾝｽﾄｰﾙ時にActivation手続きが必要となったようです。
何方か、最新版Matlab R14 SP3でのｲﾝｽﾄｰﾙ時のActivation手続きに関連して、
2台目へのｲﾝｽﾄｰﾙが出来るかどうか知りませんか?
0573デフォルトの名無しさん
2006/01/06(金) 10:46:27メーカーに聞くといいよ。
0574デフォルトの名無しさん
2006/01/06(金) 14:45:59オススメのCPUってどれだろう？
バージョンは6.5でも7でもいい
0575デフォルトの名無しさん
2006/01/06(金) 15:46:530576デフォルトの名無しさん
2006/01/07(土) 13:32:14x=????????
polar(t,x);
これで結果が楕円になるような関数xが欲しいのですが
どなたか救いの手をよろしくお願いします。
0577デフォルトの名無しさん
2006/01/07(土) 19:27:28数学の本を嫁
0578570
2006/01/08(日) 00:59:44そういえばその友人はR14SP1を使ってる。
MathworksのページにはSP3からはアクティベーションが必要みたいなことが書いてあるね。
http://www.mathworks.com/academia/student_version/faq/
ここには「同時に2台までアクティベーション可能」って書いてある。
アクティベーション情報を別マシンに移動させることもできるみたいだよ。
0579デフォルトの名無しさん
2006/01/09(月) 20:50:13どのようにやれば良いのでしょうか？
0580デフォルトの名無しさん
2006/01/12(木) 16:43:070581デフォルトの名無しさん
2006/01/12(木) 19:33:14できるけど、最初から用意されてるわけじゃないから、自分で作ればいい。
ヒルベルト変換を使うのが簡単だよ。
0582デフォルトの名無しさん
2006/01/15(日) 20:18:29なんか同じようなことを経験した人いますか？
PowerBook1.5GHz OSX tiger メモリ768Mb。30Mb程度のものをloadするのに
３０分以上かかる。windows、linux上だと１０秒ぐらいなんですが。
0583デフォルトの名無しさん
2006/01/15(日) 20:22:35http://www.geocities.jp/imananji2321/
音楽室
0584デフォルトの名無しさん
2006/01/15(日) 20:54:47http://newsreader.mathworks.com/WebX?14@@.ef13759
0585デフォルトの名無しさん
2006/01/16(月) 15:46:28読みこんだファイルに対する積分(シンプソン則等)
は出来ないのでしょうか?
初心者なものでよろしかったら
アドバイスして下さい．
0586デフォルトの名無しさん
2006/01/16(月) 23:21:55ファイルを読み込む部分と、数値積分をする部分に分けて考え、
マニュアル読んでがんばれ。
0587デフォルトの名無しさん
2006/01/17(火) 12:13:31ありがとうございます。
参考にします．
0588デフォルトの名無しさん
2006/01/23(月) 16:54:51誰か分かる人いませんか？？help me!!
あと、シミュレーション画面をGUIで表示する方法も分かる方がいれば教えてください。
0589デフォルトの名無しさん
2006/01/23(月) 22:21:55Simulinkのこと？
どのように変更したいのかなぁ。
チューナブルパラメータだと、値を変えれば反映されない？
0590デフォルトの名無しさん
2006/01/24(火) 01:48:15Simulinkのことです。
チューナブルパラメータが何だかよく分からないんだけど、要はシミュレーション結果をグラフで表示してるんだけど、途中で値を変化させてそれにともなってグラフも変化してほしいんです。
グラフをリモコン（GUI）で自在に操ってる的な感じかな。あ、だからシミュレーション途中の動いてるグラフがGUIの画面中で表示させたいんです。
無理なんですかねぇ？分かりますか？
0591デフォルトの名無しさん
2006/01/24(火) 03:18:54やりたいことが漠然としすぎててﾜｶﾗﾅｽ
もうちょいまとめてください
0592デフォルトの名無しさん
2006/01/26(木) 00:57:170593デフォルトの名無しさん
2006/01/26(木) 08:16:24チューナブルパラメータというのは、シミュレーション中に変更しても反映されるパラメータ。
リモートコントロールとGUIも結びつかないのだが、Switchみたいなダブルクリックしたら
変化するブロックなら、簡単に作れるぞ。扇風機みたいな、切,弱,中,強みたいなブロックを
作って、ダブルクリックしたらONになり他はOFFになるってのも、ちょっとがんばれば可能。
グラフが変化ってのがイメージわかん。
0594デフォルトの名無しさん
2006/01/26(木) 20:31:17for i = 1:size(A,3)
A(:, :, i) = A(:, :, i).';
end
という処理をしたいのですが、ループを使わずに、かつ、大きな
インデックス行列を作らずにすむ方法はありますか?
0595デフォルトの名無しさん
2006/01/27(金) 02:25:50パワースペクトルを求めた後,プロットすると横軸が周波数,縦軸がパワーのグラフに
なりますが、このグラフのある周波数帯域だけを積分とかすることは可能なのでしょうか？
お願いします
0596デフォルトの名無しさん
2006/01/27(金) 18:46:25ここにある課題２がどうしても分からないです。。
MATLABをつかって解くらしいんですが。。。
0597デフォルトの名無しさん
2006/01/27(金) 18:58:27先生や先輩あるいは同じ学科の友人に聞いた方が
得られるものは多いと思いますよ。今後のためにも。
0598デフォルトの名無しさん
2006/01/27(金) 18:59:49マルチ乙
0599デフォルトの名無しさん
2006/01/27(金) 21:07:27解決しました。permuteですね。
0600デフォルトの名無しさん
2006/01/28(土) 01:38:56モデルウインドウのSimulation→parameters→Workspace I/OのLoad from workspace欄のInputにチェックを入れたいんですけど、エラーが消えません。
何が悪いのかさっぱり分かりません。。。エラーもいつも同じエラーというわけでもないし…。
どうすればいいか分かりますか？簡単な例でいいので教えてくれませんか？
あと、チューナブルパラメータってS-Functionしか使えないんですか？
0601デフォルトの名無しさん
2006/01/28(土) 03:16:10んー・・・
Simulink歴はどの程度ですか？
話を聞く限りでは初めて触るような印象を受けるのですが・・・
そうだとすると説明するのが逆にﾑｽﾞｶｼｽ(´･ω･`)
ttp://www1.bbiq.jp/kougaku/matlab_intro.html
ググって引っかかったページですが，例えば【１３】でのGainブロックの値を
ｍファイルではなくブロックパラメータに直接記入し，さらにシミュレーション中に他の値に変えれば
Scopeブロックに表示される結果も変わります。
わざわざワークスペースから引っ張ってこなくても，モデルのブロックの値を変えればいいと思うのですが・・・
おっしゃっているチューナブルパラメータってのは単にこういう変数ってことでいいのかな？
0602デフォルトの名無しさん
2006/01/28(土) 04:57:31ほぼ初心者ですね。
でも、担当の教授が先月亡くなりまして…。
しかも、先輩もMATLABはあまり分からないとのことで、helpが必要なんです。。。
確かにブロックパラメータに直接数値を記入すれば、シミュレーション中にも結果は変わるようですが、ブロックパラメータに数値でなく変数をいれ、ワークスペースで変数の値を変更しても結果はシミュレーションしなおすまで変わりません…。
なぜでしょうか？不可能なんですか？
0603デフォルトの名無しさん
2006/01/28(土) 05:47:38今手元にMATLABないから確かめられないけど，
シミュレーション実行中に変えられる値と変えられない値があります。
そのせいかも知れません。
どういった値を変化させたいんですか？
0604デフォルトの名無しさん
2006/01/28(土) 09:40:32パワースペクトルをどうやって求めたのかにもよるけど、
コンピュータ上（離散領域）での積分というのは加算（Σ）と同じことなので、
ただ必要な周波数帯ごとのパワーを足し算するだけで良いんじゃないかな。
ただ、デシベルになっているときは単純に加算できないので注意ね。
0605デフォルトの名無しさん
2006/01/28(土) 10:33:00変えられる値と変えられない値とは？例えば？
う～ん。どういった値って言えばいいんでしょうか？
例えばconstantのブロックパラメータに変数aを入れて、そのままscopeにつなぐ。
で、ワークスペースかなんかで初期値をa=1とかにしとく。
んで、シミュレーションをやると、scopeにはa=1のグラフがでますよね？
そのシミュレーション途中でワークスペースでa=2とかを打つと
scopeで出てくる値もa=2に変化してくれればありがたい。
それがならない…。
ブロックパラメータに直接１と入れて、シミュレーション中にブロックパラメータの値を
２とかに変えると、それはscopeに反映されて出てくる。
違いが分からん。
0606デフォルトの名無しさん
2006/01/28(土) 15:50:56おかげでブックマークがつながらないぽ。
卒論の時期に面倒なことしないでくれよー。
0607デフォルトの名無しさん
2006/01/28(土) 16:19:05サイバネットに資料があるので読んでみるといいと思います
ttp://www.cybernet.co.jp/matlab/support/manual/r14/toolbox/simulink/index.shtml
で，変数の件ですが
ttp://www.cybernet.co.jp/matlab/support/manual/r14/toolbox/simulink/ug/how_sim8.shtml#342208
にあるようにシミュレーション実行中の変数の変更は，シミュレーションを一時停止してやらないと無理そうです。
ですから，シミュレーション実行中にワークスペースで値を変化させても，
シミュレーションを一時停止し，再実行しないと反映しないと思われます。
>>変えられる値と変えられない値とは？例えば？
変えられる値の例で言えばCostantの定数
変えられない値はSamplingtimeかな？
この資料で言うチューニング可能な変数という奴でしょう。
0608デフォルトの名無しさん
2006/01/28(土) 16:20:550609デフォルトの名無しさん
2006/01/30(月) 05:07:540610デフォルトの名無しさん
2006/01/30(月) 19:01:44変数をそんなに扱えるわけじゃないし
0611デフォルトの名無しさん
2006/02/01(水) 02:23:24>んで、シミュレーションをやると、scopeにはa=1のグラフがでますよね？
>そのシミュレーション途中でワークスペースでa=2とかを打つと
>scopeで出てくる値もa=2に変化してくれればありがたい。
ていってたから
例えだってのは分かってるけど
どの程度のことをやりたいのかがよくわからない
0612デフォルトの名無しさん
2006/02/01(水) 19:13:46A = [1 2 3 4 5 6 7 8 9];
これを3列に分けたいんですけど、何かいい方法ないでしょうか？
つまり、ベクトルAから、↓こんな行列Bをつくたいのですけど・・
B = [1 2 3; 4 5 6; 7 8 9]
環境は、MATLAB 7.0.1 です。よろしくお願いします。
0613デフォルトの名無しさん
2006/02/01(水) 23:46:350614デフォルトの名無しさん
2006/02/02(木) 05:29:29このとき、Signal Processing Toolboxにある関数pyulear.mを使用する。
ってどうやって使用するのですか？
0615デフォルトの名無しさん
2006/02/02(木) 08:48:55reshape
0616デフォルトの名無しさん
2006/02/03(金) 01:42:50伝達関数の計算をしたいんですけど、
inv(sI-A)とかやりたいんですが、
当然s*Iじゃ受け付けてくれず、
tf('s')とかした後に、s*I-Aとやると、おそらく行列ではなくなって、
逆行列とか計算できなくなるのですが･･･
いったいどうしたらいいんでしょうか？
手計算じゃきついです･･･。
0617デフォルトの名無しさん
2006/02/03(金) 01:54:46Aは適当な正方行列，ｓはスカラとしてIはそのままIとして使ってないよね？
単位行列はたしかeye関数で生成すると思ったんだけど・・・
0618デフォルトの名無しさん
2006/02/03(金) 02:24:43Iもちゃんと単位行列です。
sはラプラス変換のときのsです。
0619デフォルトの名無しさん
2006/02/03(金) 03:03:17ラプラス変換のｓ？
よくわからんけど，スカラー量でいいんだよね？
とりあえずサイバネットを見てみてくれ
MATLABについて詳しく書いてある
漏れはもう眠いから今日は寝る('A`)ﾉｼ
0620デフォルトの名無しさん
2006/02/03(金) 03:07:58答えてもらったようで、なんも答えてもらってない･･(´ д｀)
0621デフォルトの名無しさん
2006/02/03(金) 03:12:21結局ｓについて何も答えてないじゃないか
勝手にするがよろし
0622デフォルトの名無しさん
2006/02/03(金) 04:18:260623デフォルトの名無しさん
2006/02/03(金) 23:13:09tf('s')で変換するとinv()は逆行列ではくて逆システムを求める演算になります。
Mathematicaっぽい機能を使いたいときはsymbolic math機能で。
s=sym('s')
I=eye(2)
A=[ 0 1 ; 2 3]
X=inv(s*I-A)
0624623
2006/02/03(金) 23:18:49>623で書いたやり方はMatlabのインストールの仕方やネットワーク構成が複雑だと
エラーになることがあるよ(ライブラリの整合性の関係で)。
スタンドアローンで動かしているとかOSをサーバ/クライアントで統一しているとかなら問題なし。
0625デフォルトの名無しさん
2006/02/05(日) 21:54:50行列の要素で一定の値以上のものだけを
ほかの行列にしまっていきたいのですが、よい方法はないでしょうか？
できればfor文を使わない方法で教えていただきたいのですが…
よろしくお願いします
0626デフォルトの名無しさん
2006/02/05(日) 22:21:21find関数使ったらいいんじゃない?
0627デフォルトの名無しさん
2006/02/05(日) 22:26:22find関数について調べてみます。
0628625
2006/02/05(日) 22:38:15リンクからnonzerosという関数を見つけたのですが
値が0以外ではなくこちらから指定した値以上のものを抽出したいのですが…
たびたびすいません。ご助言お願いします。m(_ _；)m
0629デフォルトの名無しさん
2006/02/05(日) 23:00:12A 任意の行列
x 閾値
find(A>x) A>xの要素の"インデックス"
これで動くと思うけど。。。
0630625
2006/02/05(日) 23:11:02”インデックス”が抽出できるのはわかりました。
ただ、そこから中身をどう持ってくるのかわからないんです…
0631625
2006/02/05(日) 23:13:44ありがとうございましたm(_ _)m
0632デフォルトの名無しさん
2006/02/05(日) 23:48:38解決して何よりっす。
0633デフォルトの名無しさん
2006/02/06(月) 08:42:00こういう人って最終的にどう解決したか書かないんだよねぇ。
最近そういうの多く見てて書きたくなっただけ。
0634デフォルトの名無しさん
2006/02/06(月) 22:44:46ある行列をplotしたいんですけど、行列が長くて一回で変数に入りきらないので、
for文で回して、plotしたいのですけれど、できますか？
例えば、こんな感じです。
for j=1:10
A = fread(fReadId, [2 10000], 'real*8');
t = j:10000
plot(A(1, :), A(2,:)); % <= ここでちゃんと続きにplotできるようにしたいのですが、
end
よろしくお願いします。
0635デフォルトの名無しさん
2006/02/07(火) 02:23:08自分だったら、for文の前に
plot(0, 0, 'k.', 'visible', 'off');
hold on;
と書くかな。
end文のあとにhold offも忘れずに。
0636634
2006/02/07(火) 13:18:53なるほど、ありがとうございます。
ですが、plotしたときに、 plot(A(1, :), A(2,:), 'b.-');とか線でつなぎたいのですが、
そこまではやっぱりできないですかね？
0637デフォルトの名無しさん
2006/02/07(火) 23:26:59setでXDataとYDataを直接変更する。
動作確認してないがだいたいこんな感じ。
% 1回目
h = plot(A(1, :), A(2,:), 'b.-');
% 2回目以降
set(h, 'XData', [get(h, 'XData') A(1,:)], 'YData', [get(h, 'YData') A(2,:)])
0638デフォルトの名無しさん
2006/02/08(水) 15:33:01初期バッファが割り当てられますが、
そのバッファの中身を覗くことってできますか？
0639デフォルトの名無しさん
2006/02/08(水) 16:24:05なければ0じゃないかな
0640デフォルトの名無しさん
2006/02/09(木) 03:09:100641デフォルトの名無しさん
2006/02/09(木) 03:10:53どうすればいいんでしょうかって･･･（ ´Д｀；)
少しは自分で考えれ
話はそれからだ
0642デフォルトの名無しさん
2006/02/10(金) 00:22:25設定されているのでしょうか？分かる方いましたら教えてください
0643デフォルトの名無しさん
2006/02/10(金) 00:27:07コマンドでhelp fft
0644デフォルトの名無しさん
2006/02/10(金) 00:52:20サイバネットのページみても「お問い合わせください」だったと思うんだが。
0645デフォルトの名無しさん
2006/02/10(金) 09:15:550646デフォルトの名無しさん
2006/02/10(金) 10:04:41手元に無いので自信は無いし、やったこともないが、get_paramの
ヘルプを見ると、ヒントがあるかも。
モデルごとに覚えられるパラメータ(デフォルトフォント等)だと、
変更前と変更後でmdlファイルのdiffを取ったら一発。
0647デフォルトの名無しさん
2006/02/10(金) 22:06:3035万くらい
正確な値段は会社にカタログがあるが土日休みだ
0648デフォルトの名無しさん
2006/02/11(土) 08:21:12パーソナルライセンスだったらそんなもんだね。
ネットワークライセンスだったら、3倍くらい？
特に、個人と企業の差は無かったと思う。
たくさん買えば、少し割り引きしてくれるけど。
0649デフォルトの名無しさん
2006/02/13(月) 02:01:08シンプルにやるにはどうしたらいい？
a = randperm(1000);
b = a(1:6);
でもいいけど randperm の無駄さ加減が目に余る。
0650デフォルトの名無しさん
2006/02/13(月) 04:44:17関数の引数の初期数はどのように設定するのでしょうか？orできないのでしょうか？
function a = hoge(b, c = 1)
のようにやったらエラーがでました。
もしくは、その変数が定義されていなかったら。と制御できるのでしょうか？
function a = hoge(b , c)
if undef(c)
c = 1;
end
のような雰囲気で。どのようにしてデフォルト値を設定するのでしょう？
ライブラリ(toolbox)などでは引数を省略すると、、
ある値になるというのが普通にあったような気がするので、できるのだとは思うのです。
どなたかよろしくお願いします。
0651デフォルトの名無しさん
2006/02/13(月) 06:15:480652デフォルトの名無しさん
2006/02/13(月) 20:31:36nargin
詳しくはライブラリ関数のソース見ろ
0653デフォルトの名無しさん
2006/02/16(木) 21:53:19適当な関数が見つかりません教えてください。
0654デフォルトの名無しさん
2006/02/17(金) 01:13:060655デフォルトの名無しさん
2006/02/24(金) 22:20:18stateflow 解説 使い方 チュートリアル とかぐぐってもなかなかいいのが見つからなくて
0656デフォルトの名無しさん
2006/03/01(水) 01:39:51どなたか、どうか、この件について教えていただけないでしょうか？
よろしくお願い致します。
0657デフォルトの名無しさん
2006/03/02(木) 14:18:25w_log(:,:,t) = w;
という代入をすると、ものすごく重い。
ちなみにwのサイズは(n,n)。
w_logはサイズ(n,n,t_max)で事前に初期化済み。
どう書くのが一番速い？
0658デフォルトの名無しさん
2006/03/02(木) 20:40:27スワップしてるんじゃないか?
t_maxを物理メモリのサイズに抑えて、こまめにセーヴ・初期化
0659デフォルトの名無しさん
2006/03/03(金) 12:07:40ではないみたい
どうやら、global変数でw_logを持たせると代入も速いが、
引数渡しでw_logを持ち回ると、代入時に劇遅になるみたい。
0660デフォルトの名無しさん
2006/03/03(金) 12:23:26matlab はたしか、
関数内で変化がない場合はアドレス渡しだけやって、
変化がある場合は、値渡し（コピー）するはず。そのせいだと思う。
0661デフォルトの名無しさん
2006/03/03(金) 12:27:090662デフォルトの名無しさん
2006/03/06(月) 17:41:52このソフトをフリーでダウンロードすることってできますか？
初歩的な質問ですいません。。
0663デフォルトの名無しさん
2006/03/06(月) 22:25:22海外ならMathworks、日本国内ならサイバネット・ジャパンに問い合わせてください。
学生版なら1万円くらい、正規版はたぶん40万円くらい。
フリーウェアでMatlabに似たような文法のソフトがいくつかありますが、
なかでもOctaveはお勧めできます。
0664デフォルトの名無しさん
2006/03/10(金) 00:54:240665デフォルトの名無しさん
2006/03/10(金) 07:15:02Real-Time Workshop
0666デフォルトの名無しさん
2006/03/11(土) 14:06:01最初にDOS窓が開いてしまう。かなり見栄えが悪いです。このDOS窓を開かないよ
うにする方法を知っている方がいたら教えて下さい。
0667デフォルトの名無しさん
2006/03/11(土) 19:05:03Linux使う
0668デフォルトの名無しさん
2006/03/11(土) 20:32:04Xp/2000の場合だが
exeのｼｮ-ﾄｶｯﾄを作成して、その属性を開き、実行時の大きさを
最小化すれば殆んど目立たんよ。
0669デフォルトの名無しさん
2006/03/11(土) 22:09:54Thanks
>>668
Thanks 確かにタスクに引っ込みました。
0670デフォルトの名無しさん
2006/03/15(水) 08:57:020671デフォルトの名無しさん
2006/03/15(水) 10:21:10でもって、test1.mを（もしあれば）開けるけど、
逆に、エディタで開いているファイルを閉じるコマンド
（マウスを使わない）
は、ないでしょうか？
0672デフォルトの名無しさん
2006/03/15(水) 23:13:01すべて同じ処理させたい場合、その共通処理用のmファイルを
別に作成しなければならないのでしょうか？
GUIプログラミング作成時に自動生成されるmファイル内に、上記
共通処理用のfunctionを作成できないものなのでしょうか？
ご教授お願いいたします。
0673デフォルトの名無しさん
2006/03/16(木) 17:57:03CTRL+F4
で現在編集中のファイルが閉じる
0674デフォルトの名無しさん
2006/03/18(土) 13:05:53でもこれだと、一旦meditorに行ってからの操作になるよね。
MATLABのコマンドウィンドウから閉じる方法が知りたいんだけど。
スクリプトに書いて、最初に開いているmファイルを閉じたいんだが、
できないのかな？
0675http://www.vector.co.jp/soft/win95/util/se072729.html
2006/03/18(土) 20:14:05もしくは64bitにネイティブ対応したテキスト置換ソフトありますか？
0676672
2006/03/22(水) 23:11:430677デフォルトの名無しさん
2006/03/24(金) 08:49:38同じような症状の人はおりませんでしょうか？
0678デフォルトの名無しさん
2006/03/25(土) 09:53:04インストールできません？
こっちはインストールはできたんだけど、立ち上がらないんだよね。
ライブラリがどうのとか文句言われて。
妙な方法で解決したけど、なんだかなという感じ。
次のリリースで直ってるといいんだけど。
0679677
2006/03/25(土) 13:27:09立ち上がりませんでした。もしよかったらその解決方法を教えていただけないでしょうか？
0680デフォルトの名無しさん
2006/03/25(土) 17:29:52「だれでもわかるMATLAB即戦力ﾂ-ﾙﾌﾞｯｸ」よくできてる。
今まで、初心者向けの和文ではｻｲﾊﾞﾈｯﾄの講習用ﾃｷｽﾄに勝るものなかったからね。
0681デフォルトの名無しさん
2006/03/25(土) 18:13:49matlabを何か別のlinuxマシンに入れて、nfs経由でfedora core 5のマシンから
そのディレクトリにアクセスできるようにすると立ち上がります。
システムを入れたディスクにmatlabを入れてしまうと、
libwmlaback.soがどうのといって立ち上がりません。
nfs経由だと立ち上がるのに、そのまま/usr/localなんかにコピーしてくると立ち上がりません。
ファイルは全くいっしょのはずなのに、どうなっているんでしょう。
ちなみにバージョンはR2006aです。
0682デフォルトの名無しさん
2006/03/26(日) 18:37:52http://www.math.mcgill.ca/loisel/octave-workshop/
これのもう一つの特徴はMinGW上でコンパイルされている事です
Cygwin無しで動くのですね
0683デフォルトの名無しさん
2006/03/29(水) 23:34:50何を勉強するべきでしょうか？ C＋＋などを勉強すればいいんですか？
0684デフォルトの名無しさん
2006/03/29(水) 23:49:42いや煽りとかじゃなくて
0685デフォルトの名無しさん
2006/03/30(木) 00:37:49具体的課題があるなら習得は速いよ もち努力次第だけど
0686デフォルトの名無しさん
2006/03/30(木) 02:42:05低脳にもほどが（ｒｙ
0687デフォルトの名無しさん
2006/03/30(木) 12:34:19ありがとうございます。
>>686 すいません 「Matlab勉強しといて」って先生に言われただけで
さっぱりわかんなかったのです。低脳でした orz
0688デフォルトの名無しさん
2006/03/30(木) 12:51:57http://www.amazon.co.jp/exec/obidos/ASIN/4798008478/249-1329616-0032338
これが行列操作分かりよかったので最初の数章だけでも読むといい
0689デフォルトの名無しさん
2006/03/30(木) 13:55:45実際にやってみることだと思うがね。
とりあえず
x = 1
とコンソールに入力してみるところから始めればいいかと思うが。
0690デフォルトの名無しさん
2006/03/30(木) 14:35:20なぜまずそれをしない?
0691デフォルトの名無しさん
2006/04/03(月) 17:52:01ファイル名が２行目から
連なっているファイルfilename.txtがあります。
filename.txt
ファイルの内容はこんな感じ...
4
ichi.out
ninjin.txt
sandal.log
yatch.out
このファイルの2行目以降からの各ファイル名を
文字列として格納する文字列の配列Aを作りたいんですが、
Aはどのようにつくればよいでしょうか？
cell配列を作ればよいことは分かるのですが
Aの要素の書き換え方がわかりません。
どなたか助言願います。
0692デフォルトの名無しさん
2006/04/03(月) 18:16:26を使えばいいんですね。
事故解決しました。すみません...
0693デフォルトの名無しさん
2006/04/06(木) 13:52:300694デフォルトの名無しさん
2006/04/11(火) 15:58:10と思うんだけど、どうすれば白に固定とかできます？
0695デフォルトの名無しさん
2006/04/11(火) 21:02:310696デフォルトの名無しさん
2006/04/11(火) 23:02:19ﾊﾝﾄﾞﾙｸﾞﾗﾌｨｯｸｽ
0697デフォルトの名無しさん
2006/04/12(水) 00:04:39うん
0698デフォルトの名無しさん
2006/04/12(水) 06:05:50うーむ、
http://infoshako.sk.tsukuba.ac.jp/jdoc/MATLAB5/jhelp/techdoc/refbycat/group23.html
見当たらない。
0699デフォルトの名無しさん
2006/04/12(水) 16:52:570700デフォルトの名無しさん
2006/04/12(水) 23:20:24axisX=gca;set(axisX,'color',[1 1 1]);→バック白
0701デフォルトの名無しさん
2006/04/14(金) 00:18:43ていうかＭＡＴＬＡＢの演習問題集とかって書店で売ってる？
近場じゃ全然ＭＡＴＬＡＢ関連置いてﾅｻｽ・・・
0702デフォルトの名無しさん
2006/04/14(金) 01:57:49http://www.amazon.co.jp/exec/obidos/ASIN/478983090X/503-7139237-9691946
Amazon.co.jp：MATLABによる組み込みプログラミング入門―SimulinkとReal‐Time Workshopを使ったMeasurement & Control 計測・制御シリーズ: 本
http://www.amazon.co.jp/exec/obidos/ASIN/4789837173/qid=1144947422/sr=1-1/ref=sr_1_2_1/503-7139237-9691946
この２冊はいけっこうお勧めだと思う。
音響もしくは機械制御というそれぞれの分野にいる人間には
役立つと思う。
0703デフォルトの名無しさん
2006/04/14(金) 04:03:34一枚目は普通にインストールできてもまったく追加プロダクトが入れられない
「1. 新規および更新されたプロダクトを確認するためにThe MathWorks社のウェブサイトに接続できませんでした。」
と出てくるのだけどどうすればいいのだろう
助けて
0704デフォルトの名無しさん
2006/04/14(金) 14:55:160705デフォルトの名無しさん
2006/04/15(土) 11:33:46PLPは長ーーーくなっているはず。
0706703
2006/04/15(土) 14:52:29うちは
14-*****-*****-*****-*****
くらいの長さなんだけど、多分これでは短いのでしょう
（ウチの先生は「Toolbox全部使えるように注文したから頑張ってくれたまえHAHAHA…」といっていたのに…）
それなら追加プロダクトを使えるようにするには何の手続きをすればいいの？
0707デフォルトの名無しさん
2006/04/15(土) 16:10:59①既存ツ－ルの保守契約を結んでいないと追加のToolboxは販売してくれないから
もし、結んでないなら、まずこれを結ぶ必要あり。
②あとは追加したいToolboxの購入費用を支払うだけ。
>Toolbox全部使えるように注文したから...。
全部追加して＆保守契約を結ぶとｴﾗｲ費用だと思うよ。企業の場合はね。
ｼﾝﾎﾞﾙ計算はMathmaticaが優れているけど、数値解析、ｼｽﾃﾑ制御はMATLABの
1人勝ちだから強気の商売やってるよ。(>ｻｲﾊﾞﾈｯﾄさん)
0708デフォルトの名無しさん
2006/04/15(土) 17:31:37大変な金額が・・・。
だからOctaveに・・・。
でも目的によっては非力が玉に瑕。ごめんちゃい。
0709デフォルトの名無しさん
2006/04/16(日) 11:25:27>（ウチの先生は「Toolbox全部使えるように注文したから頑張ってくれたまえHAHAHA…」といっていたのに…）
今，Toolboxってすごくたくさんあるから，全部使えるように注文ってすごいぞ。
0710デフォルトの名無しさん
2006/04/16(日) 13:02:42ttps://web.kudpc.kyoto-u.ac.jp/lsc/modules/contents2/index.php?id=16
ttp://matlab.chase-dream.com/group/toolbox.html
0711デフォルトの名無しさん
2006/04/17(月) 09:14:550712703
2006/04/17(月) 09:51:04先生に聞いてみたところ、MATLABの注文は昔の上級生任せにしていたらしく
注文すればToolboxも自動的に全部つくと思っていたようです
道理で「安くついた」と言っていたわけだ…
現在はその上級生が変わりに注文していたMATLAB6.5の方を使っています（いくつかToolboxがついていたので）
そこで更に質問なのですが、6.5のToolboxを7.0のフォルダに移してサーチパスを通せば使えるでしょうか？
0713デフォルトの名無しさん
2006/04/17(月) 09:52:390714デフォルトの名無しさん
2006/04/17(月) 13:36:350715デフォルトの名無しさん
2006/04/17(月) 21:26:42パス通せば答えは自明でしょう～♪
0716デフォルトの名無しさん
2006/04/19(水) 20:17:39のですが、その方法が判りません。
もしできるのなら伝授下さい。お願い致します。
0717デフォルトの名無しさん
2006/04/20(木) 12:53:19% マシン:dual Op. 2.8GHz, メモリ10GBytes
% OS: WinXP X64 英語版
%
%行列の最大の大きさ
a=zeros(23170,23170);
% まで通る。これが限界。またonesはなぜか限界が低い。
% 変数名さえ変えれば、上の大きさの行列がいくつでも作れる。
% もちろんスワップ始めると、とても遅くなる。
%
%簡単ベンチマーク
clear all
a=ones(10000,10000);b=ones(10000,10000);c=zeros(10000,10000);
starttime=now*3600*24;
c=a*b;
now*3600*24-startime
%上記の実行時間が約480秒
% 結論 完成度高い。あとはindexの制限がなくなればいいのだが。
0718デフォルトの名無しさん
2006/04/20(木) 23:41:040719デフォルトの名無しさん
2006/04/21(金) 21:46:380720通りすがり１
2006/04/23(日) 10:52:27ActiveXで外部機器をmatlabで操作しています。
しかし時々（特にコンピューターに負荷がかかっているとき）
外部機器からのイベントがActiveXの３つのコールバックのうち一つ（一番左に書いてあるやつ）
をdropしてしまうようなのです。ActiveXControlはTimerObjectのようなので、
timer objectのbusymodeをqueueに設定する方法を探しています。
ちなみに外部機器の製造元はどうやって設定するのかわからないといっています。
どなたか助けてください。~
0721デフォルトの名無しさん
2006/04/23(日) 13:10:29ttp://www.cybernet.co.jp/matlab/support/quickstart/stepup_matlab/ml01/7.shtml
0722デフォルトの名無しさん
2006/04/23(日) 19:27:51u = [2+3*%i 1+4*%i];
max(u)
でuの要素の最大値を求めようと思うんですが
「maxi: Argument 1: wrong type argument expecting a real matrix」とエラーが出てしまいます
正しく最大値を求めるにはどうすればいいのでしょうか？
0723デフォルトの名無しさん
2006/04/23(日) 19:32:18なので，複素数の絶対値を取る必要があるのではないかと。
max(abs(u))
で試してみてはいかがでしょうか？
0724デフォルトの名無しさん
2006/04/23(日) 19:57:34scilabネタがあるなら使ってね
(´･ω･`)ｼｮﾎﾞｰﾝ
0725通りすがり１
2006/04/24(月) 08:12:33ここでset(timer,'BusyMode', 'queue')とっやた後、get(timer)してもBusyMode: dropがでできます。
どうすればいいのでしょうか？
0726通りすがり１
2006/04/24(月) 08:20:53変えたことにはなりませんよね？設定自体を変えるにはどうすればいいのでしょうか？
0727デフォルトの名無しさん
2006/04/24(月) 23:23:38ﾊﾝﾄﾞﾙ名の取得、そのﾌﾟﾛﾊﾟﾃｨの変更までいったじゃん。
自分で好きな設定をする.mﾌｧｲﾙ作ればいいんじゃないの？
ActiveXを使っているなら、GUIの .m or .figがあるんじゃあないの？
そこに追加してみてみて。
0728デフォルトの名無しさん
2006/04/27(木) 10:03:58WindowsかUNIXか知ることはできますか？
0729通りすがり１
2006/04/27(木) 10:31:40外部機器のハンドルをとってBusyModeを変更しようとしても、そもそもそんなプロパティ
は無いのでエラーがでます。.m fileとか.fig fileで具体的にどうすればいいのですか？
0730デフォルトの名無しさん
2006/04/27(木) 10:32:26ispcとかisunixとか、そんなんじゃなかったっけ？
0731デフォルトの名無しさん
2006/04/27(木) 11:52:26ありがとうございます。
0732デフォルトの名無しさん
2006/04/27(木) 13:17:53plot3(x,y,z)で出力した３次元のグラフィックのy,z平面にcontour(y,z)を載せる場合はどのように
したらいいですか？
hold on を使うとx,y平面に載ってしまいます。
0733デフォルトの名無しさん
2006/04/27(木) 14:21:190734デフォルトの名無しさん
2006/04/27(木) 15:59:26その方法が、わからないっす・・・。
0735デフォルトの名無しさん
2006/04/27(木) 23:24:51外部機器はUSB、LAN、GPIB、RS...etc 何のいんた－faceで制御してんの？
ToolBoXはInstrument Control ToolBox...etc 何使ってるのかな？
0736通りすがり１
2006/04/28(金) 06:47:130737デフォルトの名無しさん
2006/04/28(金) 14:40:51いっそのことyz平面がxy平面になるように座標変換すればいいじゃない
0738デフォルトの名無しさん
2006/04/28(金) 23:19:45おれっちが無知なのかもしれんけど
Instrumenｔ Control ToolBox なしでUSB機器制御できてのか？
0739デフォルトの名無しさん
2006/04/30(日) 19:28:26関数のfactorialってベクトルあるいは行列を引数にとることは
できないんですか？
factorial([1,2])=[1,2]
みたいに。というよりも、やってみたらできないので、
何か代替法があれば教えていただけると幸いです。
よろしくお願いします。
0740デフォルトの名無しさん
2006/04/30(日) 21:17:26scalar。
0741デフォルトの名無しさん
2006/04/30(日) 22:09:290742739
2006/04/30(日) 22:25:41レスありがとうございます。
しかし、scalarとは何でしょうか？
関数でもないみたいですし、ヘルプ見ても載ってないですし・・。
あるいは、ベクトルでなくスカラー上でしかfactorialは定義されない
という意味でしょうか？
0743デフォルトの名無しさん
2006/04/30(日) 22:26:350744デフォルトの名無しさん
2006/04/30(日) 22:36:20つttp://www.cybernet.co.jp/matlab/support/manual/r13/toolbox/matlab/ref/isreal.shtml
0746デフォルトの名無しさん
2006/05/01(月) 21:50:57for..end
0747デフォルトの名無しさん
2006/05/02(火) 23:57:08R14SP3 と R14SP1 はどちらがお勧めですか？
0748デフォルトの名無しさん
2006/05/03(水) 09:46:00R14SP1→R14SP3→Release 2006a(7.2)
でｱｯﾌﾟﾃﾞ-ﾄしてる。 最新はRelease 2006a(7.2)。
0749デフォルトの名無しさん
2006/05/04(木) 21:27:30それ以外は、そんなにかわらんような…
俺はR12.1を、メインで使っているような…
0750デフォルトの名無しさん
2006/05/09(火) 23:49:46scicosで作った.cosモデルをC言語もしくはC++のコードに変換して、
自分のmain関数から関数として呼び出したいんですが、
そういうことはできるんでしょうか？
やり方の書いてあるサイトや、そういうこと質問できそうなサイトがありましたら、
ぜひ教えてください。
0751デフォルトの名無しさん
2006/05/16(火) 09:17:460752デフォルトの名無しさん
2006/05/16(火) 23:41:420753デフォルトの名無しさん
2006/05/17(水) 06:04:040754デフォルトの名無しさん
2006/05/20(土) 21:52:39現状ではターミナルでしか動かないというのが現状なんですかねぇ・・・。
0755デフォルトの名無しさん
2006/05/21(日) 00:43:11有効に使うためにはどうしたらよいでしょか？
0756デフォルトの名無しさん
2006/05/21(日) 01:35:57自分レスです。
>現状ではターミナルでしか動かないというのが現状なんですかねぇ・・・。
なんか現状が2回でてきて、変な文になっています。
馬鹿な子ですいません。
0757デフォルトの名無しさん
2006/05/22(月) 21:24:42島GA
0758デフォルトの名無しさん
2006/05/22(月) 23:56:34「島GA」とは？
0759デフォルトの名無しさん
2006/05/23(火) 11:00:16/var/tmp/lm_TMW.ld: relocation error: /var/tmp/lm_TMW.ld: symbol errno,
version GLIBC_2.0 not defined in file libc.so.6 with link time reference
>>51
の方法を使ってもだめでした．解決方法をご存知のかたいらっしゃいません
でしょうか？
0760デフォルトの名無しさん
2006/05/31(水) 16:11:02if a; d = b; else; d = c; end とやるしかないんでしょうか。
0761デフォルトの名無しさん
2006/05/31(水) 22:48:310762デフォルトの名無しさん
2006/05/31(水) 23:09:55それだと
・b、cが有限でないと(NaN、±Infだと)異常な結果になる
・aがtrueのときもcが(falseのときもbが)実行される
という問題があります。
0763デフォルトの名無しさん
2006/06/05(月) 21:55:30を、どなたか教えていただけませんか？
0764デフォルトの名無しさん
2006/06/05(月) 23:29:360766デフォルトの名無しさん
2006/06/06(火) 23:57:22のように定義して、配列位置と最小値を表示させるにはどのように書いたらよいのでしょうか？
特に配列位置を表示させるのに手間取っています。
どなたか教えてくださいませんか？
0767デフォルトの名無しさん
2006/06/07(水) 00:31:000768デフォルトの名無しさん
2006/06/07(水) 00:39:46最小値のインデックスってこと？
%最小値
min(U)
%インデックス
find(U == min(U))
でどう？
0769デフォルトの名無しさん
2006/06/07(水) 00:49:29なるほど！！
findをそういう使い方をするんですね。
わかりました。ありがとうございます！！
0770デフォルトの名無しさん
2006/06/07(水) 00:50:29やってみましたが、納得いく答えが得られました。
ありがとうございます。
0771デフォルトの名無しさん
2006/06/07(水) 02:10:48maxやminの戻り値は二つあって、
[value, index] = min(U);
のように最小値の場所を得られるよ。
0772デフォルトの名無しさん
2006/06/07(水) 02:36:56わざわざありがとうございます。
やっぱり、関数を把握できるようにならないときついかな・・・。
0773デフォルトの名無しさん
2006/06/08(木) 20:17:351. 最小値の要素が複数あるとき
2. 全てNaNのとき
0774デフォルトの名無しさん
2006/06/12(月) 18:17:50できますか？
0775デフォルトの名無しさん
2006/06/12(月) 19:32:020776774
2006/06/12(月) 19:49:26あ、すいません。MATLAB使ってです。
一応、ImageAcquisitionToolBox、とImageProcessingToolboxあるんですけど
使えますかね？
0777ゆーき
2006/06/15(木) 14:40:15matlab初心者なのですが。。。
csapeで曲線を描きたくて調べたのですが、どうしてもcsape(X,Y,CONDS,VALCONDS)
のVALCONDSに何を入れたらいいのか分かりません。。。
特異値とか特異条件とかが関係してくるっぽいんですが。
どなたかご教授願います。
0778デフォルトの名無しさん
2006/06/15(木) 15:55:382本の線をプロットする際、線同士が重なってしまうところは
後からプロットしたほうが表示されますよね。
これを防ぐ方法はないのでしょうか？たとえば少しずらして表示されるとか
どなたかご教授お願いします。
0779デフォルトの名無しさん
2006/06/15(木) 20:34:04後ろの線を太くする
0780デフォルトの名無しさん
2006/06/16(金) 15:31:25scilabではdeff('')でよかったのですが同じようなことをmatlabでするにはどのようにすればよいのでしょうか？
0781778
2006/06/16(金) 22:18:18ご返信ありがとうございます！
そう言えばそれでいけますね．ありがとうございました．
0782デフォルトの名無しさん
2006/06/17(土) 02:25:00プログラムを走らせているときにでる警告文を
表示しないようにするにはどうしたらよいのでしょうか？
うっとおしくて困っています。
よろしくお願いします。
0783デフォルトの名無しさん
2006/06/17(土) 10:15:45warning off
特定の警告を抑制するには
warinig off MATLAB:***
MATLAB:***の部分は、警告が出た後lastwarnで出る
0784782
2006/06/17(土) 16:20:59ありがとうございます！
かなり快適になりました！
0785デフォルトの名無しさん
2006/06/20(火) 15:14:27matlabのプログラムで、特定の行に移動するように指定するには
どうすればよいのでしょうか
よろしくお願いいたします
0786デフォルトの名無しさん
2006/06/20(火) 17:16:09ある数値がNaNであることを判定するにはどうすればよいでしょうか？
a=NaN
とし、
a==NaN
とすると0が返ってきて、判定できません。
どなたかご教示お願いします。
0787デフォルトの名無しさん
2006/06/20(火) 18:27:36ってかマニュアル読めば良いんじゃね？
0788デフォルトの名無しさん
2006/06/20(火) 20:42:530789786
2006/06/21(水) 11:15:34早速にありがとうございました。
>787
そうですね…。安易に書き込んでしまいました。申し訳ありません。
非常に助かりました。
0790ゆーき
2006/06/21(水) 14:31:15>777
のこと分かりませんか。。。?orz
0791デフォルトの名無しさん
2006/06/21(水) 19:07:260792ゆーき
2006/06/21(水) 21:49:03感じで。。。どなたか。。。orz
0793デフォルトの名無しさん
2006/06/21(水) 23:22:52E-mail：
内容：
どなたか教えてください。
Matlabでファイルを開くからfigファイルを開こうとすると
エラー：==>open
Invalid field name'OpenFig_gifdeltinu～'
と出力され、開くことができません。
さらに、その後に普通にplotするときもなぜかfigファイルが開けず、
matlabを再起動するまで一切のfigファイルを開くことができません。
上のInvalid～のところで、ファイル名が逆（アベコベ）に表示されて
いるのも謎です。 よろしくお願いします。
0794793
2006/06/21(水) 23:24:460795デフォルトの名無しさん
2006/06/22(木) 16:53:12振幅をデータ数の半分で割っても、うまくいきません。
どうしたらいいでしょうか。
0796デフォルトの名無しさん
2006/06/22(木) 20:01:08まず目的を書こう
振幅をデータ数の半分で割るという操作にどんな意味があるかは解らないけど
もしもよくある周波数-スペクトルのグラフ等を描くためのデータを作りたいのなら
fft関数の出力をabs関数で絶対値をとるといいよ
fft関数の出力は複素数が含んでいるからね
0797デフォルトの名無しさん
2006/06/22(木) 23:07:15周波数ースペクトルのグラフを書こうとしてます。
絶対値とってるんですが、うまくいきません。
参考書などを見て、データインデックスから周波数へのの換算と、ｆｆｔデータから振幅の換算
を行っても期待する値が振幅では得られません。
0798デフォルトの名無しさん
2006/06/22(木) 23:17:04意味あるよ
もちっと勉強しようね
>>795
>796の言うとおり質問自体が？
0799デフォルトの名無しさん
2006/06/22(木) 23:58:18うまく伝えられなくてすいません。
具体的な例を出すと、sinc(t)をｆｆｔすると幅１[Hｚ]、振幅１の矩形スペクトル
になると思うんですが、振幅はそうなりません。
形は矩形になるんですけど・・・。
0800デフォルトの名無しさん
2006/06/23(金) 00:36:28そのとおり の結果になるはずだし
現実にそのとおりの 結果だよ
今晩は寝た方が...。
0801デフォルトの名無しさん
2006/06/23(金) 00:53:03そのとおりの結果でますか？
よければプログラムを・・・
0802ゆーき
2006/06/24(土) 20:17:42度々。。。
subplotで出力したfigureを保存したいのですが、print -depsやgset term postscript
enhanced を使ってもできません(subplotでなければ無問題)。
ど、どなたかぁぁぁぁぁぁ。。。。。ご教授ください。
0803デフォルトの名無しさん
2006/06/25(日) 20:22:06それはMATLABじゃなくてoctaveだよね。
Octaveスレで聞いてみたら？
0804ゆーき
2006/06/25(日) 23:16:520805デフォルトの名無しさん
2006/06/25(日) 23:27:48ないでしょ？
ここでいいんじゃね？
0806デフォルトの名無しさん
2006/06/25(日) 23:30:43総合スレいったほうがいい。
0807デフォルトの名無しさん
2006/06/25(日) 23:42:330808デフォルトの名無しさん
2006/06/25(日) 23:56:330809ゆーき
2006/06/26(月) 00:02:49ちなみにmatlabでsubplotを保存するにはどうするのですか?
あと、できれば
>>777 >>792
の質問も見て欲しいです。
お願いします。
0810デフォルトの名無しさん
2006/06/26(月) 00:03:56ﾉｼ
>>802
gsetでterminalとoutputを予め設定してからplotしてみたら。
0811デフォルトの名無しさん
2006/06/26(月) 02:36:55Universal 対応してないから使ってない
0812ゆーき
2006/06/26(月) 18:38:29gset term postscript enhanced;
gset output "a.eps";
replot;
としているのですがterminalの設定方法がわかりません。教えて下さい。
0813デフォルトの名無しさん
2006/06/26(月) 19:42:31subplotはgnuplotのmultiplotモードを使う。
gnuplotはmultiplotモードではterminalを変更できない。
従ってsubplotコマンドを使用する前に
terminalとoutputを変更すればOKのはず。
これが希望する解決策でないことはわかるけど
エレガントな解決方法はわからない。
0814ゆーき
2006/06/26(月) 20:00:48とてもありがとうございます。自分はエレガントなど求めておりませんので保存で
きればけっこうです。
ただ１つ問題が。。。
terminalとoutputの変更の仕方が分からない。。。(ハイパー素人のため)
subplotとしてではなくてグラフ１つずつ保存すると言うことでしょうか??
お願いします。
0815デフォルトの名無しさん
2006/06/26(月) 20:16:28gset terminal postscript eps
gset output "test.eps"
subplotとplot色々
gset nomultiplot
gset terminal x11
これでtest.epsに保存できる。
0816813=815
2006/06/26(月) 20:39:170817ゆーき
2006/06/26(月) 21:20:06ほんとにありがとうございました!!!!
0818デフォルトの名無しさん
2006/06/27(火) 09:18:280819ゆーき
2006/06/27(火) 13:36:05保存できました!!!!ありがとうございました!!
。。。が!!!事件です!!!!!グラフがほぼ壊滅状態に陥るほど劣化してしまい使えません。orz
やっとできたと思ったのに。。。
誰か救いの手をさしのべて下さい。
ちなみに
>>777 >>792
も見てもらえると有り難いです。
0820813=815
2006/06/27(火) 18:44:33私はエスパーではないのでry
octaveについてはUNIX板か数学板の
octaveスレに
行けば詳しい人もいるかも。
0821ゆーき
2006/06/27(火) 20:06:04ご親切ありがとうございます!!!
頑張ってみます!!
0822デフォルトの名無しさん
2006/06/29(木) 19:27:48／|
＿＿＿／ ｜ ＿＿＿＿ 0
｜ ／1
｜／
0 -1
みたいなグラフを描く場合、どういう風に記述すればいいんですか？
途中で切り替える場合が出来ません。
ifを使って、
if -1<= t < 0
a=1;
b=1;
elseif 0 <= t < 1
a=1
b=-2
end
x = a*t + b
みたいな感じで記述したらどうかなーと思ってたんですが
値が後者に固定されてしまい全く出来ません
助けてください＞＜；
0823デフォルトの名無しさん
2006/06/29(木) 20:29:45「0<=t」がまず評価され0か1に置き換えられ、
「0<1」または「1<1」が評価される。
お前の意図どおりの結果を得る式は「0<=t && t<1」。
しかしグラフを描くには、t=-2:1/256:2みたいな入力に
対応できないと困らないか？
そのためには条件分岐を使ってはいけない。
isbitween01 = 0<=t & t<1;
x(isbitween01) = t(isbitween01) + 1;
みたいに書くべし。
0824デフォルトの名無しさん
2006/06/29(木) 22:34:00hamming(sample_point)
は使えるんですけど
hanning(sample_point)
だとscript hanningを関数として使ってください
みたいなことがエラーで出ます。
初心者ですいません。
教えてください。
0825デフォルトの名無しさん
2006/07/01(土) 00:01:260826デフォルトの名無しさん
2006/07/04(火) 23:30:18でも、窓関数関数自体のスペクトルがうまく出力できないー。
0827matlabとＣ言語について
2006/07/06(木) 14:26:08フリーのscilabはlink可能だったのですが研究でmatlabを用いて数値計算をしたいと考えています．
よろしくお願いします．
0828デフォルトの名無しさん
2006/07/06(木) 15:09:450829ゆーき
2006/07/07(金) 11:43:23valcondsにいれるべき値とはどういったものかお教えください。お願いします。
0830デフォルトの名無しさん
2006/07/07(金) 12:41:50OBB(Oriented Bounding Box)で、物体頂点群から、
主成分分析をもちい(共分散行列つくって固有ベクトル)、
3つの軸を見つけることができたのですが、
直方体の物体に対する第一成分の軸が、最長の対角線上
になってしまいます。こうなるとAABBより悪い。。。
あとで考えてみると、確かに普通に適用するとそうなるのは
わかったのですが、外国の論文みてると、どうやらこれを回避
する方法があるらしいのです。
ただ、英語がよくわからないので、もう一歩が読みきれません。
あとなにかのアイデアひとつと思うですが、
ポイントを教えてください。
0831デフォルトの名無しさん
2006/07/07(金) 15:43:26スレ違い。あとマルチいくない
0832デフォルトの名無しさん
2006/07/07(金) 16:08:27すみません。撤回します。しかし、よくマルチって
わかるなぁ。全部マイナーなスレにのっけたのに。
0833デフォルトの名無しさん
2006/07/07(金) 16:12:17検出するソフトがある
0834デフォルトの名無しさん
2006/07/08(土) 13:33:34プログレスバーを移動したり、再生部分の時刻を表示する
ことは可能ですか？
0835デフォルトの名無しさん
2006/07/08(土) 13:50:24そのまでやるのか・・・キモ
0836デフォルトの名無しさん
2006/07/08(土) 14:18:46. ＼ | ノ ヽ
＼ / ● ● |
＼| ( _●_) ミ
彡､ |∪| ,/..
ヽ ヽ丿/ /⌒|
/ ＼＿＿ノﾞ＿/ / =====
〈 _ノ ====
＼ ＼＿ ＼
＼＿__） ＼ ====== (´⌒
＼ ＿＿_ ＼＿＿ (´⌒;;(´⌒;;
＼＿＿＿）＿＿＿）(´;;⌒ (´⌒;; ズザザザ
0837デフォルトの名無しさん
2006/07/08(土) 16:49:020838デフォルトの名無しさん
2006/07/10(月) 08:56:090839デフォルトの名無しさん
2006/07/14(金) 23:41:480840デフォルトの名無しさん
2006/07/15(土) 09:32:370841デフォルトの名無しさん
2006/07/15(土) 15:05:240842デフォルトの名無しさん
2006/07/16(日) 14:09:33実際のパワースペクトル密度を得たいときはどうすればいいんですか？
勉強不足ですいません。よろしくおねがいします。
0843デフォルトの名無しさん
2006/07/16(日) 23:48:210844デフォルトの名無しさん
2006/07/17(月) 20:56:13やってみます。
0845デフォルトの名無しさん
2006/07/19(水) 02:29:48連続時間フーリエと同じ結果になりました。
0846デフォルトの名無しさん
2006/07/21(金) 06:14:31●モデル：
Test.mdl
Test/Subsystem
Test/Subsystem/ConstantA
Test/Subsystem/BlockB
構成は、Subsystemのブロックパラメータで変数Aに値1を設定し、
ConstantAに変数Aを設定します。
ConstantAはBlockBにつなげてBlockBの入力にします。
コンフィグレーションはタイプ：固定ステップ、ソルバ：離散、固定ステップ：0.01
としてます。
●質問：
ここでシミュレーションを実行すると、時間ゼロではConstantAの変数A=0(NaN?)
がBlockBに入力され、次のステップでA=1がBlockBに入力されているようなので
すが、これを時間ゼロからBlockBにA=1が入力されるようにすることはできるの
でしょうか。
BlockBはS-Function(Level2 Mファイル)で、変数Aで割り算しているところがあ
るのですが、当初動作確認のためBlockBのブロックパラメータで変数Aを設定
して問題なく動作していたため、上記のように時間ゼロ時に変数Aで割るとゼロ割
になってしまいます。BlockBで工夫すればゼロ割は避けられますが、プログラ
ムの都合上時間ゼロでConstantAの変数Aが確定している方が良いです。
なお、実際作成しているモデルはブロックパラメータが複数有り、複数のConstant
出力をMuxしてBlockBに入力しています。
0847デフォルトの名無しさん
2006/07/21(金) 15:32:08今回、エクセルのデータ（ｘとｙの値が５０個ずつくらい）をMATLABでプロットしたいんですが、
データの読み込み方もそれをプロットするやり方もわかりません＞＜
webで調べても難しいことばっか言っててわかんなかったです…
こんな私にどなたかご教授ください～
0848デフォルトの名無しさん
2006/07/21(金) 18:23:320849デフォルトの名無しさん
2006/07/21(金) 21:37:09氏ね
0850デフォルトの名無しさん
2006/07/21(金) 22:52:46うまくいくかわかりませんが，IC（初期値を設定するブロック）を
ConstantとS-funの間に挿入してはどうでしょうか。
http://www.cybernet.co.jp/matlab/support/manual/r14/toolbox/simulink/slref/ic.shtml
>>847
つ ささやかなヒント
http://www.cybernet.co.jp/matlab/support/manual/r14/toolbox/matlab/matlab_prog/ch_impo6.shtml#173379
0851デフォルトの名無しさん
2006/07/21(金) 23:44:32皆に酷評されてるけど それくらい自分で努力せにゃあかんよ
最低だぞー
0852847
2006/07/24(月) 11:46:38さっぱり何を言っているのか泣
今回はしかたなく、エクセルのデータをテキストデータにして、
a=dlmread('データ')
x=a(:,1)
y=a(:,2)
plot(x,y)
ってやりました。
問題は、
x1 y1 y2 … yn
x2 y1 y2 … yn
x3 y1 y2 … yn
：
：
っていうデータで、x1とx2の間のxは等間隔でy1～ynと対応しています。
これをきちんと、xとyのグラフにプロットするにはどういうふうにしたらいいんでしょうか。
バカかと思われてもしかたありませんが解決策を教えてください…
0853デフォルトの名無しさん
2006/07/24(月) 12:06:02x1 y1
x2 y2
x3 y3
.
.
xn yn
じゃねーの？
0854デフォルトの名無しさん
2006/07/24(月) 12:13:58x1とx2の間が等間隔なんですが値がなく、それとy1～ynが対応してるんです…
XYDATA=(X++(Y..Y))
っていう…。それをきちんとxとyのグラフにプロットがしたいんです
0855デフォルトの名無しさん
2006/07/25(火) 22:59:11回答ありがとうございます。私が作成したMファイルに問題があったようで、
直接フィードスルーがfalseになったままだったのが原因でした。
trueにしたところ無事t=0で入力端子からの値が確定しています。
ただ、t=0の一回のみ実行のファンクションInitialConditionsやStartでは
入力端子からの読み込みはできないようで、結局Outputsファンクションで
block.CurrentTime==0のときだけ入力値から初期化するようにせざるを得
ないようです。Startで入力値で初期化できたらスマートだったのですが。
なお、ICはt=0のときのみ値を変えたい場合に利用できるようで、ひとつ勉
強になりました。
0856デフォルトの名無しさん
2006/07/26(水) 23:30:18xnにones()を掛ければ、ynと同じサイズ になるよｗ
あとは行ごとに列ベクトルに変換して、 cat()で Endかな？
0857デフォルトの名無しさん
2006/07/27(木) 16:36:23プログラムどうやって作るの？
0858デフォルトの名無しさん
2006/07/27(木) 20:24:42^
0859デフォルトの名無しさん
2006/07/28(金) 01:49:52subplotで別のウィンドウに表示するには
どうすればいいか教えてください…
0860ど素人
2006/07/28(金) 01:59:56ぜひお願いします。
A=new Array(3, 8, 19, 20, 22, 24, 45);
i=0;
j=A.length-1;
a=19;
k=Math.floor((i+j)/2);
while(i<=j){
if(A[k] == a){
break;
}
else if(A[k] > a){
j = k - 1;
}
else{
i = k + 1;
}
k = Math.floor((i+j)/2);
}
if(i>j) k=-1;
document.write(a," ",k);
0861デフォルトの名無しさん
2006/07/28(金) 02:42:56なぜ matlab スレ？
0862デフォルトの名無しさん
2006/07/28(金) 02:44:55マルチ＆板違い
死ね
0863デフォルトの名無しさん
2006/07/28(金) 08:45:141000000 * 1.05^10
0864デフォルトの名無しさん
2006/07/28(金) 10:31:200865デフォルトの名無しさん
2006/07/28(金) 15:31:480866デフォルトの名無しさん
2006/07/28(金) 23:12:24悩み意味不明。
0867デフォルトの名無しさん
2006/07/29(土) 15:21:30止まるというプログラムを書きたいのですが、
その場合、どの部分でも止めたらいいわけでなく、排他制御というか
ある箇所までは絶対に処理を行って欲しいのですが、
このようなプログラムはMATLABで可能ですか？
0868デフォルトの名無しさん
2006/07/30(日) 05:06:17ボタン押したときにフラグ立てて
終了したいところで判定すればいいんじゃない？
0869デフォルトの名無しさん
2006/07/30(日) 05:09:170870859
2006/07/30(日) 13:20:36すみません、元の信号と、フィルタ処理をした後の2種類の信号があって、
周波数を横軸にしたmagnitudeのグラフを2種類並べたfigure1と、phaseのグラフを2種類並べたfigure2を別々に表示させたいんです。
1つの信号の周波数特性自体はfreq関数で表示できるのですが、上記のようにして利得と位相どうしを比較するfigureを作りたいです。
0871デフォルトの名無しさん
2006/07/30(日) 14:38:02hold と subplot の組合せで
0872859
2006/07/31(月) 19:05:09holdだと重ねて表示されてしまいますね、縦に並べて表示したかったんですが…。
でもありがとうございました！
0873デフォルトの名無しさん
2006/08/01(火) 00:29:360874デフォルトの名無しさん
2006/08/01(火) 14:53:42MATLAB R-linkでRの関数は使えるようになったんですが。
0875デフォルトの名無しさん
2006/08/01(火) 15:07:580876デフォルトの名無しさん
2006/08/01(火) 16:56:55仕様です。
画面表示と印刷（含EPS）では違う描画ルーチンを使っているようで、
グラフのグリッド幅なんかも違ってきますよね。
0877デフォルトの名無しさん
2006/08/02(水) 20:49:02uをsymbolとして、
if u >= 1, a=1
else a=0
end
みたいなことをやりたいんですが。。。
こんな感じでやると「'>='は'sym'classには定義されてない」って怒られます。
0878デフォルトの名無しさん
2006/08/05(土) 21:42:28forル－プ＋textreadで、20Mのテキストデ－タを読ませたら
完了するまで約5Hかかった。(tic toc 計測)
即刻、forル－プはずして、textreadのオプション使用したら
約2minで完了。
MATLABのforル－プは怖っ。
0879デフォルトの名無しさん
2006/08/08(火) 12:41:370880デフォルトの名無しさん
2006/08/08(火) 20:22:38正規のライセンスでもそうな(ry
0881デフォルトの名無しさん
2006/08/08(火) 21:07:05for...endはBigSizeのデ－タ処理には向かない。
PCのパフォ－マンスで改善はされるけど。
使わないのが賢明。
工夫。工夫。 さ。
0882デフォルトの名無しさん
2006/08/11(金) 23:38:160883デフォルトの名無しさん
2006/08/15(火) 23:29:00ttp://wwwcs23.cs.kobe-u.ac.jp/~mori/Matlab/mexmemo.html
0884デフォルトの名無しさん
2006/08/17(木) 11:42:19mcc -m xxxx.m
とコンパイルしたファイル群をMATLABのないPCで実行しようとすると
「MCLMCRRT74.DLLが見つかりません」
と言われるんですが、どうしたらいいんでしょう。
問題のDLLは見つかったんですが、同じディレクトリに入れて実行すると
エラーメッセージさえ出さずに異常終了します。
MNTALABは7.2.0.232 (R2006a)、MCCはV4.4 (R2006a)です。
0885デフォルトの名無しさん
2006/08/17(木) 11:45:230886デフォルトの名無しさん
2006/08/17(木) 21:52:46ttp://www.cybernet.co.jp/matlab/support/techkwdb/search_view.php?id=1399
0887デフォルトの名無しさん
2006/08/29(火) 15:04:23Mファイルを人からもらってそのソースコードを見たいのですが、Editorで開いても文字化けして見ることができません。誰か教えて下さい。よろしくお願いします。
0888デフォルトの名無しさん
2006/08/29(火) 19:51:520889デフォルトの名無しさん
2006/08/29(火) 21:02:58*.mファイルは、暗号化されている。
0890デフォルトの名無しさん
2006/09/02(土) 16:49:54(s^2 + Xs + 1)/(s^2 + Ys + 1)みたいな伝達関数から
B=[なんたらかんたら]
A=[うんたらちんたらら]
みたいに一発でもとめたいんですけど。
0891デフォルトの名無しさん
2006/09/02(土) 22:45:32マルチポスト乙
0892デフォルトの名無しさん
2006/09/03(日) 01:01:00例えばインクリメントの++をMATLABでも使えるようにするとか…
すでに存在するの2項演算子の「+」ならplus関数をオーバーライドすればいいみたいですが、
存在しない演算子を新たに定義することは可能なんでしょうか。
0893デフォルトの名無しさん
2006/09/03(日) 18:56:03答えられない香具師は書かなくていいです。わかりやすく書ける人だけ書いてください。
0894デフォルトの名無しさん
2006/09/04(月) 01:13:17何かいい方法はないでしょうか？
よろしくお願いします。
0896デフォルトの名無しさん
2006/09/04(月) 01:38:23お前自身も「答えられない香具師」じゃないか。ここに書かなくていいです。
0897デフォルトの名無しさん
2006/09/04(月) 01:48:30ヒント：
ttp://dl.cybernet.co.jp/matlab/support/manual/r14/toolbox/matlab/math/funfun14.shtml#381254
2重積分のところも読むといいかも。
0898デフォルトの名無しさん
2006/09/04(月) 04:50:07ありがとうございます。
でも、なんで２重なんですか？？
0899デフォルトの名無しさん
2006/09/04(月) 14:26:290900デフォルトの名無しさん
2006/09/04(月) 22:02:32関数の作り方が書いてある
0901デフォルトの名無しさん
2006/09/05(火) 00:34:29なんとかできました。
あと、数式じゃなくてfftして得られたグラフやデータを積分することもできるんでしょうか？
0902デフォルトの名無しさん
2006/09/05(火) 01:32:48台形公式、ガウス積分，シンプソン法なんかでググれば出てくると思います。
0903デフォルトの名無しさん
2006/09/05(火) 02:33:09パラメータの刻み幅はどこで決まるんですか？
0904デフォルトの名無しさん
2006/09/05(火) 02:34:47パラメータの刻み幅はどこで決まるんですか？
0905デフォルトの名無しさん
2006/09/05(火) 03:47:37刻み幅を指定して積分を行いたいなら，自分で積分のルーチンを作るしかないと思います。
0906デフォルトの名無しさん
2006/09/05(火) 03:58:130907デフォルトの名無しさん
2006/09/05(火) 04:01:060908デフォルトの名無しさん
2006/09/05(火) 22:57:020909デフォルトの名無しさん
2006/09/07(木) 12:30:37>>1
0910デフォルトの名無しさん
2006/09/09(土) 03:07:250911デフォルトの名無しさん
2006/09/09(土) 09:31:49@bartlett - Bartlett window.
@barthannwin - Modified Bartlett-Hanning window.
@blackman - Blackman window.
@blackmanharris - Minimum 4-term Blackman-Harris window.
@bohmanwin - Bohman window.
@chebwin - Chebyshev window.
@gausswin - Gaussian window.
@hamming - Hamming window.
@hann - Hann window.
@kaiser - Kaiser window.
@nuttallwin - Nuttall defined minimum 4-term Blackman-Harris window.
@rectwin - Rectangular window.
@tukeywin - Tukey window.
@triang - Triangular window.
0912デフォルトの名無しさん
2006/09/11(月) 22:58:59orz...。
0913デフォルトの名無しさん
2006/09/12(火) 00:44:51@hanning - Hanning window
も仲間にいれてあげて。
@gausswin - Gaussian window.
@blackmanharris - Minimum 4-term Blackman-Harris window
これあったんだな。自分で実装して使ってた。
ところで、上のほうにある
startup.mでデフォルトの線の太さを指定って需要ある？
必要なら今度学校から書き込むが。
0914デフォルトの名無しさん
2006/09/12(火) 14:16:56ライセンスマネージャーを起動しようとしても
---------------------------------------------------------------------
Checking license file for local hostname and local hostid . . .
Taking down any existing license manager daemons . . .
sort: オプションが違います -- 1
詳しくは `sort --help' を実行して下さい.
/usr/local/matlab/etc/lmdown: line 160: echo: write error: パイプが切断されまし た
Killing license manager daemon processes . . . (localhost.localdomain)
process id = 3898 killed . . . (/bin/sh)
Starting license manager . . .
Debug logfile = /var/tmp/lm_TMW.log
Waiting 300 secs for MATLAB vendor daemon to come up . . .
Type your interrupt character (usually CTRL-C) to quit.
Time = 3 secs : still waiting . . .
Time = 30 secs : still waiting . . .
--------------------------------------------------------------------
のようなメッセージが出てきて起動できません。
何が原因なのでしょうか？
よろしくお願いします。
0915デフォルトの名無しさん
2006/09/12(火) 22:54:32ライセンス形態とかの情報は？
0916デフォルトの名無しさん
2006/09/13(水) 14:24:39何とか解決できました。
すみませんでした。
0917デフォルトの名無しさん
2006/09/13(水) 14:55:58つまりたとえば任意のTrueTypeフォントで「Ａ」という文字を300x300ドットの
ビットマップ画像（行列）にするとかいうことをしたいのですが・・・
0918デフォルトの名無しさん
2006/09/26(火) 11:55:23A=['a' 'b' 'c']
B=[1 2 3]
という二つの変数があった時、
Bの名前を、Aの要素に変えることはできますか?
似ているんだけど、
または、forループを使って、
変数Aの要素が名前の新変数、a,b,cを3つ作ることはできますか?
0919デフォルトの名無しさん
2006/09/27(水) 11:05:40例えば、
y=[-2 -4 -6 -4];
bar(y);
としたとき、
ゼロからのプロット（軸の上から下に向かって表示）ではなく、
負の無限大？から（軸の下側から上に向かって表示）プロットしたいとき
どのようにすればよいか
分かる方いましたら教えてください。
0920デフォルトの名無しさん
2006/09/27(水) 18:06:27○とか×でデータ点をプロットするときに，
○，×の大きさを変更することはできるのでしょうか？
'FontSize'だとエラーではじかれるし，
'LineWidth'だと輪郭が太くなるだけなので困っています．
フィギュアプロパティでフォントサイズを変更すれば，できることは
できるのですが，プログラムで書きたいのです．
どなたかご教示をお願いいたします
0921デフォルトの名無しさん
2006/09/28(木) 00:21:50'MarkseSize'
0922デフォルトの名無しさん
2006/09/28(木) 00:22:43'MarkerSize'
0923920
2006/09/28(木) 16:58:35ありがとうございました！助かりました。
0924デフォルトの名無しさん
2006/10/04(水) 07:22:430925デフォルトの名無しさん
2006/10/04(水) 12:27:320926デフォルトの名無しさん
2006/10/04(水) 22:53:06matlabでポリガンマ関数って扱えますか？
よろしくお願いします
0927926
2006/10/04(水) 23:15:45psiですね
0928デフォルトの名無しさん
2006/10/05(木) 05:16:09specgramの結果にスプラインをかけて、分解能ブロックごとに色づけではなく
滑らかに表示することはできないのでしょうか？ どなたか行っている方が
いらっしゃいましたら、ご教授お願いいたします。
0929デフォルトの名無しさん
2006/10/05(木) 15:22:03matlabでfunctionのなかにfunctionを定義することってできますか？
0930デフォルトの名無しさん
2006/10/05(木) 15:34:00できますよ。
一行目とそれより下に二つ以上書いて、試してみて。
Linuxで学校のmatlabを使ってるんですけど、
Signal Processing Toolboxのライセンスがあるとか、って
どうやって確かめられるんでしょうか。
学校に聞くのはなしで。
0931デフォルトの名無しさん
2006/10/05(木) 15:34:33一行目とそれより下、合わせて2つ以上です。ごめん
0932デフォルトの名無しさん
2006/10/05(木) 16:38:53普通にsmoothingフィルタかけてみたら？
0933デフォルトの名無しさん
2006/10/05(木) 22:43:23HELPでSignal Processing ToolboxのHELPが出れば、めでたしめでたしだー。
0935デフォルトの名無しさん
2006/10/05(木) 23:49:01そうなのか？
忘れたが、verだったかversionだったかで出るのではないか？
そのものズバリのlicenseなんかもあったような…
0936デフォルトの名無しさん
2006/10/06(金) 12:04:21license.datの中身にどのツールボックスが入ってたか書いてる気がする。
INCREMENT Signal_Toolbox
とかだったかな？記憶がはっきりしない。
0937困ってます
2006/10/07(土) 08:32:53自分のＭＡＴＬＡＢでfevalを使用すると
>feval(f,2)
??? 入力引数 'x' が定義されていません.
と返ってきます。ｘ＝２じゃないのですか？
それでいつもfevalを回避してプログラム作る羽目になっているのですが
これは一体何が原因かわかりません。どなたかアドバイスください。
0938デフォルトの名無しさん
2006/10/07(土) 14:55:48でうごいたけど？
0939困ってます
2006/10/07(土) 23:34:10レスありがとうございます。
>feval(f,2) ではなく
>feval('f',2) としなきゃいけなかったんですね。
馬鹿らしいけどマニュアルにはfeval(functionname, x)と書いているのが
多くて・・・。多謝です！
0940デフォルトの名無しさん
2006/10/08(日) 01:36:57多くの例はたとえば
f = 'abs';
feval(f, 2);
という状況を意味しているのでしょう
0941困ってます
2006/10/08(日) 02:13:21難しいプログラムも結構くだらないところでつまずいているものですね。
というか自分がどじなだけだろうけど。
ご示唆ありがとうございました。
0942デフォルトの名無しさん
2006/10/09(月) 23:27:34そのtoolboxは使えるようにならない？
0943デフォルトの名無しさん
2006/10/10(火) 08:46:51なるわけない。
ライセンス(License.dat)でインストールしているtoolboxは管理されてる。
0944デフォルトの名無しさん
2006/10/11(水) 18:14:230945デフォルトの名無しさん
2006/10/12(木) 01:05:400946デフォルトの名無しさん
2006/10/12(木) 16:32:220947デフォルトの名無しさん
2006/10/15(日) 16:20:310948デフォルトの名無しさん
2006/10/16(月) 13:14:580949デフォルトの名無しさん
2006/10/16(月) 20:21:290950デフォルトの名無しさん
2006/10/19(木) 01:14:28で描いたx1軸とx2軸の
両方を日付目盛りにする
方法を教えて下さい
お願いします
0951デフォルトの名無しさん
2006/10/21(土) 01:18:52引数のサイズがでかいとかいわれるんだけど…
なんじゃそりゃー！
0952デフォルトの名無しさん
2006/10/21(土) 01:57:14日記ならご自分ブログでどうぞ
0953デフォルトの名無しさん
2006/10/21(土) 01:59:10fminsearchの初期値つまり求めたい変数に数百×数百（千だったかも）の行列を突っ込んだ．
0954デフォルトの名無しさん
2006/10/21(土) 12:25:51試してないけど fminsearch なら確かに妥当っぽい動作な気もする。
初期値からちょっと変えて計算、ちょっと変えて計算ってやるんだろうし。
0955デフォルトの名無しさん
2006/10/21(土) 15:42:27optimization toolbox持ってないとダルイのね・・・
0956デフォルトの名無しさん
2006/10/21(土) 15:59:150957デフォルトの名無しさん
2006/10/22(日) 20:16:58octaveだとあるのだが。
UNIX版matlabだと system('sleep 10'); なんてやってますが、
windows版でできる方法はないでしょうか。
0958デフォルトの名無しさん
2006/10/23(月) 02:58:580959デフォルトの名無しさん
2006/10/23(月) 05:04:13次のように出てしまいます。どこがおかしいのかわかりません。
どなたかアドバイスいただけませんか。
>>f=@(x,y) x*y
??? f=@(x,y) x*y
|
エラー: "identifier"が見つかりません。"("が見つかりました。
0960デフォルトの名無しさん
2006/10/27(金) 05:43:27デフォルトだと、青、緑、赤茶、みたいな設定になっていて、それを変えたいんです。
教えてください！！！お願いします！！！
0961デフォルトの名無しさん
2006/10/27(金) 06:57:20棒グラフの軸の設定はどうやったらできるかわかりますか。
例えば、yの値を0から10までにしたい・・など・・。
0962デフォルトの名無しさん
2006/10/29(日) 02:54:56やりたいことは
0:pi/4:pi
みたいなことなのですが、pi を含んでしまうもので。
0963デフォルトの名無しさん
2006/10/29(日) 03:14:55じゃ駄目？
0964デフォルトの名無しさん
2006/10/29(日) 04:57:25ありがとう
0965デフォルトの名無しさん
2006/10/29(日) 22:43:30非負実数は黒く
負実数は白く
みたいなイメージのだしかたってありますか？
0966デフォルトの名無しさん
2006/10/29(日) 23:05:52a=rand(1,10)-0.5; %正負のある行列を用意
plot(find(a>0),a(find(a>0)),'*k',find(a>0),a(find(a<=0)),'*w') %正の数だったら黒く、負の数だったら白く*でプロット
0967デフォルトの名無しさん
2006/10/29(日) 23:08:25に訂正。ほんとにバカだ。ﾊﾊﾊ
0968デフォルトの名無しさん
2006/10/29(日) 23:16:310969デフォルトの名無しさん
2006/10/29(日) 23:19:24にお礼言うのを忘れてた．ありがとう．
0970デフォルトの名無しさん
2006/10/30(月) 07:54:48imreadすると自動的に縦×横×ＲＧＢの３次元配列になってしまいます。
これに１次元さらに加えて、それが何枚目の画像かを分かるようにするにはどうすればいいでしょうか？
0971デフォルトの名無しさん
2006/10/30(月) 12:04:50そのままその固有ベクトルを求めようとすると、
「??? エラー: ==> eig
EIG への入力は、NaN または Inf を含まなければなりません。」
と出ます。
どうしたらエラーを出さずに固有ベクトルを求めることができるのでしょうか？
0972デフォルトの名無しさん
2006/10/30(月) 23:03:31Aをその行列として，
imagesc(A, [min(A(:), max(A(:)))]);
colormap(gray);
0973デフォルトの名無しさん
2006/10/30(月) 23:04:40は
magesc(A, [min(A(:)), max(A(:))]);
に訂正
0974デフォルトの名無しさん
2006/10/31(火) 09:29:26色を全てを自由にって訳には行かないだろうが、
インスペクタでいじればOK.縦軸を0から10にするには
set(gca,'ylim',[0 10]);
カラーを直接いじる方法もあったがちょっと思い出せね。
>>970
4次元行列にしてもよし、セルを使ってもよし。
画像サイズが異なることを想定してセルにしたほうがいいな。
3次元データx,yがあった場合、
z{1} = x;
z{2} = y;
みたいな感じでOKよ。
>>971
相関行列にゼロ割りによるInfが含まれている。
相関行列計算式は省くが、分母が0になっていることが原因。
「EIG への入力は・・・」は言葉通り受け取ればよし。
つまり相関係数に0が含まれている公算が大だ。
回避する方法としては、入力行列に影響が出ない程度微細なノイズを
加えることが考えられる。詳しいデータ無いからこれ以上はわかんね。
0975デフォルトの名無しさん
2006/10/31(火) 22:05:18a='method'
と打ってみてびびった。
0976デフォルトの名無しさん
2006/11/01(水) 18:24:37読み込まないと思います。
・・その％の使い方ではなく、文字として「％」を表示させたいのですが、
どうしたいいのでしょうか。。。
0977デフォルトの名無しさん
2006/11/03(金) 11:04:360978デフォルトの名無しさん
2006/11/03(金) 23:26:19ありがとうございます。はい、￥％とか使ってしてみたんですが、うまく表示されなくて。。。
title('5%');
の場合はできるんですが、
title(sprintf('%d \%',j*5),'FontSize',12);
のように使うとできないんです・・・。
使っているVers.は7.0.1なんですが。。。
0979デフォルトの名無しさん
2006/11/05(日) 12:34:40みたいな等比数列を作りたいのですが、
2:3:19
的に楽にできる方法ってありますか？
0980デフォルトの名無しさん
2006/11/05(日) 14:13:400981デフォルトの名無しさん
2006/11/05(日) 14:14:172 .^ (1:10)
0982977
2006/11/05(日) 14:15:07だから「％％」だってば。
0983デフォルトの名無しさん
2006/11/06(月) 02:26:42なるほろけっと。とんくす
0984デフォルトの名無しさん
2006/11/06(月) 02:27:39もちろん win の話です。
0985デフォルトの名無しさん
2006/11/06(月) 02:28:160986デフォルトの名無しさん
2006/11/06(月) 02:39:24とかどう？今手元にMatlabないから試してないけど。
0987デフォルトの名無しさん
2006/11/06(月) 04:37:54ほぇー。そんな手があるとは。動いた。ありがとう。
0988デフォルトの名無しさん
2006/11/06(月) 04:40:25!explorer .
でできた。なんにしてもありがとう。
0989デフォルトの名無しさん
2006/11/07(火) 00:36:06先日ＭＡＴＬＡＢを使い始めまして
ＭＡＴＬＡＢでどうにか公共交通機関のシミュレーションを行いたいと思っています。
線の上で四角が移動するというようなプログラムを作りたいのですがどうすればよろしいでしょうか？
よろしくお願いします。
0990デフォルトの名無しさん
2006/11/08(水) 02:20:060991デフォルトの名無しさん
2006/11/08(水) 14:06:18先日MATLABを使い始めて、取り込んだ筆跡データをSOMでマップ化（分類）することには成功したのですが・・・
この先、どうやって認識させればいいのでしょうか？
誰か、アドバイスを宜しくお願いします。
0992デフォルトの名無しさん
2006/11/08(水) 15:52:52ってか matlab スレ。
0993デフォルトの名無しさん
2006/11/08(水) 16:07:59レスありがとうござます。
どこで聞けばいいのか分からなかったので、すいませんでした。
classificationってなんでしょうか？
あと、認識→認証の間違いでした。
0994デフォルトの名無しさん
2006/11/08(水) 22:42:48認識、認証は対した問題じゃないけど
Classification(分類)とClustering(クラスタリング）はかなり意味が違うこの分野の用語
>>991 でマップ化（分類）といっているけど、SOM は Unsupervised(教師なし）のはずで、
たぶん教師なしで、”クラスタリング”したんじゃないかと思う。
そこを説明なしで区別したかったので、分類という言葉を使いたくなかった。
0995デフォルトの名無しさん
2006/11/08(水) 23:28:50レスありがとうございます。
用語が全然分かってませんでした。
あと、質問文が間違ってました。
筆跡データをMATLABのSOM TOOLBOXというものを使い、クラスタリングしてマップ化したものから、
新たに筆跡データに入力して一番類似しているデータを探し出すのはどうしたらいいのでしょうか？
0996デフォルトの名無しさん
2006/11/09(木) 02:56:37ってか matlab スレ。
0997デフォルトの名無しさん
2006/11/09(木) 02:57:200998デフォルトの名無しさん
2006/11/09(木) 10:19:420999デフォルトの名無しさん
2006/11/09(木) 12:06:531000(´･д･`)DOn'03 ◇JcpHA18Tqs
2006/11/09(木) 12:07:391001１００１
Over 1000Threadもう書けないので、新しいスレッドを立ててくださいです。。。
レス数が1000を超えています。これ以上書き込みはできません。