function [At,Ph,A2]=EnvelopPhase(X,f0,fs) HX=imag(hilbert(X)); [M,N]=size(X); t=0:1/fs:((N-1)/fs); Ac=X.cos(2*pi*f0*t)+HX.sin(2*pi*f0*t); As=HX.cos(2*pi*f0*t)-X.sin(2*pi*f0*t); Ph=atan(As./Ac); A2=Ac.*Ac+As.*As; At=sqrt(A2);这段代码正确吗

这段代码可能存在一些问题。

首先，`X.cos`，`X.sin`，`HX.cos` 和 `HX.sin` 都不是有效的 Matlab 语法。如果 `X` 和 `HX` 都是矩阵，则需要使用 `.*` 而不是 `*` 来进行逐元素的乘法和除法。

其次，`Ac` 和 `As` 的计算公式似乎有误，应该是：

Ac = X.*cos(2*pi*f0*t) + HX.*sin(2*pi*f0*t);
As = HX.*cos(2*pi*f0*t) - X.*sin(2*pi*f0*t);

最后，对于 `Ph`，需要考虑分母为零的情况，可以使用 `atan2` 函数来避免这种情况：

Ph = atan2(As, Ac);

综上所述，修正后的代码如下：

function [At, Ph, A2] = EnvelopPhase(X, f0, fs)
    HX = imag(hilbert(X));
    [M, N] = size(X);
    t = 0:1/fs:((N-1)/fs);
    Ac = X.*cos(2*pi*f0*t) + HX.*sin(2*pi*f0*t);
    As = HX.*cos(2*pi*f0*t) - X.*sin(2*pi*f0*t);
    Ph = atan2(As, Ac);
    A2 = Ac.*Ac + As.*As;
    At = sqrt(A2);
end