integration using set of matrix equation

1 次查看（过去 30 天）
areej abdulshaheed
评论： Walter Roberson 2020-3-2
if I have the accelartion
and I would like to calculate the theta and its derivative angles
L=1.6; n=16; l=0.1; m=0.1; I=0.0001 g=9.8; kn=2;
alpha=pi/6; mu=0.3; a=0.0625; T=32;
xl_0=0; zl_0=0; xD1_0=0; zD1_0=0; pD1_0=0; G=m*n*g;
theta(1:n-1)=-2*alpha*sin(kn*pi/n)*sin(2*kn*pi*s/L+2*kn*pi*(1:n-1)/n+kn*pi/n);
thetaD(1:n-1)=(-4*kn*pi*alpha/L)*sin(kn*pi/n)*cos(2*kn*pi*s/L+2*kn*pi*(1:n-1)/n+kn*pi/n)*sd;
phi(1:n)=alpha;
phi(2:n)=phi(1:n-1)+theta(1:n-1);
how can calculate them
can anyone help me please?

请先登录，再回答此问题。

回答（1 个）

Walter Roberson
Walter Roberson 2020-3-2
In https://www.mathworks.com/matlabcentral/answers/506724-how-to-solve-this-equation-for-laplace-transform-with-matlab#comment_802772 I show how to use heaviside() to implement piecewise calculations in a way that can be integrated (or fourier transformed, or laplace transformed.) See in particular the helper function I named R. The first parameter is the variable name; the second parameter is the lower bound that applies; the third parameter is the upper bound that applies; the last parameter is the value that is to apply inside that range.
  4 个评论
显示 2更早的评论
areej abdulshaheed
areej abdulshaheed 2020-3-2
it is the- s- in theta equation
sin(2*kn*pi*s/L+2*kn*pi*(1:n-1)/n+kn*pi/n);
Walter Roberson
Walter Roberson 2020-3-2
Is the same as the s in that equation, or is the the second derivative of the s in that equation?
The is effectively a function of t but in the equation for theta you are treating s as a variable. I do not see any t in the theta equation ? Is the idea that at each different angle, you would want a different time-based formula ?

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Calculus 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by