※ 引述《flockbody ()》之銘言：
: t=0:0.001:8;
: y=(square(2*pi*t/8)+1)/2;
: x=fft(y);
x1=abs(x); %將振幅大小計算出來
magx1=fftshift(x1); % 將頻率為零的資料移到頻譜圖中心
N = length(t); % 點數
fs = 1/0.001; % 取樣頻率
freqStep = fs/N; % 頻域的頻率的解析度
freq = freqStep*(-(N-1)/2:(N-1)/2); % 頻域的頻率刻度
plot(freq,magx1);
: plot(y)
用plot(y)畫出來並不是振福大小，且對應之頻率也不對
: 這是一個震幅為0~1頻率為1/8HZ的方波
: 作傅立葉轉換
: 出來的圖形 (y) 應該要是sinc的波形
: 可是我run出來的圖形卻很奇怪
: 請問一下各位大大
: 是我哪邊寫錯了嘛
: 謝謝
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.123.121.56