matlabで日付と時刻データをシリアル日付に変更したい

Question

前提・実現したいこと

１日～３０日までの一時間ごとのデータ、つまり３０×２４個の時系列データについて解析を行いたいと思い、datenumを用いて日付と時刻の二つの値からシリアル日付を生成することにしました。

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

エラー: datenum (line 188)
DATENUM は失敗しました。

原因:
エラー: datenum (line 167)
DATENUM への入力引数は、サイズが同じであるか、スカラーでなければなりません。

該当のソースコード

yr = repmat((0:23),24,1);　　%一時間ごとの時刻
mo = repmat((1:30)',1,30);　 %日付
time = datestr(datenum(yr(:),mo(:),1));

試したこと

サイズが同じでなければならない、というのは今回の場合は、時系列データの縦と横の数が等しくなければならないということでしょうか。つまりいまはyrは24個、moは30個要素がありますがこの値が等しくなければならない？

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

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

Answer 1

間に抜けがない1時間ごとの時刻、が欲しいだけであれば、日付と時刻の組み合わせを作る必要はありません。1日は24時間と決まっていますので。

time = datestr(datenum(2020,4,1)+(0:30*24-1)/24);

（先頭の日付は適当ですのでご自分で合わせてください）

---

あなたが元にしたソースコードでは、年×月の行列のまま timeseries にしているのではなく、ベクトルに並べ直して使っていますね？

ts = timeseries(y(:),time,'name','AirlinePassengers');

つまり、日付×時刻の行列に並べる必要はないということです。

Comments

[notes]

ご回答いただきありがとうございます。
時刻が欲しいというより、一か月分の一時間ごとに測定された時系列データを扱うために、行を日付、列を時刻として、その組み合わせを作成しようとしています。
示して頂いたコードは、4月1日から一か月分、一時間ごとの時刻をシリアル日付に変換している、という理解でよろしいでしょうか?私は作成したシリアル日付：timeを、時系列データ値：tsにtimeseriesを用いて対応づけようと考えています。(参考サイト:https://jp.mathworks.com/help/stats/examples/time-series-regression-of-airline-passenger-data.html)
この場合、示して頂いたコードでも対応付けて時系列データオブジェクトを作成することは可能なのでしょうか。

[yuki23]

元にしたソースコードがあるなら先に示していただきたかったです

[yuki23]

> 時刻が欲しいというより、一か月分の一時間ごとに測定された時系列データを扱うために、行を日付、列を時刻として、その組み合わせを作成しようとしています。
手段と目的を混同されていますね。目的は1ヶ月分の1時間毎の time を作ることで、「行を日付、列を時刻として、その組み合わせを作成」は手段に過ぎないはずです。
そのような遠回りな手段を取らなくても目的は達成できるはずです。

> この場合、示して頂いたコードでも対応付けて時系列データオブジェクトを作成することは可能なのでしょうか。
聞く前にご自分で試してみていただけないでしょうか

[yuki23]

ちなみに、参考サイトがそのような遠回りな手段をとっているのは、月によって日数が異なるからだと思われます

[notes]

お返事いただきありがとうございます。
自分の必要とする答えを求めるあまり、提示すべき情報や考えるべきことを放棄してしまっていました。
今後このようなことがないよう、十分に調査して考えてから、必要な情報を載せることを心がけます。
matlabを全く触っていなかったため、参考にしたサイトと同じようなコードを書くことができれば自分の目的を達成することにつながると勘違いしていました。
yuki23さん、およびBearded-Ockhamさんからいただいた回答をもとにもう一度自分で考えてみようと思います。そのうえで、どうしても分からないようなことが出てきた場合はコメントにてお聞きする形を取らせていただければと考えています。よろしくお願いいたします。

Answer 2

まず、datenumに引数を3個渡した場合は、datenum(年,月,日)と解釈されるので時刻部分は生成されません (そもそもyrはyear、moはmonthの略です)。時刻まで生成したいときは、datenum(年,月,日,時,分,秒)の6個の引数をすべて指定する必要があります。

次に一度に複数のシリアル日付値を作りたいときは、すべての日付・時刻の組み合わせを指定しなければいけません。例えばdatenum(2020, 6, 1:2, 0:1, 0, 0)とすれば2020/6/1 00:00:00 と2020/6/2 01:00:00のふたつの値が返ってきます。勝手にすべての組み合わせを求めてくれるわけでは、ありません。なので、日付の引数は[1が24個 2が24個 ... 30が24個]、時の引数は[0:23が30個]のベクトルを渡す必要があります。

---

引数に列ベクトルと行ベクトルを混ぜると変な結果が出ますが、バグでしょうね。例えば月に5:6の行ベクトル、日に(11:12)'の列ベクトルを渡すと

>> datenum(2020, 5:6, (11:12)')
ans =
      737922      737922
      737922      737922

つまり、2020/5/11が2x2に並んだ行列が返ってきます。これは明らかに期待する動作ではありません。

Comments

[notes]

ご回答いただきありがとうございます。
引数にyr、moを入れれば勝手に全ての組み合わせを作成してくれると思っていましたが違うんですね。
全ての組み合わせを引数とする書き方は、日付は[(1:30)＊24]、時は[(0:23)＊30]であっていますでしょうか。
また、私はyr、moを別に定義しているのですが、これを引数として用いる場合は、日付の引数は[mo(:)＊24]、時の引数は[yr(:)＊30]とする、という理解でよろしいでしょうか。
引数に行と列ベクトルを用いると結果にバグが生じるとおっしゃっていますが、この結果を用いて計算や作図を行なっても正常な結果が得られるのでしょうか？