Using piecewise-constant functions in ode Simulink

Hi,

I want to integrate the multiplication of a cosine with a
piecewise-constant function (e.g. 1 from 0 to 1/1000, -1 from 1/1000 to
2/1000, ...)

The time axis is equally distributed among the time (i.e., transient
change every 1/1000) and the values are stored in an array in MATLAB.

I'm struggling arround for 2 days on this issue :( So help is greatly
appreciated!!


What I did was (among other tries): Create a pulse block with rising
edge every 1/1000 which triggers a "Triggered Subsystem". This calls a
MATLAB function which uses find_system/set_param to change the value of
a constant block. This constant is then multiplied by the cosine and
followed by an integrator.

However, this seems to work for some signal but for some it introduces
huge errors.

For example

p = 1 for [0,1/1000], -1 for [1/1000,2/1000] and 1 for [2/1000,3/1000]
x = cos(100*pi*t);

y = (p*x)-> 1/s block

y at t=3/1000 should be 8.00487e-4 if I solve analytically but Simulink
has 8.29766161000261e-4 :( :(

Thanks a lot!
Peter