- 1、本文档共15页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数据采集与信号处理作业基于FFT的功率谱分析程序设计与应用
数据采集与信号处理作业
一、基本内容: 基于FFT的功率谱分析程序设计与应用
1.基本要求
1) 对一个人为产生的信号进行采用FFT变换方法进行功率谱分析。
已知信号x(n)=120.0*COS(2*3.14*SF*n/FS)
式中: n=0,1,2 ……N-1
SF---信号频率
FS---采样频率
其FFT变换结果X(k)可用下面提供的FFT子程序求出,计算功率谱的公式为:
W(k)=2(XR(k)2 +XI(k)2)/N
式中: k=0,1,2 ……N/2-1
XR(k)--- X(k)的实部
XI(k)--- X(k)的虚部
请用VB,VC或C++Builder编译器编程,或采用MATLAB计算,或采用高级语言调用MATLAB计算。处理结果为采用窗口显示时域波形和频域波形。
基于MATLAB6.5的源程序如下:
clf;
clc;
Fs=300; %采样频率
n=0:1:300;N=128;
Sf=20; %信号频率
xn=120*cos(2*pi*Sf*n/Fs); %产生波形序列
cxn=xcorr(xn,unbiased); %计算序列的自相关函数
figure(1);plot(n,xn);grid on;title(余弦曲线的时频图);
nfft=1024;
CXk=fft(cxn,nfft);
Pxx=abs(CXk); %求序列的PSD;
index=0:round(nfft/2-1);
k=index*Fs/nfft;
plot_Pxx=Pxx(index+1);
figure(2);plot(k,plot_Pxx);title(余弦曲线的功率谱曲线图);
xlabel(Frequency(Hz));
ylabel(Magnitude);
grid on;
2)对实验所采集的转子振动信号进行频谱分析
SF=1000; %采样频率1024/508HZ
fid = fopen(D:\fanbo_45HZ_1024_1000HZ\sanjiao_45HZ_1024_1000HZ);%转子信号
[a,N]= fscanf(fid,%f);
fclose(fid);
y=fft(a,N); %FFT运算
Pyy =sqrt(y.*conj(y))*2.0/N; %取功率普密度
f=(0:length(Pyy)-1)*SF/length(Pyy);
LPyy=20*log10(Pyy);
plot(f(1:N/2),Pyy(1:N/2),black); %输出FS/2点幅频谱图
grid;
2.讨论
1)信号经过均值化处理或不经过均值化处理的结果比较
%含直流分量而未均值化的信号
clc;clf;
Fs=300; %采样频率
n=0:1:300;N=128;
Sf=20; %信号频率
xn=120*cos(2*pi*Sf*n/Fs)+60; %产生波形序列
window=boxcar(length(xn)); %矩形窗
nfft=512; %采样点数
[Pxx,f]=periodogram(xn,window,nfft,Fs); %直接法
subplot(2,1,1);
plot(f,Pxx);title(含有直流分量的余弦曲线未均值化的功率谱波形图);
xlabel(Frequency(Hz));
ylabel(Magnitude);
%含直流分量且均值化的信号
clc;clf;
Fs=300; %采样频率
n=0:1:300;N=128;
Sf=20; %信号频率
xn=120*cos(2*pi*Sf*n/Fs)+60;%产生波形序列
z=mean(xn);
h=xn-z;
window=boxcar(length(xn)); %矩形窗
nfft=512;%采样点数
[Pxx,f]=periodogram(h,window,nfft,Fs); %直接法
subplot(2,1,2);
plot(f,Pxx);title(含有直流分量的余弦曲线均值化后的功率谱波形图);
xlabel(Frequency(Hz));
ylabel(Magnitude);
余弦函数信号x的直流分量即x的均值,它的存在将导致余弦函数在f0 (即图形的波峰处的频率)处的功率降低,因此,必须对带有直流分量的信号进行均值化处理,以消除直流分量的干扰。由图3可以看出,带有直流分量的余弦函数不经均值化处理时的功率谱在低频段含有一个很大的功率分量,导致余弦函数在f0处的功率(或主频率分量)降低。而经均值化处理后,余弦函数在f0处的功率明显增加,这
您可能关注的文档
- 广场车展策划部活动方案.doc
- 建筑涂料施工质量管理措施.doc
- 建筑节能及绿色建筑的发展运用策略探讨.doc
- 建筑施工测量方法研究.doc
- 建筑装饰工程公司管理手册.doc
- 建筑施工图教学课件PPT.ppt
- 建设工程有限公司重大安全隐患整改方案.doc
- 建设工程竣工验收报告空白模板.doc
- 建设工程招投标与合同管理FIDIC施工合同条件教学PPT.ppt
- 彩妆化妆品市场营销策划方案.doc
- 四川省德阳市罗江中学2025届高三考前热身化学试卷含解析.doc
- 山东省枣庄现代实验学校2025届高三下学期第五次调研考试化学试题含解析.doc
- 吉林省长春市十一高中等九校教育联盟2025届高三一诊考试生物试卷含解析.doc
- 2025届江苏省盐城市伍佑中学高考仿真模拟化学试卷含解析.doc
- 2025届广西贺州中学高考冲刺押题(最后一卷)生物试卷含解析.doc
- 安徽省池州市贵池区2025届高三第一次模拟考试生物试卷含解析.doc
- 宁夏银川一中2025届高三(最后冲刺)化学试卷含解析.doc
- 广东省广州市增城区四校联考2025届高考压轴卷化学试卷含解析.doc
- 2025届邯郸市第一中学高考生物必刷试卷含解析.doc
- 2025届安徽省安庆市石化第一中学高考仿真卷化学试卷含解析.doc
最近下载
- 山东省泰安市泰山区2023-2024学年上学期期中考试七年级英语试题.docx VIP
- 财务管理实务 财务管理实务 实训指导书.doc
- 2024年新人教版九年级上册化学全册教案.docx
- 上海通用汽车GVDP整车开发流程.pdf VIP
- 苏教版四年级上册同步奥数培优 第十四讲 智巧问题.pdf VIP
- 2024年江苏省苏州高新区招聘“两新”组织党建专职党务工作者6人历年【综合基础知识500题】高频考点模拟试题及参考答案解析.docx VIP
- 2024年武汉市东湖生态旅游风景区管理委员会公开招聘机关派遣制工作人员9人笔试备考试题及答案解析.docx VIP
- 人教版普通高中地理选择性必修2区域发展.pdf
- 人体解剖学常考重要知识点.doc
- 2023年江苏苏州相城区招考聘用两新组织党建专职党务工作者23人笔试历年难易错点考题荟萃附带答案详解.docx VIP
文档评论(0)