Question

matlab中使用FFT进行谐波分析的程序

本人毕业设计，需要用matlab进行谐波分析，希望高手给出应用FFT，加窗函数，以及插值三种方法在处理同一谐波数据时的对比，可参考如下数据
谐波次数 幅值（A） 相位（°）
1 380.4 0
3 2.035 60
5 3.156 135
7 1.042 157.5
9 4.291 0
11 16.13 60
13 2.016 0

根据表设信号表达式（谐波次数记为n，n=2k+1）

1. 用MATLAB进行仿真谐波电流信号图

2. 当取采样频率为3000HZ，数据长度为1024采样点时，利用三种方法编写MATLAB语言程序得到信号频谱，得出的谐波幅度谱

如果运行结果通过的话，高分悬赏，如有任何问题请私信联系我，如果可以发给我程序，可以邮箱

Answer 1

参考
%二.FFT应用举例

%例1：x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t)。采样频率fs=100Hz，分别绘制N=128、1024点幅频图。

clf;
fs=100;N=128; %采样频率和数据点数
n=0:N-1;
t=n/fs; %时间序列
x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t); %信号
y=fft(x,N); %对信号进行快速Fourier变换
mag=abs(y); %求得Fourier变换后的振幅
f=n*fs/N; %频率序列
subplot(2,2,1),plot(f,mag); %绘出随频率变化的振幅
xlabel('频率/Hz');
ylabel('振幅');title('N=128');grid on;
subplot(2,2,2),plot(f(1:N/2),mag(1:N/2)); %绘出Nyquist频率之前随频率变化的振幅
xlabel('频率/Hz');
ylabel('振幅');title('N=128');grid on;

%对信号采样数据为1024点的处理
fs=100;N=1024;n=0:N-1;t=n/fs;
x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t); %信号
y=fft(x,N); %对信号进行快速Fourier变换
mag=abs(y); %求取Fourier变换的振幅
f=n*fs/N;
subplot(2,2,3),plot(f,mag); %绘出随频率变化的振幅
xlabel('频率/Hz');
ylabel('振幅');title('N=1024');grid on;
subplot(2,2,4)
plot(f(1:N/2),mag(1:N/2)); %绘出Nyquist频率之前随频率变化的振幅
xlabel('频率/Hz');
ylabel('振幅');title('N=1024');grid on;

%运行结果：

%fs=100Hz，Nyquist频率为fs/2=50Hz。整个频谱图是以Nyquist频率为对称轴的。并且可以明显识别出信号中含有两种频率成
%分：15Hz和40Hz。由此可以知道FFT变换数据的对称性。因此用FFT对信号做谱分析，只需考察0~Nyquist频率范围内的福频特性。若没有给
%出采样频率和采样间隔，则分析通常对归一化频率0~1进行。另外，振幅的大小与所用采样点数有关，采用128点和1024点的相同频率的振幅是有不同的表
%现值，但在同一幅图中，40Hz与15Hz振动幅值之比均为4：1，与真实振幅0.5：2是一致的。为了与真实振幅对应，需要将变换后结果乘以2除以N。

追问

追答

自己改吧！

Answer 2

很多人问multilab的问题，，。可我就是不会啊