FFT谐波分析

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）;