| 查看: 937 | 回复: 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 |
回复此楼» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
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
100