clc clear all close all f = xlsread('C:\Users\Lenovo\Desktop\幅频-相频特性曲线.xlsx','B2:B56'); % 读取原始频率 A = xlsread('C:\Users\Lenovo\Desktop\幅频-相频特性曲线.xlsx','C2:C56'); % 读取原始幅值 phy = xlsread('C:\Users\Lenovo\Desktop\幅频-相频特性曲线.xlsx','D2:D56'); % 读取原始相位 %% 绘制幅频特性曲线 figure(1); loglog(f ,A,'b'); %x2=linspace(min(f),max(f),500); %yy = spline(f, A, x2); %semilogy(x2, yy, '-r'); xlabel('频率（Hz)') ylabel('幅值（mV/nT)') title('幅频特性曲线') %% 绘制相频特性曲线 figure(2); loglog(f ,phy,'b'); xlabel('频率（Hz)') ylabel('相位（°)') title('相频特性曲线') %% 根据幅频以及相频特性曲线绘制频率响应 %Z=A.*exp(-j*phy); Z = tfest(idfrd(A.*exp(1j*phy*pi/180), 2*pi*f), 2, 1);

这段代码的作用是读取一个 Excel 文件中的数据，并绘制幅频特性曲线和相频特性曲线。同时，它还使用 `tfest` 函数基于幅频和相频特性曲线来估计系统的频率响应。

其中，第一部分是读取数据：

f = xlsread('C:\Users\Lenovo\Desktop\幅频-相频特性曲线.xlsx','B2:B56'); % 读取原始频率
A = xlsread('C:\Users\Lenovo\Desktop\幅频-相频特性曲线.xlsx','C2:C56'); % 读取原始幅值
phy = xlsread('C:\Users\Lenovo\Desktop\幅频-相频特性曲线.xlsx','D2:D56'); % 读取原始相位

这个代码块使用 `xlsread` 函数从 Excel 文件中读取三列数据，分别是频率、幅值和相位。

第二部分是绘制幅频和相频特性曲线：

figure(1);
loglog(f ,A,'b');
xlabel('频率（Hz)')
ylabel('幅值（mV/nT)')
title('幅频特性曲线')

figure(2);
loglog(f ,phy,'b');
xlabel('频率（Hz)')
ylabel('相位（°)')
title('相频特性曲线')

这个代码块使用 `loglog` 函数绘制幅频和相频特性曲线，并为图像添加标签和标题。

第三部分是使用 `tfest` 函数估计系统的频率响应：

Z = tfest(idfrd(A.*exp(1j*phy*pi/180), 2*pi*f), 2, 1);

这个代码块使用 `idfrd` 函数创建一个频率响应模型，然后使用 `tfest` 函数基于该模型来估计系统的频率响应。其中，`2` 表示模型的阶数，`1` 表示模型的采样时间。最终的结果保存在变量 `Z` 中。