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
|
function calculate_cont_frm(filename, num, den)
f_temp = mopen(filename, 'wt');
H = cont_frm(num, den);
[A, B, C, D] = abcd(H);
cont_frm_write(A, B, C, D, f_temp);
mclose(f_temp);
endfunction
function cont_frm_write(varargin)
loop = argn(2);
f_temp = varargin(loop);
mfprintf(f_temp, '[')
for k = 1:(loop-1)
variable = varargin(k)
[m, n] = size(variable)
for y = 1:n
for z = 1:m
if (k == loop-1) & (y == n) & (z == m) then
mfprintf(f_temp, '%.17g', variable(z, y))
else
mfprintf(f_temp, '%.17g,', variable(z, y))
end
end
end
end
mfprintf(f_temp, ']')
endfunction
|
