1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
|
function varargout = iddata(varargin)
// Stores plant input-output data
//
// Calling Sequence
// plantData = iddata(yData,uData,Ts)
// plantData = iddata([],uData,Ts)
// plantData = iddata(yData,[],Ts)
// plantData = iddata(yData,uData)
// plantData = iddata(yData)
//
// Parameters
// uData : nx1 matrix of the plant input data
// yData : nx1 matrix of the plant output data
// Ts : non-negative real number
// plantData : iddata type module
//
// Description
// It is a iddata type module that stores the plant input and output data with sampling time Ts. The time unite is in second.
//
// Examples
// uData = idinput(1024,'PRBS',[0 1/20],[-1 1])
// yData = rand(1024,1)
// Ts = 0.1
// plantData1 = iddata(yData,uData,Ts)
// plantData2 = iddata(yData,[],Ts)
// plantData3 = iddata([],uData,Ts)
//
// Authors
// Ashutosh Kumar Bhargava
[lhs,rhs] = argn(0)
if rhs == 2 | rhs == 1 then
Ts = 1
elseif rhs == 3 then
Ts = varargin(3)
else
error(msprintf(gettext("%s:Incorrect number of input arguments.\n"),"iddata"))
end
if Ts <= 0 || typeof(Ts) <> 'constant' || ~size(Ts,'*') then
error(msprintf(gettext("%s:Inconsist sampling time ""Ts"", ""non negative real number"" expected.\n"),"iddata"))
end
if rhs == 1 then
OutputData = varargin(1)
InputData = []
elseif rhs == 2 | rhs == 3 then
OutputData = varargin(1)
InputData = varargin(2)
if isrow(InputData) then
InputData = InputData'
end
end
if size(OutputData,'*') & size(InputData,'*') then
if size(OutputData,'r') ~= size(InputData,'r') then
error(msprintf(gettext("%s:The numbers of the plant out datas must be equal to the numbers of the plant input datas.\n"),"iddata"))
end
end
t = tlist(['iddata','OutputData','InputData','Ts','TimeUnit'],OutputData,InputData,Ts,'seconds')
varargout(1) = t
endfunction
function %iddata_p(mytlist)
f = fieldnames(mytlist)
if ~size(mytlist(f(1)),'*') & ~size(mytlist(f(2)),'*') then
mprintf(' Empty sample data.\n')
else
outputSize = size(mytlist(f(1)))
inputSize = size(mytlist(f(2)))
if prod(outputSize) then
sampleSize = max(outputSize)
elseif prod(inputSize) then
sampleSize = max(inputSize)
end
mprintf(' Time domain sample data having %d samples.',sampleSize)
if round(mytlist(f(3)))-mytlist(f(3)) == 0 then
mprintf('\n Sampling Time = %d',mytlist(f(3)))
else
mprintf('\n Sampling Time = %f',mytlist(f(3)))
end
mprintf(' %s',mytlist(f(4)))
mprintf('\n')
if prod(outputSize) then
mprintf('\n Output channel \n')
for ii = 1:min(outputSize)
yString = 'y' + string(ii)
mprintf(' %s\n',yString)
end
end
if prod(inputSize) then
mprintf('\n Input channel \n')
for ii = 1:min(inputSize)
uString = 'u' + string(ii)
mprintf(' %s\n',uString)
end
end
end
mprintf('\n')
endfunction
|
