clc
clear?all???%清除內存所有變量
close?all????%關閉所有打開的圖形窗口
N=512;??%?信號長度（變量@@@@@@@）
f1=50.3;
fs=2500;??%?采樣頻率
dt=1/fs;??%?采樣間隔
t=[0:N-1]*dt;??%?時間序列
xn=0;
A=[207.25.54.63.52.31.61.11.30.90.60.40.50.70.850.640.510.390.40.20.1];
theta=[756042263266582784654662157243231867351948];
for?h=1:21
xn=xn+A（h）*cos（2*pi*f1*h*t+theta（h）*pi/180）;??
end
xn=[xnzeros（1N-512）];??%?原始信號的值序列
%subplot（321）;
figure（1）;
plot（txn）
xlabel（‘時間/s‘）ylabel（‘電壓/V‘）title（‘光伏發電系統電壓波形‘）??%?變量@@@@@@@
NN=2^nextpow2（N）;??%?執行512點FFT
XN=fft（xnNN）/NN*2;??%?共軛復數，具有對稱性
f0=1/（dt*NN）;?