wang1.m

一些关于系统辩士的matlab程序

  -1.7840062e+00
  -1.8738079e+00
  -1.6435210e+00
  -1.9205835e+00
  -1.8983514e+00
  -2.0115416e+00
  -2.0274915e+00
  -2.1844445e+00
  -2.6165120e+00
  -2.3547626e+00
  -2.2840804e+00
  -2.4405111e+00
  -3.1366299e+00
  -2.7885653e+00
  -2.8688423e+00
  -2.8088904e+00
  -3.0407631e+00
  -3.0028091e+00
  -2.7428312e+00
  -2.9542300e+00
  -3.1276963e+00
  -3.4215410e+00
  -3.3044819e+00
  -3.3123991e+00
  -3.5030936e+00
  -3.5068629e+00
  -3.4025207e+00
  -3.2182942e+00
  -4.0613031e+00
  -3.8330139e+00
  -3.4281714e+00
  -3.5741378e+00
  -3.7810054e+00
  -3.6787014e+00
  -3.9904974e+00
  -4.0312600e+00
  -3.7941451e+00
  -3.5839069e+00
  -3.7654623e+00
  -4.0439278e+00
  -3.7912284e+00
  -4.0427111e+00
  -4.1620908e+00
  -4.1285820e+00
  -4.1839122e+00
  -4.1239757e+00
  -3.9896314e+00
  -4.3338925e+00
  -4.0690042e+00
  -3.9301366e+00
  -4.1433849e+00
  -4.0486658e+00
  -4.2149817e+00
  -3.8153278e+00
  -4.0480738e+00
  -4.6089852e+00
  -4.2518282e+00
  -4.6324850e+00
  -4.0768172e+00
  -4.2145554e+00
  -4.1761131e+00
  -4.1391910e+00
  -3.9731339e+00
  -4.3967663e+00
  -4.3689611e+00
  -4.2266453e+00
  -4.1621576e+00
  -4.5883764e+00
  -4.4825311e+00
  -4.3714347e+00
  -4.8825451e+00
  -4.1605049e+00
  -4.6034589e+00
  -4.1731938e+00
  -4.5666045e+00
  -4.2875192e+00
  -4.2553240e+00
  -4.4593045e+00
  -4.3160767e+00
  -4.4583212e+00
  -4.3260581e+00
  -4.6602291e+00
  -4.5907050e+00
  -4.1670080e+00
  -4.5753402e+00
  -4.4070516e+00
  -4.7748831e+00
  -4.5406984e+00
  -4.6172899e+00
  -4.5183061e+00
  -4.2729654e+00
  -4.2750742e+00
  -4.2388931e+00
  -4.4431868e+00
  -4.3771932e+00
  -4.6901222e+00
  -4.5710519e+00
  -4.2285942e+00
  -4.3046126e+00
  -4.3982107e+00
  -4.7369437e+00
  -4.4239909e+00
  -4.4122896e+00
  -4.5726440e+00
  -4.4344024e+00
  -3.9808718e+00
  -4.5867479e+00
  -4.8184261e+00
  -5.0311892e+00
  -4.6501977e+00
  -4.6341291e+00
  -4.7344039e+00
  -4.0943190e+00
  -4.3080065e+00
  -4.3684031e+00
  -4.4191541e+00
  -5.0784201e+00
  -4.1673480e+00
  -4.4995591e+00
  -4.6861409e+00
  -4.9169481e+00
  -4.4890137e+00
  -4.0899194e+00
  -4.4347988e+00
  -4.3390610e+00
  -4.5456268e+00
  -4.3029920e+00
  -4.8150547e+00
  -4.3169182e+00
  -4.4050309e+00
  -4.5474941e+00
  -4.4711698e+00
  -4.1931483e+00
  -4.5050874e+00
  -4.3339708e+00
  -4.7718981e+00
  -4.2758041e+00
  -4.1490631e+00
  -4.6296448e+00
  -4.7151349e+00
  -4.7053049e+00
  -4.3751684e+00
  -4.6653093e+00
  -4.5331563e+00
  -4.7953193e+00
  -4.3079263e+00
  -4.6662151e+00
  -4.3395587e+00
  -4.3622955e+00
  -4.5478357e+00
  -4.4545931e+00
  -4.5793856e+00
  -4.5886861e+00
  -4.2855309e+00
  -4.5069249e+00
  -4.7349168e+00
  -4.2020533e+00
  -4.7251694e+00
  -4.4310223e+00
  -4.6845282e+00
  -4.5780457e+00
  -4.4362284e+00
  -4.3059005e+00
  -4.5840807e+00
  -4.4487229e+00
  -4.5603121e+00
  -4.5127649e+00
  -4.7546192e+00
  -4.5425257e+00
  -4.2187387e+00
  -4.6836525e+00
  -4.7765272e+00
  -4.2015448e+00
  -4.2832016e+00
  -4.3591602e+00
  -4.8088085e+00
  -4.2056206e+00
  -4.6360521e+00
  -4.7586884e+00
  -4.6571369e+00
  -4.4351441e+00
  -4.3353324e+00
  -4.3554408e+00
  -4.3607034e+00
  -3.8364583e+00
  -4.1969950e+00
  -3.9028654e+00
  -3.9121475e+00
  -3.5170036e+00
  -3.8112111e+00
  -3.1095031e+00
  -3.2728568e+00
  -2.8792903e+00
  -2.8737436e+00
  -2.7103977e+00
  -2.6401986e+00
  -2.2563457e+00
  -1.8188415e+00
  -1.5632302e+00
  -1.1547342e+00
  -1.1395372e+00
  -1.5050926e+00
  -1.0910093e+00
  -7.2662407e-01
  -1.0202397e+00
  -3.9784067e-01
  -2.7826012e-01
   7.8534037e-02
   2.9122158e-01
   5.2249192e-01
   4.6101347e-01
   6.0869412e-01
   1.0952164e+00
   1.3177961e+00
   1.0148722e+00
   1.2589352e+00
   1.1670662e+00
   1.4954873e+00
   1.7642891e+00
   1.6115416e+00
   1.8716055e+00
   1.9636091e+00
   1.9387463e+00
   2.0589137e+00
   2.6422523e+00
   2.1611821e+00
   2.3545491e+00
   2.2417379e+00
   2.6086979e+00
   2.7889940e+00
   2.6828759e+00
   2.9327515e+00
   2.7917040e+00
   3.2596092e+00
   3.0665101e+00
   3.2924637e+00
   3.6343485e+00
   3.3447045e+00
   3.5802468e+00
   3.4178732e+00
   3.3573632e+00
   3.2844637e+00
   3.4396055e+00
   3.7936616e+00
   3.7780366e+00
   3.4402281e+00
   3.6484621e+00
   3.4651451e+00
   3.8228362e+00
   3.9901437e+00
   3.6801754e+00
   3.3022551e+00
   4.0312151e+00
   4.5623050e+00
   4.0558815e+00
   4.0348827e+00
   4.0573219e+00
   4.2166252e+00
   4.4652967e+00
   4.3747532e+00
   4.0518020e+00
   4.1143846e+00
   4.1793615e+00
   3.9774798e+00
   4.3973773e+00
   4.2148437e+00
   4.5097490e+00
   4.1747681e+00
   4.5581766e+00
   3.8619303e+00
   4.6193205e+00
   3.9636506e+00
   4.0690536e+00
   4.3006301e+00
   4.0105380e+00
   4.4051321e+00
   4.4307601e+00
   4.5951230e+00
   4.2591045e+00
   4.4322193e+00
   4.4061149e+00
   4.5052025e+00
   4.4981371e+00
   4.1802961e+00
   4.7387960e+00
   4.3154014e+00
   4.2156017e+00
   4.2044442e+00
   4.6481906e+00
   4.3825095e+00
   4.9601908e+00
   4.3844333e+00
   4.4003725e+00
   4.5152103e+00
   4.2012414e+00
   4.3853978e+00
   4.1284543e+00
   4.6580027e+00
   4.5120989e+00
   4.2094613e+00
   4.4717279e+00
   4.3814036e+00
   4.7412269e+00
   4.7389752e+00
   4.6002193e+00
   4.5692285e+00
   4.6333742e+00
   4.3219255e+00
   4.5503106e+00
   4.3997918e+00
   4.6540836e+00
   4.4474670e+00
   4.8595438e+00
   4.5207431e+00
   4.7421293e+00
   4.4821553e+00
   4.1903315e+00
   4.5649780e+00
   4.4534136e+00
   4.2163827e+00
   4.5413163e+00
   4.5507575e+00
   4.4242638e+00
   4.7738946e+00
   4.3399152e+00
   4.3129478e+00
   4.5006405e+00
   4.5811349e+00
   4.4009308e+00
   4.4497270e+00
   4.2876358e+00
   4.9075368e+00
   4.8941308e+00
   4.2016030e+00
   4.8943281e+00
   4.3278245e+00
   4.4466632e+00
   4.6190228e+00
   4.4435767e+00
   4.7766205e+00
   4.7537276e+00
   4.1904622e+00
   4.2964826e+00
   4.3583973e+00
   4.5949140e+00
   4.4243602e+00
   4.6525880e+00
   4.5571968e+00
   4.3085007e+00
   4.5401313e+00
   4.3627887e+00
   4.3868465e+00
   4.7918310e+00
   4.1835390e+00
   4.1102819e+00
   4.6118143e+00
   4.2026582e+00
   4.4857047e+00
   4.5629800e+00
   4.4606561e+00
   4.4002476e+00
   4.4769847e+00
   4.7562839e+00
   4.3869294e+00
   4.6359298e+00
   4.3900410e+00
   4.6294924e+00
   4.3702689e+00
   4.7536892e+00
   4.5526936e+00
   4.8344215e+00
   4.5797802e+00
   4.6889677e+00
   4.7740451e+00
   4.6026119e+00
   4.4385773e+00
   4.4933719e+00
   4.3230807e+00
   4.1045566e+00
   4.5460626e+00
   4.6143246e+00
   4.5162540e+00
   4.3366638e+00
   4.2240618e+00
   4.4238893e+00
   4.1462637e+00
   4.2454286e+00
   3.9331598e+00
   3.8948664e+00
   3.1488792e+00
   3.1892210e+00
   3.2510162e+00
   3.2252252e+00
   2.3183373e+00
   2.4798201e+00
   2.6532480e+00
   2.0936932e+00
   1.7353119e+00
   1.3739493e+00
   1.4506562e+00
   1.3430154e+00
   1.1187083e+00
   8.4918247e-01
   7.0328399e-01
   6.4316143e-01
   5.1971886e-01
  -1.5162068e-01
   1.2052025e-01
  -3.0147388e-01
  -2.2958373e-01
  -6.2346013e-01
  -4.0932390e-01
  -6.0647495e-01
  -1.1113595e+00
  -1.2097869e+00
  -8.8969501e-01
  -1.1694061e+00
  -1.8738243e+00
  -1.5087019e+00
  -1.6517797e+00
  -2.0267079e+00
  -2.2115641e+00
  -1.9944317e+00
  -2.4059726e+00
  -2.1328746e+00
  -2.3755621e+00
  -2.8012190e+00
  -2.3813939e+00
  -2.8107919e+00
  -2.8124929e+00
  -2.4090865e+00
  -2.9145211e+00
  -2.9713041e+00
  -2.8102469e+00
  -2.8148803e+00
  -3.2429671e+00
  -3.3368753e+00
  -3.4291307e+00
  -3.2077998e+00
  -3.4132530e+00
  -3.3441272e+00
  -3.2247759e+00
  -3.2969245e+00
  -3.7326487e+00
  -3.7863968e+00
  -3.5074536e+00
  -3.5904946e+00
  -3.5627671e+00
  -3.6959563e+00
  -4.1137436e+00
  -4.0138097e+00
  -3.9843282e+00
  -3.9588160e+00
  -3.9048668e+00
  -3.9582027e+00
  -4.0981508e+00
  -3.9524514e+00
  -3.9996989e+00
  -4.0925084e+00
  -4.1480339e+00
  -4.3672718e+00
  -3.8956823e+00
  -3.8762080e+00
  -4.0509304e+00
  -4.3571024e+00
  -4.0921073e+00
  -4.4490981e+00
  -4.1621959e+00
  -4.1436762e+00
  -4.0332194e+00
  -4.0129965e+00
  -4.1111476e+00
  -4.5152523e+00
  -4.0903584e+00
  -4.5064728e+00
  -4.5231960e+00
  -3.7441318e+00
  -4.3846116e+00
  -4.4377148e+00
  -4.2646448e+00
  -4.5051972e+00
  -4.2498010e+00
  -4.4004881e+00
  -4.4101399e+00
  -4.5078027e+00
  -4.3717441e+00
  -4.6204701e+00
  -4.2122300e+00
  -4.5247412e+00
  -4.5156286e+00
  -4.6160148e+00
  -4.5115794e+00
  -4.1005094e+00
  -4.3268850e+00
  -4.8615304e+00
  -4.6110299e+00
  -4.5294853e+00
  -4.6541188e+00
  -4.5213695e+00
  -4.5894155e+00
  -4.2598391e+00
  -4.5761665e+00
  -4.7191693e+00
  -4.4098142e+00
  -4.7386481e+00
  -4.4523945e+00
  -4.1951887e+00
  -4.4447258e+00
  -4.3171511e+00
  -4.4388393e+00
  -4.5655256e+00
  -4.2392130e+00
  -4.8349669e+00
  -3.8971542e+00
  -4.4505600e+00
  -4.3557899e+00
  -4.5765930e+00
  -4.7212205e+00
  -4.8678072e+00
  -4.2127125e+00
  -4.7907355e+00
  -4.6583950e+00
  -4.5369147e+00
  -4.7240429e+00
  -4.4634371e+00
  -4.4943035e+00
  -4.6433472e+00
  -4.4105378e+00
  -4.7181180e+00
  -4.6401107e+00
  -4.3414086e+00
  -4.4773080e+00
  -4.4245672e+00
  -4.2481645e+00
  -4.4163959e+00
  -4.5309209e+00
  -4.4475187e+00
  -4.3793514e+00
  -4.8636198e+00
  -4.5491264e+00
  -4.7565469e+00
  -4.3944206e+00
  -4.3233200e+00
  -4.2841032e+00
  -4.6919912e+00
  -4.6141109e+00
  -4.2750346e+00
  -4.4715192e+00
  -4.5847653e+00
  -4.5479786e+00
  -4.4786214e+00
  -4.3533311e+00
  -4.6057674e+00
  -4.6759671e+00
  -4.4998048e+00
  -4.4766410e+00
  -4.5248562e+00
  -4.3789613e+00
  -4.0024824e+00
  -4.2787140e+00
  -4.7740045e+00
  -4.5964611e+00
  -4.8840294e+00
  -4.6945127e+00
  -4.2801625e+00
  -4.1925471e+00
  -4.0514049e+00
  -4.6780844e+00
  -4.3446364e+00
  -4.6301453e+00
  -4.2428413e+00
  -4.2858060e+00
  -4.5913603e+00
  -4.9917715e+00
  -4.5320754e+00
  -4.7605224e+00
  -4.7326928e+00
  -4.6040778e+00
  -4.8204755e+00
  -4.1836010e+00 ];
t=[0:0.1:60];
N=length(y);
for i=1:4
clear ftk ;clear f;
for k=i+1:N
f(:,k)=[-y((k-1):-1:(k-i));u((k-1):-1:(k-i))];
end
ftk=f';
clear F;clear Y;
 for k=i+1:N
F(k-1,:)=[ftk(k,:)];
Y(k-1,:)=[y(k)];
end
clear q1;clear q;
q1=inv(F'*F)*F';
q=q1*Y;
J(i)=0;
for k=i+1:N
 J(i)=J(i)+(y(k)-ftk(k,:)*q)^2;
 J(i)=J(i)/N;
end
Ym=ftk*q;
% figure(i)
 subplot(2,2,i)
plot(t,y,'m',t,Ym,'g')
grid on
axis([0 60 -5 +5])
hold on
%title('compare the real output with the model output')
hold on
subplot(2,2,i)
e=y-Ym
plot(t,e,'c')
grid on
title('error of the real output with the model output&compare the real output with the model output')
end