用你写好的fft函数编写一个chirp z变换的c代码，记住要写全

Chirp Z变换是一种在数字信号处理中用于分析随时间变化频率（即调制）信号的方法，它是傅里叶变换（FFT）和Z变换的结合。以下是一个简单的C语言示例，使用自定义的FFT函数来计算Chirp Z变换。这里假设你已经有了一个名为`fft()`的快速傅立叶变换函数。请注意，实际项目中可能需要根据具体的硬件平台或库进行优化。

#include <stdio.h>
#include <math.h>

// 假设你有一个自定义的FFT函数
void fft(double complex *in, double complex *out, int n) {
    // 这里省略了FFT的具体实现，你可以根据实际情况提供
}

// Chirp信号的参数
#define SAMPLING_PERIOD 0.001
#define START_FREQ 1000
#define END_FREQ 5000
#define LENGTH 10000

double chirpz(double t) {
    return START_FREQ + (END_FREQ - START_FREQ) * pow(t / SAMPLING_PERIOD, 2);
}

int main() {
    double time[LENGTH];
    double freq_domain_data[LENGTH/2 + 1]; // 只取一半结果，因为正负频是互为镜像的

    for (int i = 0; i < LENGTH; i++) {
        time[i] = i * SAMPLING_PERIOD;
        freq_domain_data[i] = 0;
    }

    // 计算Chirp信号
    for (int i = 0; i < LENGTH; i++) {
        freq_domain_data[i] = chirpz(time[i]) * sin(2 * M_PI * chirpz(time[i]) * time[i]);
    }

    // 应用FFT
    fft(freq_domain_data, freq_domain_data, LENGTH); // 注意调整长度为2的幂

    // 打印结果
    printf("Chirp Z Transform:\n");
    for (int i = 0; i <= LENGTH/2; i++) {
        printf("%.4f + j%.4f\n", creal(freq_domain_data[i]), cimag(freq_domain_data[i]));
    }

    return 0;
}