相位偏折术matlab代码

相位偏折是一种信号处理技术，通常用于对信号进行滤波或者变换操作。下面是一个简单的MATLAB实现相位偏折的代码示例：

% 设置输入信号
fs = 1000;  % 采样频率
t = 0:1/fs:1-1/fs;  % 时间序列
f1 = 10;  % 输入信号频率
f2 = 100;  % 滤波器中心频率
x = sin(2*pi*f1*t);  % 输入信号

% 设计滤波器
N = 50;  % 滤波器阶数
fc = f2/fs;  % 归一化截止频率
h = fir1(N, fc);  % 线性相位滤波器

% 相位偏折
y = filter(h, 1, x);  % 线性相位滤波
delay = mean(grpdelay(h));  % 滤波器延迟
y = y(delay+1:end);  % 去除延迟
y = real(y .* exp(-1j*2*pi*f2*t));  % 相位偏折

% 绘图
figure;
subplot(2,1,1);
plot(t, x);
title('输入信号');
xlabel('时间(s)');
ylabel('幅值');
subplot(2,1,2);
plot(t, y);
title('相位偏折后的信号');
xlabel('时间(s)');
ylabel('幅值');

该代码首先生成一个10Hz正弦波作为输入信号，然后设计一个中心频率为100Hz的FIR滤波器，使用filter函数进行线性相位滤波，并使用grpdelay函数求出滤波器的延迟。接着，对滤波后的信号进行相位偏折操作，最后将结果绘制成图表。