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
|
c---------------------------------------------------------------------------
c---------------------------------------------------------------------------
subroutine qzcel(t,x,nx,z,nz,u,nu,rpar,nrpar,ipar,nipar,nclock,
& out,nout,flag)
double precision t,x(*),u(*),rpar(*),out(*),z(*)
integer ipar(*),flag
c QZCEL, Alvaro:9-6-95
c Continous block, MIMO
c Gives quantized signal by ceiling method
c rpar(i) quantization step used for i input
c out=rpar*(floor(u/rpar))
c
integer i
common /dbcos/ idb
c
if(idb.eq.1) then
write(6,'(''qzcel t='',e10.3,'' flag='',i1)') t,flag
endif
goto(10,20,20) flag
c flag=1
10 continue
do 15 i=1,nu
out(i)=rpar(i)*ANINT(u(i)/rpar(i)-0.5)
15 continue
return
c flag=2 OR 3
20 continue
c ERROR
write(*,'(''ERROR; block ceil flag:'',i2,'' t '',e10.3)') flag,t
return
end
|
