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=idfrd(A.*exp(1j*phy*pi/180), f*2*pi,0.1) sys = tfest(Z,2,1);

你的代码看起来是用 MATLAB 绘制幅频特性曲线和相频特性曲线，并且根据这些特性曲线绘制频率响应。你的代码中使用了 xlsread 函数读取 Excel 文件中的数据，然后使用 loglog 函数绘制幅频特性曲线和相频特性曲线。最后，你使用 idfrd 函数创建一个频率响应对象 Z，并使用 tfest 函数估计系统传递函数。请问你的代码有没有报错？