0001 function [sys, VERSION] = model_NSF_TM_SC
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023
0024
0025
0026
0027 VERSION ='$Id:$';
0028 utils.helper.msg(utils.const.msg.MNAME, ['running ', mfilename]);
0029
0030
0031 sys =ssm();
0032
0033
0034 sys.paramnames = ...
0035 {'p_p1_N_TM_SC_x' 'p_p2_N_TM_SC_x' 'p_p3_N_TM_SC_x' 'p_p4_N_TM_SC_x' 'p_p5_N_TM_SC_x' ...
0036 'p_g1_N_TM_SC_x' 'p_g2_N_TM_SC_x' 'p_g3_N_TM_SC_x' 'p_g4_N_TM_SC_x' 'p_g5_N_TM_SC_x' 'p_g6_N_TM_SC_x' ...
0037 'p_p1_N_TM_SC_y' 'p_p2_N_TM_SC_y' 'p_p3_N_TM_SC_y' 'p_p4_N_TM_SC_y' 'p_p5_N_TM_SC_y' ...
0038 'p_g1_N_TM_SC_y' 'p_g2_N_TM_SC_y' 'p_g3_N_TM_SC_y' 'p_g4_N_TM_SC_y' 'p_g5_N_TM_SC_y' 'p_g6_N_TM_SC_y' ...
0039 'p_p1_N_TM_SC_z' 'p_p2_N_TM_SC_z' 'p_p3_N_TM_SC_z' 'p_p4_N_TM_SC_z' 'p_p5_N_TM_SC_z' ...
0040 'p_g1_N_TM_SC_z' 'p_g2_N_TM_SC_z' 'p_g3_N_TM_SC_z' 'p_g4_N_TM_SC_z' 'p_g5_N_TM_SC_z' 'p_g6_N_TM_SC_z' ...
0041 ...
0042 'p_p1_N_TM_SC_t' 'p_p2_N_TM_SC_t' 'p_p3_N_TM_SC_t' 'p_p4_N_TM_SC_t' ...
0043 'p_g1_N_TM_SC_t' 'p_g2_N_TM_SC_t' 'p_g3_N_TM_SC_t' 'p_g4_N_TM_SC_t' 'p_g5_N_TM_SC_t' ...
0044 'p_p1_N_TM_SC_e' 'p_p2_N_TM_SC_e' 'p_p3_N_TM_SC_e' 'p_p4_N_TM_SC_e' ...
0045 'p_g1_N_TM_SC_e' 'p_g2_N_TM_SC_e' 'p_g3_N_TM_SC_e' 'p_g4_N_TM_SC_e' 'p_g5_N_TM_SC_e' ...
0046 'p_p1_N_TM_SC_p' 'p_p2_N_TM_SC_p' 'p_p3_N_TM_SC_p' 'p_p4_N_TM_SC_p' ...
0047 'p_g1_N_TM_SC_p' 'p_g2_N_TM_SC_p' 'p_g3_N_TM_SC_p' 'p_g4_N_TM_SC_p' 'p_g5_N_TM_SC_p' ...
0048 ...
0049 'p_gDC_N_TM_SC_x' 'p_gDC_N_TM_SC_y' 'p_gDC_N_TM_SC_z' 'p_gDC_N_TM_SC_t' 'p_gDC_N_TM_SC_e' 'p_gDC_N_TM_SC_p' ...
0050 };
0051 sys.paramvalues = ...
0052 [0 0 0 0 0 ...
0053 0 0 0 0 0 0 ...
0054 0 0 0 0 0 ...
0055 0 0 0 0 0 0 ...
0056 0 0 0 0 0 ...
0057 0 0 0 0 0 0 ...
0058 ...
0059 0 0 0 0 ...
0060 0 0 0 0 0 ...
0061 0 0 0 0 ...
0062 0 0 0 0 0 ...
0063 0 0 0 0 ...
0064 0 0 0 0 0 ...
0065 ...
0066 0 0 0 0 0 0 ...
0067 ];
0068 sys.paramsigmas = ...
0069 [0 0 0 0 0 ...
0070 0 0 0 0 0 0 ...
0071 0 0 0 0 0 ...
0072 0 0 0 0 0 0 ...
0073 0 0 0 0 0 ...
0074 0 0 0 0 0 0 ...
0075 ...
0076 0 0 0 0 ...
0077 0 0 0 0 0 ...
0078 0 0 0 0 ...
0079 0 0 0 0 0 ...
0080 0 0 0 0 ...
0081 0 0 0 0 0 ...
0082 ...
0083 0 0 0 0 0 0 ...
0084 ];
0085
0086
0087
0088 for i_params=1:length(sys.paramnames)
0089 cmd = [sys.paramnames{i_params}, '=sym(''', sys.paramnames{i_params},''');'];
0090 eval(cmd);
0091 end
0092
0093
0094
0095
0096 f = 0.9*1/sqrt(2);
0097 f1 = f* 25/20.35*1.8164e-14;
0098 A1 =[-5e-5+p_p1_N_TM_SC_x 1 0 0 0 ;...
0099 0 -5e-5+p_p2_N_TM_SC_x 1 0 0 ;...
0100 0 0 -5e-5+p_p3_N_TM_SC_x 0 0 ;...
0101 0 0 0 -4.5e-5+p_p4_N_TM_SC_x 0 ;...
0102 0 0 0 0 -1e-4+p_p5_N_TM_SC_x];
0103 B1 =[ 0; 0; 1; 1; 1];
0104 C1 =f1*[3.409875000e-9+p_g1_N_TM_SC_x, 0.0004114874890+p_g2_N_TM_SC_x, 91.66094784+p_g3_N_TM_SC_x, -90.69787011+p_g4_N_TM_SC_x, -0.9584727300+p_g5_N_TM_SC_x];
0105 D1 =f1*(1+p_g6_N_TM_SC_x) ;
0106
0107 f2 = f* 25/20.35*1.8106e-14;
0108 A2 =[-5e-5+p_p1_N_TM_SC_y 1 0 0 0 ;...
0109 0 -5e-5+p_p2_N_TM_SC_y 1 0 0 ;...
0110 0 0 -5e-5+p_p3_N_TM_SC_y 0 0 ;...
0111 0 0 0 -4.5e-5+p_p4_N_TM_SC_y 0 ;...
0112 0 0 0 0 -1e-4+p_p5_N_TM_SC_y];
0113 B2 =[ 0; 0; 1; 1; 1];
0114 C2 =f2*[3.968625000e-9+p_g1_N_TM_SC_y, 0.0004786624900+p_g2_N_TM_SC_y, 106.6504479+p_g3_N_TM_SC_y, -105.5263929+p_g4_N_TM_SC_y, -1.119200000+p_g5_N_TM_SC_y];
0115 D2 =f2*(1+p_g6_N_TM_SC_y) ;
0116
0117 f3 = f* 25/20.35*1.8262e-13;
0118 A3 =[-5e-5+p_p1_N_TM_SC_z 1 0 0 0 ;...
0119 0 -5e-5+p_p2_N_TM_SC_z 1 0 0 ;...
0120 0 0 -5e-5+p_p3_N_TM_SC_z 0 0 ;...
0121 0 0 0 -4.5e-5+p_p4_N_TM_SC_z 0 ;...
0122 0 0 0 0 -1e-4+p_p5_N_TM_SC_z];
0123 B3 =[ 0; 0; 1; 1; 1];
0124 C3 =f3*[2.322875000e-9+p_g1_N_TM_SC_z, 0.0002807874940+p_g2_N_TM_SC_z, 62.49814869+p_g3_N_TM_SC_z, -61.84802733+p_g4_N_TM_SC_z, -0.6460363600+p_g5_N_TM_SC_z];
0125 D3 =f3*(1+p_g6_N_TM_SC_z) ;
0126
0127 f4 = f* 25/20.35*1.2985e-11;
0128 A4 =[-9.5e-5+p_p1_N_TM_SC_t 1 0 0 ;...
0129 0 -9.5e-5+p_p2_N_TM_SC_t 0 0 ;...
0130 0 0 -9e-5+p_p3_N_TM_SC_t 0 ;...
0131 0 0 0 -1e-4+p_p4_N_TM_SC_t];
0132 B4 =[ 0; 1; 1; 1 ];
0133 C4 =f4*[-0.001894924902+p_g1_N_TM_SC_t, 38.36974950+p_g2_N_TM_SC_t, 170.7544552+p_g3_N_TM_SC_t, -209.1199847+p_g4_N_TM_SC_t ];
0134 D4 =f4*(1+p_g5_N_TM_SC_t) ;
0135
0136 f5 = f* 25/20.35*1.2932e-11;
0137 A5 =[-9.5e-5+p_p1_N_TM_SC_e 1 0 0 ;...
0138 0 -9.5e-5+p_p2_N_TM_SC_e 1 0 ;...
0139 0 0 -9.5e-5+p_p3_N_TM_SC_e 0 ;...
0140 0 0 0 -1e-4+p_p4_N_TM_SC_e];
0141 B5 =[ 0; 0; 1; 1 ];
0142 C5 =f5*[8.897160875e-9+p_g1_N_TM_SC_e, -0.001959277150+p_g2_N_TM_SC_e, 392.6841350+p_g3_N_TM_SC_e, -392.6800000+p_g4_N_TM_SC_e ];
0143 D5 =f5*(1+p_g5_N_TM_SC_e) ;
0144
0145 f6 = f* 25/20.35*2.5628e-12;
0146 A6 =[-9.5e-5+p_p1_N_TM_SC_p 1 0 0 ;...
0147 0 -9.5e-5+p_p2_N_TM_SC_p 1 0 ;...
0148 0 0 -9.5e-5+p_p3_N_TM_SC_p 0 ;...
0149 0 0 0 -1e-4+p_p4_N_TM_SC_p];
0150 B6 =[ 0; 0; 1; 1 ];
0151 C6 =f6*[1.560165288e-8+p_g1_N_TM_SC_p, -0.003439734350+p_g2_N_TM_SC_p, 689.4850150+p_g3_N_TM_SC_p, -689.4800000+p_g4_N_TM_SC_p ];
0152 D6 =f6*(1+p_g5_N_TM_SC_p) ;
0153
0154 A =blkdiag(A1, A2, A3, A4, A5, A6, A1, A2, A3, A4, A5, A6);
0155 B_TM_SC_noise =blkdiag(B1, B2, B3, B4, B5, B6, B1, B2, B3, B4, B5, B6);
0156 C =blkdiag(C1, C2, C3, C4, C5, C6, C1, C2, C3, C4, C5, C6);
0157 D_TM_SC_noise =blkdiag(D1, D2, D3, D4, D5, D6, D1, D2, D3, D4, D5, D6);
0158
0159
0160 B_TM_SC_DC =zeros(54,12);
0161 D_TM_SC_DC =eye(12) + diag([p_gDC_N_TM_SC_x, p_gDC_N_TM_SC_y, p_gDC_N_TM_SC_z, p_gDC_N_TM_SC_t, p_gDC_N_TM_SC_e, p_gDC_N_TM_SC_p, ...
0162 p_gDC_N_TM_SC_x, p_gDC_N_TM_SC_y, p_gDC_N_TM_SC_z, p_gDC_N_TM_SC_t, p_gDC_N_TM_SC_e, p_gDC_N_TM_SC_p ]);
0163
0164
0165
0166
0167
0168
0169
0170 sys.name = 'NSF_TM_SC';
0171
0172 sys.amats = {A};
0173 sys.bmats = {B_TM_SC_noise B_TM_SC_DC};
0174 sys.cmats = {C};
0175 sys.dmats = {D_TM_SC_noise D_TM_SC_DC};
0176 sys.timestep = 0;
0177
0178 sys.ssnames = {'NSF_TM_SC'};
0179
0180 sys.inputnames = {'Noise_F_TM_SC' 'DC_F_TM_SC'};
0181 sys.inputvarnames = ...
0182 {{'n_fx1_TM_SC' 'n_fy1_TM_SC' 'n_fz1_TM_SC' 'n_tx1_TM_SC' 'n_ty1_TM_SC' 'n_tz1_TM_SC' ...
0183 'n_fx2_TM_SC' 'n_fy2_TM_SC' 'n_fz2_TM_SC' 'n_tx2_TM_SC' 'n_ty2_TM_SC' 'n_tz2_TM_SC'}...
0184 {'DC_fx1_TM_SC' 'DC_fy1_TM_SC' 'DC_fz1_TM_SC' 'DC_tx1_TM_SC' 'DC_ty1_TM_SC' 'DC_tz1_TM_SC' ...
0185 'DC_fx2_TM_SC' 'DC_fy2_TM_SC' 'DC_fz2_TM_SC' 'DC_tx2_TM_SC' 'DC_ty2_TM_SC' 'DC_tz2_TM_SC'}...
0186 };
0187
0188 sys.outputnames = {'Dist_F_TM_SC'};
0189 sys.outputvarnames = ...
0190 {{'Dist_fx1_TM_SC' 'Dist_fy1_TM_SC' 'Dist_fz1_TM_SC' 'Dist_tx1_TM_SC' 'Dist_ty1_TM_SC' 'Dist_tz1_TM_SC' ...
0191 'Dist_fx2_TM_SC' 'Dist_fy2_TM_SC' 'Dist_fz2_TM_SC' 'Dist_tx2_TM_SC' 'Dist_ty2_TM_SC' 'Dist_tz2_TM_SC'}};
0192
0193 end