最初に、Matlab FFT に慣れていました。たとえば、次のようにします。
clc
fs = 4;
t = 0:1/fs:1-1/fs;
x1 = sin((2*pi)*t);
y1 = fft(x1,n);
FT_power1 = abs(y1(1:floor(n/2))).^2;
FT_phase1= (((angle(y1(1:floor(n/2))))) * 180/pi);
[c1,i1] = max(FT_power1);
phase(1) = FT_phase1(i1);
そして、ナイキストとサンプリング周波数パラメータを使用する理由については十分に理解していると思います。ただし、出力の配列に対して同じ手順を実行することを想定しています (つまり、私が持っている別のサブルーチンからの正弦波です)。これは、出力配列に 2 つの列があることを意味します。最初の列は時間ベクトル用で、2 番目の列は実際の正弦応答。では、fs パラメータを組み込むにはどうすればよいでしょうか。私は突然混乱しています。詳細が必要な場合はお知らせください。ありがとうございました。