Numerical computation using Matlab

Question

Muhammad Bilal 2020-7-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/569373-numerical-computation-using-matlab

关闭： MATLAB Answer Bot 2021-8-20

在 MATLAB Online 中打开

Hi,

I am trying to solve the expression for theta with T as an input:

Input: T

Output: theta

How can i solve to find the solution for a range of T?

Ls = 0.090;
r1 = 0.035;
Ks = 445e3;
xp = 0.1403;
func = Ks .*((sqrt((r1^2 + xp^2) - (2.*r1.*xp.*cos(theta))))...
    - Ls).*r1.*(xp./(sqrt((r1^2 + xp^2) - (2.*r1.*xp.*cos(theta))))).*sin(theta) - T;

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

此问题已关闭。

Answer 1

Alan Stevens 2020-7-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/569373-numerical-computation-using-matlab#answer_469695

在 MATLAB Online 中打开

Here's a (not very elegant!) way:

Ls = 0.090;
r1 = 0.035;
Ks = 445e3;
xp = 0.1403;
theta0 = 0;
TT = 0:50:500;
Theta = zeros(size(TT));
for i = 1:length(TT)
        
    T = TT(i);
    func = @(theta) Ks.*((sqrt((r1^2 + xp^2) - (2.*r1.*xp.*cos(theta))))...
    - Ls).*r1.*(xp./(sqrt((r1^2 + xp^2) - (2.*r1.*xp.*cos(theta))))).*sin(theta) - T;
    Theta(i) = fzero(func, theta0);
        
end
plot(TT,Theta*180/pi),grid
xlabel('T'), ylabel('\theta [degrees]') 