连续扫频-FFT-系统辨识

阿W呀

于 2024-12-30 16:55:29 发布

阅读量1k

点赞数 8

分类专栏：系统辨识以及Matlab相关文章标签： matlab

本文链接：https://blog.csdn.net/catzhaojia/article/details/144829028

版权

系统辨识以及Matlab相关专栏收录该内容

6 篇文章

订阅专栏

利用连续扫频的方式进行系统模型辨识。

连续扫频相对于单点扫频来说，操作更简单，利用simulink中的“Chirp Signal”模块作为激励信号，设置随着时间变化的“变频”曲线对被控对象进行扫频，然后用FFT对输入信号与输出信号进行分析，得到各自的幅值、相位，然后求得输入、输出信号的幅值比、相位差，有了幅值比、相位差，同样就可以采用Simulink自带的模型辨识工具箱去获得被控对象的传递函数。

具体步骤：

01. Simulink模型

建立扫频模型如上图所示，并将系统步长设置为定步长0.0001，使用龙格库塔求解器（ode4）。后运行仿真，输入信号和输出信号用“to Workspace”模块将数据实时保存至“Workspace”中。

02. 编写 fft.m代码

对导入Workspace里的u=out.u.Data;outy=out.y.Data;数据进行FFT分析，并得到幅值比、相位差。

u=out.u.Data;
outy=out.y.Data;

close all;
fs = 10000; % 对应采样时间0.0001s。
N = length(u);
n = 0: (N - 1);
fn  = n * fs / N;
%% 
y = fft(u,N);
mag  = abs(y);
pha = angle(y);
%% 
y1 = fft(outy,N);
magl  = abs(y1);
phal = angle(y1);
%% 
%%计算幅值比，相位差
for i = 1 : N / 2
    mag_el(i) = magl(i)/mag(i);
    w(i) = fn(i) * 2 * pi;
    ph_el(i) = (phal(i) - pha(i)) * 180 / pi;
    mag_log(i) = 20 * log10(mag_el(i));
end
%%对相位差的数据进行预处理
for i = 1: length(ph_el)
    if ph_el(i) > 0
        ph_el(i) =  ph_el(i) - 360;
    end
end
%% 画Bode图
subplot(211);semilogx((w),(mag_log));grid on;
subplot(212);semilogx((w),(ph_el));grid on;