MATLABプログラミング　質問箱　その2

1 名前：デフォルトの名無しさん [2006/11/10(金) 11:28:22 ] MATLAB　日本語版 www.cybernet.co.jp/matlab/ MATLAB 本家 www.mathworks.com/ Yahoo MATLAB dir.yahoo.co.jp/Science/Mathematics/Software/MATLAB/ MATLAB入門 feedback.mech.fukui-u.ac.jp/matlab/matlab0.htm MATLAB チュートリアル www-isys.eei.metro-u.ac.jp/~sho/pub/ttmatlab/Index.html 365 名前：デフォルトの名無しさん [2008/01/19(土) 12:47:42 ] polyfit使えば一発。ただ原点に返って正規方程式を逆行列で解くのもありだな。 366 名前：デフォルトの名無しさん [2008/01/19(土) 14:03:31 ] >>364 MATLAB以前の問題なのか… >>365 polyfitって、フィットさせる式に制限があるからね。 367 名前：デフォルトの名無しさん [2008/01/19(土) 21:23:54 ] >>365 ployfitって単なる多項式しか出来ないんですよね >>366 MATLAB以前の問題かもしれないんですけど、 やるとしたらデジタルで出てきたデータをz変換してそれを行列使って最小二乗してやるということになるのでしょうか？ 368 名前：デフォルトの名無しさん [2008/01/20(日) 15:18:25 ] だから単純なたとえばy = a*x^2 + …とか>>360に書いてあることはわかるんだが、 伝達関数の場合のプロセスが分からないっていうことなのだが教えてもらえないのか？ なにを調べたらいいか分からないからここに来たんだが 369 名前：デフォルトの名無しさん [2008/01/20(日) 18:50:24 ] >>368 あなたはどういうデータを持っていますか？ 理論的には、どういう操作をしますか どういう式にあわせて、係数を推定したいですか？ 370 名前：デフォルトの名無しさん mailto:sage [2008/01/20(日) 19:08:25 ] わかった具体的に質問しよう まず>>357のような伝達関数で出力されたデータを>>360のような行列表記するとどうなるのか そもそもこのやり方以外に方法があるのならどういったものがあるのか ということをよければ教えてくれ 371 名前：デフォルトの名無しさん [2008/01/20(日) 20:11:17 ] 実空間で解いたらいいんじゃないの？ ２次の遅れ系なら解は２つの指数関数の和だろ？ それに入力パルス（今の場合e^(-τs) なのかな？）との叩き込み積分が掛かっているだけ。 それをコスト関数にしてfminsearchにかければいい。 行列表記ができるのはpolyfit関数なんかが使える線形な式のみじゃないの？ 372 名前：デフォルトの名無しさん [2008/01/21(月) 13:58:02 ] 実空間でfminsearch使ったら出来そうだけど、ζの大きさによって式が変わるから、 せっかくz空間で来まった式で表現できるG(z)がもったいないなぁと思ったんですが、そんなものですか？ あと、e^(-τ*s)っていうのはむだ時間を考えたときのです。いるかいらないかは実験データを見てから決めようと思うんですけど。 今functionの使い方を勉強してる段階ですが、fminsearchでやってみます。 373 名前：デフォルトの名無しさん mailto:sage [2008/01/21(月) 23:06:45 ] 370==372 なのか？ なにか教えてもらえると思うとすぐに下手に出るのか。 374 名前：372 [2008/01/21(月) 23:24:44 ] というよりイライラしていたら言葉が悪くなるんだよ 375 名前：デフォルトの名無しさん [2008/01/22(火) 21:13:48 ] 問題をあなたが理解していないのでは、と思わざるを得ない。 なぜなら、あなたが何をしたいのかが、伝わらないからである。 伝達関数の求め方はいろいろあるんだよ。 376 名前：デフォルトの名無しさん [2008/01/22(火) 23:54:33 ] 何をしたいか：伝達関数を出したい 教えて欲しいこと：伝達関数の求め方 求め方として最小二乗かなと思って>>353と書き込んだんだがどこが伝わらなかったんでしょうか？ 377 名前：デフォルトの名無しさん [2008/01/23(水) 09:50:12 ] >>376 それで伝わっていると思うのなら、 あんたの理解が甘い。 先生に教えてもらえ。 378 名前：デフォルトの名無しさん [2008/01/23(水) 13:29:48 ] linux(openSuSE 10.3)でmatlabを使おうとするとjavaのAssertionに引っかかって 死にます。java1.5だとそうなるらしいのですが、-nojvmというオプションを渡すと GUIなしで起動するという話もみました。でも_JAVA_OPTIONという環境変数に-nojvmを 渡しても、Matlabが起動しません。(matlabのactivationはまだできてない)。どなたか同様の問題を乗り越えた方はいませんか？ 379 名前：デフォルトの名無しさん mailto:sage [2008/01/23(水) 18:24:49 ] polarでMarkerSizeを変更する方法を教えてください 380 名前：デフォルトの名無しさん mailto:sage [2008/01/23(水) 20:08:24 ] ふつうのグラフで軸をスケーリングするには AXIS([XMIN XMAX YMIN YMAX]) とすればできますが、 subplotを使ったグラフの軸をスケーリングするには どうしたらいのでしょうか？ 381 名前：デフォルトの名無しさん mailto:sage [2008/01/23(水) 21:17:09 ] >>380 同じ 382 名前：デフォルトの名無しさん mailto:sage [2008/01/23(水) 21:51:14 ] >>381 同じでできるんですか？ AXIS([XMIN XMAX YMIN YMAX])としてもできなかったので 聞きにきたんですが。もう一度試してみます。 もしこれでできない場合、他に方法はないでしょうか？ 383 名前：デフォルトの名無しさん [2008/01/23(水) 23:35:12 ] >>382 subplot(211) plot(x,y) axis([xmin xmax ymin ymax]) % subplor(212) plot(x,y) axis([xmin xmax ymin ymax]) 384 名前：デフォルトの名無しさん mailto:sage [2008/01/24(木) 09:38:12 ] >>382 subplot(2,1,1); plot(x, y); set(gca, 'xlim', [xmin xmax]); set(gca, 'ylim', [ymin ymax]); でもいけるのでは？ 385 名前：デフォルトの名無しさん mailto:sage [2008/01/24(木) 13:55:38 ] >>381>>383-384 できました。ありがとうございました！ 自分が馬鹿でした。 axis設定してからplotしてました。 plotしてからaxisなんですね。 386 名前：デフォルトの名無しさん [2008/01/27(日) 13:14:09 ] >>384 set(gca, 'xlim', [xmin xmax], 'ylim', [ymin ymax]); で桶では？ 387 名前：デフォルトの名無しさん [2008/01/27(日) 18:20:12 ] 逆ラドン変換iradon()を使うとデフォルトでRampフィルタがかかってしまうのですが フィルタ無しの方法はありませんか？ 海外のMATLABのヘルプでは値でNo filteringにできるそうなんですがバージョンが 古いのかバージョン7.2では出来ません。 388 名前：デフォルトの名無しさん [2008/01/27(日) 20:46:38 ] 逆投影なんかツールボックスを使わなくても数行でかけるのに(´･ω･`) 389 名前：387 [2008/01/28(月) 04:32:33 ] なんとかヒントでもよければ教えてもらえないでしょうか？ 専門外のことなのでどうにもこうにもツールに頼ってしまってorz 390 名前：デフォルトの名無しさん mailto:sage [2008/01/28(月) 09:50:41 ] edit その関数 で、中身見てフィルタ部分カットしたのを書くのはどうですか？ 391 名前：デフォルトの名無しさん [2008/01/28(月) 12:49:16 ] まあ、簡単に考えると配列の方をirotateで回転させ、そこに加算していけば２行で書けるんだけどね。 それで角度でループさせるの。 iradonの中を触るんだったら、Hというのがフィルター関数なので、この部分を取り去るか、1にしてしまうか。 勉強のために１から作りたいんだったら、Ｃ言語で出ている教科書を買って、ネチネチと移植してもいい。 おれは院の先輩の薦めで「SPECT画像再構成の基礎」（医療科学社）という本の中に出ているフーリエ変換法のプログラムを移植したよ。 392 名前：デフォルトの名無しさん mailto:sage [2008/01/29(火) 00:08:17 ] グラフの大きさを変えるにはどうすればいいのでしょうか？ いつもマウスでウィンドウのサイズを変えて 縦横比を見やすく調整しているんですか コマンドでできないでしょうか？ 393 名前：デフォルトの名無しさん mailto:sage [2008/01/29(火) 02:07:23 ] position 394 名前：デフォルトの名無しさん mailto:sage [2008/01/29(火) 12:54:18 ] set(2,'position',[30 246 1120 420] って感じですね？ できました！ ありがとうございます 395 名前：デフォルトの名無しさん mailto:sage [2008/01/29(火) 15:41:21 ] figure(1); print -deps 〜.eps のようにグラフをepsファイルに出力してるんですが その出力したepsファイルのサイズ（縦横比）を変えることはできないでしょうか？ set 'position'でできるかと思ったらできませんでした。 396 名前：デフォルトの名無しさん [2008/01/29(火) 16:20:06 ] 自由落下した物体が床で跳ね返る時、速度の符号が 変わりますよね。その時ゼロクロッシングとかいう エラーでちゃいます。ゼロクロッシングの検知を 無効にするとうまくシミュレーションできません。 何でですか？ 397 名前：デフォルトの名無しさん mailto:sage [2008/01/29(火) 20:45:50 ] 483 名前：名無しさん＠５周年[] 投稿日：2008/01/29(火) 16:04:56 ID: 自由落下した物体が床で跳ね返る時、速度の符号が 変わりますよね。その時ゼロクロッシングとかいう エラーでちゃいます。ゼロクロッシングの検知を 無効にするとうまくシミュレーションできません。 何でですか？ 398 名前：デフォルトの名無しさん [2008/01/30(水) 16:50:48 ] set(gcf,'PaperPositionMode','auto')かな？ 399 名前：395 mailto:sage [2008/01/30(水) 20:04:51 ] >>398 printに-zbufferとか-openglオプションってのを つけてみてもできなくて困っていたのですが ずばりそれでできました！ありがとうございます 400 名前：デフォルトの名無しさん [2008/02/05(火) 13:58:03 ] MEXをVisualC++2008ExpressEditionで 開発するときの設定の仕方ってどこか にありませんでしょうか？m(_ _)m 401 名前：デフォルトの名無しさん [2008/02/05(火) 17:04:41 ] x=[0:1000]に対応するf(x)が手元にあって 各xについてx-x/10<=x <= x+x/10の範囲でｆの移動平均をとって平滑化をしたいんですが ループを使わない方法はないでしょうか？ 402 名前：デフォルトの名無しさん mailto:sage [2008/02/06(水) 21:10:33 ] 100*100*100のxyz空間を用意して,各点に対して同じ100*100*100のAという行列の値を カラーマップに対応する色で描画するにはどのようにすればいいのでしょうか？ X = 1:100;　Y = 1:100;　Z = 1:100; [XX,YY,ZZ] = meshgrid(X,Y,Z); によってxyz空間とそれぞれの点に値を持つ100*100*100のA行列は用意できている状態です。 例えばカラーマップの範囲が0(青)〜10(赤)でありA(50,50,50)=10の場合 出力されるfigureの(50,50,50)の部分を赤に色付けしたいです 403 名前：デフォルトの名無しさん [2008/02/08(金) 03:53:49 ] Control System Toolboxを使われている方に質問です。 ゼロ次ホールドの連続時間伝達関数である、 H = (1-exp(-sT))/s ;　　T:サンプリング周期[sec] という表現をmatlab上で再現するにはどのようにしたらよいのでしょうか？ s = tf('s'); を使って直接表現しようと試みたのですが、上手くいきませんでした。 連続時間システムの制御器の設計及びコンピュータへの実装の際に 使用したいと考えています。 よろしくお願いします。 404 名前：デフォルトの名無しさん [2008/02/08(金) 12:02:16 ] 構造体にセットしたデータをfwriteとかでバイナリファイルに書き込むのはどうしたら良いのでしょうか？ 405 名前：デフォルトの名無しさん [2008/02/12(火) 21:18:32 ] 素人質問ですが教えていただければ幸いです。 10個の自動計測されたデータが列方向に並んでいるcsvファイルを csvreadで一気に読み込んだ後に(例えばAとする)、 10個の列要素に小分け(A1〜A10)しようと考えています。 各要素の1行目は説明なので2行目から最終行までを選択して取り込むには どうすればよいのでしょうか? 思いつきでA1=A(:,2:inf)とかしようとしたらエラーがでてしまいました。 多分while文とかでNaNがでるまで読み込み続けるとかすると思ったのですが 力尽きてしまいました。 やさしい人。お願い。 406 名前：デフォルトの名無しさん [2008/02/12(火) 21:48:36 ] infじゃなくってendじゃないの？ 407 名前：デフォルトの名無しさん mailto:sage [2008/02/12(火) 22:15:32 ] >>405 よく読まずに書いてしまうが、infではなくendを使ってみたらどう 408 名前：デフォルトの名無しさん mailto:sage [2008/02/12(火) 23:12:03 ] >>２　請けた仕事について、２ｃｈで無料で意見伺う無神経さｗ というか綱渡りで仕事してるからけつに火がついてしまうのに懲りてないな。 ベンチャー高専隊 409 名前：デフォルトの名無しさん mailto:sage [2008/02/12(火) 23:47:16 ] A(:,1)=[] で１列消せるはず 410 名前：405 [2008/02/13(水) 00:47:05 ] 407さん、ありがとうございます！おっしゃる通りで解決しました。 さらに色々触ってみて疑問が沸いたので厚かましく質問させてください。 （長文ですみません、誰も呼んでくれないだろうなぁ・・・） 例えば、 ----------------- time　data1 data2 1 10 100 2 11 99 3 10 102 4 9 98 | | | ----------------- というデータをcsvファイルから読み込んで、 data1とdata2を1/1000倍してから csvファイルとかExcelファイルに上と同じ書式でE書き込もうとして 下記のコードを書いてみました。 結果、xlswriteではうまく動くのですが、csvwriteではエラーが出ました。 セル配列を使えなかったようです。 私が自宅で使っているMATLABはバージョンが古く(ver.6.5)、xlswriteが使えません。 できれば古いバージョンでも動くファイル処理関数で対応できないものかと思っています。 411 名前：405 [2008/02/13(水) 00:48:21 ] (続き）書いてみて失敗したコードです。ご指導お願いします！！ ------------ clear all Alldata = csvread('csvsample.csv', 1, 0); Time = Alldata(2:end, 1); Data1 = Alldata(2:end, 2); Data2 = Alldata(2:end, 3); Data1 = Data1/1000; Data2 = Data2/1000; Dataset = [Time Data1 Data2]; Dataset_cell = num2cell(Dataset) Label = {'Time', 'Data1', 'Data2'}; Outdata = {Label ; Dataset_cell}; %xlswrite('newxls', Outdata) csvwrite('newcsv.csv',Outdata) 412 名前：デフォルトの名無しさん mailto:sage [2008/02/13(水) 01:28:43 ] みなさんMATLABのfigureをWordとかPowerPointにはるのにどうしてますか？ 今まではBMPで保存してXnViewの自動トリミングで余白を削ってはってました。 ただ、BMPだとWordでグラフを拡大縮小したりするとフォントがかすれたりして見た目がかなり汚くなります ベクターフォーマットのemfファイルで保存するのが一番きれいなのですが、グラフの余白を簡単に削る方法が見つからないです XnViewでもemfファイルを読み込めるんですが、XnViewで自動トリミングしてemfファイルに保存すると、 ファイルサイズが10倍以上になって、しかもビットマップデータみたいに拡大縮小すると汚くなります グラフの余白を簡単に削れてきれいに貼り付けられるうまい方法とかないですか？ 413 名前：405 [2008/02/13(水) 02:19:36 ] 自己解決しました。スレ汚しすみませんでした。 泥臭く下記のように書いてみました。 データ整理するのにxlsreadとかxlswriteって便利だなぁと思いました。 --------- clear all Alldata = csvread('csvsample.csv', 1, 0); Time = Alldata(2:end, 1); Data1 = Alldata(2:end, 2); Data2 = Alldata(2:end, 3); Data1 = Data1/1000; Data2 = Data2/1000; Fid=fopen('newcsv3.csv', 'wt'); Label = {'Time', 'Data1', 'Data2'}; for i=1:length(Label) fprintf(Fid, '%s,', Label{i}); end fprintf(Fid,'\n'); Dataset = [Time, Data1, Data2]; [Rownum, Colnum] = size(Dataset); for i=1:Rownum for j=1:Colnum fprintf(Fid, '%d,', Dataset(i,j)); end fprintf(Fid, '\n'); end fclose(Fid); 414 名前：デフォルトの名無しさん [2008/02/13(水) 07:14:12 ] すいません、皆さんみたいにハッカーになるためにはどんなことが必要ですか？ 415 名前：デフォルトの名無しさん [2008/02/15(金) 12:54:18 ] >>412 グラフの大きさを指定して、拡大縮小せずに使えるようにする 416 名前：デフォルトの名無しさん [2008/02/15(金) 15:21:13 ] MATLABでドロネー分割をつかって、３角形のパッチを利用した形状モデルを作っているのですが、パッチにテクスチャを貼り付けるにはどうすればいいでしょうか？ 417 名前：デフォルトの名無しさん mailto:sage [2008/02/15(金) 16:11:34 ] 3次元の立体形状から有限要素法用のメッシュを作り出す方法を教えてください 418 名前：デフォルトの名無しさん [2008/02/15(金) 17:19:00 ] テクスチャマッピングの仕方を教えて欲しいです。 わかりやすいサイトとかがあればお願いします。 よろしくお願いします。 419 名前：デフォルトの名無しさん [2008/02/15(金) 18:09:31 ] 要求が高いなぁｗ まず、www.google.co.jp。日本語以外のサイトも調べよう。 420 名前：デフォルトの名無しさん [2008/02/18(月) 22:58:56 ] ２次元のグレースケールの画像を、Z軸がピクセル値で３次元表示させたいんですけど、どうすればよろしいですかね？ 421 名前：デフォルトの名無しさん mailto:sage [2008/02/19(火) 04:52:22 ] >>420 surf 422 名前：デフォルトの名無しさん mailto:sage [2008/02/20(水) 16:17:32 ] すいません， a = 'Hello World' を b[1] = 'Hello' b[2] = 'World' に分割するコマンドはありませんでしょうか． 423 名前：デフォルトの名無しさん mailto:sage [2008/02/20(水) 16:40:36 ] doc regexp 424 名前：デフォルトの名無しさん [2008/02/20(水) 16:55:39 ] >>422 [ str,rem ] = strtok( a,' ' ); これで、最初のスペースまでの単語がstrに入り、残りがremに入る。 remの中をさらに切り分けたければこれを繰り返す。 str系の関数はC言語とほぼ同じだな。 425 名前：デフォルトの名無しさん [2008/02/21(木) 23:23:05 ] WAVファイルをパワースペクトル表示する方法を教えてください お願いします 426 名前：デフォルトの名無しさん [2008/02/22(金) 01:14:46 ] ぐぐれ 427 名前：デフォルトの名無しさん mailto:sage [2008/02/23(土) 23:13:43 ] ググっても見つからなかったんで質問させてください。 h1 = plot(x1) hold on h2 = plot(x2) h3 = plot(x3) ... といくつかのプロットを重ねた後、一つだけ(例えばx2のプロット)だけ消したい、または変更したいという場合、どうしたらいいですか？ 全部プロットしなおせばいいんですが、ちょっと重い上に何度もやるプロットなんで、全部やり直すというのは避けたいんです。 428 名前：427 mailto:sage [2008/02/23(土) 23:15:57 ] 書き忘れてました。プロットツールは使わず、関数で出来る方法で探しています。 429 名前：デフォルトの名無しさん mailto:sage [2008/02/24(日) 10:19:58 ] www.google.co.jp/search?hl=ja&rls=GGGL%2CGGGL%3A2006-17%2CGGGL%3Aja&q=matlab+line&btnG=%E6%A4%9C%E7%B4%A2&lr= 430 名前：427 mailto:sage [2008/02/24(日) 22:00:13 ] >>429 ありがとうございます！ 431 名前：デフォルトの名無しさん [2008/02/25(月) 11:05:16 ] グラフの凡例中で改行するにはどうすればよいでしょうか？ plot(X1) hold on plot(X2) plot(X3) legend({'X1';'test'},'X2','X3') とタイトルや座標軸中の改行と同じように入力しても無理でした。 432 名前：デフォルトの名無しさん [2008/02/25(月) 13:33:55 ] 100枚の画像あってサイズが縦240横352なんですが 横の左50画素、右46画素をカットして横幅を256画素に、 30×30の位置の色と同じものを上に縦16×横240くっつけて 256×256の画像にしたいのですが、どうすればいいでしょうか？ 433 名前：デフォルトの名無しさん mailto:sage [2008/02/25(月) 13:40:08 ] つ[ImageMagick] なんでもかんでもMATLABでやりたいならしらね。 つーか、MATLABエキスパートに任せた。 434 名前：422 mailto:sage [2008/02/25(月) 19:24:02 ] >>423 >>424 ありがとうございました！ 6.5なので423での'split'は無理みたいでした． 435 名前：デフォルトの名無しさん [2008/02/25(月) 20:47:37 ] 質問です！ ある画像の一部分（オブジェクト）の画素値の平均をとりたいのですが、 ラベル付けして、オブジェクトの全ピクセルの座標は取ったのですが、 for i=r(min):r(max) for j=c(min):c(max) heikin=mean2(Img(i,j,:)); だと、オブジェクトの形が四角ではないのでオブジェクトの外側の画素値まで 計算してしまいます。 何かいい方法がありましたら教えてください。 436 名前：デフォルトの名無しさん [2008/02/25(月) 21:04:19 ] >>435 オブジェクトがなんなのかわかりませんが、その内側と外側で違う特徴的なものをfindで引っかけてやればいいのでは？ たとえばオブジェクトの外側が-100より必ず小さいのであれば m = find( Img( : ) < -100 ); heikin = mean( Img( m )); の２行でできるはず。 Imgの範囲を絞りたいのだったら、 i = r(min):r(max); j = c(min):c(max); m = find( Img( i : j ) < -100 ); heikin = mean( Img( m )); forループは極力使うのは止めましょう。 437 名前：435 [2008/02/25(月) 22:08:15 ] ご回答ありがとうございます。 バッチリ上手く行きました。 差しでがましいですが もう一つ質問が… 二値化の画像にラベル付けをしてあるのですが、 ある座標を与えて、その座標を含むオブジェクトのラベル番号、 もしくはそのオブジェクトを構成している全座標を入手したいです。 [r,c] = find(bwselect(Img2,i,j)); これが僕の絞り出した答えでしたが無理でした。。 いいアイディアがありましたら教えてください。 438 名前：435 [2008/02/25(月) 22:08:54 ] ご回答ありがとうございます。 バッチリ上手く行きました。 差しでがましいですが もう一つ質問が… 二値化の画像にラベル付けをしてあるのですが、 ある座標を与えて、その座標を含むオブジェクトのラベル番号、 もしくはそのオブジェクトを構成している全座標を入手したいです。 [r,c] = find(bwselect(Img2,i,j)); これが僕の絞り出した答えでしたが無理でした。。 いいアイディアがありましたら教えてください。 439 名前：デフォルトの名無しさん [2008/02/26(火) 21:13:42 ] >>432 画像の読み込みをまずやってみよう 440 名前：デフォルトの名無しさん mailto:sage [2008/02/27(水) 11:35:03 ] フルカラーを使ってのムービーを作成する際のコーデックは何がおすすめですか？ 441 名前：デフォルトの名無しさん mailto:sage [2008/02/27(水) 13:22:39 ] 無圧縮AVI 後で別のソフトで変換 442 名前：デフォルトの名無しさん [2008/02/28(木) 22:32:06 ] 質問させてください！ 統計的に（最小２乗法などを用いて）２次近似モデルを式として作成する関数を探しています。 例えば実験にて２つの入力を調整し、１つの出力が得られるものがあるとします。 そのときの計測条件を[x,y]（２入力、ともにベクトル）、計測値を[z]（１出力、ベクトル）としたときの 計測値Zに対するX,Yの２次近似モデルを作成したいのです。 アウトプットイメージとしては Z=係数*X^2+係数*Y^2+係数*X*Y+係数*X+係数*Y+係数 のような感じです。 プロの皆さん、どうか教えてください〜 443 名前：デフォルトの名無しさん [2008/02/29(金) 03:55:17 ] >>442 z=[ z1,z2,......,zn ] '、x=[ x1,x2,......,xn ] ' , y = [ y1,y2,......,yn ] ' とし z = A*x^2 + B*y^2 + C*x*y + D * x + E*y + F とすれば a = [ x.^2 y.^2 x.*y x y ones( length( x ),1) ]; x = [ A B C D E F] ' ; とおいて x = a\z; で解けるはず。未知数が６つなので、データ数が少なくとも６つ以上ないと無理だけどね。 重み係数を使う必要がある場合は、ノルムを最小にするようにしないといけないけどね。 444 名前：デフォルトの名無しさん [2008/02/29(金) 03:56:57 ] >>437 find関数の中は条件文じゃないと機能しないのでは？ 445 名前：デフォルトの名無しさん [2008/02/29(金) 21:13:31 ] >>443 ありがとうございますー、やってみます！ 446 名前：デフォルトの名無しさん [2008/03/02(日) 12:22:49 ] >重み係数を使う必要がある場合は、ノルムを最小にするようにしないといけないけどね。 この意味が理解できないｗ 447 名前：デフォルトの名無しさん mailto:sage [2008/03/03(月) 17:50:26 ] >>446 ノルムというのは距離のこと 448 名前：デフォルトの名無しさん [2008/03/03(月) 18:43:02 ] 重み付き最初二乗法はどうやって記述できるの？ 449 名前：デフォルトの名無しさん mailto:sage [2008/03/03(月) 21:14:23 ] なぜmatlabスレで数式をきくんだ？ matlab への翻訳がわからないのなら数式は記述しておくとか、 数式が記述してあるページとかを示しておくべきなんじゃないのか？ 数式がわからないというのはスレ違いだろ？ 450 名前：デフォルトの名無しさん [2008/03/03(月) 21:32:41 ] >>449 ウルサイ！ 451 名前：デフォルトの名無しさん mailto:sage [2008/03/03(月) 22:27:27 ] シミュ板の scilab のスレッドが過疎で レスがいただけないようなので、こちらで質問させてください。 MATLAB で ( x^2 ) + ( x * y ) + ( y^2 ) = 1; 主軸が斜めになっている楕円は どのように plot2d すればいいのでしょうか？ 452 名前：デフォルトの名無しさん [2008/03/03(月) 23:41:30 ] >>447 いや、そもそも変なこと書いてるのではないか？ 453 名前：デフォルトの名無しさん [2008/03/04(火) 00:26:46 ] あるwaveファイルにQ値、ゲインをある特定の周波数に持たせたい （EQをかけたい）のですが、どうすればよいのでしょうか？ フィルタを使う（fdatool）ではダメですよね？ 宜しくお願いします 454 名前：デフォルトの名無しさん mailto:sage [2008/03/05(水) 23:19:31 ] >>453 ここを参考にすべし www.musicdsp.org/files/Audio-EQ-Cookbook.txt 455 名前：デフォルトの名無しさん mailto:sage [2008/03/06(木) 19:29:46 ] dl.cybernet.co.jp/matlab/support/manual/r2006/toolbox/matlab/visualize/?/matlab/support/manual/r2006/toolbox/matlab/visualize/chvolv19.shtml このような3次元の物体の可視化を行っているのですが、物体の表面で金属のような質感を出すにはどのように 設定したらよいのでしょうか？ 456 名前：デフォルトの名無しさん [2008/03/12(水) 16:33:14 ] 質問させていただきます。 「MATLABによる制御系設計」という本で勉強しているんですが、 MATLABの制御系コマンドを使えるフリーソフトを探しています。 検索すると無料でScilab,Octave,Freemat,有料でMATXというのがあるようですが、 制御系コマンドを使う時、互換性が最も高いなどお勧めの ものがあれば教えてください。 457 名前：デフォルトの名無しさん mailto:sage [2008/03/13(木) 00:42:00 ] 有料だけどMATLABっていうのがあるらしいよ 458 名前：デフォルトの名無しさん [2008/03/14(金) 11:53:20 ] >>457 456 もちろんMATLABは知ってますが・・・・ 研究でMATLABは使ってきたのですが、これからは 家でも勉強したいのであのような質問をしたのですが。 拙い質問と感じられたのでしたら申し訳。 459 名前：デフォルトの名無しさん mailto:sage [2008/03/14(金) 12:07:50 ] Scilabでいいんじゃね 460 名前：デフォルトの名無しさん mailto:sage [2008/03/14(金) 15:25:43 ] Octaveでいいんじゃね 461 名前：デフォルトの名無しさん mailto:sage [2008/03/14(金) 15:54:57 ] Freematでいいんじゃね 462 名前：デフォルトの名無しさん mailto:sage [2008/03/14(金) 21:24:08 ] とりあえず見つけたやつ全部自分で試して見ればいいんじゃね 463 名前：デフォルトの名無しさん [2008/03/15(土) 00:59:47 ] 学生ならMATLAB買え 464 名前：デフォルトの名無しさん [2008/03/15(土) 20:20:24 ] >462 ごもっともだと思います。 知り合いにもそういう意見をもらったんで、 そうしてみます。とりあえず、Scilabインストールしてみました。 金銭的に厳しいのでフリーのソフトから自分なりに試して みようと思います。 >463 それもそうですね。 MATLABをインストールできるだけお金が貯まったそうしようと思います。 ４月から社会人なので。 意見ありがとうございます。 465 名前：デフォルトの名無しさん mailto:sage [2008/03/21(金) 11:58:18 ] 前の職場でmatlab使ってたけど今のところでmatlab買えない pyrhon.numpy と maxima どっちがお勧め？

---

---

[ 続きを読む ] / [ 携帯版 ]

read.cgi ver5.27 [feat.BBS2 +1.6] / e.0.2 (02/09/03) / eucaly.net products.
担当:undef