Request for Matlab text output

classic Classic list List threaded Threaded
3 messages Options
Reply | Threaded
Open this post in threaded view
|

Request for Matlab text output

Lukas Reichlin-4
Dear Octave community,

For the control package, I'm planning to increase Matlab compatibility by supporting additional properties for LTI models.
Now I wonder how Matlab displays the models defined in the following script file:

http://sourceforge.net/p/octave/control/ci/default/tree/devel/units.m

Could someone with access to a somewhat recent Matlab version (plus Control System Toolbox) please run the script 'units.m' from the link above and send me its output? (copy-pasted text, not a screenshot if possible)

Thanks in advance,
Lukas


Reply | Threaded
Open this post in threaded view
|

Re: Request for Matlab text output

nrjank
Lukas Reichlin-4 wrote
Dear Octave community,

For the control package, I'm planning to increase Matlab compatibility by supporting additional properties for LTI models.
Now I wonder how Matlab displays the models defined in the following script file:

http://sourceforge.net/p/octave/control/ci/default/tree/devel/units.m

Could someone with access to a somewhat recent Matlab version (plus Control System Toolbox) please run the script 'units.m' from the link above and send me its output? (copy-pasted text, not a screenshot if possible)

Thanks in advance,
Lukas
output from v2015a
workspace attached.

ran a ver, then your units.m script

>> ver
----------------------------------------------------------------------------------------------------
MATLAB Version: 8.5.0.197613 (R2015a)
MATLAB License Number: 186008
Operating System: Microsoft Windows Server 2012 Standard Version 6.2 (Build 9200)
Java Version: Java 1.7.0_60-b19 with Oracle Corporation Java HotSpot(TM) 64-Bit Server VM mixed mode
----------------------------------------------------------------------------------------------------
MATLAB                                                Version 8.5        (R2015a)
Simulink                                              Version 8.5        (R2015a)
Aerospace Blockset                                    Version 3.15       (R2015a)
Aerospace Toolbox                                     Version 2.15       (R2015a)
Bioinformatics Toolbox                                Version 4.5.1      (R2015a)
Communications System Toolbox                         Version 6.0        (R2015a)
Computer Vision System Toolbox                        Version 6.2        (R2015a)
Control System Toolbox                                Version 9.9        (R2015a)
Curve Fitting Toolbox                                 Version 3.5.1      (R2015a)
DSP System Toolbox                                    Version 9.0        (R2015a)
Data Acquisition Toolbox                              Version 3.7        (R2015a)
Database Toolbox                                      Version 5.2.1      (R2015a)
Datafeed Toolbox                                      Version 5.1        (R2015a)
Econometrics Toolbox                                  Version 3.2        (R2015a)
Financial Instruments Toolbox                         Version 2.1        (R2015a)
Financial Toolbox                                     Version 5.5        (R2015a)
Fixed-Point Designer                                  Version 5.0        (R2015a)
Fuzzy Logic Toolbox                                   Version 2.2.21     (R2015a)
Global Optimization Toolbox                           Version 3.3.1      (R2015a)
Image Acquisition Toolbox                             Version 4.9        (R2015a)
Image Processing Toolbox                              Version 9.2        (R2015a)
Instrument Control Toolbox                            Version 3.7        (R2015a)
MATLAB Coder                                          Version 2.8        (R2015a)
MATLAB Compiler                                       Version 6.0        (R2015a)
MATLAB Compiler SDK                                   Version 6.0        (R2015a)
MATLAB Report Generator                               Version 4.1        (R2015a)
Mapping Toolbox                                       Version 4.1        (R2015a)
Model Predictive Control Toolbox                      Version 5.0.1      (R2015a)
Neural Network Toolbox                                Version 8.3        (R2015a)
OPC Toolbox                                           Version 3.3.3      (R2015a)
Optimization Toolbox                                  Version 7.2        (R2015a)
Parallel Computing Toolbox                            Version 6.6        (R2015a)
Partial Differential Equation Toolbox                 Version 2.0        (R2015a)
Robust Control Toolbox                                Version 5.3        (R2015a)
Signal Processing Toolbox                             Version 7.0        (R2015a)
SimEvents                                             Version 4.4        (R2015a)
SimMechanics                                          Version 4.6        (R2015a)
Simscape                                              Version 3.13       (R2015a)
Simulink 3D Animation                                 Version 7.3        (R2015a)
Simulink Coder                                        Version 8.8        (R2015a)
Simulink Control Design                               Version 4.2        (R2015a)
Simulink Design Optimization                          Version 2.7        (R2015a)
Spreadsheet Link EX                                   Version 3.2.3      (R2015a)
Stateflow                                             Version 8.5        (R2015a)
Statistics and Machine Learning Toolbox               Version 10.0       (R2015a)
Symbolic Math Toolbox                                 Version 6.2        (R2015a)
System Identification Toolbox                         Version 9.2        (R2015a)
Wavelet Toolbox                                       Version 4.14.1     (R2015a)
>> units

sys =
 
  a =
               x1          x2          x3
   x1      0.3325     -0.5981     -0.1902
   x2     -0.2371      0.3206     -0.5838
   x3       0.581      0.1985  -0.0008821
 
  b =
           u1      u2
   x1   2.585       0
   x2       0       0
   x3  0.1873       0
 
  c =
            x1       x2       x3
   y1   -1.795  -0.5445   0.7394
   y2   0.8404        0    1.712
   y3   -0.888  -0.6003        0
   y4        0     0.49   -2.138
 
  d =
            u1       u2
   y1        0    0.124
   y2    1.355        0
   y3        0   -1.961
   y4        0  -0.1977
 
Sample time: unspecified
Discrete-time state-space model.


sys_ss =
 
  a =
                 power     voltage     current
   power        0.3325     -0.5981     -0.1902
   voltage     -0.2371      0.3206     -0.5838
   current       0.581      0.1985  -0.0008821
 
  b =
             angle  length
   power     2.585       0
   voltage       0       0
   current  0.1873       0
 
  c =
                power  voltage  current
   energy      -1.795  -0.5445   0.7394
   mass        0.8404        0    1.712
   force       -0.888  -0.6003        0
   frequency        0     0.49   -2.138
 
  d =
                angle   length
   energy           0    0.124
   mass         1.355        0
   force            0   -1.961
   frequency        0  -0.1977
 
Sample time: 1 seconds
Discrete-time state-space model.


sys_tf =
 
  From input "angle" to output...
               -4.502 z^2 + 2.961 z - 0.6717
   energy:  ------------------------------------
            z^3 - 0.6523 z^2 + 0.1906 z - 0.2858
 
          1.355 z^3 + 1.61 z^2 + 1.896 z - 1.116
   mass:  --------------------------------------
           z^3 - 0.6523 z^2 + 0.1906 z - 0.2858
 
               -2.296 z^2 + 1.2 z + 0.1663
   force:  ------------------------------------
           z^3 - 0.6523 z^2 + 0.1906 z - 0.2858
 
                  -0.4006 z^2 - 3.305 z + 0.8963
   frequency:  ------------------------------------
               z^3 - 0.6523 z^2 + 0.1906 z - 0.2858
 
  From input "length" to output...
   energy:  0.124
 
   mass:  0
 
   force:  -1.961
 
   frequency:  -0.1977
 
Sample time: 1 seconds
Discrete-time transfer function.


sys_zpk =
 
  From input "angle" to output...
             -4.5016 (z^2 - 0.6578z + 0.1492)
   energy:  -----------------------------------
            (z-0.8344) (z^2 + 0.1821z + 0.3425)
 
          1.3546 (z-0.4034) (z^2 + 1.592z + 2.042)
   mass:  ----------------------------------------
            (z-0.8344) (z^2 + 0.1821z + 0.3425)
 
              -2.296 (z-0.6363) (z+0.1139)
   force:  -----------------------------------
           (z-0.8344) (z^2 + 0.1821z + 0.3425)
 
                  -0.40058 (z+8.513) (z-0.2628)
   frequency:  -----------------------------------
               (z-0.8344) (z^2 + 0.1821z + 0.3425)
 
  From input "length" to output...
   energy:  0.12405
 
   mass:  0
 
   force:  -1.9609
 
   frequency:  -0.1977
 
Sample time: 1 seconds
Discrete-time zero/pole/gain model.


sys_frd =

  From input "angle" to:

    Frequency(rad/s)          energy              mass              force      
    ----------------          ------              ----              -----      
         0.0010         -8.7603 + 0.0415i  14.8297 - 0.0744i  -3.6837 + 0.0141i
         0.0017         -8.7599 + 0.0692i  14.8290 - 0.1241i  -3.6836 + 0.0235i
         0.0028         -8.7587 + 0.1155i  14.8270 - 0.2070i  -3.6832 + 0.0392i
         0.0046         -8.7556 + 0.1925i  14.8214 - 0.3452i  -3.6822 + 0.0653i
         0.0077         -8.7470 + 0.3208i  14.8059 - 0.5752i  -3.6794 + 0.1088i
         0.0129         -8.7230 + 0.5334i  14.7632 - 0.9564i  -3.6717 + 0.1810i
         0.0215         -8.6572 + 0.8819i  14.6455 - 1.5814i  -3.6504 + 0.2994i
         0.0359         -8.4801 + 1.4359i  14.3290 - 2.5752i  -3.5931 + 0.4880i
         0.0599         -8.0299 + 2.2456i  13.5239 - 4.0296i  -3.4475 + 0.7657i
         0.1000         -7.0331 + 3.1932i  11.7387 - 5.7383i  -3.1250 + 1.0988i

  From input "angle" to:

Frequency(rad/s)         frequency    
----------------         ---------    
     0.0010         -11.1245 + 0.0668i
     0.0017         -11.1239 + 0.1115i
     0.0028         -11.1221 + 0.1859i
     0.0046         -11.1171 + 0.3100i
     0.0077         -11.1031 + 0.5166i
     0.0129         -11.0645 + 0.8589i
     0.0215         -10.9583 + 1.4202i
     0.0359         -10.6726 + 2.3124i
     0.0599          -9.9458 + 3.6175i
     0.1000          -8.3336 + 5.1476i

  From input "length" to:

    Frequency(rad/s)    energy  mass   force   frequency
    ----------------    ------  ----   -----   ---------
         0.0010         0.1240    0   -1.9609   -0.1977
         0.0017         0.1240    0   -1.9609   -0.1977
         0.0028         0.1240    0   -1.9609   -0.1977
         0.0046         0.1240    0   -1.9609   -0.1977
         0.0077         0.1240    0   -1.9609   -0.1977
         0.0129         0.1240    0   -1.9609   -0.1977
         0.0215         0.1240    0   -1.9609   -0.1977
         0.0359         0.1240    0   -1.9609   -0.1977
         0.0599         0.1240    0   -1.9609   -0.1977
         0.1000         0.1240    0   -1.9609   -0.1977
 
Sample time: 1 seconds
Discrete-time frequency response.
units_workspace.mat
Reply | Threaded
Open this post in threaded view
|

Re: Request for Matlab text output

Lukas Reichlin-4

On 20.10.2015, at 14:35, NJank <[hidden email]> wrote:

> Lukas Reichlin-4 wrote
>> Dear Octave community,
>>
>> For the control package, I'm planning to increase Matlab compatibility by
>> supporting additional properties for LTI models.
>> Now I wonder how Matlab displays the models defined in the following
>> script file:
>>
>> http://sourceforge.net/p/octave/control/ci/default/tree/devel/units.m
>>
>> Could someone with access to a somewhat recent Matlab version (plus
>> Control System Toolbox) please run the script 'units.m' from the link
>> above and send me its output? (copy-pasted text, not a screenshot if
>> possible)
>>
>> Thanks in advance,
>> Lukas
>
> output from v2015a
> workspace attached.
>
> ran a ver, then your units.m script
>
>>> ver
> ----------------------------------------------------------------------------------------------------
> MATLAB Version: 8.5.0.197613 (R2015a)
> MATLAB License Number: 186008
> Operating System: Microsoft Windows Server 2012 Standard Version 6.2 (Build
> 9200)
> Java Version: Java 1.7.0_60-b19 with Oracle Corporation Java HotSpot(TM)
> 64-Bit Server VM mixed mode
> ----------------------------------------------------------------------------------------------------
> MATLAB                                                Version 8.5      
> (R2015a)
> Simulink                                              Version 8.5      
> (R2015a)
> Aerospace Blockset                                    Version 3.15      
> (R2015a)
> Aerospace Toolbox                                     Version 2.15      
> (R2015a)
> Bioinformatics Toolbox                                Version 4.5.1    
> (R2015a)
> Communications System Toolbox                         Version 6.0      
> (R2015a)
> Computer Vision System Toolbox                        Version 6.2      
> (R2015a)
> Control System Toolbox                                Version 9.9      
> (R2015a)
> Curve Fitting Toolbox                                 Version 3.5.1    
> (R2015a)
> DSP System Toolbox                                    Version 9.0      
> (R2015a)
> Data Acquisition Toolbox                              Version 3.7      
> (R2015a)
> Database Toolbox                                      Version 5.2.1    
> (R2015a)
> Datafeed Toolbox                                      Version 5.1      
> (R2015a)
> Econometrics Toolbox                                  Version 3.2      
> (R2015a)
> Financial Instruments Toolbox                         Version 2.1      
> (R2015a)
> Financial Toolbox                                     Version 5.5      
> (R2015a)
> Fixed-Point Designer                                  Version 5.0      
> (R2015a)
> Fuzzy Logic Toolbox                                   Version 2.2.21    
> (R2015a)
> Global Optimization Toolbox                           Version 3.3.1    
> (R2015a)
> Image Acquisition Toolbox                             Version 4.9      
> (R2015a)
> Image Processing Toolbox                              Version 9.2      
> (R2015a)
> Instrument Control Toolbox                            Version 3.7      
> (R2015a)
> MATLAB Coder                                          Version 2.8      
> (R2015a)
> MATLAB Compiler                                       Version 6.0      
> (R2015a)
> MATLAB Compiler SDK                                   Version 6.0      
> (R2015a)
> MATLAB Report Generator                               Version 4.1      
> (R2015a)
> Mapping Toolbox                                       Version 4.1      
> (R2015a)
> Model Predictive Control Toolbox                      Version 5.0.1    
> (R2015a)
> Neural Network Toolbox                                Version 8.3      
> (R2015a)
> OPC Toolbox                                           Version 3.3.3    
> (R2015a)
> Optimization Toolbox                                  Version 7.2      
> (R2015a)
> Parallel Computing Toolbox                            Version 6.6      
> (R2015a)
> Partial Differential Equation Toolbox                 Version 2.0      
> (R2015a)
> Robust Control Toolbox                                Version 5.3      
> (R2015a)
> Signal Processing Toolbox                             Version 7.0      
> (R2015a)
> SimEvents                                             Version 4.4      
> (R2015a)
> SimMechanics                                          Version 4.6      
> (R2015a)
> Simscape                                              Version 3.13      
> (R2015a)
> Simulink 3D Animation                                 Version 7.3      
> (R2015a)
> Simulink Coder                                        Version 8.8      
> (R2015a)
> Simulink Control Design                               Version 4.2      
> (R2015a)
> Simulink Design Optimization                          Version 2.7      
> (R2015a)
> Spreadsheet Link EX                                   Version 3.2.3    
> (R2015a)
> Stateflow                                             Version 8.5      
> (R2015a)
> Statistics and Machine Learning Toolbox               Version 10.0      
> (R2015a)
> Symbolic Math Toolbox                                 Version 6.2      
> (R2015a)
> System Identification Toolbox                         Version 9.2      
> (R2015a)
> Wavelet Toolbox                                       Version 4.14.1    
> (R2015a)
>>> units
>
> sys =
>
>  a =
>               x1          x2          x3
>   x1      0.3325     -0.5981     -0.1902
>   x2     -0.2371      0.3206     -0.5838
>   x3       0.581      0.1985  -0.0008821
>
>  b =
>           u1      u2
>   x1   2.585       0
>   x2       0       0
>   x3  0.1873       0
>
>  c =
>            x1       x2       x3
>   y1   -1.795  -0.5445   0.7394
>   y2   0.8404        0    1.712
>   y3   -0.888  -0.6003        0
>   y4        0     0.49   -2.138
>
>  d =
>            u1       u2
>   y1        0    0.124
>   y2    1.355        0
>   y3        0   -1.961
>   y4        0  -0.1977
>
> Sample time: unspecified
> Discrete-time state-space model.
>
>
> sys_ss =
>
>  a =
>                 power     voltage     current
>   power        0.3325     -0.5981     -0.1902
>   voltage     -0.2371      0.3206     -0.5838
>   current       0.581      0.1985  -0.0008821
>
>  b =
>             angle  length
>   power     2.585       0
>   voltage       0       0
>   current  0.1873       0
>
>  c =
>                power  voltage  current
>   energy      -1.795  -0.5445   0.7394
>   mass        0.8404        0    1.712
>   force       -0.888  -0.6003        0
>   frequency        0     0.49   -2.138
>
>  d =
>                angle   length
>   energy           0    0.124
>   mass         1.355        0
>   force            0   -1.961
>   frequency        0  -0.1977
>
> Sample time: 1 seconds
> Discrete-time state-space model.
>
>
> sys_tf =
>
>  From input "angle" to output...
>               -4.502 z^2 + 2.961 z - 0.6717
>   energy:  ------------------------------------
>            z^3 - 0.6523 z^2 + 0.1906 z - 0.2858
>
>          1.355 z^3 + 1.61 z^2 + 1.896 z - 1.116
>   mass:  --------------------------------------
>           z^3 - 0.6523 z^2 + 0.1906 z - 0.2858
>
>               -2.296 z^2 + 1.2 z + 0.1663
>   force:  ------------------------------------
>           z^3 - 0.6523 z^2 + 0.1906 z - 0.2858
>
>                  -0.4006 z^2 - 3.305 z + 0.8963
>   frequency:  ------------------------------------
>               z^3 - 0.6523 z^2 + 0.1906 z - 0.2858
>
>  From input "length" to output...
>   energy:  0.124
>
>   mass:  0
>
>   force:  -1.961
>
>   frequency:  -0.1977
>
> Sample time: 1 seconds
> Discrete-time transfer function.
>
>
> sys_zpk =
>
>  From input "angle" to output...
>             -4.5016 (z^2 - 0.6578z + 0.1492)
>   energy:  -----------------------------------
>            (z-0.8344) (z^2 + 0.1821z + 0.3425)
>
>          1.3546 (z-0.4034) (z^2 + 1.592z + 2.042)
>   mass:  ----------------------------------------
>            (z-0.8344) (z^2 + 0.1821z + 0.3425)
>
>              -2.296 (z-0.6363) (z+0.1139)
>   force:  -----------------------------------
>           (z-0.8344) (z^2 + 0.1821z + 0.3425)
>
>                  -0.40058 (z+8.513) (z-0.2628)
>   frequency:  -----------------------------------
>               (z-0.8344) (z^2 + 0.1821z + 0.3425)
>
>  From input "length" to output...
>   energy:  0.12405
>
>   mass:  0
>
>   force:  -1.9609
>
>   frequency:  -0.1977
>
> Sample time: 1 seconds
> Discrete-time zero/pole/gain model.
>
>
> sys_frd =
>
>  From input "angle" to:
>
>    Frequency(rad/s)          energy              mass              force      
>    ----------------          ------              ----              -----      
>         0.0010         -8.7603 + 0.0415i  14.8297 - 0.0744i  -3.6837 +
> 0.0141i
>         0.0017         -8.7599 + 0.0692i  14.8290 - 0.1241i  -3.6836 +
> 0.0235i
>         0.0028         -8.7587 + 0.1155i  14.8270 - 0.2070i  -3.6832 +
> 0.0392i
>         0.0046         -8.7556 + 0.1925i  14.8214 - 0.3452i  -3.6822 +
> 0.0653i
>         0.0077         -8.7470 + 0.3208i  14.8059 - 0.5752i  -3.6794 +
> 0.1088i
>         0.0129         -8.7230 + 0.5334i  14.7632 - 0.9564i  -3.6717 +
> 0.1810i
>         0.0215         -8.6572 + 0.8819i  14.6455 - 1.5814i  -3.6504 +
> 0.2994i
>         0.0359         -8.4801 + 1.4359i  14.3290 - 2.5752i  -3.5931 +
> 0.4880i
>         0.0599         -8.0299 + 2.2456i  13.5239 - 4.0296i  -3.4475 +
> 0.7657i
>         0.1000         -7.0331 + 3.1932i  11.7387 - 5.7383i  -3.1250 +
> 1.0988i
>
>  From input "angle" to:
>
> Frequency(rad/s)         frequency    
> ----------------         ---------    
>     0.0010         -11.1245 + 0.0668i
>     0.0017         -11.1239 + 0.1115i
>     0.0028         -11.1221 + 0.1859i
>     0.0046         -11.1171 + 0.3100i
>     0.0077         -11.1031 + 0.5166i
>     0.0129         -11.0645 + 0.8589i
>     0.0215         -10.9583 + 1.4202i
>     0.0359         -10.6726 + 2.3124i
>     0.0599          -9.9458 + 3.6175i
>     0.1000          -8.3336 + 5.1476i
>
>  From input "length" to:
>
>    Frequency(rad/s)    energy  mass   force   frequency
>    ----------------    ------  ----   -----   ---------
>         0.0010         0.1240    0   -1.9609   -0.1977
>         0.0017         0.1240    0   -1.9609   -0.1977
>         0.0028         0.1240    0   -1.9609   -0.1977
>         0.0046         0.1240    0   -1.9609   -0.1977
>         0.0077         0.1240    0   -1.9609   -0.1977
>         0.0129         0.1240    0   -1.9609   -0.1977
>         0.0215         0.1240    0   -1.9609   -0.1977
>         0.0359         0.1240    0   -1.9609   -0.1977
>         0.0599         0.1240    0   -1.9609   -0.1977
>         0.1000         0.1240    0   -1.9609   -0.1977
>
> Sample time: 1 seconds
> Discrete-time frequency response.
> units_workspace.mat
> <http://octave.1599824.n4.nabble.com/file/n4673019/units_workspace.mat>  

Thank you very much for your quick response!

I'm surprised that Matlab doesn't display the input/output/state units at all. I've expected that Matlab would display the units together with the corresponding input/output/state names, or at least in an extra line similar to input/output groups.

Best regards,
Lukas