example2_2.m

一个黑盒对输入输出矩阵进行辨识模型并输出系数

clear all;  hold off;
num=[1  3.6];                                         %传递函数的分子
den=[1 5.8  4.8];                                       %传递函数的分母
w=[0.01, 0.02,0.03,0.04,0.05,0.06,0.07,0.08,0.1,0.13,0.2].*2.*pi;             %角频率

H=freqs(num,den, w);                                    %频率响应特性
na=2;    nb=1;
P=real(H);  Q=imag(H);  M=length(w);                     %频率响应特性的实部和虚部

Am=sqrt(P.*P+Q.*Q);
Phi=atan(Q./P).*57.3;

Am =[ 0.75, 0.74, 0.74, 0.73,  0.72,  0.70 , 0.69,  0.68,  0.64 ,  0.59 , 0.48];                   %幅频特性
Phi=[ -0.06,  -0.12,  -0.17,  -0.23,  -0.28, -0.33,  -0.38,  -0.45,  -0.52, -0.63,  -0.82];        %相频特性

P=Am.*cos(Phi);                                          %频率响应特性的实部
Q=Am.*sin(Phi);                                          %频率响应特性的虚部

H=P+i.*Q;                                                %测得的频率响应特性


[n,d]=freq2levy(H,w,2,1);                                 %  Levy法频率特性拟合

[n1,d1]=invfreqs(H,w,1,2);                               % 工具箱中的频率特性拟合函数

plot(log(w),Phi);