如何找到非线性信号的总谐波失真。使用matlab

Question

如何找到非线性信号的总谐波失真。例如，Forced Van der pol Oscillator，代码如下所示。我已经尝试过这个＆＃39; thd＆＃39;在matlab中运行，但我想我错过了一些东西。 这是等式    的 X＆＃39;＆＃39; -mu（1-x ^ 2（t））的 X＆＃39;（T）+ X（t）= P COS（W * t）的

function vdpo()
t=0:0.001:10
mu=2
x0=-2;
v0=2;
p=10; w=7;
 [t,x]= ode45(@f, t, [x0,v0])
function dxdt=f(t,x)
dxdt1=x(2); dxdt2= mu(1-x(1)^2)*x(1)+p*cos(w*t); 
dxdt=[dxdt1 ;dxdt2];
end
end

Answer 1

尝试下面的代码，其中函数f(t,x)是我们的ODE方程，我们调用函数ode45来使用Runge-Kutta方法来解决它。

function [x]=vdpo()
  t=0:0.001:10
  mu=2
  x0=-2;
  v0=2;
  p=10; w=7;
  [t,x]= ode45(@f, t, [x0,v0])
  function dxdt=f(t,x)
    dxdt1=-x(2)-x(1)+(x(1)^3)/3; dxdt2=-x(1)+p*cos(w*t); 
    dxdt=[dxdt1 ;dxdt2];
  end
end

然而，它实际上是一个数学问题而不是编程问题。我们要做的第一件事是通过定义y=x'+((x^3)/3-x)*mu将方程转换为更方便的形式，然后我们有2 First Order Ordinary Differential Equations所以我们可以调用ode45来解决它。我在here(page2)看了他们。

致电

X=vdpo();
x=X(:,1);
thd(x) 

我们可以得到如下答案：

P.S。不是关于THD部分。