系统辨识理论及应用(课后题答案第三章3233)国防工业.docVIP

下载本文档

35
0
约2.61千字
约 5页
2017-06-26 发布于重庆
举报
版权申诉

系统辨识理论及应用(课后题答案第三章3233)国防工业.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

系统辨识理论及应用(课后题答案第三章3233)国防工业

1、系统辨识——连续系统传递函数——脉冲传递函数 function h=Continuous_system_transferFcn(N,G,dt) % N——系统阶数 % G——采样数据（个数大于等于2N+1） % G为一维行向量 % dt——采样间隔 if nargin3 errordlg(not enough input varibles,error hint); else g_NN=zeros(N,N); for i=1:N g_NN(i,:)=G(i+1:i+1+N-1); end g_N=-G(1:N); a=inv(g_NN)*g_N; %% x的求解 syms x for i=1:N X(i)=x^i; end f=X*a+1; x=double(solve(f)); %%极点的求解 p=log(x)/dt; c_NN=zeros(N,N); for i=1:N c_NN(i,:)=x.^(i-1); end c_N=G(1:N);k=inv(c_NN)*c_N; p k z=zeros(1,N); p=p; k=k; Continuous_TransferFcn=0; for i=1:N Continuous_TransferFcn=Continuous_TransferFcn+zpk(z(i),p(i),k(i)); end Continuous_TransferFcn end end 例题 3.1(P32) G=[0 0.1924 0.2122 0.1762]; N=2; dt=1; Continuous_system_transferFcn(N,G,dt) p = -0.4934 -0.7085 k = 1.6280 -1.6280 Continuous_TransferFcn = 0.35024 s --------------------- (s+0.4934) (s+0.7085) Continuous-time zero/pole/gain model. 习题3.2(P34) G=[0 0.196 0.443 0.624 0.748 0.831]; N=3; dt=0.2; Continuous_system_transferFcn(N,G,dt) p = -0.0633 -1.7846 -11.1860 k = 1.1249 -1.3399 0.2150 Continuous_TransferFcn = -0.08507 s (s-253.1) ------------------------------- (s+0.06329) (s+1.785) (s+11.19) Continuous-time zero/pole/gain model. 2 系统辨识——离散系统传递函数——脉冲传递函数 function h=Discrete_system_transferFcn(N,G,dt) % N——系统阶数 % G——采样数据（个数大于等于2N+1） % G为一维行向量 % dt——采样间隔 if nargin3 errordlg(not enough input varibles,error hint); else g_NN=zeros(N,N); for i=1:N g_NN(i,:)=G(i+1:i+1+N-1); end g_N=-G(N+2:2*N+1); a1=inv(g_NN)*g_N; a=zeros(N,1); for j=1:N a(j,1)=a1(N+1-j,1); end B=zeros(N+1,N+1); B=diag(linspace(1,1,N+1)); for i=1:N+1 for j=1:N+1 if (i==j)(iN+1)(jN+1) B(i+1:N+1,j)=a(1:N+1-i,1); end end end g__N=G(1:N+1); b=B*g__N; a b num=b; den=[1 a]; Discrete_TransferFcn=tf(num,den,dt); Discrete_TransferFcn end end 例题 3.2(P33) G=[0 7.157039 9.491077 8.563839 5.930506 2.845972 0.144611]; N=3; dt=0.05; Discrete_system_transferFcn(N,G,dt) a = -