计算机仿真技术及CAD.rar

clear B=[0 1 0.5]; %建立仿真模型 A=[1 -1.5 0.7]; C=[1 -1 0.2]; th0=poly2th(A,B,C); u=sign(randn(350,1));%产生输入信号 e=randn(350,1);%产生干扰信号 y=idsim([u,e],th0);%计算系统输入 z=[y u]; %生成输入输出数据组 idplot(z); pause [Go,No]=spa(z);%进行数据谱估计分析 disp('Find ARX-Model by LS')%估计ARX-Model(LS法) A2=arx(z,[2 2 1]); present(A2) disp('Find ARX-Model by IV')%估计ARX-Model(IV法) I2=iv4(z,[2 2 1]); present(I2) e1=resid(z,A2);pause%计算比较两种估计模型残差 e2=resid(z,I2);pause plot([e1,e2]); title('ARX-Model by LS and one by IV') grid pause yarxls=idsim(u,A2);%利用函数idsim进行模型仿真 subplot(221) %比较两种估计模型的仿真结果 plot([y,yarxls]); title('Ral output and ARX-LS model output') grid subplot(222) yarxiv=idsim(u,I2); plot([y,yarxiv]); title('Real output and ARX-IV model output') subplot(223) compare(z,A2);%利用函数compare比较两种模型仿真结果 title('Real output and ARX-IV model output') grid subplot(224) compare(z,I2); title('Real output and ARX-IV model output') grid pause disp('Compare transfer function') %比较仿真模型和两种估计模型的频率特性 [Garxls,Narxls]=th2ff(A2); [Garxiv,Narxiv]=th2ff(I2); bodeplot([Go,Garxls,Garxiv]); grid pause bodeplot(No,Narxls,Narxiv); grid