Matlab freqz 函数实现

最新推荐文章于 2024-10-24 07:30:00 发布

Thomas会写字

最新推荐文章于 2024-10-24 07:30:00 发布

阅读量842

点赞数 1

分类专栏：数学计算文章标签：数据结构 java 算法

共享知识，欢迎转载

本文链接：https://blog.youkuaiyun.com/kakiebu/article/details/133385919

版权

数学计算专栏收录该内容

16 篇文章

订阅专栏

Matlab freqz 函数代码说明

%               jw               -jw              -jmw
%        jw  B(e)    b(1) + b(2)e + .... + b(m+1)e
%     H(e) = ---- = ------------------------------------
%               jw               -jw              -jnw
%            A(e)    a(1) + a(2)e + .... + a(n+1)e

代码实现

static public Complex[] Freqz(double[] B, double[] A, double[] X)
{
    if (A == null || B == null || X == null)
        return null;
    Complex[] dst = new Complex[X.Length];

    for (int i = 0; i < X.Length; i++)
    {
        Complex cb = new Complex(0, 0);
        Complex ca = new Complex(0, 0);

        for (int j = 0; j < B.Length; j++)
        {
            double m = j * X[i];
            Complex c = new Complex(Math.Cos(m), -Math.Sin(m));
            cb += B[j] * c;
        }

        for (int j = 0; j < A.Length; j++)
        {
            double m = j * X[i];
            Complex c = new Complex(Math.Cos(m), -Math.Sin(m));
            ca += A[j] * c;
        }

        dst[i] = cb / ca;
    }
    return dst;
}