%%%做系统识别很重要,方法上完全符合系统识别最基础的理论 function [sun]=main(n) fplot(,],'r'); x=ones(n+,); :n+ x(j+)=cos(pi*(n+-j)/(n+)); end first=ones(n+,); f=./(x+); %原函数 last=first; :n+ last(j)=(-)*last(j-); end A=ones(n+,n+); A(:,)=first; A(:,n+)=last; :n+ :j A(:,j)=x.*A…