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
|
function iddataplot(varargin)
// Plot the iddata class
//
// Calling Sequence
// iddataplot(plantData)
// Parameters
// plantData : iddata type module
//
// Description
// Takes the iddata type input. Plot the respective input and output data
//
// Examples
// uData = idinput(1024,'PRBS',[0 1/20],[-1 1]);
// yData = rand(1024,1);
// Ts = 0.1;
// plantData1 = iddata(yData,uData,Ts)
// iddataplot(plantData1)
// plantData2 = iddata(yData,[],Ts)
// iddataplot(plantData2)
// plantData3 = iddata([],uData,Ts)
// iddataplot(plantData3)
//
// Authors
// Ashutosh Kumar Bhargava, Bhushan Manjarekar
[lhs rhs] = argn(0)
if rhs <> 1 then
error(msprintf(gettext("%s: Wrong number of input arguments."),"iddataplot"))
end
iddataData = varargin(1)
if typeof(iddataData) <> 'iddata' then
error(msprintf(gettext("%s:Wrong type for input argument %d: ""iddata"" expected."),"iddataplot",1))
end
// figure()
// xtitle('Plant i/o Data','Time('+iddataData.TimeUnit+')','Amplitude')
uData = iddataData.InputData
yData = iddataData.OutputData
timeLength = max(size(uData,'r'),size(yData,'r'))
timeData = ((1:1:timeLength)*iddataData.Ts)'
// timeData = timeData(1:length(iddataData)-1)
if size(uData,'*') && size(yData,'*') then
firstIndex = 2
else
firstIndex = 1
end
// ploting y data
h = gcf()
if h.figure_name == "Plant Input Output Data"; then
clf()
end
if size(yData,'*') then
secondIndex = size(yData,'c')
if secondIndex == 1 then
outputString = 'y'
else
outputString = []
for ii = 1:secondIndex
outputString = [outputString 'y'+string(ii)]
end
end
// disp(outputString)
for ii = 1:secondIndex
subplot(firstIndex,secondIndex,ii);plot(timeData,yData(:,ii))
xtitle(outputString(ii))
end
end
if size(uData,'*') then
secondIndex = size(uData,'c')
// disp(secondIndex)
if secondIndex == 1 then
outputString = 'u'
else
outputString = []
for ii = 1:secondIndex
outputString = [outputString 'u'+string(ii)]
end
end
// disp(outputString)
for ii = 1:secondIndex
if size(yData,'*') then
temp = 1
else
temp = 0
end
subplot(firstIndex,secondIndex,ii+secondIndex*temp);plot2d(timeData,uData(:,ii),2)
// pause
xtitle(outputString(ii))
end
end
h = gcf()
// disp(h)
h.figure_name= "Plant Input Output Data";
endfunction
|
