- 1、本文档共5页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
FFT谐波频谱相谱分析Matlab程序二则
FFT谐波频谱相谱分析Matlab程序二则
by ggihhimm
程序一:从Excel表中读取数据,然后分析其频谱、相谱和总谐波畸变率,并画出茎状图(stem),图示如下(测试数据含直流分量):
%该程序从Excel表中读取数据,然后分析其频谱、相谱和总谐波畸变率
%一组数据+采样率 显示频谱 下面的阀值应适当调大 注意I1rms的定义
clear all;
close all;
clc;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
y=xlsread(HCwithDI_2.xls); %一列数据:excel表格的第一列要知道数据的采样率,为
%保证结果的准确性,最好为整数倍周期
L=length(y);
mdflag=1; %是否标点,1是,0否
del=1e-1; %删除的阀值,根据数据幅值大小来设置
Fs=50*L; %采样率
f0=50; %信号基波频率
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
f=0:L-1;
f=f*Fs/L/f0;
yfft=fft(y);
yfft(L/2+1:L)=[];
f(L/2+1:L)=[];
mag=abs(yfft)/L*2;
ang=angle(yfft)*180/pi;
mag(1)=mag(1)*0.5;
mag1=mag;
mag1(mag=del)=[];
ang1=ang;
ang1(mag=del)=[];
f1=f;
f1(mag=del)=[];
L1=length(f1);
subplot(211);
hdl1=stem(f1(1:L1),mag1(1:L1),b:o);
set(hdl1,MarkerFaceColor,red,Marker,square);
ylabel(各次谐波幅值/A);
tem=get(gca,XLim);
xwidth=tem(end)-tem(1);
x=zeros(1,L1);
zhi=zeros(2,L1);
dotstr=cell(2,L1);
xright=1.5;
x(1)=f1(1)+xwidth/100*xright;
if f1(L1)xwidth-1
x(L1)=f1(L1)+xwidth/100*xright;
else
x(L1)=f1(L1)-xwidth/100*10;
end
zhi(1,1)=mag1(1);
zhi(2,1)=ang1(1);
zhi(1,L1)=mag1(L1);
zhi(2,L1)=ang1(L1);
dotstr{1,1}=num2str(mag1(1),%8.3f);
dotstr{2,1}=num2str(ang1(1),%8.3f);
dotstr{1,L1}=num2str(mag1(L1),%8.3f);
dotstr{2,L1}=num2str(ang1(L1),%8.3f);
for i=2:L1-1
x(i)=f1(i)+xwidth/100*xright;
zhi(1,i)=mag1(i);
zhi(2,i)=ang1(i);
zhi(1,i)=mag1(i);
zhi(2,i)=ang1(i);
dotstr{1,i}=num2str(mag1(i),%8.3f);
dotstr{2,i}=num2str(ang1(i),%8.3f);
end
%计算Irms和THD
I1rms=max(mag1(abs(f1-1)0.5))/sqrt(2);
Irms2=0.5*sum(mag1(2:L1).^2);
if f1(1)==0
Irms2=Irms2+mag1(1)^2;
dotstr{2,1}=;
else
Irms2=Irms2+0.5*mag1(1)^2;
end
Irms=sqrt(Irms2)
THD=sqrt(Irms2/I1rms^2-1)*100
if mdflag
for i=1:L1
text(x(i),zhi(1,i),dotstr{1,i});
end
end
subplot(212);
hdl2=stem(f1(1:L1),ang1(1:L1),m:o);
set(hdl2,MarkerFaceColor,red,Marker,
您可能关注的文档
- dvd光驱不读盘及笔记本光驱不读盘的处理办法.docx
- d杆菌对温度制粒工艺及pH值的耐受性.pdf
- E. coli 谷氨酸脱羧酶高产菌株选育及发酵条件研究.pdf
- EAST中性束注入器用高压缓冲器分析.pdf
- Eclipse+MyEclipse的优化.doc
- ECR 用于大面积薄膜研究.pdf
- EDTA标准溶液的标定及水硬度的测定.ppt
- eetop.cnGPS 接收机的灵敏度分析.pdf
- eetop.cn环路参数计算与噪声研究.pdf
- EGFRTKIs治疗非小细胞肺癌进展.pdf
- 中国国家标准 GB/T 18233.4-2024信息技术 用户建筑群通用布缆 第4部分:住宅.pdf
- GB/T 18233.4-2024信息技术 用户建筑群通用布缆 第4部分:住宅.pdf
- GB/T 18978.210-2024人-系统交互工效学 第210部分:以人为中心的交互系统设计.pdf
- 《GB/T 18978.210-2024人-系统交互工效学 第210部分:以人为中心的交互系统设计》.pdf
- 中国国家标准 GB/T 18978.210-2024人-系统交互工效学 第210部分:以人为中心的交互系统设计.pdf
- GB/T 16649.2-2024识别卡 集成电路卡 第2部分:带触点的卡 触点的尺寸和位置.pdf
- 《GB/T 16649.2-2024识别卡 集成电路卡 第2部分:带触点的卡 触点的尺寸和位置》.pdf
- 中国国家标准 GB/T 16649.2-2024识别卡 集成电路卡 第2部分:带触点的卡 触点的尺寸和位置.pdf
- GB/T 17889.4-2024梯子 第4部分:铰链梯.pdf
- 《GB/T 17889.4-2024梯子 第4部分:铰链梯》.pdf
文档评论(0)