MATLABでカーブフィットしたデータをテキスト形式で保存する方法を教えてください| OKWAVE

MATLABでカーブフィットしたデータを、関数「xlswrite」「save」を用いて、xls形式・txt形式で保存する方法を教えてください

下記は、周波数応答関数を入力し、カーブフィットを行うプログラムです。カーブフィットのデータをエクセル形式、テキスト形式で保存をする場合、関数[save]または[xlswrite]を使用し，どのように書けばいいのでしょうか。どうか何卒ご教授お願い致します。 ---------------------- clear variables global;comgui('close all') fm='wagawara_ryou10_'; n1=101; n2=142; n=0; for i=n1:n2 fnm=strcat(fm,int2str(i),'_h.txt') [IIw,Mag,Phi,Q,Q]=textread(fnm,'%f %f %f %f %f','headerlines',3); n=n+1; IIxf(:,n)=complex(Mag.*cos(Phi*pi/180),Mag.*sin(Phi*pi/180)); end %---初期固有振動数データと減衰比データの読み込み．2005.9.13--- fn='wagawara_ryou10_mif_pk_dm.txt' [peak,zeta]=textread(fn,'%f %f','headerlines',2); pkze=[peak,zeta]; %---処理開始--- iiplot %伝達関数のプロット XF(5).po=pkze %初期固有振動数と減衰比を代入。 idcom('e 15 371'); %idcom('e i w') i:帯域幅，w:中心周波数 idcom('est'); idcom('eup .05 .002 -10'); %idcom('eup dstep fstep num i') out=id_rm(XF(5),[1 1 1 1]); XF(3)=res2xf(out,IIw); iicom('IIxhOn'); [som2,ga2,pbs2,cps2]=res2nor(IIres,IIpo,IDopt); ※注意 'wagawara_ryou10_'←入力データ strcat←組み込み関数 textread←組み込み関数 complex←組み込み関数 'wagawara_ryou10_mif_pk_dm.txt'←入力データ textread←組み込み関数 idcom←組み込み関数 id_rm←組み込み関数 res2xf←組み込み関数 iicom←組み込み関数 res2nor←組み込み関数 ---------------------- どうぞよろしくお願いします。

カーブフィットを行うプログラムについて

カーブフィットを行うプログラムについて質問があります。下記は、周波数応答関数を入力し、カーブフィットを行うプログラムですが、プログラムの下４行の意味が全く分かりません。（下４行とは、■マークが先頭に書いてあるものです）どうか何卒ご教授お願い致します。わかりにくい質問の仕方をしてしまい、申し訳ないです。 ---------------------- clear variables global;comgui('close all') fm='wagawara_ryou10_'; n1=101; n2=142; n=0; for i=n1:n2 fnm=strcat(fm,int2str(i),'_h.txt') [IIw,Mag,Phi,Q,Q]=textread(fnm,'%f %f %f %f %f','headerlines',3); n=n+1; IIxf(:,n)=complex(Mag.*cos(Phi*pi/180),Mag.*sin(Phi*pi/180)); end %---初期固有振動数データと減衰比データの読み込み．2005.9.13--- fn='wagawara_ryou10_mif_pk_dm.txt' [peak,zeta]=textread(fn,'%f %f','headerlines',2); pkze=[peak,zeta]; %---処理開始--- iiplot %伝達関数のプロット XF(5).po=pkze %初期固有振動数と減衰比を代入。 idcom('e 15 371'); %idcom('e i w') i:帯域幅，w:中心周波数 idcom('est'); idcom('eup .05 .002 -10'); %idcom('eup dstep fstep num i') ■out=id_rm(XF(5),[1 1 1 1]); ■XF(3)=res2xf(out,IIw); ■iicom('IIxhOn'); ■[som2,ga2,pbs2,cps2]=res2nor(IIres,IIpo,IDopt); ---------------------- よろしくお願い致します。失礼します。

周波数応答関数の計算

打撃実験を行って解析するのですが下記のプログラムが分かりません。 [IIw,Mag,Phi,Q,Q]=textread(fnm,'%f %f %f %f %f','headerlines',3); IIxf(:,n)=complex(Mag.*cos(Phi*pi/180),Mag.*sin(Phi*pi/180)); 読み込んだファイルの内容が Np Nc / Freq Hm Ha Coh XPwr 6401 5 0 11.8941 0 0.938239 1.78076E-006 1 27.7519 130.402 0.915975 2.26936E-008 2 3.95322 -38.8943 0.973454 1.64854E-008 . . 一列目が周波数と思います。２列と３列が何か分かりません。また、IIxfは周波数応答関数だと思うのですがどの公式を用いているのでしょうか?

テキストへの保存方法

N個(毎回変動します)のテキストファイルがあります。その中に360行２列のデータがあります。このN個のテキストを1つのテキストファイルにまとめて、 360行(２×N)列のデータにして保存したいです。 '------------------保存するプログラム------------------------- SaveName = Dir2.Path & "\" & Nam & ".txt" Open SaveName For Output As #1 For i Write #1, V1(1,i),VRMS(1,i),V1(2,i)・・・V1(N,i),VRMS(N,i)　 next i Close #1 '-------------------------------------- 上のプログラムの４行目 Write #1, V1(1,i),VRMS(1,i),V1(2,i)・・・V1(N,i),VRMS(N,i)　ですが毎回変動するため何列になるかわかりません。行の場合でしたらFor NEXTを使えばよいと思うのですが、横に増えてい場合どのようにすればよいのでしょうか？ VB6で作成しています。アドバイスよろしくおねがいします。

エクセルのデータを大量のテキストファイルに書き込みたい。

エクセルの縦のセルのデータ(下のCells( 4 * j + i, 4))を繰り返しで、400個のテキストファイルに書き込みたいのですが、下の""内のjは繰り返しせずに、j.txtとなって入力されてしまいます。なにかいい方法がありませんでしょうか。ご協力お願いします。（ちなみにvbaは初心者です。） Private Sub CommandButton2_Click() Dim n As Long n = FreeFile Dim j As Integer For j = 0 To 399 For i = 1 To 4 a = Cells( 4 * j + i, 4) f2 = ActiveWorkbook.Path & "\j.txt" Open f2 For Append As #n Print #n, a Close #n Next i Next j End Sub 例. 1 2 3 4 5 6 7 8 9 10 12 13 1,2,3,4→1.txt 5,6,7,8→2.txt 9,10,11,12→3.txt

保存データの￥nを に変換

http://okwave.jp/qa/q7387963.html 以前こちらの質問をさせていただきまして、 htmlspecialcharsはHTMLとしてクライアントに返すときに使用することがわかりました。しかし、ここで問題が起こりました。フォームから入力されたデータをdata.txtに一行一データとして保存していくプログラムなのですが、以前のやり方は、保存時にhtmlspecialchars変換してしまう方法で以下のようにやっていました。 $honbun=htmlspecialchars($_POST["honbun"],ENT_QUOTES); $honbun=str_replace("\n"," ",$honbun); $file=fopen("data.txt",a"); fputs($file,$honbun."\n"); //行末に\nをつけて保存 fclose($line); こうすることで、data.txtには --------------------------------------------- こんにちは、山田です。 元気ですか？\n こんばんは、鈴木です。 元気です。\n --------------------------------------------- このように行の最後以外は改行を にできていました。しかし、最初に書いた通り、 htmlspecialcharsは保存時ではなく、出力時に使う。ということでしたのでこのやり方は正しくないのでこんな感じに書き直してみました。 $honbun=$_POST["honbun"];　//まだこの時点ではhtmlspecialcharsしない $honbun=str_replace("\n"," ",$honbun);; $file=fopen("data.txt",a"); fputs($file,$honbun."\n"); //行末に\nをつけて保存 fclose($line); そして、htmlとして出力時に $f=file("data.txt"); foreach($f as $key=>$value){ print(htmlspecialchars($value));　//ここでhtmlspecialcharsするよう変更 print(" "); 　 } こうやってみました。が、しかし、ここで問題が発生。 html出力結果が、こんにちは、山田です。 元気ですか？こんばんは、鈴木です。 元気です。このように、<br　/>で改行して表示したいにもかかわらず、 を文字列として表示してしまいました。いったいどうやったら解決できるでしょうか？やはり、以前のやりかたに戻すべきでしょうか？わかる方、ご回答よろしくお願い致します。

fortran90

fortran90 現在 fortran90 で計算プログラムを作っております。プログラムはかけたのですが、実行画面にエラーが出ます。解決方法が分かりません。どなたか分かる方いらっしゃいますか。プログラム implicit real*8 (a-h,o-z) parameter(n=900) open (1,file='004.txt',status='old') do 10 i=1,n read (1,*,end=11)xx,yy,zz,ww theta=acos(zz) phi1=acos(xx/sqrt(1-zz*zz)) phi2=asin(yy/sqrt(1-zz*zz)) write(*,*)theta,phi1,phi2 10 continue 11 end プログラム内容あるフォルダーに保存してあるテキストファイルから４×４行列のデータを読み込み、そのデータからθの値とφ１、φ２のアタを求めるプログラムです。実行分のエラー画面には、 acos: DOMAIN error です。

教えてください

a(複素平面)を正則関数ｆの孤立特異点とする。この時aがｎ位の零点ならばaはf '/fの一位の極で、Res{(f '/f),a}=n をしめせ

数学の問題がわかりません＞＜

数学の問題がわかりません！よければ教えてください＞＜ I = [0,∞)とおく。 f , f_n ; I → R はI上で有界な関数とし、関数列{f_n}[n=1,∞]は関数 f に I 上で一様収束するとする。 (1) 各n∈Nに対してlim[x→∞]f_n(x) = a_n ∈Rが成り立つならば、数列{a_n}[n=1,∞]はCauchy列であることを示せ。 (2) (1)と同じ条件の下でlim[n→∞]a_n = A とおくとき、lim[x→∞]f(x) = Aであることを示せ。回答よろしくお願いします！

phpでのソートについて

phpのソートについて教えてください。以下のようなカンマ区切りのログファイルlog.txtがあるとしまして、 100,200,a,b,c, 300,100,d,e,f, 500,60,g,h,i, 50,300,j,k,l, 1000,60,m,n,o, このデータから１番目のデータ（数値）が２番目のデータ（数値）より大きいデーターのみを対象にして ※つまりは以下のみ対象 300,100,d,e,f, 500,60,g,h,i, 1000,60,m,n,o, ここから１番目のデータ（数値）から２番目のデータ（数値）を引いた数が大きい順にソートしてファイルに保存させたいのですが、 200,d,e,f,（300-100なので200） 440,g,h,i,（500-60なので440） 940,m,n,o,（1000-60なので940） ↓ ※最終的にこの順番で新しいログファイルnewlog.txtへ保存させたい 940,m,n,o, 440,g,h,i, 200,d,e,f, これらの処理を効率よく１度で行う方法はございますでしょうか。わかりにくい説明で申し訳ありません。お忙しい中恐縮ですがご教授いただけましたら幸いです。何卒宜しくお願い致します。

Perl5で同時刻のデータを統合したい

perl初心者です。下記のような同時刻の2つのデータを1つのデータに統合させたいのですがうまく出来ず困っています。どうかお知恵を貸していただけないでしょうか。データは時刻(時:分:秒), 値1, 値2になっています。 test1.txtの同時刻の後ろにtest2.txtの値1と値2を入れ、欠測値には-999を入れるプログラムを作っています。厄介なのは、開始時刻がtest1.txtよりtest2.txtが早い場合や終了時刻がtest2.txtよりtest1.txtが遅い場合がある事です。 test1.txt 10:13:14, 3.1, 0.1 10:13:15, 6.1, 0.3 10:13:16, 8.7, 0.2 10:13:17, 12.8, 0.3 10:13:18, 13.4, 0.5 10:13:19, 15.2, 0.4 test2.txt 10:13:16, 32.5, 0.01 10:13:17, 33.1, 0.03 10:13:18, 36.2, 0.02 10:13:19, 34.3, 0.01 10:13:20, 33.8, 0.04 10:13:21, 32.6, 0.09 10:13:22, 32.1, 0.08 希望結果 test.txt 10:13:14, 3.1, 0.1, -999.0, -999.00 10:13:15, 6.1, 0.3, -999.0, -999.00 10:13:16, 8.7, 0.2, 32.5, 0.01 10:13:17, 12.8, 0.3, 33.1, 0.03 10:13:18, 13.4, 0.5, 36.2, 0.02 10:13:19, 15.2, 0.4, 34.3, 0.01 10:13:20, -999.0, -999.0, 33.8, 0.04 10:13:21, -999.0, -999.0, 32.6, 0.09 10:13:22, -999.0, -999.0, 32.1, 0.08 以下が自分が作成したプログラムです。 open IN_1, "test1.txt"; open IN_2, "test2.txt"; open OUT, ">test.txt"; while ($input1 = <IN_1>) { # test1.txtの処理 chomp $input1; # 改行削除 @input1 = split(/,/, $input1); # カンマ区切り $n = $n + 1; $time1[$n] = @input1[0]; @time1 = split(/:/, $time1[$n]); # 時刻をコロン区切り @hours1[$n] = @time_l[0]; @min1[$n] = @time_l[1]; @sec1[$n] = @time_l[2]; @a[$n] = @input1[1]; # 値1 @b[$n] = @input1[2]; # 値2 } while ($input2 = <IN_2>) { # test2.txtの処理 chomp $input2; @input2 = split(/,/, $input2); $m = $m + 1; $time2[$m] = @input2[0]; @time2 = split(/:/, $time2[$m]); @hours2[$m] = @time2[0]; @min2[$m] = @time2[1]; @sec2[$m] = @time2[2]; @c[$m] = @input2[1]; # 値1 @d[$m] = @input2[2]; # 値2 } # 開始と終了時刻の計算 # test1.txtの時刻 $hours1_S = @hours1[1]; # 開始時 $hours1_E = @hours1[$n]; # 終了時 $min1_S = @min1[1]; # 開始分 $min1_E = @min1[$n]; # 終了分 $sec1_S = @sec1[1]; # 開始秒 $sec1_E = @sec1[$n]; # 終了秒 $time1_S = $hours1_S*3600 + $min1_S*60 + $sec1_S; # 開始時刻を秒に計算 $time1_E = $hours1_E*3600 + $min1_E*60 + $sec1_E; # 終了時刻を秒に計算 # test2.txtの時刻 $hours2_S = @hours2[1]; $hours2_E = @hours2[$m]; $min2_S = @min2[1]; $min2_E = @min2[$m]; $sec2_S = @sec2[1]; $sec2_E = @sec2[$m]; $time2_S = $hours2_S*3600 + $min2_S*60 + $sec2_S; $time2_E = $hours2_E*3600 + $min2_E*60 + $sec2_E; if($time1_S <= $time2_S){ # 開始時刻の比較 $starttime = $time1_S; }else{ $starttime = $time2_S; } if($time1_E <= $time2_E){ # 終了時刻の比較 $endtime = $time2_E; }else{ $endtime = $time1_E; } $j = $endtime - $starttime; # 全体のデータ個数 # test.txtへ出力 for($i=1; $i<=$j; $i++){ if(@hours1[$i] != @hours2[$i] && @min1[$i] != @min2[$i] && @sec1[$i] != @sec2[$i]){ printf OUT "%2d:%2d:%2d %5.1f% 5.1f %5.1f %7.2f\n", @hours1[$i], @min1[$i], @sec1[$i], @a[$i], @b[$i], -999, -999; }elsif(@hours1[$i] == @hours2[$i] && @min1[$i] == @min2[$i] && @sec1[$i] == @sec2[$i]){ printf OUT "%2d:%2d:%2d %5.1f% 5.1f %5.1f %7.2f\n", @hours1[$i], @min1[$i], @sec1[$i], @a[$i], @b[$i], @c[$i], @d[$i]; }else{ printf OUT "%2d:%2d:%2d %5.1f% 5.1f %5.1f %7.2f\n", @hours2[$i], @min2[$i], @sec2[$i], -999@, -999, @c[$i], @d[$i]; } } close IN_1; close IN_1; close OUT; とても汚いプログラムになってしまいました… このプログラム以外でも構いませんのでどうかよろしくお願い致します。

複素解析の極と留数を求める問題

複素解析の極と留数を求める問題 f(z)=z/sinh(z)の留数はなぜ(-1)^n nπiなのですか？ f(z)の極はz=nπiだということは理解できるのですが、その後の留数を求めるときの式変形ができません。 Res[f(nπi)]=lim[z→nπi](z-nπi)*z/sinh(z) これ以降の式変形、導出ができません。

テキストデータExcel取込時の文字化け、その他

VBA初心者です。やりたいことは以下の通りです。カンマ区切りの文字が記載された.txtデータをExcelファイルに取り込みたいと思っています。 (1)自動でフォルダを開いて手動で任意のテキストデータを選ぶ (2)データを文字化けなくカンマ区切りでセルごとにExcelに表示 (3)元データの.txtはファイルによって空白行が1行のときもあれば2行のときもありまちまちなので自動で空白行をスキップして呼び込みたい (4)呼び込んだ.txtの内容を元に読み込先のExcelに関数（I列に来るべき.txtのデータを参照してJ列「収支」K列「勝率」が出る）が仕込んであるが、ファイルによってデータの行数が違うため最終行を取得したい特定の.txtファイルではなくこちらで選べるようにしたいので以下に示したソースの通りGetOpenFilenameで自動でファイルを開いて手動で選べるようにしています。その際、一部別の漢字に変換されたり、カンマがひとつだけ・に代わっていたり文字化けがあります。呼び込みたい.txtは空白行が1行あったり、2行あったりするのでそこをスキップして純粋に文字のある行から取り込みたい。読み込み先のExcelには事前に項目を作っているため、空白行をスキップできれば自動的に項目の下にデータが表示されるようになっています。※予め空白のないテキストを呼び込んでテストすると（文字化けは別として）上手くいきました。 (4)は難しくてもせめて(1)～(3)までは自力でできないかと調べていますがうまくいきません。UTF-8をANTI形式に保存しても直ったり直らなかったり、またVBAのADODB.Streamオブジェクト（Microsoft ActiveX Data Objects x.x Library）を試そうとするもファイルパスやターゲットで特定のファイルを指定する部分があり、それを指定しないで使うにはどうすればいいのかがわかりません。ご査収いただけましたら嬉しいです。下記はネットのソースを一部直して作成。これを利用してできないでしょうか。 Option Explicit '=================================================================================================== Private Const g_cnsTitle As String = "テキストファイル読み込み" Private Const g_cnsFilter As String = "全てのファイル (*.*),*.*" Sub Txt読み込み() Dim varFileName As Variant Dim intFree As Integer Dim strRec As String Dim strSplit() As String Dim i As Long, j As Long varFileName = Application.GetOpenFilename(FileFilter:="txtファイル(*.txt),*.txt", _ Title:="txtファイルの選択") If varFileName = False Then Exit Sub End If intFree = FreeFile '空番号を取得 Open varFileName For Input As #intFree 'txtファイルをオープン i = 0 Do Until EOF(intFree) Line Input #intFree, strRec '1行読み込み i = i + 1 strSplit = Split(strRec, ",") 'カンマ区切りで配列へ For j = 0 To UBound(strSplit) Cells(i, j + 1) = strSplit(j) Next '配列をそのまま入れる方法も、ただし全て文字列として入力される 'Range(Cells(i, 1), Cells(i, UBound(strSplit) + 1)) = strSplit Loop Close #intFree End Sub

留数定理による実定積分の計算について

留数定理による実定積分の計算について現在複素積分について勉強中のものです。 ∫^{+∞}_{-∞}f(x)exp{itx}dxという形の積分の計算なのですが t>0については ∫^{+∞}_{-∞}f(x)exp{itx}dx=2πiΣ^{m}_{k=1}Res{f(z)exp{itz}} となりf(z)が偶関数のとき ∫^{+∞}_{0}f(x)cos(tx)dx=πiΣ^{m}_{k=1}Res{f(z)exp{itz}} となりf(z)が奇関数のとき ∫^{+∞}_{0}f(x)sin(tx)dx=πΣ^{m}_{k=1}Res{f(z)exp{itz}} となりますが t<0のときはどうなるのでしょうか。マイナスになるだけでしょうか。よろしくお願いします。

再帰関数とユークリッドの互除法を用いて最大公約数を求める

学校の課題で、再帰関数とユークリッドの互除法を用いて１０～１０００００までの最大公約数を求めるという問題が出て自分でプログラムを作ってみたのですが無限ループに入ったり、関数が使えてないみたいでできません。プログラムを見て頂いて、どこを改善したらいいかを教えてください。 #include <stdio.h> /* 正整数 n, m の最大公約数を計算する */ int gcd(int n, int m) { int res; res = n % m; if (res == 0) return m; /* 最大公約数が求まった */ return gcd(m, res); /* 再帰呼び出し */ } int main(void) { int i,j; for(i=10000;i>=10;i--){ for(j=10;j=10000;j++){ printf("%d to %d no saidaikouyakusuu ha %d \n", i, j, gcd(i, j)); } } return (0); } です。期限が今日の夜までで、ぎりぎりなんですがよろしくお願いします。

配列の受け渡し

全く原因が分からないので質問させていただきます．次のようなプログラムがあります．メイン関数で配列を宣言して，これを関数funcに渡して，func内で処理を行う．その後，メイン関数でprintf()を使って配列の内容を表示させる．私の書いたソースは次のようになっています． int main(void) { double data[]={1.0, 2.0, 3.0, 4.0}; for(i=0; i<4; i++) { printf("%f\n", data[i]); } func(data, 4); for(i=0; i<4; i++) { printf("%f\n", data[i]); } } void func(double *f, int n) { 配列に処理を行う } メイン関数の中で，一つ目のprintf()はちゃんと配列dataの初期値を表示します．ところが，二つ目のprintf()の手前で，「例外」とエラーがでてきます．（コンパイルは通っています．）これが何故なのか，さっぱり分かりません．関数func()の最後にprintf()を書いてみたところ，ちゃんと処理を行った後の正しい値が表示されます．なので，関数func()自体にバグはないように思うのです．原因が分かる方，いらっしゃいましたら，どうぞよろしくお願いします．

留数定理による実定積分の計算について

留数定理による実定積分の計算について現在複素積分について勉強中のものです。 ∫^{+∞}_{-∞}f(x)exp{itx}dxという形の積分の計算なのですが tを実数とし,kはΣの添え字,mは極の個数,iは虚数とします. このときtがt<0のとき ∫^{+∞}_{-∞}f(x)exp{itx}dx=-2πiΣ^{m}_{k=1}Res{f(z)exp{itz}} となりf(z)が偶関数のとき ∫^{+∞}_{0}f(x)cos(tx)dx=πiΣ^{m}_{k=1}Res{f(z)exp{itz}} となりf(z)が奇関数のとき ∫^{+∞}_{0}f(x)sin(tx)dx=-πΣ^{m}_{k=1}Res{f(z)exp{itz}} となるこれで合っていますでしょうか？よろしくお願いします。

関数左辺値参照返り値　

こんにちは。宜しくお願い致します。 > たとえば，下のサンプルプログラムでは関数fが参照を返すようになっているため，この関数fを代入の左辺においた，f(a, n) = 10;のような式が許される． > > #include <iostream> > using namespace std; > > int& f(int* a, int n) > { (2)～ return a[n]; > } > > int main(void) > { > const int n = 5; > int arr[n] = { 1, 2, 3, 4, 5 }; (1)～ f(arr, 3) = 10; > for (int i = 0; i < n; ++i) > cout << arr[i] << endl; > return 0; > } 上のサンプルプログラムで左辺値に参照を返す関数が(1)～あるのですが動作が分りません。関数内部の(2)～return a[n];というところで、関数の返り値int&の参照はa[n]という変数になるので、(1)～の結果としてa[n]に10が代入されるのでしょうか？教えてください。

複数ファイルの出力

質問があります。1600個のデータを読み込んで（16個ずつの組で100組）、1番目、3番目、5番目、7番目のデータをss[1].txtとして出力し、同様に2番目、4番目、6番目、8番目のデータをsk[1].txt、9番目、11番目、13番目、15番目のデータをks[1].txt、10番目、12番目、14番目、16番目のデータをkk[1].txtとして出力させます。以下、ss[2].txt,sk[2].txt,ks[2].txt,kk[2].txt・・・に関してはそれぞれ直前のファイルに16を足したデータ（例えばss[2].txtなら、17番目、19番目、21番目、23番目のデータを出力）を出力させるというプログラムです。以下に自分で作ったプログラムを掲載します。 void ecg_rr(fp,data_max) FILE *fp; { int i,j=0; int c=1; char fname[64]; data[0][i]=trend_data[0][i]; for(i,j=0;i<100,j<100;i++,j++) { sprintf(fname,"ss[%d].txt",1+i); fp = fopen(fname,"w"); fprintf(fp,"%4d\n",c); fprintf(fp,"%8.8f\n",data[0][0+j*16]); fprintf(fp,"%8.8f\n",data[0][2+j*16]); fprintf(fp,"%8.8f\n",data[0][4+j*16]); fprintf(fp,"%8.8f\n",data[0][6+j*16]); sprintf(fname,"sk[%d].txt",1+i); fp = fopen(fname,"w"); fprintf(fp,"%4d\n",c); fprintf(fp,"%8.8f\n",data[0][1+j*16]); fprintf(fp,"%8.8f\n",data[0][3+j*16]); fprintf(fp,"%8.8f\n",data[0][5+j*16]); fprintf(fp,"%8.8f\n",data[0][9+j*16]); fclose(fp); /*ks,kkに関しても同様*/ } } これだと100個ずつファイルを出力できるのですが、データが全て0として出力されてしまいます。

NHPP尤度関数の導出過程を教えて頂きたい

独立同一確率分布に従う n 個の確率変数 X_i (i=1,2,...,n) はある確率過程の到着時間間隔を表すとし，f (x_i) を X_i の密度関数とします．この場合，尤度関数はX_i の同時確率になるので， L(x_1,...,x_n)=Π_{i=1}^{n} f(x_i) となるのはよく本などで見ますが，NHPPの到着時刻を表す確率変数 S_k = Σ_{i=1}^{k} X_i を考える場合，尤度関数は L(s_1,...,s_n)=exp{-Λ(S_n)} Π_{k=1}^{n} λ(s_k) となるのはなぜでしょうか．Λ(・)とλ(・)はそれぞれNHPPの平均値関数と強度関数です．なんらかのNHPPの性質または f(・)とλ(・)の関係を使って導出したのでしょうか．詳しい導出過程を教えて頂ければとても助かります．よろしくお願い致します．

MATLABでカーブフィットしたデータをテキスト形式で保存する方法を教えてください

みんなの回答

お礼 2007/02/17 08:59

お礼 2007/02/17 08:59

関連するQ&A

MATLABでカーブフィットしたデータを、関数「xlswrite」「save」を用いて、xls形式・txt形式で保存する方法を教えてください

カーブフィットを行うプログラムについて

周波数応答関数の計算

テキストへの保存方法

エクセルのデータを大量のテキストファイルに書き込みたい。

保存データの￥nを<br />に変換

fortran90

教えてください

数学の問題がわかりません＞＜

phpでのソートについて

Perl5で同時刻のデータを統合したい

複素解析の極と留数を求める問題

テキストデータExcel取込時の文字化け、その他

留数定理による実定積分の計算について

再帰関数とユークリッドの互除法を用いて最大公約数を求める

配列の受け渡し

留数定理による実定積分の計算について

関数左辺値参照返り値

複数ファイルの出力

NHPP尤度関数の導出過程を教えて頂きたい

注目のQ&A

カテゴリ

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

専門家に質問してみよう

MATLABでカーブフィットしたデータをテキスト形式で保存する方法を教えてください

みんなの回答

お礼 2007/02/17 08:59

お礼 2007/02/17 08:59

関連するQ&A

MATLABでカーブフィットしたデータを、関数「xlswrite」「save」を用いて、xls形式・txt形式で保存する方法を教えてください

カーブフィットを行うプログラムについて

周波数応答関数の計算

テキストへの保存方法

エクセルのデータを大量のテキストファイルに書き込みたい。

保存データの￥nを<br />に変換

fortran90

教えてください

数学の問題がわかりません＞＜

phpでのソートについて

Perl5で同時刻のデータを統合したい

複素解析の極と留数を求める問題

テキストデータExcel取込時の文字化け、その他

留数定理による実定積分の計算について

再帰関数とユークリッドの互除法を用いて最大公約数を求める

配列の受け渡し

留数定理による実定積分の計算について

関数 左辺値 参照 返り値

複数ファイルの出力

NHPP尤度関数の導出過程を教えて頂きたい

注目のQ&A

カテゴリ

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

専門家に質問してみよう

関数左辺値参照返り値　