Definte Integration of a Received Sequence

Hello all,

I have the sequence v of length K-by-1. I want to do integration of v over nn*nsamp to (nn+1)*nsamp for nn=0, ..., N-1 and K=nsamp*N. I did this in MATLAB as following:

y=zeros(N,1);
for nn=0:N-1
y(nn)=sum(v(nn*nsamp:(n+1)*nsamp))
end

Is this right?

[Mahdi]

Unlike some other programming languages, the index for the first element in a matrix/vector is 1 in MATLAB. So for example, if A=[5 6 7], then A(1) gives 5. While A(0) gives an error. So in your case, you want to define nn from 1 to N-1:

for nn=0:N-1

You might also want to look into the trapz built-in function to integrate. What is the range of x-values you're integrating over?

[S. David]

Yes you are right, I should have written

...
y(nn+1)=sum(v(nn*nsamp+1:(n+1)*nsamp))
...

Actually, v corresponds to the continuous time signal v(t) where

t=0:Ta:(K-1)*Ta

I want to do the integration:

yn=int_{n*Ts}^{(nn+1)*Ts} v(t) dt % n=0,...,N-1

where Ta=Ts/nsamp.