1fft 是离散fourier变换后的数据,也就是谱但它是用快速传立叶算法来实现的离散传立叶变换
2这个说明低频能量高,n对应的是基频的位置,但不是频率本身,因为离散传立叶变换的频率也是离散的,而且跟采样率有关,但不是个简单的序收
3直接用 fft,你就得到谱了,然后,把它fftshift一下,使基频位於中心,这样中心向两边的谱值就(近似)对应不同的谐波成分
disp(输出函数 -------------=输出值;
如:要输出 2z^2+3z -------------=a,
则:disp('2z^2+3z -------------=a')
顺祝你2015幸福快乐。
clc;%清理窗口显示内容
clear %清除工作间变量
close all;%关闭所有图形
n=input('n=');%命令窗中手动输入,2,3,………谐波次数
T0=001;%周期
A=1;%振幅
t=0:00001:005;%时间
z=square(2pi100t);%方波脉冲
x=0;%均值,0
pi=31415926;
w0=2pi/T0;%周期数,频率
for i=1:n
fw(i)=(2i-1)w0;%每次谐波频率
a(i)=4A/(pi(2i-1));%每次谐波的振幅
y(i,:)=a(i)sin(fw(i)t);%信号
fai(i)=0;%相位,没用到,可去掉
x=x+y(i,:);%合成信号
figure(1)%生成第一图窗
a=1:100;%定义无效,可去掉
subplot(n,1,i);%生成子图
plot(t,y(i,:));%画子图
xlabel('时间/s');ylabel('振幅')%x,y轴名称
title(['方波',num2str(a(i)),'次谐波'])%图形名称
end;
figure(2);%生成第二图形窗
subplot(2,1,1);%生成子图
plot(t,z);%画子图
xlabel('时间/s');%x轴名称
axis([-0,005,-2,2]);%设置轴范围
ylabel('振幅')%y轴名称
title('方波')%子图名称
subplot(2,1,2);%生成第二子图
plot(t,x);%画子图
xlabel('时间/s');%x轴名称
ylabel('振幅');%y轴名称
title([num2str(a(i)),'次谐波合成图'])%子图名称
以上就是关于matlab 含有谐波的傅立叶变换全部的内容,包括:matlab 含有谐波的傅立叶变换、如何在matlab实现5次谐波提取、求matlab大神解释一下这个程序每一步的含义,越详细越好,跪谢(悬赏可再加)等相关内容解答,如果想了解更多相关内容,可以关注我们,你们的支持是我们更新的动力!
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)