%clear?all
close?all
%load?sin.mat?%使用自己的數據

n1=1:21;%取一個周期

N=length（n1）;
t=ScopeData（n11）;
data=ScopeData（n12）;
%t=sinData（n11）;
%data=sinData（n12）;
figureplot（tdata）;%原數據圖
title（‘data‘）;
%DFT

%%
%基函數
W=exp（-j*2*pi/N）;

%%
%DFT變換
X=zeros（1N）;
for?k=1:N
????for?i=1:N
?????????X（k）=X（k）+data（i）*W^（（i-1）*（k-1））;
????end
end
figure（2）plot（n1-1X）;%頻譜
title（‘頻譜‘）
%%
%DFT反變換
x=zeros（1N）;
for?i=1:N
????for?k=1:N
????????x（i）=x（i）+X（k）*W^（-（k-1）*（i-1））;
????end
end
x=x/N;
figure（3）plot（tx）;%恢復數據
title（‘recover?data‘）;
%result=int（X1‘w‘-pipi）/（2*pi）;
%?w=-pi;
%?result1=eval（result）;
%?w=pi;
%?result2=eval（result）;
%?（result2-result1）/（2*pi）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????201??2020-02-17?17:47??code\ReadMe.txt
?????文件?????????733??2017-12-04?21:15??code\myDFT.m
?????文件???????15524??2017-12-04?20:14??code\signal.slx
?????文件?????????759??2017-12-04?18:47??code\sin.mat
?????目錄???????????0??2017-12-04?21:14??code\