f1.m

多变量极点配置

%将非循环矩阵A转化为循环矩阵A1
function [a,k1]=f1(b)
global k1
k1=fix(10*rand(n2,n1));
a=A-b*k1
return