• 締切済み

matlabで価格時系列

matlabで価格の時系列を発生させたいです。 現在はこんな感じです。 x1(1)=100 for t = 1:1000 x1(t+1) = x1(t) + randn; end これをfor文を使わないで書くことは可能でしょうか。 計算時間を短縮したいです。

みんなの回答

回答No.1

filter 関数使えばできます。 t=1:1000; delta = [0 randn(1,length(t)-1)]; x = filter(1,[1 -1],delta,100); でできないかな？ matlab手元に無いんで未確認です。

  1. カテゴリ
  2. [技術者向] コンピューター
  3. プログラミング・開発
  4. その他（プログラミング・開発）

関連するQ&A

  • MATLABの計算精度

    いま， 　A=[0.803, -0.803, 1.2124, 0.803, -0.803, -0.803, 1.2124, 0.803]; とします． このとき，次の２通りの演算をします． (i) a=0; for I=1:8 　if( A(I)＞＝0 ) 　　a=a+A(I)*0.25; 　else 　　a=a+A(I)*0.75; 　end end (ii) Bi=[0.25, 0.75, 0.25, 0.25, 0.75, 0.75, 0.25, 0.25]; b=sum(Bi.*A); (i)と(ii)によって得られる和の値をMATLABのコマンド上で表示すると， 同じ値として表示されるのですが，a-b　と計算すると，0 ではなく， 8.04e-15のようになります．（現環境ではMATLABが使えないので，確認できませんが，もっと大きいデータでやったとき場合には，上記のように誤差が生じました）． 困っているのは，現在組んでいるプログラムでは計算した値の正負によって 異なった計算をする必要があり，計算精度が非常に重要になってきます． そのため計算結果が例えば　+1.4e-15となるか　-6.4e-16　となるかによってその後の計算結果が大きく変わってしまいます． MATLABでは，上記(i)(ii)のような誤差が生じるためか，+1.4e-15となってほしいところが，なぜか　-6.4e-16となってしまいうまくプログラムが動きません． MATLABでは，for文によって和をとるのではなく，一回でまとめて計算しなければ誤差が生じるのでしょうか？ また，MATLABで計算精度を変更できないのでしょうか？ （formatコマンドによって，「format short g」「format double g」としましたが有効な小数点が増減するだけで結果は同じでした．） よりよい精度で計算するにはどうしたらよいのでしょうか？ 回答よろしくお願いします．

  • matlabのfor文について

    matlabについて基本的な質問があります． たとえば a1 = 11 a2 = 32 a3 = 53 a4 = 63 for i = 1:length(a1) ○○○ end for i = 1:length(a2) 　○○○ end for i = 1:length(a3) 　○○○ end for i = 1:length(a4) 　○○○ end の時に，for文を一つにまとめるにはどうすればいいですか？

  • MATLABでの行列の全要素の和

    MATLABでの行列の全要素の和を計算するのに下の例のようにfor文を使うのとsum関数を使うのとで結果が違います。何故ですか？ 計算結果の差をとってみると極々微小な差が生じてしまいます(e-10位)。 シミュレーションの中でこのような処理を行っているのでより正確な方が知りたいです。 例 A=rand(10); plus=0; sum=0; for i=1:10 　for j=1:10 　　plus=plus+A(i,j); 　end end sum=sum(A(:));

  • matlabでの長時間の計算について

    すみません、現在matlabで数値解析を行っております。 質問を以下にまとめたので、よろしければアドバイスをお願いします。 -------------------------------------------------------------- matlabで数十分かかると考えられる計算を行おうとしています。 その場合、計算後の大容量のデータは一般的にどのように保存すべきでしょうか？ これまで、パラメータを50個用意してそれぞれ計算を行い、 後に見比べるために、目的の変数の値をプロットし、画像(角度などを指定してjpegに)として保存していました。 （for文を使って一定量パラメータを変化させる毎に、saveasを使ってグラフを保存していた。） 今まではパラメータ１つあたり数十秒程度で済む計算だったのでこれでも良かったのですが、 今後はパラメータ１個あたりに対しかかる計算時間が数十分になります。。 そこで、今までと同じように保存していくやり方では、もしグラフの保存方法（角度など）を間違えてしまった場合に かなりの損失になってしまうと思ったので質問しました。。 何か効果的なデータの保存の仕方はありませんでしょうか？ 以後もパラメータをfor文を使って一定量変化させ、50個用意して計算していこうと思うのですが・・ 一応現状で考えているのは、saveを使ってworkspace内の変数（全て）をそれぞれ保存していくというものです。 何かその他に効果的なデータの保存方法などありましたら是非教えてください。 -------------------------------------------------------------- 上記のように、 計算時間のかかる計算をmatlabを使って行おうと考えているのですが、 計算機を走らせている間、動作が著しく重くなり、ほとんど他の作業が出来なくなります。 （これは避けられませんよね・・？） 私はmatlabを学校で使っているので、帰宅する際、計算機を走らせて帰ろうと考えています。 （節電家でなくてすみません。。） 問題は途中でメモリ不足などで計算がストップしてしまわないかということです・・ 一応、パラメータを変える毎に初期化すべき変数はzerosを使って初期化しているのですが、 パラメータ毎の目的の変数(6つ)は、多いとそれぞれ10^6×200ほどの配列となる可能性があります。 せっかく手間暇かけて走らせたのに、途中でメモリ不足などと表示されたら辛いなぁと思い・・ 現在使っている計算機はubuntu11.04で、メモリは992.8MiBなのですが、やはりこの計算機だと際どいでしょうか・・？ -------------------------------------------------------------- コードに関する質問なのですが・・ uが50000000*200の配列で、それぞれの要素に値が入っているとします。 これをそのまま1行目～50000000行目まで、順にグラフにプロットすると、かなりの時間がかかるので、 1行目から、一定数（例えば100）毎の行の配列を グラフにプロットさせたいのですが、やり方がわかりません・・ x = -10: 10: 200 ; t = 0: 100: 50000000 ; waterfall(x,t,u) ; とするとエラーが出ます。 どなたか訂正出来る方、お願いします。。 waterfall以外でも全然構いません。 -------------------------------------------------------------- 以上長々と質問すみませんでした。。 どなたかご指導頂けると助かります・・ あと、分類というのがよくわからなかったので、その他となっています・・ ごめんなさい。。

  • プログラムについて(matlab)

    matlabで以下のプログラムを実行しているのですが、 グラフが表示されません。 何故でしょうか、、。 ----------------------------------------- t=10 for r=400:1:2000 %波長 if 400<= r <700 y= 18*t.^(0.75); elseif 700<= r <1050 y= 18.*t.^(0.75).*10.^(0.002.*(r-700)); elseif 1050<= r <1150 y= 90*t.^(0.75)*1*1; elseif 1150<= r <1200 y= 90.*t.^(0.75).*1.*10.^(0.018.*(r-1150)); elseif 1200<= r <1400 y= 90*t.^(0.75)*1*8; elseif 1400<= r <1500 y= 5600*t.^(0.25); end plot(r,y) end ------------------------------------------------------ お詳しい方よろしくお願いします.

  • MATLABのプログラミングで

    MATLABでプログラミングを行っているのですが、forを使い for trt=1:rety(1); xt2(trt,:)=xt1(trt):xt1(trt)+(x-1); yt2(trt,:)=yt1(trt):yt1(trt)+(x-1); end for t3t=1:re/ysuu torie(:,:,t33(:,t3t))=saweq(xt2(gjir1(trt),:),yt2(gjir1(trt),:),:); end この様な変数を作製しています。 上のアルゴリズムは　a:b　(a=1:10) (b=a+5) で複数の位置データを記録する物で、 下のアルゴリズムは画像データを重ね合わせる物です。 これらのアルゴリズムでforを使わずにより高速化できる方法があれば連絡をお願いします。

  • MATLABのデータ抽出についておたずねしたいことがあります。

    MATLABのデータ抽出についておたずねしたいことがあります。 MATLAB初心者ですが、以下のようなプログラムを現在書いています。 f = input('file name? ---> ','s'); m = csvread(f,66,2,[66,2,1089,3]); time=(0:0.002:2.046); t=reshape(time,1024,1); current=m(:,1)*1000; voltage=m(:,2)*10; m2=[t current voltage]; a=m2(1,2); b=m2(2:40,2); if a>max(b); g=a; end for I=2:1024; for J=1:1024; for K=39+I; if K>1024; K=1024; end c=m2(I,2); d=m2(J:I-1,2); e=m2(I+1:K,2); if c>max(e)&&c>max(d); g=c; else end end end end plotyy(t,current,t,voltage) ylim([0,700]) ifの条件を満たしたときに、m2の2列目の値だけでなく条件を満たしたその行の3列全てのデータを抽出し、条件を満たした行だけの新たな行列を作りたいのですが、どのような命令を書けばよいでしょうか。 お手数をおかけしますが、ご教授よろしくお願いいたします。

  • MATLABのアップグレード問題

    MATLAB初心者です。 MATLAB6.1（おそらく）で開発されたコードがMATLAB7.01で動きません。 アドバイスを頂ければと思います。 1: out = {}; 2: for i=1:length(files), 3: out = [out, helper([d '/' files{i}], pat(2:end))]; 4: end 上記のコードの3行目で「Conversion to cell from char is not possible.」エラーが発生します。 MATLAB6.1から6.5のアップグレード問題で、空のセル配列の連結が原因のようです。 そこで下記のように修正しました。 1: out = {}; 2: for i=1:length(files), 3: if i==1 out = helper([d '/' files{i}], pat(2:end)); 4: else out = [out, helper([d '/' files{i}], pat(2:end))]; 5: end 6: end 1回目のループの場合は連結せず、2回目以降は連結するように修正しています。 このように修正しても意図通り動いていません。 修正方法は間違っていますか。

  • MATLABで外れ値をプロットする方法を教えてくだ

    　MATLABで外れ値をプロットする方法を教えてください。 　具体的には下記のとおりです。 　MATLABで、中央値＋中央値から±20％離れた時系列データをプロットするにはどうすればよいでしょうか？ 冗長なprogramしか作成できないため、よい関数か、方法があれば教えていただきたく思い、 質問しました。 まず、時系列データですが、数十万点数からなる信号になっています。 この信号の中央値をだし、その中央値から±20%の範囲でひっかかるデータをplotしたく思っています。 中央値算出は、関数を知っていますので処理が早いのですが中央値から±20%の範囲でひっかかるデータをプロットするのにいまは、以下のような考えで行っています。 行データ、size(ｘ（：））　ans＝：1000000　1とおきます。 まず、 １）中央値を求める。 ２）median+20%、median-20%を求める。 ３） hold on; 　for i=1:1:data-length 　　if（x(i)<=median-20% & x(i)>=median+20%） plot x(i) endif 　end というようにしています。しかし、これではデータ一点一点検討する必要があるので冗長で、処理に時間が かかります。 そのため、for文を使わずに、 medi-20(:)=function(x(:),median-20) medi+20(:)=function(x(:),median+20) で計算できるような、関数functionを教えていただけると幸いです。 よろしくお願いします。

  • matlabで複数条件のif文を行列計算で行うには？

    matlabで複数条件のif文を行列計算で行うには？ matlabではfor文などは行列計算で回したほうが高速化できると聞きました。 for文でのやりかたはわかったのですが、複数条件でのif文ではどのように書けばよいのでしょうか？ 自分のプログラムは以下のようなものです for i=1:100 for j=1:100 theta_kari=atan(abs(B(j,i))/abs(A(j,i))); if A(j,i)>=0&&B(j,i)>=0 theta(j,i)=theta_kari; elseif A(j,i)<0&&B(j,i)>=0 theta(j,i)=theta_kari+(pi/2); elseif A(j,i)<0&&B(j,i)<0 theta(j,i)=theta_kari+pi; elseif A(j,i)>=0&&B(j,i)<0 theta(j,i)=theta_kari+(pi/2*3); end end end よろしくお願いします

注目のQ&A

カテゴリ

専門家に質問してみよう

職業から探して質問する

あなたにピッタリな商品が見つかる！ OKWAVEセレクト

質問する