?clear?all;?close?all;?clc;?

fs=1000;
N=1024;
t=0:1/fs:N/fs;

x=2*cos（2*pi*40*t）;

degree2=90;
Dfreq=degree2/360;%fs/4?變頻
Dphase=0;%0.1;???????
for?k=1:length（x）
????Dphase=Dpha