生成一段Chirp信号进行分数阶傅里叶变换并画出三维图的MATLAB代码

以下是生成并绘制分数阶Chirp信号的MATLAB代码：

% 生成分数阶Chirp信号
alpha = 0.8; % 分数阶阶数
t = 0:0.001:1; % 时间范围
f0 = 10; % 起始频率
f1 = 100; % 终止频率
y = chirp(t, f0, 1, f1, 'quadratic'); % 生成Chirp信号
yf = fracFourier(y, alpha, t); % 计算分数阶傅里叶变换

% 绘制三维图
figure;
fsurf(real(yf), imag(yf), abs(yf));
xlabel('Re');
ylabel('Im');
zlabel('Amplitude');
title(['Fractional Fourier transform of Chirp signal, \alpha = ' num2str(alpha)]);

其中，`fracFourier`函数是进行分数阶傅里叶变换的自定义函数，其代码如下：

function yf = fracFourier(y, alpha, t)
    N = length(y);
    yf = zeros(size(y));
    for k = 1:N
        tk = t(k);
        for n = 1:N
            tn = t(n);
            yf(k) = yf(k) + y(n)*exp(-1i*2*pi*tk*tn/(N^(1-alpha)));
        end
        yf(k) = yf(k) * (t(2)-t(1))^alpha;
    end
end

运行代码后，将会生成一个分数阶Chirp信号的三维图，其中X轴表示实部，Y轴表示虚部，Z轴表示振幅。