svdnew.m

PMSM直接转矩控制系统的空间辨识,辨识出系统介次,和状态空间模型,从而调节控制器参数

function [u1,s1,v1]=svdnew(x)
[uu,ss,vv]=svd(x);
u=[-1,1,-1,1,-1,1,-1,-1,1,-1,1,1,-1,1,1,1,-1,1,-1,-1,1,1,1,1,-1,-1,1,1,1,-1,-1,-1,1,-1,1,-1,1,1,1,-1,1,1,1,-1,-1,1,1,-1,1,1,1,1,-1,1,1,1,1,1,-1,-1,-1,1,1,-1,1,-1,1,1,-1,-1,1,1,-1,-1,-1,-1,-1,-1,1,1,-1,1,1,-1,-1,1,-1,1,1,1,-1,-1,-1,-1,1,-1,1,1,1,1,1,1,1,1,-1,1,-1,1,1,1,1,-1,-1,-1,-1,-1,1,-1,-1,1,1,1,-1,1,1,-1,1,1,-1,1,-1,1,-1,1,1,-1,1,-1,-1,1,-1,-1,-1,1,-1,1,1,-1,-1,-1,-1,1,1,-1,-1,-1,1,1,1,-1,1,-1,1,-1,-1,-1,1,-1,-1,-1,1,1,-1,-1,1,-1,-1,-1,-1,1,-1,-1,-1,-1,-1,1,1,1,-1,-1,1,-1,1,-1,-1,1,1,-1,-1,1,1,1,1,1,1,-1,-1,1,-1,-1,1,1,-1,1,-1,-1,-1,1,1,1,1,-1,1,-1,-1,-1,-1,-1,-1,-1,-1,1,-1,1,-1,-1,-1,-1,1,1,1,1,1,-1,1,1,-1,-1,-1,1,-1,-1,1,-1,
-1,1,-1,1,-1,1,-1,-1,1,-1,1,1,-1,1,1,1,-1,1,-1,-1,1,1,1,1,-1,-1,1,1,1,-1,-1,-1,1,-1,1,-1,1,1,1,-1,1,1,1,-1,-1,1,1,-1,1,1,1,1,-1,1,1,1,1,1,-1,-1,-1,1,1,-1,1,-1,1,1,-1,-1,1,1,-1,-1,-1,-1,-1,-1,1,1,-1,1,1,-1,-1,1,-1,1,1,1,-1,-1,-1,-1,1,-1,1,1,1,1,1,1,1,1,-1,1,-1,1,1,1,1,-1,-1,-1,-1,-1,1,-1,-1,1,1,1,-1,1,1,-1,1,1,-1,1,-1,1,-1,1,1,-1,1,-1,-1,1,-1,-1,-1,1,-1,1,1,-1,-1,-1,-1,1,1,-1,-1,-1,1,1,1,-1,1,-1,1,-1,-1,-1,1,-1,-1,-1,1,1,-1,-1,1,-1,-1,-1,-1,1,-1,-1,-1,-1,-1,1,1,1,-1,-1,1,-1,1,-1,-1,1,1,-1,-1,1,1,1,1,1,1,-1,-1,1,-1,-1,1,1,-1,1,-1,-1,-1,1,1,1,1,-1,1,-1,-1,-1,-1,-1,-1,-1,-1,1,-1,1,-1,-1,-1,-1,1,1,1,1,1,-1,1,1,-1,-1,-1,1,-1,-1,1,-1];

y=[16.37732,16.29774,16.21455,16.13498,16.04456,15.95776,15.8818,15.79499,15.71542,15.64308,15.5599,15.46947,15.39352,15.30671,15.22714,15.14757,15.06438,14.97396,14.91247,14.84737,14.78226,14.70992,14.63397,14.54716,14.50738,14.4893,14.48568,14.47483,14.44227,14.42419,14.46036,14.4893,14.51099,14.50376,14.46759,14.44951,14.41334,14.40249,14.36994,14.36632,14.41334,14.46759,14.50738,14.56163,14.68822,14.81481,14.90162,15.01374,15.14395,15.24884,15.37905,15.53819,15.69734,15.84563,15.97584,16.09881,16.26157,16.42072,16.57624,16.72092,16.85836,16.96325,17.03559,17.11155,17.22729,17.33579,17.45153,17.53472,17.65408,17.79876,17.93258,18.03386,18.13874,18.28704,18.41363,18.50767,18.57277,18.59086,18.60171,18.60894,18.59809,18.56554,18.55469,18.52575,18.54384,18.57277,18.57639,18.56192,18.54384,18.51852,18.50767,18.55469,18.60894,18.63426,18.63064,18.59809,18.52937,18.50043,18.43895,18.41001,18.41001,18.4498,18.51129,18.63788,18.7934,18.93446,19.06105,19.18041,19.27807,19.39381,19.55295,19.70124,19.84592,19.96889,20.05932,20.10995,20.14612,20.15336,20.15336,20.17144,20.15697,20.11357,20.08102,20.0991,20.13165,20.14612,20.18591,20.22931,20.27271,20.3342,20.42462,20.52951,20.63802,20.72844,20.81525,20.90929,21.04311,21.15886,21.26736,21.38672,21.46991,21.56756,21.65075,21.71224,21.73032,21.75926,21.77734,21.81351,21.8533,21.90755,21.92564,21.92202,21.9184,21.86415,21.81713,21.79543,21.74841,21.68692,21.5965,21.52054,21.48799,21.46991,21.43012,21.39757,21.35055,21.32885,21.27821,21.23481,21.12992,21.03588,20.94907,20.8695,20.76461,20.67419,20.591,20.50058,20.41739,20.3125,20.22569,20.13889,20.05932,19.96889,19.864,19.77358,19.69763,19.6072,19.51678,19.42636,19.3504,19.24913,19.16233,19.07914,18.94531,18.89106,18.84404,18.71745,18.5366,18.5583,18.47511,18.40278,18.28704,18.20023,18.10981,18.023,17.95067,17.90364,17.8928,17.92173,17.92173,17.90726,17.90364,17.88556,17.83492,17.81322,17.78429,17.7662,17.73365,17.70472,17.63238,17.56366,17.46239,17.4226,17.4226,17.44068,17.45153,17.46962,17.46962,17.44068,17.37196,17.29601,17.20558,17.12963,17.02836,16.95602,16.87283,16.78964,16.70284,16.62326,16.54008,16.44604,16.36646,16.28328,16.19647,16.10605,16.05903,16.04094,15.99754,15.99392,16.00477,15.99754,15.95414,15.89627,15.85286,15.78053,15.69734,15.61777,15.53819,15.46586,15.38267,15.3031,15.21629,15.14034,15.05353,14.97396,14.89077,14.8112,14.72801,14.64844,14.56163,14.48568,14.39887,14.35909,14.341,14.2976,14.29036,14.23973,14.18909,14.14931,14.15292,14.19633,14.2542,14.2976,14.31207,14.35909,14.42419,14.50376,14.56525,14.61589,14.64844,14.68822,14.72439,14.76056,14.77865,14.83652,14.91247,14.99566,15.10417,15.23799,15.39714,15.51649,15.62862,15.73712,15.88542,16.01563,16.14222,16.30136,16.4605,16.6088,16.739,16.89091,17.02112,17.19112,17.34664,17.49494,17.62514,17.71195,17.7662,17.83492,17.93982,18.03386,18.12428,18.21108,18.30512,18.43171,18.55469,18.64511,18.73553,18.86574,18.98148,19.06829,19.11531,19.13701,19.14424,19.13701,19.10807,19.06829,19.06105,19.01765,19.0068,19.03573,19.02489,18.99233,18.97063,18.94531,18.95978,18.99595,19.03212,19.06105,19.06829,19.04659,18.97786,18.92361,18.85127,18.80064,18.7934,18.81872,18.87297,18.95255,19.10446,19.24913,19.37572,19.53487,19.65061,19.79528,19.92188,20.05932,20.21123,20.3342,20.42101,20.47164,20.50058,20.50058,20.5259,20.5259,20.51143,20.44271,20.39931,20.39931,20.42101,20.42462,20.46079,20.52951,20.56207,20.62717,20.69589,20.761,20.84418,20.92376,20.97078,21.0612,21.18417,21.2746,21.38672,21.46629,21.53139,21.57841,21.65437,21.67245,21.67607,21.69054,21.70862,21.70862,21.72671,21.75926,21.77734,21.76649,21.70501,21.6182,21.56033,21.51331,21.45906,21.38672,21.28545,21.18779,21.1263,21.08652,21.0395,21.00333,20.97078,20.93099,20.87312,20.8044,20.71398,20.61994,20.52951,20.44271,20.34867,20.27271,20.17144,20.0774,19.9906,19.89656,19.80614,19.71933,19.64337,19.55295,19.45168,19.36849,19.27807,19.19488,19.09722,19.01403,18.91638,18.84404,18.74638,18.67043,18.58362,18.4932,18.39916,18.30874,18.22917,18.13874,18.07002,17.97598,17.88918,17.79876,17.71195,17.63238,17.53472,17.4877,17.46962,17.45877,17.44068,17.4009,17.36835,17.31409,17.22729,17.18027,17.15495,17.10793,17.07538,17.02836,16.9741,16.8656,16.78602,16.74624,16.73539,16.76432,16.76432,16.78964,16.78241,16.739,16.68114,16.59433,16.51114,16.43157,16.34838,16.26157,16.17477,16.08796,16.01563,15.92159,15.84201,15.76968,15.67925,15.59606,15.50926,15.42969,15.38267,15.35735,15.32118,15.29224,15.29224,15.27416,15.23076,15.16204,15.10778,15.02821,
14.92332,14.84013,14.76418,14.69184,14.61227,14.5327,14.44951,14.36632,14.29036,14.20356,14.1276,14.05165,13.96484,13.88889,13.79123,13.73336,13.65379,13.5706,13.51997,13.44401,13.37167,13.29572,13.21976,13.13296,13.06062,13.02807,13.0136,12.97381,12.92318,12.88701,12.88701,12.90509,12.90148,12.86531,12.80744,12.76042,12.68808,12.64468,12.60851,12.55787,12.57957,12.60851,12.60489,12.65914,12.72063,12.79297,12.82552,12.90148,12.97743,13.07509,13.17636,13.30295,13.47295,13.6357,13.74783,13.91782,14.07697,14.22164,14.38079,14.54716,14.67376,14.77503,14.84737,14.91247,15.03545,15.13672,15.23799,15.34288,15.44415,15.59968,15.72627,15.83478,15.93967,16.0952,16.20732,16.31221,16.3737,16.40263,16.42433,16.43157,16.42072,16.3954,16.40625,16.3954,16.42072,16.45689,16.4605,16.44604,16.44242,16.41348,16.42072,16.46412,16.50391,16.53284,16.55093,16.51114,16.4605,16.43157,16.38455,16.35923,16.3737,16.43157,16.47859,16.61965,16.76432,16.92708,17.06091,17.19835,17.31047,17.45515,17.59983,17.74812,17.90726,18.03024,18.12066,18.16406,18.20385,18.21108,18.22917,18.25087,18.25448,18.2364,18.22193,18.25448,18.30874,18.35938,18.43895,18.54022,18.65596,18.80064,18.92723,19.0719,19.1985,19.307,19.45891,19.59636,19.74826,19.87847,20.0123,20.13527,20.22569,20.33058,20.41739,20.48249,20.51143,20.5476,20.58739,20.61994,20.68866,20.75376,20.7827,20.78631,20.78631,20.7827,20.76461,20.74653,20.72483,20.65611,20.5693,20.49696,20.47888,20.46441,20.44271,20.44271,20.41739,20.40654,20.37037,20.3125,20.24378,20.15697,20.06655,19.97613,19.89294,19.80252,19.71933,19.63252,19.55295,19.45168,19.36849,19.27807,19.19126,19.10446,19.01765,18.92361,18.84404,18.75723,18.66319,18.57639,18.48235,18.39916,18.30874,18.2364,18.14598,18.06641,17.96875,17.87833,17.79514,17.70833,17.636,17.53472,17.45153,17.36473,17.27792,17.1875,17.10069,17.05729,17.02836,17.02836,17.00665,16.97049,16.9307,16.89453,16.81134,16.75347,16.7173,16.6739,16.63773,16.57624,16.51476,16.42433,16.34114,16.27604,16.25796,16.24349,16.21094,16.20732,16.17477,16.11328,16.04456,15.95414,15.87095,15.78776,15.71542,15.62862,15.55266,15.46947,15.39714,15.30671,15.21629,15.13672,15.04991,14.97758,14.88715,14.8112,14.76418,14.73163,14.70269,14.68822,14.68461,14.65929,14.61227,14.56163,14.49653,14.43142,14.35185,14.26866,14.19994,14.1059,14.03356,13.95399,13.8708,13.78761,13.70805,13.6357,13.614,13.49103,13.40061,13.31742,13.2146,13.15828,1,13.00275,12.95211,12.88339,12.82552,12.7568,12.67723,12.59766,12.55787,12.53979,12.55787,12.55787,12.52532,12.51447,12.53617,12.58681,12.62297,12.61212,12.59404,12.57957,12.57234,12.57234,12.53979,12.57234,12.61574,12.65914,12.72063,12.8219,12.97381,13.12211,13.20168,13.34997,13.49826,13.63932,13.79123,13.94676,14.1059,14.26143,14.3844,14.54355,14.70631,14.84013,14.99204,15.16204,15.28501,15.39714,15.48032,15.53819,15.66117,15.76968,15.8818,15.96137,16.10243,16.23987,16.36646,16.48582,16.59433,16.74624,16.87645,16.9741,17.02474,17.05367,17.07538,17.08261,17.06453,17.02112,17.01027,16.97049,16.9741,16.99219,16.98495,16.96687,16.9524,16.91623,16.9307,16.97049,17.01751,17.04282,17.05729,17.01751,16.96325,16.909,16.8656,16.82943,16.82219,16.87283,16.91985,17.04644,17.19112,17.35388,17.48047,17.63238,17.76259,17.86748,18.03747,18.193,18.33767,18.4715,18.56554,18.61979,18.65234,18.66681,18.67405,18.68489,18.68851,18.63426,18.61617,18.63426,18.68128,18.71745,18.7717,18.80787,18.84404,18.90914,18.99595,19.08275,19.15509,19.23105,19.29977,19.38657,19.51317,19.61444,19.71571,19.80252,19.87486,19.93996,20.00145,20.03038,20.03761,20.04123,20.05932,20.03038,20.01953,20.04123,20.01953,19.98336,19.90017,19.80975,19.71571,19.66869,19.6072,19.5204,19.43359,19.3504,19.25998,19.21296,19.16956,19.13339,19.08275,19.0502,18.99957,18.93446,18.84404,18.7717,18.67043,18.58001,18.4932,18.40278,18.31959,18.22917,18.1496,18.08449,17.97598,17.8928,17.80237,17.7228,17.62514,17.54919,17.46962,17.37196,17.29239,17.2092,17.11878,17.03559,16.9524,16.8656,16.78602,16.69922,16.61965,16.53284,16.44604,16.35923,16.28328,16.2037,16.10605,16.02648,15.93967,15.84925,15.77691,15.73351,15.74074,15.76244,15.75521,15.71904,15.70457,15.67925,15.625,15.5816,15.57436,15.56351,15.55266,15.52011,15.46586,15.38267,15.3031,15.27055,15.27416,15.3031,15.33565,15.36097,15.37182,15.34288,15.28863,15.21629,15.14034,15.0463,14.95949,14.8763,14.80035,14.72801,14.65929,14.5761,14.4893,14.40249,14.32292,14.23973,14.15292,14.08058,14.04442,14.02633,14.00825,14.00825,14.03356,14.0408,14.00825,13.96123,13.92867,13.86357];

p=1;
m=1;
r=18;
t=10;
N=200;
s=r;
s11=0.5*s;
s22=0.5*s;

G1=getg(u,y,p,m,r,t,N,s,s11,s22);
n=getn(G1);
for i=1:n
   for j=1:n
        s1(i,j)=ss(i,j);
    end
end

for i=1:p*r
    for j=1:n
        u1(i,j)=uu(i,j);
    end 
end 
[c,d]=size(x);
for i=1:n
    for j=1:d
        v1(i,j)=vv(i,j);
    end 
end