frd matlab,频率响应转换成状态空间，frd算法实现

已知测试或者仿真的snp，也就是一个系统的频率响应H(f)。如频率f=0.010:0.01:20(GHz)，每个频点对应一个n*n*2的向量(幅度Mag，相位Angel)矩阵，如下为10MHz和20MHz的数据

!                  MAG       ANG               MAG       ANG               MAG       ANG               MAG       ANG

10.00000000       0.006484    52.120192       0.999115    -0.904316       0.003935    89.507948       0.001519   -81.671603

0.999598    -0.865390       0.010093   111.563606       0.001550   -80.601228       0.003930    89.485880

0.003992    89.056528       0.001495   -79.770792       0.002632   109.538554       0.999102    -0.852198

0.001491   -80.735664       0.003966    88.602613       0.999311    -0.852870       0.002605   108.782551

20.00000000       0.012766    81.674063       0.997753    -1.684988       0.007881    87.911435       0.003112   -87.769232

0.999070    -1.680573       0.012187    99.972446       0.003087   -86.348349