How to plot a Damped and driven oscillation

I try to compute a simple damped and driven oscillation but it just won't work out. Is there anybody who can help me? This is my code:

    q = [0:10];
q = q(:);
    for(x=1:10)
        q(x);
        for(t=0:0.1:1)
            q(x)= exp(-t)*(exp(2*t)-exp(-2*t));
        end
    end

 plot(q)

[suman shah]

for dot product use .* together

Answer 1

[Chad Greene]

Hi Lukas,

For starters, it looks like you're doing a lot of unnecessary steps. I think the code above can be rewritten simply as

t = 0:0.1:1; 
q = exp(-t).*(exp(2*t)-exp(-2*t));
plot(t,q)

However, that's not going to give you a damped oscillation. Maybe you want something more like this?

t = 0:0.1:100; 
q = exp(-t/10).*sin(2*t);
plot(t,q)

[Lukas Goldschmied]

Thanks for your help. This is not exactly what i was looking for. It should look like this in the end: https://lp.uni-goettingen.de/get/image/2285

[Chad Greene]

What's the mathematical expression for the signal you're trying to plot?

[Lukas Goldschmied]

a*exp(-gt)*sinh(a*t) with and g are const.

[Chad Greene]

This should do the trick:

a = 5; 
g = 7; 

t = 0:0.01:1; 
q = a*exp(-g*t).*sinh(a*t); 

plot(t,q)