| 查看: 818 | 回复: 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 |
回复此楼» 猜你喜欢
想换工作。大多数高校都是 评职称时 认可5年内在原单位取得的成果吗？
已经有7人回复
-
需要合成515-64-0，50g，能接单的留言
已经有4人回复
-
自荐读博
已经有4人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有4人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
带资进组求博导收留
已经有10人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
-
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）
已经有3人回复
-
临港实验室与上科大联培博士招生1名
已经有8人回复
送红花一朵 
5楼2017-11-12 21:24:08
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
- 专业: 电力系统
4楼2017-11-12 19:43:43