资源描述：
《电磁铁，数据资料.doc》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、升磁过程降磁过程输入电压传感器输出输入电压传感器输出00.002303.82710.021293.74620.029283.71530.043273.64740.062263.61850.097253.38960.108243.21370.197233.03380.277222.79190.325212.611100.424202.386110.579192.214120.721182.026130.834171.793140.886161.624151.081151.513161.223141.299171.374131.217181.52712
2、1.044191.722110.922201.834100.797211.95390.624222.18980.515232.38570.394242.67260.331252.81150.206263.04640.189273.19730.105283.42620.097293.58810.035303.82700.002电磁铁加电压和去电压时的曲线及模型：loadrise.txtloaddown.txtu1=rise(:,1);f1=rise(:,2);u2=down(:,1);f2=down(:,2);使用曲线拟合的方法辨识电磁铁模型：a=p
3、olyfit(u1,f1,4)f11=polyval(a,u1);a=-0.0000-0.00010.0068-0.01820.0218b=polyfit(u2,f2,4)f21=polyval(b,u2);b=-0.00000.0005-0.00170.0566-0.0272曲线比较：拟合后的电磁铁模型曲线：使用系统辨识的方法的到的模型：升磁曲线：na=2;nb=2;nc=1;nk=0;m1=armax(z1,[na,nb,nc,nk])Discrete-timeIDPOLYmodel:A(q)y(t)=B(q)u(t)+C(q)e(t)A(q)
4、=1-1.409(+-0.2812)q^-1+0.4393(+-0.2656)q^-2B(q)=-0.002376(+-0.009807)+0.01037(+-0.01396)q^-1C(q)=1-1.01(+-0.3598)q^-1EstimatedusingARMAXfromdatasetz1Lossfunction0.000951157andFPE0.00142674Samplinginterval:1去磁曲线：Discrete-timeIDPOLYmodel:A(q)y(t)=B(q)u(t)+C(q)e(t)A(q)=1-1.921(+-
5、0.138)q^-1+0.9663(+-0.107)q^-2B(q)=-0.009032(+-0.02145)+0.01512(+-0.02722)q^-1C(q)=1-1.01(+-0.2539)q^-1EstimatedusingARMAXfromdatasetz2Lossfunction0.00179829andFPE0.00254758Samplinginterval:1Matlab程序：closeallclearallclcloadrise.txtloaddown.txtu1=rise(:,1);f1=rise(:,2);u2=down(
6、:,1);f2=down(:,2);figure(1);set(figure(1),'name','采集得到的数据');plot(u1,f1,'g',u2,f2,'r')holdona=polyfit(u1,f1,4)f11=polyval(a,u1);plot(u1,f11,'b')b=polyfit(u2,f2,4)f21=polyval(b,u2);plot(u2,f21,'k')figure(2)set(figure(2),'name','拟合得到的曲线');plot(u1,f11,'b')holdonplot(u2,f21,'r')z1=
7、[f1,u1];z2=[f2,u2];m=2;switchmcase1m1=n4sid(z1,2);m2=n4sid(z2,2);case2na=2;nb=2;nc=1;nk=0;m1=armax(z1,[na,nb,nc,nk])m2=armax(z2,[na,nb,nc,nk])case3m1=arx(z1,[221]);m2=arx(z2,[221]);case4m1=pem(z1);m2=pem(z2);endfigure(3);set(figure(3),'name','升磁过程辨识结果比较');compare(z1,m1)figure(
8、4);set(figure(4),'name','降磁过程辨识结果比较');compare(z2,m2)figure(5)
