获取真实信号的谐波(Matlab代码实现）

 💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

---

💥1 概述

• 要获取真实信号的谐波，可以使用信号处理技术，如傅里叶变换。下面是一些步骤：

  1. 采集信号：使用传感器或仪器采集真实信号的样本数据。确保采样频率足够高，以捕捉到信号的高频成分。

  2. 预处理信号：对采集到的信号进行预处理，例如去除噪声、滤波等。这可以提高后续分析的准确性。

  3. 进行傅里叶变换：将预处理后的信号应用傅里叶变换，将信号从时域转换到频域。傅里叶变换将信号分解为一系列频率成分，包括基频和谐波。

  4. 分析频谱：分析傅里叶变换后的频谱图，查找主要频率成分和对应的幅度。谐波通常是基频的整数倍，可以通过查找频谱图中的峰值来确定谐波的存在。

  5. 提取谐波：根据频谱分析的结果，提取出谐波成分的频率和幅度信息。可以根据需要选择提取特定谐波或提取所有谐波。

  请注意，以上步骤是一种常见的方法，具体的步骤可能因应用和信号特性而有所不同。此外，信号的采样频率和长度也会影响谐波的分析结果。

📚2 运行结果

 

部分代码：

function [wc,w0,a0,ak,bk,c0,ck]=get_harmonics(y,pas)
%[wc,w0,a0,ak,bk,c0,ck]=get_harmonics(y,pas)
% given a signal x(t),y=x(t)
% the function get_harmoniques returns the pulse w0(rd/s) of the signal y(t)
% the coefficients ak (a1,a2,...)  and bk (b1,b2,...)  of the trigonomical
% representation  of the  signal y(t):
% y(t)= a0/2+ a1.cos(w0t)+b1.sin(w0t)+a2.cos(2w0t) +b2*sin(t*w0)+...
% it also returns the coefficients c0 and ck of the complex representation
% y(t)= ...+c(-1).exp(-jw0t)+c(0)+c(1).exp(jw0t)+c(2).exp(2jw0t)+...
% wc is the interval pulse;
% to plot the amplitude spectrum: stem(wc,abs(ck))
%[wc,w0,a0,ak,bk,c0,ck]=get_harmoniques(y,pas)
% be a signal x(t).
% For more theorical details look at www.azzimalik.com
N=length(y);
Y=fft(y);
if size(Y,2)==1;
    Y=Y';
end
w=2*(0:(N-1))*pi/(N*pas);
w0=w(2);
nc2=fix(N/2);
if nc2~=N/2;
    nc3=nc2+1;
    wc=[fliplr(-w(2:nc2+1)) w(1:nc2+1)];
    Yc=[Y(nc2+2:N) Y(1:nc2+1)];Yc=Yc/N;
else
    nc3=nc2;
    wc=[fliplr(-w(2:nc2+1)) w(1:nc2+1)];
    Yc=[Y(nc2+1:N) Y(1:nc2+1)];Yc=Yc/N;
end 

🎉3 参考文献

部分理论来源于网络，如有侵权请联系删除。

1. Bracewell, R. N. (2000). The Fourier Transform and Its Applications. McGraw-Hill Education.

2. Oppenheim, A. V., & Schafer, R. W. (2010). Discrete-Time Signal Processing. Pearson Education.

🌈4 Matlab代码实现