Solving for theta at different time steps

3 次查看（过去 30 天）
Lee Quan
Lee Quan 2017-7-5
编辑： James Tursa 2017-7-5
I have the following equation:
47.9618*theta + 15.9809*sin(2*theta) = 20.780136*t
Given, at time t = 0, theta = 0 radians. I want to find theta value at t = 1 second, 2 seconds, etc.
Could this be automated on matlab?

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

采纳的回答

James Tursa
James Tursa 2017-7-5
编辑：James Tursa 2017-7-5
E.g.,
>> t = 0;
>> f = @(theta) 47.9618*theta + 15.9809*sin(2*theta) - 20.780136*t;
>> fzero(f,0)
ans =
0
>> t = 1;
>> f = @(theta) 47.9618*theta + 15.9809*sin(2*theta) - 20.780136*t;
>> fzero(f,0)
ans =
0.2649
>> t = 2;
>> f = @(theta) 47.9618*theta + 15.9809*sin(2*theta) - 20.780136*t;
>> fzero(f,0)
ans =
0.5651
  3 个评论
显示 1更早的评论
Lee Quan
Lee Quan 2017-7-5
Could you point me to some resources where the '@' notation is used so that I can familiarize myself with it? The @ part is something I can't find in matlab documentation online.
James Tursa
James Tursa 2017-7-5
编辑：James Tursa 2017-7-5
https://www.mathworks.com/help/matlab/matlab_prog/anonymous-functions.html?searchHighlight=anonymous%20function&s_tid=doc_srchtitle

请先登录，再进行评论。

更多回答（1 个）

Walter Roberson
Walter Roberson 2017-7-5
编辑：Walter Roberson 2017-7-5
There is no closed form solution for that. Use fzero or fsolve:
eqn = @(theta,t) 47.9618*theta + 15.9809*sin(2*theta) - 20.780136*t;
for t = 1 : 5
theta0 = randn()
fzero( @(theta) eqn(theta,t), theta0)
end

类别

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

标签

Community Treasure Hunt

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

Start Hunting!

Translated by