相位梯度进行积分 matlab代码

以下是使用Matlab实现相位梯度积分的示例代码：

% 生成信号
f0 = 10; % 基频
fs = 1000; % 采样率
t = 0:1/fs:1-1/fs; % 时间向量
x = sin(2*pi*f0*t); % 正弦波信号

% 相位梯度积分
phi = unwrap(angle(hilbert(x))); % 求解相位
dphi = diff(phi); % 相位梯度
dphi = [dphi(1) dphi]; % 补齐长度

% 绘图
subplot(2,1,1)
plot(t,x)
xlabel('时间 (s)')
ylabel('幅值')
title('原始信号')

subplot(2,1,2)
plot(t,dphi)
xlabel('时间 (s)')
ylabel('相位梯度')
title('相位梯度积分')

代码中首先生成一个10Hz的正弦波信号，然后通过Hilbert变换求解相位，再对相位进行差分得到相位梯度，最后绘制相位梯度随时间的变化图。