MATLABプログラミング 質問箱 その2
at TECH
1:デフォルトの名無しさん
06/11/10 11:28:22
MATLAB 日本語版
URLﾘﾝｸ(www.cybernet.co.jp)
MATLAB 本家
URLﾘﾝｸ(www.mathworks.com)
Yahoo MATLAB
URLﾘﾝｸ(dir.yahoo.co.jp)
MATLAB入門
URLﾘﾝｸ(feedback.mech.fukui-u.ac.jp)
MATLAB チュートリアル
URLﾘﾝｸ(www-isys.eei.metro-u.ac.jp)
2:デフォルトの名無しさん
06/11/10 17:20:58
480*640の大きさの画像がn枚連続であって、画像の各画素における平均値、標準偏差、
標準偏差/平均値を求めるんですが、これであってるでしょうか？
tmpにはn枚目の画像の画素値(640*480分)が入っています。
分散は最終的にbunに入れてます。
for x=1:480
for y=1:640
heikin(x,y)=mean(tmp(x,y,:));
for i=1:n
bunsan(1,j)=(tmp(x,y,i)-heikin(x,y)).^2;
end
bun(x,y)=1/n*sum(bunsan);
end
end
cont=sqrt(bun)./heikin;
std=sqrt(bun);
3:デフォルトの名無しさん
06/11/12 10:51:58
なんかつっこみ所は多いけど、まぁいいんじゃね？
あってるかどうかの確認ぐらい自分でできるっしょ？
2x3x4 ぐらいのデータ作ってためせっつー話。
4:デフォルトの名無しさん
06/11/13 15:57:25
HSLでの赤から青までを使った、カラースケールの簡単な作り方はありませんか？
5:デフォルトの名無しさん
06/11/13 19:51:23
画像データを読み込むときに、ある範囲だけをいきなり抽出することはできるのでしょうか？
6:デフォルトの名無しさん
06/11/19 13:19:43
読み込んでから削りゃいいだろ
7:デフォルトの名無しさん
06/11/22 04:12:03
>>4
カラースケールマトリックスを作成の後、不要部分を削除
8:デフォルトの名無しさん
06/11/26 00:40:39
save関数で出力するディレクトリを指定することはできますか？
カレントディレクトリにしか出力できないのでしょうか
9:デフォルトの名無しさん
06/11/29 09:36:18
C++ で matlab チックに行列計算できるライブラリってある？
A(行列) - B(行列) とか他にも関数名がほぼ同じとか。
ちょっと急に思ったもので。
10:デフォルトの名無しさん
06/11/29 12:04:29
つ MATLAB API
11:デフォルトの名無しさん
06/11/30 18:44:51
sumっていう関数でベクトル要素の和を計算しようとしてるんですが、
何項目から何項目までの和みたいな指定はできるのでしょうか？
12:デフォルトの名無しさん
06/11/30 19:08:24
sum(v(1:10))
13:デフォルトの名無しさん
06/11/30 19:21:39
なぜかそれでできません。。
14:デフォルトの名無しさん
06/11/30 20:09:54
すいません。解決しました
15:デフォルトの名無しさん
06/12/02 02:27:46
たとえば
URLﾘﾝｸ(www.cfar.umd.edu)
の Figure 1 みたいな図ってなんという関数を使うのでしょう？
surf ではないようだし。
16:デフォルトの名無しさん
06/12/02 04:07:30
unix上のmatlabでmovieでアニメーションを作成しています。
movie2aviでは非圧縮しかないのですが、圧縮動画を造る方法はないでしょうか？
17:デフォルトの名無しさん
06/12/04 11:01:19
aviではなくquicktimeのmovを作成する方法はありませんか？
18:デフォルトの名無しさん
06/12/06 01:31:06
matlabで以下のような部分を書いたのですが
if 1<j
passenger_out(j,k)=trm_on(j,k);
stp_out(1,j)=stp_out(1,j)+trm_on(j,k);
trm_on(j,k)=0;
end
なぜか、trm_on(j,k)=0;を入れるとpassenger_out(j,k)まで０になってしって
passenger_outがゼロ行列になってしまうのですがなぜでしょうか。
よろしくお願いします。
19:１８
06/12/06 02:19:23
すみません。
自己解決しました。
簡単なミスを犯してました。
お騒がせしました。
20:デフォルトの名無しさん
06/12/07 10:15:07
責任とってｹｺｰﾝしる
21:デフォルトの名無しさん
06/12/07 10:21:27
俺のMATLAB7.01
mファイル入力中に突然エディタに何も入力することが出来なくなり
ワークスペース見たら
java.lang.IndexOutOfBoundsException
at java.util.ArrayList.RangeCheck(Unknown Source)
at java.util.ArrayList.get(Unknown Source)
（以下数十行に渡りエラー文）
こんな症状に覚えがある人いますかね？
22:デフォルトの名無しさん
06/12/07 19:50:56
秀丸使ったら？
23:デフォルトの名無しさん
06/12/07 20:07:46
トイレシミュレータ実装したやつ誰だよ
腹筋がおかしくなっちゃったじゃないか
24:デフォルトの名無しさん
06/12/08 21:21:38
MATLAB初心者です。
dw1dtoolで、左上のファイルから信号を読み込み、画像を表示するように、
specgramdemoでも同じように、ファイルから信号を読み込めるように
したいのですがどうすればいいのでしょうか....
なにせ初心者ですのでどこを参考にすればいいのかも分かりません。
よろしくお願いします。
25:デフォルトの名無しさん
06/12/09 01:30:30
>>21
SAKURA が断然お薦めです。
26:デフォルトの名無しさん
06/12/11 21:13:33
matlabR2006aです。
sys = tf([1 2], [1 2 3]);
感じで伝達関数が
伝達関数:
s + 2
-------------
s^2 + 2 s + 3
tf object ってデータ型みたいなんですが
guiでedittextに表示するためには
なんていう変換コマンドを使えばいいですか？(たとえばstr2doubleみたいな）
27:26
06/12/11 22:31:34
無理やり
diary on
sys=tf・・・・
diary off
で、string型のデータにすることで
自己解決しますた！
28:デフォルトの名無しさん
06/12/12 20:40:00
通信系の研究室で信号処理を学んでいる、matlab初心者です。
送信した信号を受信機でエネルギー検波(二乗検波)するプログラムを作りたいのですが、二乗検波の部分をどのように書いたらいいのかが分かりません…。
ちなみに、信号はUWBのパルス信号です。
どなたか、ヒントだけでもいいので教えて下さる方いらっしゃいませんか？
29:蟹男
06/12/13 00:17:39
matlabR12をwindowsXPにproductCDとDocumentCDとでインストールし、いざ起動しようとすると、
matlab6のロゴが一瞬画面に現れて終了してしまいます。
どなたか上のような症状を経験した人はいませんか？
OS：windowsXP Professional
matlabのVer.：matlabR12
ちなみに他のXPでは動作確認をしております。
なにか知っていたらお聞かせください。
30:デフォルトの名無しさん
06/12/17 18:59:52
MATLABって個人で買うと幾ら位するもんなの？
一応うちの研究室にもアカウントあるんだが
自分の家のPCでもやりたいと思ってて。
まぁあんまり高いようだと諦めますが。
31:デフォルトの名無しさん
06/12/17 19:21:07
>>30
うちから研究室のマシンにログインして使えばいいじゃん。
32:デフォルトの名無しさん
06/12/17 19:40:27
>>31
いやまぁ現状ではそうなんだけどね。
学割利くなら学生のうちに買っておけば
就職してからも使えるからいいかなぁと思って。
33:デフォルトの名無しさん
06/12/17 19:53:53
>>32
最近は、最低年に2回はバージョンが上がるぞ。
特に、2006bは変更が大きいと思うが。
そこまで考えなければ、そんなに変わらないが…
値段はしらね。
34:デフォルトの名無しさん
06/12/18 12:25:22
すみませんお聞きしたいのですが、
Matlabのプロットのプログラムをscilabのプログラムにコンバートしたいのですが、
どう書き換えれば良いでしょうか
q = [a0 a0*a1 a0*a2];
p = [1.0 b1 b2];
s = 0:0.01:pi;
subplot(2,2,2);semilogx(s, 20.0*log10(abs(freqz(q,p,s))));
axis([0 pi -20 20]);grid;
xlabel('Frequency \omega [rad]');
ylabel('|H(e~{j\omega})| [dB]');
これがMatlabのプロットのプログラムです。
よろしくお願いいたします。
35:デフォルトの名無しさん
06/12/20 05:36:41
matlab で 3D ポイントをプロットした後、
ぐりぐり回して好きなところでキャプチャできるようなやつって
あります？
matlab でなくてもいいです、この際
36:デフォルトの名無しさん
06/12/20 06:47:26
やつって何だよ バカか?
関数ならviewとprint使え
37:35
06/12/20 07:54:49
>>36
ばーか。plot3 のデフォルトでできるっつーの orz
38:デフォルトの名無しさん
06/12/20 22:09:14
>>34
亀レスですが・・・。
おまいMATLABもScilabも使うな！！
おまいにとっては猫に大判。
おまいさぁ〜
自助努力ってもんが ねぇーのかよ！！！！！
39:デフォルトの名無しさん
06/12/22 15:22:12
質問に答えられない人は書かなくていいです。答えられる人だけ書いてください。
40:デフォルトの名無しさん
06/12/22 19:25:15
>>38
質問に答えるくらいしてから罵ろうよ・・・。
41:デフォルトの名無しさん
06/12/22 20:34:51
>>34
m2sci使ってみた？
42:デフォルトの名無しさん
06/12/23 15:28:20
>>39
ここはmatlabの板
板ちゃうｗ
43:デフォルトの名無しさん
06/12/23 20:03:45
だったら誘導くらいしようよ。
罵るだけじゃなくてさ。
44:デフォルトの名無しさん
06/12/25 11:48:24
line線を書くことができますが、その際に、先端を矢印にしたいのですが
どのようにすればよいのでしょうか？
45:デフォルトの名無しさん
06/12/29 07:53:36
box offを標準の設定にする方法はありませんか？
46:デフォルトの名無しさん
07/01/12 08:32:07
このスレって全然機能してないな
47:デフォルトの名無しさん
07/01/12 08:54:07
ヽ / /⌒＼
／ヽヽ|/⌒＼ii|＼
／ ／ヾゞ//／＼＼|
|／ |;;;;;;|/ハ ＼|
|;;;;/／⌒ヽ
|;/（ ＾ω＾） >>1おっおっおっ乙枯ー
. |｛ ∪ ∪
|;;ヾ.,____,ノ
|;;; |
|;;;;;|
|;;;;;|
48:デフォルトの名無しさん
07/01/12 14:49:36
>>29 とりあえず同じ現象にあたった.。コマンドプロンプトで
matlabのディレクトリに移動してmatlab -nojvmで立ち上げると起動は出来る。
javaのバージョンによっておかしくなるみたいだけどよく判らない。
49:デフォルトの名無しさん
07/01/14 12:53:04
>2
劇遅レスだけど、
heikin = mean(tmp, 3);
hensa = std(tmp, 3);
cont = hensa ./ heikin;
でいいんじゃないか？
50:デフォルトの名無しさん
07/01/16 14:11:01
wavelet toolboxはどうやってインストールするんですか？
51:デフォルトの名無しさん
07/01/18 00:36:44
とある画像のピクセル色情報をRGB色空間にプロットってどうすればいいでしょうか？
52:デフォルトの名無しさん
07/01/20 01:36:45
大学の先生が作っているページって希にすっげー良いのがあるだろ？
そういうのをまとめてテンプレにしないか？
53:デフォルトの名無しさん
07/01/20 03:08:02
まず言いだしっぺのおまえからだ。
54:デフォルトの名無しさん
07/01/20 16:44:01
こんにちは
早速ですが質問させてください
関数内でプロットする場合に、入力変数で
xlabel('テキスト')のテキストを指定することってできますか
よろしくお願いいたします
55:デフォルトの名無しさん
07/01/20 20:40:12
function [result]=ploter(labelname)
label(labelname);
ploter('テキスト');
でいいんじゃないか？
とりあえずやってみれ
56:54
07/01/20 21:32:22
>>55
できました！！！
ありがとうございました。大変助かりました。
57:デフォルトの名無しさん
07/01/22 14:57:12
function [sys,x0]=boxp(t,x,u,flag,ts);
global Boxp
if flag == 2,
if any(get(0,'Children')==Boxp),
if strcmp(get(Boxp,'Name'),'daisya1'),
set(0,'currentfigure',Boxp);
H=get(gca,'UserData');
x=[u(1) u(1)+6*sin(u(2))];%振子
y=[-3 -3-6*cos(u(2))];%振子
set(H(1),'XData',x,'YData',y);
set(H(2),'XData',u(1),'YData',-2);%台車の位置
set(H(3),'XData',x,'YData',y);
drawnow;
end
end
58:デフォルトの名無しさん
07/01/22 14:58:27
axis([-10 10 -10 10]);
x=[0 0];
y=[-3 -9];
%H(1)=plot(x,y,'color','r','LineWidth',5,'EraseMode','background');
H(1)=plot(x,y,'LineWidth',5,'EraseMode','background');
H(2)=plot(0,-2,'s','color','r','MarkerSize',50);
H(3)=plot(x,y,'s','color','b','MarkerSize',20);
set(gca,'UserData',H);
plot([-7 7],[0.-1 0.-1],'color','k','LineWidth',15);
sys=[0 0 0 2 0 0];
x0=[];
end
end
実行速度を実際の実行スピードにしたいのですがどのようにすればいいでしょうか？
59:デフォルトの名無しさん
07/01/23 23:16:06
こんばんは
質問させてください
Y=X'*B*Xの対角成分が欲しい場合には、
diag(X'*B*X)ですが，対角成分しかいらないのに
X'*B*Xという計算をするのは不経済なので、
なんとかfor文を使わずに対角成分だけを計算する方法って
ありませんかね？
よろしくお願いします
60:デフォルトの名無しさん
07/01/24 06:52:51
T=1/fs;
[B2,F2,T2]=specgram(B,1024,1/T);
figure
f30=find(F2<30);f=[min(F2),max(F2(f30))];
t=[min(T2),max(T2)];
imagesc(t,f,abs(B2(f30,:)),[0,100]);colorbar;
specgramだと、色の値を各計測ごとの相対値でしか表せないから
（例えば一回目の計測の最大値が１００としたら赤が１００に対応し、
二度目の計測の最大値が６０としたら 今度は赤が６０に対応して表示されてしまう。）
imagescを用いてみたのですが、specgramの左下に原点のあるのグラフから
imagescのグラフは左上となってしまいます。
どうにかして元に戻す方法はありませんか？
61:デフォルトの名無しさん
07/01/24 07:58:07
MATLABのGUIってバグがありまくりじゃない？
エディットの数値を
handles.XXXXXX=str2double(get(hObject,'String'));
guidata(hObject, handles);
でセーブしても値が渡されない。
62:デフォルトの名無しさん
07/01/24 10:10:21
61です。 できました。
63:デフォルトの名無しさん
07/01/24 15:55:08
こんにちは。
等高線contourでライン幅を指定することってできますか？
よろしくお願いします
64:63
07/01/26 16:46:32
63です。
申し訳ありませんが，もう少し詳しく質問させてください．
等高線プロットcontourでラインスペックを指定するときに，
ヘルプでは
「contour(...,LineSpec) は、LineSpecで指定したラインタイプと
色を使ってコンターを描画します。」
と書かれていますが，
いざ，contour(X,Y,Z,'LineWidth',2)
とやると
??? エラー: ==> contours
入力引数が多すぎます
エラー: ==> C:\MATLAB6p5\toolbox\matlab\specgraph\contour3.m
行番号: 80 ==>[c,msg] = contours(varargin{1:nin});
エラー: ==> C:\MATLAB6p5\toolbox\matlab\specgraph\contour.m
行番号: 62 ==>[c,h,msg] = contour3(varargin{:});
エラー: ==> C:\MATLAB6p5\work〜.m
行番号: ** ==> contour(X1,Y1,Z1,'LineWidth',2);hold on
と出てしまいます．
どうすればラインスペックを指定できるのでしょうか．
どなたかご教示のほどよろしくお願いいたします．
65:デフォルトの名無しさん
07/01/26 18:02:44
その書式はMATLAB7以降でないと受け付けないようだ
6.5の場合はhハンドルとsetって
↓これでいけるはず。
[c,h]=contour(X,Y,Z);
set(h,'LineWidth',2);
66:デフォルトの名無しさん
07/01/27 22:56:50
質問させてください．
f(x,y)=0のような陰関数をプロットするにはどうすれば良いでしょうか．
よろしくお願いいたします．
67:64
07/01/27 22:59:59
>>65
ありがとうございました．
大変助かりました．
68:デフォルトの名無しさん
07/01/31 02:32:41
>>66
これでどうっすか？
for x=1:1000
y=find(f(x)<0.2 & f(x)>-0.2); %幅は適当に変更
plot(x,y,'.')
hold on
end
ところで、私も質問があります。
時系列の離散データをz変換したいのですが、どうすればいいんでしょうか？
ztransは関数で入力しなければならないみたいだし。
知ってる人、教えてください。
69:デフォルトの名無しさん
07/02/01 21:28:18
質問させてください
↓をベースにしてコッホ曲線を描くためにはどうすればよいのでしょうか。
%koch曲線表示プログラムkoch.m
%(x0,y0)から(x1,y1)までの線分についてのkoch曲線を描く
function y = koch(x0, y0, x1, y1, level)
if level == 0
% level が0
%(x0, y0)と(x1,y1)を結ぶ
%このとき、hold onとしておく必要がある→ 別の線分と重ね合わせ
else
%level が0でないとき
%(x0,y0)と(x1,y1)を3等分する→ (x2,y2)と(x3,y3)を増やす
%(x2,y2)(x3,y3)を結ぶ線分を60度回転し、(x4,y4)とする
% 生成した4辺(x0,y0) (x2,y2), (x2,y2) (x4, y4), (x4, y4) (x3, y3)
% (x3,y3)(x1,y1)についてlevel = level – 1としてkoch関数を呼ぶ
end
end
70:６９
07/02/05 14:26:18
質問させてください。
sys=[];
elseif flag == 4
ns =t/ts;
sys = (1 + floor(ns + 1e-1*(1+ns)))*ts;
elseif flag ==0;
animinit('daisya1');
[flag,Boxp] = figflag('daisya1');
axis([-10 10 -10 10]);
hold on;
このプログラムで、可変ステップ時はうごくのですが、
固定ステップにするとエラーが出てしまいます。
どのようにすればいいでしょうか？
71:初心者
07/02/11 16:27:50
MATLAB初心者です。
MATLABで下記のような離散時間状態空間モデルを
予測誤差法を使いパラメータ推定をするには
どうプログラムすればよいのでしょうか？
いままで、紙と鉛筆での計算はやってきましたが
プログラムは苦手なので、困っています。
離散時間状態空間モデル
x(k+1)=Ax(k)+Bu(k)+Ke(k)
y(k)=Cx(k)+Du(k)+e(k)
---------------------------
u(k)：入力，x(k)：状態
y(k)：出力，e(k)：誤差
初期条件として
u(k)=[2 2 3 3 5 5 3 3];
y(k)=[1 1 1 2 1 2 1 1];
です。
72:初心者
07/02/11 16:35:44
71番のMATLAB初心者です。
e(k)：誤差と書きましたが
e(k)：正しくはノイズです。
失礼しました。
MATLABで下記のような離散時間状態空間モデルを
予測誤差法を使いパラメータ推定をするには
どうプログラムすればよいのでしょうか？
いままで、紙と鉛筆での計算はやってきましたが
プログラムは苦手なので、困っています。
離散時間状態空間モデル
x(k+1)=Ax(k)+Bu(k)+Ke(k)
y(k)=Cx(k)+Du(k)+e(k)
---------------------------
u(k)：入力，x(k)：状態
y(k)：出力，e(k)：ノイズ
初期条件として
u(k)=[2 2 3 3 5 5 3 3];
y(k)=[1 1 1 2 1 2 1 1];
です。
73:７３
07/02/14 15:07:52
質問させてください。
実行してFigureに表示された動画をAVI形式に圧縮する方法は
どのようにすればよろしいでしょうか？
74:デフォルトの名無しさん
07/02/14 17:49:18
簡単ですよ。ビデオカメラで撮影してキャプチャーボードで取り込んでエンコーダーで
エンコードしてAVIにすれば完成です。
75:デフォルトの名無しさん
07/02/14 20:41:36
aviは圧縮形式じゃありません
76:デフォルトの名無しさん
07/02/15 08:30:47
圧縮形式のAVIを作りたいといういみでしょ? 圧縮用のコーデックを
好みで選んで圧縮すりゃいいんだよ。
77:デフォルトの名無しさん
07/02/16 01:26:29
>>73
MATLABムービーにしてmovie2aviか
URLﾘﾝｸ(dl.cybernet.co.jp)
ここ参照する
78:デフォルトの名無しさん
07/02/17 02:32:01
pythonからmatlabを呼ぶソフトがwindowsだとインストールでこける
pymatをpython2.4 か 2.5にコンパイルしたものをだれか持ってないですか？
79:デフォルトの名無しさん
07/02/17 18:42:56
URLﾘﾝｸ(claymore.engineer.gvsu.edu)
コンパイルには成功したけど
pythonの最新版なのが悪いのかmatlabとの通信がまったく機能してない
80:78
07/02/18 22:27:18
pymatだけど動いた
dctを使ったsampleプログラムが動かないだけだった
dctはtool boxないと動かないみたいだね
81:デフォルトの名無しさん
07/03/02 11:52:25
LineWidthのデフォルト値を設定することはできませんか？
82:デフォルトの名無しさん
07/03/04 23:08:10
こんばんわ
早速ですが質問させてください
function のなかで function を定義するのはどうすれば良いのでしょうか。
83:デフォルトの名無しさん
07/03/05 00:16:37
hoge.m
function hoge
foo
end
function foo
end
84:82
07/03/05 09:18:41
>>83
成程、関数の終了時にendが必要なのですね。
早速のご回答、ありがとうございました。
85:デフォルトの名無しさん
07/03/06 19:51:30
なくても動いたと思うが。
86:デフォルトの名無しさん
07/03/17 04:32:55
リモートデスクトップ環境ではmatlabは実行できないのでしょうか？？？？
87:デフォルトの名無しさん
07/03/18 06:57:16
おはようございます
早速ですが，1*pの行ベクトルxがあるときに，
xを縦にn行並べたn*p行列はfor文を使わないで生成できますでしょうか．
よろしくお願いします．
88:デフォルトの名無しさん
07/03/19 23:16:17
教えてください
Optimization Toolboxのfminuncで
関数y=func(x,z)をxについて最小化したいんですが，
どうすればよいでしょうか．zはこちらで与えたい定数です．
ヘルプには
X = fminunc(FUN,X0,OPTIONS,P1,P2,...) は、
問題に関連したパラメータ P1,P2,... を、直接、関数
FUN に渡します。
と書いてあったんで
x=fminunc(@func,x0,[],z)
というように，zを与えようとしたのですが，うまく行きません．
ご教示いただきますと幸いです．
よろしくお願いいたします．
89:デフォルトの名無しさん
07/03/20 08:41:43
>>87
repmat
90:デフォルトの名無しさん
07/03/20 20:26:35
>>87
ones(n, 1) * x;
91:デフォルトの名無しさん
07/03/28 17:20:31
複素成分の符号を逆（+なら-、-なら+）にする関数ってありますか？
92:91
07/03/28 17:22:06
間違えました。複素数の虚数部でした。
93:デフォルトの名無しさん
07/03/28 17:49:47
conj
94:デフォルトの名無しさん
07/03/28 17:52:01
最近は共役って言葉を数学で教えなくなったんだな
95:デフォルトの名無しさん
07/03/30 15:40:00
円周率は大体３、という教育に期待するほうが無駄
96:デフォルトの名無しさん
07/03/30 15:59:13
わけも分からず３．１４と教えられる方が問題あると思うよ
現代の小学生は正方形の中に
円をかいて４より小さいって自分で調べて３をみつけだしているからね
97:デフォルトの名無しさん
07/03/30 19:48:43
ご自分のときはどうだったのかな？
私のときは方眼紙に缶とかを利用して円を描き、
その面積を計算したりしたよ。
98:デフォルトの名無しさん
07/03/31 07:54:43
円を描いて32等分。
並べ替えて平行四辺形みたいな形状をつくり、縦ｘ横。
更に班内で平均。
99:デフォルトの名無しさん
07/04/10 17:24:07
構造体を与えると固有のパスコードが返されるような関数ってありますか？
例えば
a.b = 1;
a.c = 'test';
なときに
x = somefunc(a);
で
x = 'ZSGADRERASVFDSAD'
みたいなのが返ってくるような。
入力する構造体の構造と値が一致した場合のみ戻り値xも一致して、
そうでない場合は全然違う値が出てくる、というようなものを探して
います。何か良い方法を知っている方がいましたら教えてください。
100:デフォルトの名無しさん
07/04/13 18:33:46
MATLABやScilabで式長の制限はあるんでしょうか？
Scilabでかなり長い式を定義しようとしたら，
invalid factorと言われて困っております．．
101:デフォルトの名無しさん
07/04/14 09:23:29
>>99
そのものずばりなものは見たことないよ。
fieldnames()とかを使ってフィールド名を出してきたり、struct2cellで値を出してきたりして、あとはMD5なんかに突っ込むしかないんじゃないかなー。
102:デフォルトの名無しさん
07/04/18 13:22:38
初心者のしつもんです
simulinkで
y(0)=0
y'=1-xy
の微分方程式を作るにはどうやるのか教えてください
やってみたけどわかりませんでした
よろしくお願いします
103:デフォルトの名無しさん
07/04/18 13:49:42
すみませんスレ間違いました
出直してきます
104:デフォルトの名無しさん
07/04/23 18:07:27
ソースにIsim、placeってあるんですがどういう意味ですか？
105:デフォルトの名無しさん
07/04/24 09:19:10
あげ
106:デフォルトの名無しさん
07/04/28 14:31:42
isocolorsで3次元上の物体の表面に色をつけることができますが、同様に、alpha値を指定することはできないでしょうか？
107:デフォルトの名無しさん
07/05/01 16:54:04
一からMATLAB/simulinkを勉強しようと思っているのですが、
どういったところから始めるべきでしょうか？
解説書を見ながら、何かのデータをいじくっていけばいいのでしょうか？
108:デフォルトの名無しさん
07/05/02 18:10:44
clearで変数をクリアするとき、指定した変数以外をクリアする方法が知りたいです。
-regexpで正規表現を使うってのはわかるんですが
a ab abcって変数があるとき
clear -regexp [^\<ab\>]
とするとabcだけがクリアされちゃいます。
109:１０８
07/05/02 18:13:52
すいません変数abだけクリアしたいのです。
110:１０８
07/05/02 18:15:05
またまたすいません。
変数ab以外です。
111:デフォルトの名無しさん
07/05/04 15:19:36
１０８ですけど誰か答えられる人いないんですか?
112:デフォルトの名無しさん
07/05/04 16:50:58
clear allですべて削除した後、必要なものだけ再作成。
113:デフォルトの名無しさん
07/05/05 04:10:57
正規表現の話だろ？もはや matlab も関係ないよ。正規表現を勉強しろ、とだけ。
114:デフォルトの名無しさん
07/05/05 12:32:09
>>108
まあそれだとabcにマッチするよな当然。
そしてbはclearされないけどxはクリアされると。
115:デフォルトの名無しさん
07/05/06 04:50:08
リモートデスクトップでmatlabは使えないのでしょうか？
116:デフォルトの名無しさん
07/05/06 12:07:48
割れものだと使えません。
117:デフォルトの名無しさん
07/05/06 12:41:06
自宅で試行錯誤するだけならFreeMATでもいいんじゃないかな
URLﾘﾝｸ(freemat.sourceforge.net)
場合によっては octave-forge から持ってくるにしても
関数群が貧弱であれだけど
118:デフォルトの名無しさん
07/05/06 15:05:37
>>116
そうなの？
119:デフォルトの名無しさん
07/05/10 03:07:47
初心者です。質問させてください。
imhistでプロットしたヒストグラムの数値を
計算に使うにはどうしたらいいんでしょうか？
どこかに格納されてたりするんでしょうか。
お願いします。
120:デフォルトの名無しさん
07/05/10 03:42:55
ヘルプをみれば返り値として受け取れると書いてある
121:119
07/05/10 10:14:03
ありがとうございます〜。
122:デフォルトの名無しさん
07/05/26 08:23:39
3Dボリュームの可視化で、表面の色を液体金属のようにしたいのですが、
どのようなパラメータを設定すればよいのでしょうか？
123:デフォルトの名無しさん
07/05/29 17:47:05
gzipで圧縮されているデータを解凍してから表示しているのですが、
解凍作業をmatlab上から行えないでしょうか？
124:デフォルトの名無しさん
07/05/30 00:08:44
!でコマンドラインよべなかったっけ？
125:デフォルトの名無しさん
07/05/30 08:03:46
>>124
dosとかunix
126:デフォルトの名無しさん
07/06/01 02:50:30
matlabで3次元グラフを表示させてその表示した三次元グラフの表面積って計算可能ですか？
127:デフォルトの名無しさん
07/06/05 00:50:58
表面積の計算は自分でやってください。
128:デフォルトの名無しさん
07/06/08 23:27:58
関数外のグローバルな変数を関数内で用いるには、変数として代入するしかないのでしょうか？
function Y = f(X,Za,Zb,･･･)
if X=1
Y=sum(Za)
elseif X=2
Y=sum(Zb)
elseif...
というように必要なグローバル変数Z*が非常に多いので、代入以外に方法があればと考えています。
Z*を一つの変数にまとめることも考えましたが、すでにZ*それぞれが３次元変数なので、出来れば避けたいです。
129:デフォルトの名無しさん
07/06/09 05:32:17
>>128
matlab global でぐぐるくらいしようよ。一発だよ。
130:デフォルトの名無しさん
07/06/16 11:35:13
x -y
y x
の逆行列を求めるにはどうやるの？
131:デフォルトの名無しさん
07/06/16 19:30:47
くだらない質問でごめんなさい。
色々調べてみてもイマイチ理解が出来ていないんですが･･･。
matlabでローパスフィルタを作製しようと思ったときに参考にできるサイト
ってありますかね？
132:デフォルトの名無しさん
07/06/16 23:18:28
>>130
inv([x -y ; y x])
133:デフォルトの名無しさん
07/06/16 23:19:09
>>131
ローパスフィルタを作ること自体はfilter関数に何を入れてやればいいかだけ。
信号処理論の勉強をすべし。
134:デフォルトの名無しさん
07/06/19 15:53:58
hann関数を使用すると
未定義の関数、またはメソッド 'hann' (タイプ'double' の入力引数)。
と表示されます
w = hann(1024);
と使っているのですがうまく動きません
だれかご教授を
135:デフォルトの名無しさん
07/06/19 20:05:39
Signal Processing Toolboxが入ってないんじゃない?
136:134
07/06/19 21:31:13
多分それだと思います
すぐに確認できないんで明日確認して報告します
137:134
07/06/20 16:04:45
確認したらSignal Processing Toolboxがありませんでしたので
練習ついでにhanning関数作りました。
ありがとうございました。
138:デフォルトの名無しさん
07/06/20 21:02:08
hannning関数をつくるなんてすごいなぁ
139:デフォルトの名無しさん
07/06/21 08:27:30
その辺の関数ならOctaveのがそのまま使えるんじゃないの?
140:デフォルトの名無しさん
07/06/27 02:11:43
2次元データを平滑化する関数はありませんか？
141:デフォルトの名無しさん
07/06/29 14:48:56
forループで時間を進めて
plotで表示しているのですが、最後のplotしか表示されません。
途中経過も表示したいのですがどのようにしたらよいのか教えてください
142:デフォルトの名無しさん
07/06/29 20:39:13
figure
hold on
for i=1:n
plot(なんとか)
end
hold off
こんな感じ？
forの前にhold onしてforのあとにhold offかな？
143:デフォルトの名無しさん
07/06/29 20:41:16
あー、もしかしてこっちかな？
figure
for i=1:n
plot(なんとか)
drawnow
end
drawnowですぐ結果表示だからね。
質問をはっきりして欲しい。
144:デフォルトの名無しさん
07/07/01 07:16:55
>>141
ソースも書かずになんで原因が分かると思ったんだ?
>>142,143 はひまでしょうがないから
たまたまあてずっぽうで答えてくれただけだからな。
145:デフォルトの名無しさん
07/07/01 08:38:58
>>143でうまくいきました。
どうもありがとう。
146:デフォルトの名無しさん
07/07/01 12:06:46
Linux版で動画を作成したいのですが、movie2aviは非圧縮しか対応していません。
圧縮フォーマットは何でもよいのですが、何かよい方法はないでしょうか？
147:デフォルトの名無しさん
07/07/01 12:35:54
>>146
mplayer とか ffmpeg とかいろいろあるじゃん。
148:デフォルトの名無しさん
07/07/03 13:23:05
教授から、作成したmファイルがコンパイルできないで困っていると相談を受けました
ソースを解読しながら修正することになったのですが、如何がコンパイルできない、どの関数が実行できないなどの情報はどこかに集まっていないのでしょうか？
149:デフォルトの名無しさん
07/07/03 14:18:59
教授のくせに頭悪いな
150:デフォルトの名無しさん
07/07/03 17:41:19
せめてどんなエラーが出るか書くとかしてね。
151:デフォルトの名無しさん
07/07/03 21:56:27
ふーむ m ファイルをコンパイルか・・・
152:デフォルトの名無しさん
07/07/03 22:25:59
「移動平均」ならぬ「移動最大値」（or「移動最小値」）
を計算できるような方法はございますでしょうか？
（自分でforとmax使って作ると遅いので・・・）
153:148
07/07/04 11:49:11
>>149
どうもすいません。確かに前提のお話が足りていませんね。
>>150
言いたいことは分かりますが、ただ教授も神様じゃないんです。
一応教授から聞いた症状ですが
・使用したコンパイラはmcc
・環境はmatlab6.5でコンパイルできた
・matlab7.1にしたらコンパイルできなくなっていた（これは共同研究者の環境ら
しいです
と、ここまでは聞いた話です。忘れてください。
それで、自分で試した限りでは
・matlab6.5ではコンパイルできず
・matlab7.0ではexeファイルを生成できるけど、実行できない
となっていました。
プログラムは用意したmp3(またはwav)ファイルを読み込んで、その波形から特徴
のあるところを探すものらしいですが、400行ちょっとあって、どこが原因なのか
探すのが一苦労です。
今後もmファイルのコンパイルにはお世話になる可能性が高いので、何が障害にな
るか傾向なんかがあるとやりやすいのですが…
154:デフォルトの名無しさん
07/07/05 07:42:43
たったの400行で何言ってんだ……
155:デフォルトの名無しさん
07/07/05 08:50:00
おまえソースコードもってんだろ？エラーメッセージぐらいなんでかかんの？
156:デフォルトの名無しさん
07/07/05 18:35:14
なんで教授が２ｃｈできくわけ？
教授なら外国の掲示板で聞くと思うが
157:デフォルトの名無しさん
07/07/05 19:59:17
いつ教授が2chできいたんだ？
それに教授だからとなぜ2chで聞くのがおかしいんだ？
教授は崇高な存在だとでも思っているのか？
158:デフォルトの名無しさん
07/07/05 21:22:59
スカートの中を手鏡で覗く教授に似てるなと思っただけ
159:デフォルトの名無しさん
07/07/07 23:02:18
Final
Fantasy
Tactics
160:デフォルトの名無しさん
07/07/08 14:31:49
以前のhelpだと
inv は、LINPACK のサブルーチン ZGEDI と ZGEFA を使っています。
とありますが、
最近のは
invは、LAPACKルーチンを使って、逆行列を計算します。
行列 関数
実数 DLANGE, DGETRF, DGECON, DGETRI
複素数 ZLANGE, ZGETRF, ZGECON, ZGETRI
となっています。具体的には、何が違うのでしょうか？
161:デフォルトの名無しさん
07/07/08 15:13:50
対称スペクトルをifftしても虚部の時間信号で出てくるのはなぜでしょうか？
162:デフォルトの名無しさん
07/07/19 05:08:44
Matlabのアカデミックパックの価格って誰かご存じないでしょうか？
163:デフォルトの名無しさん
07/07/19 13:53:10
>>162
基本システムは，去年購入した価格で１０万円強だったが，価格改定されてるかも．
サイバネットシステムのサイトからメールで資料請求しなされ．
すぐにカタログと価格表を送ってくれるよ．
164:デフォルトの名無しさん
07/07/20 01:05:28
うむ。基本パックが11万円くらい。各ツールボックスが4万5000円くらい。
165:デフォルトの名無しさん
07/07/22 01:23:23
MATLABにLISPのmap関数の様なものはありますでしょうか?
166:デフォルトの名無しさん
07/07/23 01:09:35
>>165
なんでそんなものが必要だと思うのかを書かないの?
167:デフォルトの名無しさん
07/07/23 01:32:22
>>166
for文を使わずに，
'A.*B'や'A./B'の様に書きたいからです．
168:デフォルトの名無しさん
07/07/23 09:53:55
「何を」for文を使わずに'A.*B'や'A./B'の様にの様に書きたいのかを書いたら?
それは本当にlispのmap関数のようなものがないとできないのか?
169:デフォルトの名無しさん
07/07/23 10:17:18
>>168
おまえはまず興味があるのなら、Yes か No で答えてあげたらどうだ？
疑問文に疑問文で答えろと学校で教えているのか？ってやつだ
170:デフォルトの名無しさん
07/07/23 12:05:52
さいきんなぜだか知らんがこういうやつ多いよ
171:デフォルトの名無しさん
07/07/23 12:06:51
質問するなら何がやりたくて何が分からないのか全部最初に書く。
小出しにして質問者に質問させるような質問に親切に答えてくれる人は少ない。
>>169はそういう親切な人みたいだから、君が答えてあげたら?
172:デフォルトの名無しさん
07/07/23 12:11:23
>>171
あんたが質問の意味も理解できない馬鹿なだけだろ
173:デフォルトの名無しさん
07/07/23 12:13:49
mapでやるようなことをmatlabでは違う方法でやるのが普通。
>>172はちゃかしてるだけで結局答えることも
どう聞くべきかを教えてやることもしないんだろ?
174:デフォルトの名無しさん
07/07/23 13:10:04
普通じゃない方法はやるなと主張するわけだな？
175:デフォルトの名無しさん
07/07/23 14:57:43
>>174
そんなこと一言も言ってないのになんで勝手に決めるんだ?
でお前は質問に答えられもしないし
どう聞けばより答が得られやすいかを教えてあげるわけでもないんだろ?
176:初心者
07/07/23 21:36:26
MATLABのツールボックスを使わずに、
TSPの10都市問題をエネルギー関数を用いて解かせるプログラムを組みたいのですが、
全くわかりません。。
参考にできるサイトなどありましたらお教えいただけませんでしょうか。
177:デフォルトの名無しさん
07/07/23 21:41:48
pythonならmapなかった？
178:177
07/07/23 21:44:44
ちなみに
python matlab compatible
で検索すると幸せになれるかも
179:デフォルトの名無しさん
07/07/24 00:09:10
>>176
検索したら巡回セールスマン問題の解法アルゴリズムは
たくさん出てくると思うんだけど。
それをプログラムにしていく過程でまた具体的な質問が出てきたら
ここで質問したらいいと思うが、その質問だと全部ゼロから書いてくれ
って言ってるようなもんだぞ。
コツとしてはMATLABは忘れてpseudo codeでいいから求める手順を
まずしっかり書き出してみることだな。
>>177,178
matlabとmapで検索しまくって出てきた答?
180:デフォルトの名無しさん
07/07/25 02:09:20
>>176
traveling salesman annealingで検索したらそのまんまなのが見つかるな。
さすがに宿題でこのコード出すわけにはいかないだろうけど(笑)。
181:初心者
07/07/25 22:41:54
デフォルトの名無しさん、丁寧なご回答ありがとうございます。
参考にさせていただきます。
182:デフォルトの名無しさん
07/08/03 18:20:14
3つの列ベクトルx(i),y(j),z(k)をそれぞれx軸、y軸、z軸に設定してグラフを書きたいです。数式によって高さを与えるz平面のプロットは可能ですが実験で得られたデータを高さに与えて平面を書き出すことは可能ですか？
183:デフォルトの名無しさん
07/08/03 22:12:48
高さがあるのに平面とは…
184:デフォルトの名無しさん
07/08/05 04:35:53
mov2aviでPowerPointで再生できる動画を作成したいのですが、何を指定すればよいのでしょうか？
185:デフォルトの名無しさん
07/08/06 12:46:32
>>184
help mov2avi
何が出来て何が出来なかったのか言ってくれないと判らない
186:デフォルトの名無しさん
07/08/10 21:15:59
画像に、フィルターをかけるとき、フィルタのパラメータを調整したものを複数生成して
それらを一覧で表示したいのですが、
cat関数で横一列に繋げると、枚数が増えたとき表示のバランスが取れません。
何枚生成するかは、そのときそのときで違う場合で、よい方法はないでしょうか？
187:デフォルトの名無しさん
07/08/10 22:16:01
3次元のグラフというか立体図面を書きたい。
x,yの平面に高さzで、高さを面で現したような。
n行3列の行列があって(例えばdat)、行毎に[x,y,z]のデータの組が順不同に入っている。
例えば dat = [ 1 2 0.5; 1 3 0.6; 3 8 1.5; 1 1.5 0.6];
plot3(dat(:,1),dat(:,2),dat(:,3),'.') では上手く行かないし、surfaceを使えば良いのか、検討が付かない。
何か手掛かりだけでも、誰か教えてください。
188:デフォルトの名無しさん
07/08/11 02:38:41
>>186
で、そのフィルターが1個のときと1000個のときで
具体的にどう表示して欲しいんだ?
>>187
よく分からんけどsurfは? とりあえずhelp graph3dしろ。
189:186
07/08/11 18:33:27
>>188
えっと、ある程度表示がよくなればまったく問題ないのですが
画像の右上に、パラメーターの数値を書き込んで
それらの画像をおおよそウィンドウ内部に
おおよそ正方形くらいになるように並べられらばと思っています。
190:デフォルトの名無しさん
07/08/12 00:37:31
1個のときの表示、2個のときの表示、・・・と全部具体的にイメージできてる?
それが出きてればそれをプログラムに書けばいいだけ。文章から推測するに、
君はやろうとしていることを具体的に細分化して理解しようとしていない。
「おおまかにこんな感じ」ではプログラムにならない。
コンピュータはそんなにかしこくないんです。
191:デフォルトの名無しさん
07/08/24 10:49:42
すごい初歩的なことなのですが少数点以下をもつ数の累乗根を出すような関数やコマンドは何でしょうか？
192:デフォルトの名無しさん
07/08/24 11:14:08
ヘルプは当然見たんだろうけど、どうやって調べた?
193:デフォルトの名無しさん
07/08/25 03:05:55
>>192
helpでコマンド一覧みたいなのを見たのですが整数の累乗根を出すものしか見つけられなかったので質問させてもらいました
194:デフォルトの名無しさん
07/08/25 04:42:51
>>193
うちでは整数のみのそのコマンドが見つからないです。
195:デフォルトの名無しさん
07/09/12 18:12:42
等間隔で行を抽出するにはどうすればいいですか？
例えば、1行目、11行目、21行目…という風にです。
196:デフォルトの名無しさん
07/09/13 01:12:14
1:10:end
これじゃダメ？
197:デフォルトの名無しさん
07/09/13 18:58:31
matlabで箱ヒゲ図を描きたいです。
ヒゲは、errorbarを使えば何とかなるのですが、
箱を描く方法が思いつきません。
どなたか、ご教授してください。
URLﾘﾝｸ(www.sci.kagoshima-u.ac.jp)
こういう図が描きたいのです。
198:デフォルトの名無しさん
07/09/13 19:20:26
株価のローソク足に似てるなーと思って「MATLAB ローソク足」でググったら一番上にそれらしきものがでてきたよ
Financial Time Series Toolboxがいるらしい
そこでCandle Plotを描くだけならだれかが自作のmファイルを後悔してるかも知れないと思って
「MATLAB Candle Plot」でググったらCandle.mを無性公開しているサイトを発見したぜ
199:デフォルトの名無しさん
07/09/13 19:20:30
>>195
fid = fopen(filename);
while 1
line = fgetl(fid);
if length(line)==0
break;
end
// 任意の操作をここに。
for ii = 1:10 // 10行を読み飛ばす。
fgetl(fid);
end
end
fclose(fid);
でOK？
200:197
07/09/13 19:57:27
>>198
多謝。それです。
「ローソク足」というものがあるのですね。知りませんでした。
見つけていただいたコードを参考にしたいと思います。
以下、チラ裏。
plot( [x1 x2 x2 x1 x1], [uh uh lh lh uh] ));
で、箱が描けるとは。
matlab暦5ヶ月目で、新鮮な発見。
検索スキルとmatlab力を磨いて出直してきます。
201:195
07/09/15 16:14:20
>>196>>199
できました、ありがとうございます。
202:デフォルトの名無しさん
07/09/23 16:27:08
MATLABには行列内部の最小の値を取る関数はあるけど
最小の値のある場所のインデックスを返す関数ってないよね
203:デフォルトの名無しさん
07/09/23 16:40:45
同じ関数
204:デフォルトの名無しさん
07/09/23 19:07:22
help その関数
で、わかりますよ＾＾
205:デフォルトの名無しさん
07/09/25 11:19:11
ygridの位置を指定するにはどうすればいいのでしょうか？
206:デフォルトの名無しさん
07/09/26 00:01:34
もっと詳しく書いて欲しいところだが、
set(gca, 'YTick', [ ... ]);
とかでどう？
207:デフォルトの名無しさん
07/09/26 16:05:17
matlabのグラフで，プロットした波形にマウスを当てると，
その時点での縦軸と横軸の値が出ますが，
それを，コマンドプロンプトでやりたいんですが，
関数かなにかありますか？
検索しても見つからなかったので，
やはり自分でプログラムを組むべきなのでしょうか？
208:デフォルトの名無しさん
07/09/27 15:34:32
>>207
つ[gnuplot]
209:207
07/09/27 23:44:56
>>208
ご教授ありがとうございます．
gnuplotは，さわりだけしか使ったことないんで，
どうなるか分かりませんが，やってみます．
210:207
07/09/27 23:47:12
ということは，
そういったことができるmatlabコマンドはないんですね・・・
211:デフォルトの名無しさん
07/09/28 14:52:23
画像処理でアフィン変換できる関数が見当たらないのですが気のせい？
結構基礎的な関数の気がするんだけど
212:デフォルトの名無しさん
07/09/28 22:45:39
そのくらい自分でやれ、ってことじゃないの？
213:デフォルトの名無しさん
07/09/29 11:00:49
まとめサイトはありませんか？
214:デフォルトの名無しさん
07/09/30 15:08:48
v=1*50000 vector & A=50000*50000 matrix
v*AをやろうとするとOut of memoryになっていまいます。
メモリーを追加しないで、効率よく計算するにはどうしたらいいですか？
よいBuilt-in関数とかないっすか？
215:デフォルトの名無しさん
07/09/30 19:10:11
>>202
ヘルプ見なさい
216:デフォルトの名無しさん
07/09/30 22:31:40
>203,204,215
helpみました
吊ってきます…
217:デフォルトの名無しさん
07/10/01 17:27:37
Axes上に複数のLineオブジェクトを作成しました。
このLineをマウス操作で移動したり、リサイズしたりしたいと思ってます…。
Axes上にLineを引くときは、function axes1_ButtonDownFcn( )でいけました。
しかし、function axes1_ButtonDownFcn( )の中にselectmoverisizeを書くと、
Axes自体が移動・リサイズされてしまいました。
そもそも、Lineオブジェクトをselectmoveresizeで操作するのは無理なんでしょうか？
218:デフォルトの名無しさん
07/10/02 15:25:07
3次元プロットで、（例えばsurfなどで）値がaxisで指定した範囲外にあってもプロットされてしまうのですが、
axisで指定した範囲内のみプロットするようにするにはどうしたらよいのでしょうか？
219:デフォルトの名無しさん
07/10/02 16:06:06
関数を定義する場合必ず別ファイルにしないといけないのですか？
同一ファイル内で定義することはできないものでしょうか？
220:デフォルトの名無しさん
07/10/02 18:45:41
>>218
aってのを3次元の行列だとして
surf(a(xmin:xmax,ymin:ymax,zmin:zmax))
じゃだめ？
意味取り違えてたらごめんね。
>>219
ｽﾚﾘﾝｸ(tech板:82-番)
同じスレなので出来れば読んでね。
221:デフォルトの名無しさん
07/10/04 03:19:13
matlabってリーもとデスクトップじゃつかえないのだね
222:デフォルトの名無しさん
07/10/11 14:00:14
bata版のNAG ToolBox for MATLAB
無料評価版
URLﾘﾝｸ(www.nag-j.co.jp)
223:デフォルトの名無しさん
07/10/11 14:23:25
pythonで大部分の数値計算できるけどmatlabってやっぱりまだ需要あるの？
224:デフォルトの名無しさん
07/10/11 16:29:25
>>223
C で python のできることは全部できるけど python って需要あんの？
225:デフォルトの名無しさん
07/10/12 00:19:45
俺紙とエンピツで全部やってるけどコンピュータって需要(ry
226:デフォルトの名無しさん
07/10/12 05:57:48
さすがにMATLABのスレでそりゃないだろ。
手作業でやったら一生かかっても終わらん計算してるだろうに。
227:デフォルトの名無しさん
07/10/23 14:25:52
質問なのですが、下記のプログラムをsimulinkでブロック図にて表すにはどうすればよろしいでしょうか？
clear allsyms z n p T
syms z n p T
h=0.1;
A=[exp(-h) int(exp(-T),0,h);-0.5 1];
b=[0;0.5];
c=[exp(-p) int(exp(-T),0,p)];
G=c/(z*eye(2)-A)*b;
R=z/(z-1);
Y=R*G;
y=iztrans(Y);
t=1;
for m=0:9,
v(1,m+1)=subs(y,n,m);
for l=0:9,
w(1,l+t)=subs(v(1,m+1),p,l*1e-2);
end
t=t+10;
end
figure
plot(w)
grid on
次ページ最新レス表示スレッドの検索類似スレ一覧話題のニュースおまかせリスト▼オプションを表示暇つぶし2ch
4096日前に更新/230 KB
担当:undef