"Matthew " <firstname.lastname@example.org> wrote in message news:email@example.com... > So what I'm trying to do is simulate different damping ratios (labeled > "C") and simply use those in a bunch of equations for outputting a plot of > the displacement. However I'm having matrix agreement issues and I can't > figure out why. My code is: > > t=0:0.1:5; > C=[.01 .2 .6 .1 .4 .8]; > Wn=2; > Wd=Wn*sqrt(1-C.^2); > A=sqrt(((C.*Wn*1).^2 + Wd.^2)./(Wd.^2)); > phi=atan(Wd./(C.*Wn*1)); > x=A.*exp(-C.*Wn*t).*sin(Wd.*t + phi); > plot(x,t) > > I've found that the issue takes place in my "x" equation between Wn*t and > I can't figure out what is wrong. Excuse all the dot assignments, as I'm > not sure if all them were needed or not.
.* is elementwise multiplication. The two quantities being multiplied must be EXACTLY the same size (in which case the result is that same size) or one must be scalar (in which case the result is the size of the other.)
* is matrix multiplication. The number of columns in the first quantity being multiplied must be equal to the number of rows in the second (in which case the result has the same number of rows as the first and the same number of columns as the second), or one must be scalar (in which case the result is the size of the other.)
Since Wn is scalar the multiplication C.*Wn is valid and results in a vector the same size as C. Then (C.*Wn)*t is attempting to perform matrix multiplication on two vectors that do NOT satisfy the requirements, and that's the cause of the problem. I'm guessing you want to perform elementwise multiplication there, but C and t are not the same size so that won't work. You may want to change your t to something like: