| 查看: 956 | 回复: 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 |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有11人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
送红花一朵 
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
25