FFT

2017-01-31 12 views
0

と私はオクターブの練習cos(x)にプロットし、私はプロットFFT

enter image description here

ことによってそれを補間する必要がオクターブ/ Matlabの、プロットのcos(x)の中に近似のcos(x)の

fplot("[cos(x)]", [0, 2*pi]) 

N私は

x = zeros(n,1); 
for i=1:n 
    x(i,1)= (-1) + (i-1/2)*(2/n); 
end 
によって計算 [0, 2*pi]に等距離支持点であります

用語をプロットして近似するにはどうすればよいですか?

+0

期間1と2は、あなたのポストでも同じです。 – Adriaan

+0

ここに質問がありません.... – Andy

答えて

0

実際に質問が何かについても混乱しています。私はあなたがその式をプロットしたいと思います。

行ベクトルk = 0:(n/2 - 1)nを偶数とします)があります。

次に、同じ長さの行ベクトルとして係数をdする必要があります。

を次にベクトルx(列ベクトルを転置によって行ベクトルから作られ、例えばx = x.')あなたの列を定義

は、関数の

左項がある (あなたががそれらを得るところから、私は知りません)その後、

leftterm = sum(d .* exp(i * x * k), 2) 

右用語:

rightterm = sum(fliplr(d) .* exp(- i * x * (k.+1)), 2) 

、彼らが与えるAlltogether:

f = sum(d .* exp(i * x * k) + fliplr(d) .* exp(- i * x * (k .+ 1)), 2) 

とあなたがそれをプロットします

plot(x,f)