🎄teratailクリスマスプレゼントキャンペーン2024🎄』開催中!

\teratail特別グッズやAmazonギフトカード最大2,000円分が当たる!/

詳細はこちら
MATLAB

MATLABはMathWorksで開発された数値計算や数値の視覚化のための高水準の対話型プログラミング環境です。

Q&A

解決済

1回答

2175閲覧

MATLAB 変数の最適化を行うプログラムについて(for文使用中のエラー)

MT1111

総合スコア6

MATLAB

MATLABはMathWorksで開発された数値計算や数値の視覚化のための高水準の対話型プログラミング環境です。

0グッド

0クリップ

投稿2019/12/11 16:08

前提・実現したいこと

プログラミング初学者です.MATLABを用いて変数の最適化を行うプログラムを組んでいます.
最適化を行う関数ltは,αとβの2変数で構成されており,一つ前のタイムステップ値を積算する関数になっております.変数αとβについては,100*1の乱数行列を作成しました.for文を使用して,時系列データと変数ごとの値をそれぞれ行と列に表示させることを目標にしておりますが,以下のエラーメッセージが発生しました.

発生している問題・エラーメッセージ

位置 2 のインデックスが配列範囲を超えています (1 を超えてはならない)。 random_demo (line 16) lt(n,m) = lt(n-1,m) + rt_input(n,m);

該当のソースコード

MATLAB

1 clear all 2 3sm_2010_day=xlsread('snowmelt_2010.xlsx','C1:C365');%日データ読み込み 4sm_1hr=sm_2010_day/24;%日データ/24 5sm=repelem(sm_1hr,24);%時間データへの換算 6 7alfa=-0.001+(1+0.001)*rand(100,1);%変数α 乱数発生  8beta=0.001+(1-0.001)*rand(100,1);%変数β 乱数発生 9 10n_end=100; 11 12for n=2:n_end 13 for m=1:100; 14 rrt(n) = sm(n); 15 rt_input(n,m) = rrt(n) * exp( alfa(m) * rrt(n) ); 16 lt(n,m) = lt(n-1,m) + rt_input(n,m); 17 qt(n,m) = sqrt(lt(n,m))*beta(m); 18 lt(n,m) = lt(n,m) - qt(n,m); 19 end 20end 21

試したこと

ここに問題に対して試したことを記載してください。

補足情報(FW/ツールのバージョンなど)

ここにより詳細な情報を記載してください。

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答1

0

ベストアンサー

MATLAB

1lt(n,m) = lt(n-1,m) + rt_input(n,m);

のところですが、変数ltが作られていないのに、いきなりlt(n-1,m)で中身を見ています。細かいことを省いていえば、存在しないものを見ようとしているのでエラーになっています。

ltの初期値がすべて0でいいのなら、forの前にlt = zeros(n_end, 100);を入れればいいです。

MATLAB

1lt = zeros(n_end, 100); 2for n=2:n_end 3 for m=1:100; 4 rrt(n) = sm(n); 5 % 以下、略

ちなみにfor m=1:100;の最後のセミコロンは不要ですね。

投稿2019/12/11 16:59

Bearded-Ockham

総合スコア430

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.36%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問