| 查看: 940 | 回复: 4 | ||
[求助]
fft频谱分析
|
|
求助，大家帮帮忙，最近在学fft频谱分析，原始信号中含有衰减直流分量、衰减交流负分量、稳态交流分量、间谐波，去除衰减直流分量是减去其均值，但这这是用于周期的，具有衰减分量的提取不准，不知道怎么样较好的去除衰减直流分量和衰减交流分量，图片是生成正弦信号的（幅值不太准），用仿真数据做fft频谱分析的时候完全提不准，下面是我的程序 close all; %先关闭所有图片 N=1024; %采样点数 Fs=1024; %采样频率(Hz) t=[0:1/Fs:N/Fs]; %采样时刻 %信号 S1=100*sin(2*pi*50*t)+20*exp(-t/2)+50*sin(2*pi*35*t).*exp(-t/5)+50*sin(4*pi*50*t); %生成正弦信号 S=S1-mean(S1); plot(S); title('原始信号'); figure; Y = fft(S,N); %做FFT变换 Ayy = (abs(Y)); %取模 plot(Ayy(1:N)); %显示原始的FFT模值结果 title('FFT 模值'); figure; Ayy=Ayy/(N/2); %换算成实际的幅度 Ayy(1)=Ayy(1)/2; F=([1:N]-1)*Fs/N; %换算成实际的频率值 plot(F(1:N/2),Ayy(1:N/2)); %显示换算后的FFT模值结果 title('幅度-频率曲线图'); figure; Pyy=[1:N/2]; for i=1:N/2 Pyy(i)=phase(Y(i)); %计算相位 Pyy(i)=Pyy(i)*180/pi; %换算为角度 end; plot(F(1:N/2),Pyy(1:N/2)); %显示相位图 title('相位-频率曲线图');  $7$)JCXUCS9X~6V{)RLISI0.png |
回复此楼» 猜你喜欢
应该是93bebmhtak前后十一个字符比较关键
已经有33人回复
-
FileCode能看出啥？
已经有19人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有5人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有13人回复
-
是这周出结果还是下周出结果？
已经有3人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有9人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有4人回复
-
帮忙看看fileCode
已经有13人回复
mophez
木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 3707.4
- 散金: 175
- 红花: 11
- 沙发: 2
- 帖子: 1260
- 在线: 101.3小时
- 虫号: 4205858
- 注册: 2015-11-08
- 性别: GG
- 专业: 电力系统
2楼2017-10-11 08:06:59
3楼2017-11-12 15:59:15
mophez
木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 3707.4
- 散金: 175
- 红花: 11
- 沙发: 2
- 帖子: 1260
- 在线: 101.3小时
- 虫号: 4205858
- 注册: 2015-11-08
- 性别: GG
- 专业: 电力系统
ssddjj4楼2017-11-12 19:43:43
5楼2017-11-12 21:24:08
