代码：
% x = [-0.9, 0.81]; [y, L, B] = QCoeff(x, 3) % Unquantized parameters r = 0.9; theta = pi/3; a1 = -2*r*cos(theta); a2 = r*r; p1 = r*exp(j*theta); p2 = p1';% Quantized parameters: N = 3; [ahat, L, B] = QCoeff([a1, a2], 3); rhat = sqrt(ahat(2)); thetahat = acos(-ahat(1)/(2*rhat)); p1hat = rhat*exp(j*thetahat); p2 = p1';% Changes in pole locations Dp1 = abs(p1 - p1hat) Dp1 = 0.0344
运行结果：
转载于:https://www.cnblogs.com/ky027wh-sx/p/6601448.html
【描述】 公元19881231年，一颗巨大的陨石坠落在世界的政治文化中心cs。于是，灾难降临了，地球上出现了一系列反常的现象。当人们焦急万分的时候，一支由cs科学家组成的考察队赶到了出事地点。经过一番侦察，科学家们发现陨石上刻有若干行密文，每一行都包含5个整数： 1 1 1 1 6 0 0 6 3 57 8 0 11 3 284...
点击打开链接 //求SUM（gcd(i,n), 1<=i<=n） /*g(n)=gcd(i,n),根据积性定义g(mn)=g(m)*g(n)(gcd(m,n)==1)所以gcd(i,n)是积性的，所以f(n)=sum(gcd(i,n))是积性的，f(n)=f(p1^a1*p2^a2*...*pn^an)=f(p1^a1)*f(p...