方程式は次のとおりです。
ここで、、a
およびb
は定数であり、両方とも等しい0.0130
これは私が使用しているコードです:
% Solving the equation for zero.
f = @(theta) ((a+b).*theta)./((theta.^2)-(a.*b)) - tan(theta); % Notice the dots (.)
% Now plot it to get an idea of where the zeros are.
theta = 0:1:100;
for i=1:length(theta)
hold on
plot(theta(i),f(theta(i)),'-o') % Look for the zeros
end
% Now find the roots.
cnt = 1;
for ii = [0,2,50] % This vector has the guesses.
rt(cnt) = fzero(f,ii); % Pass each guess to FZERO.
cnt = cnt + 1;
end
これは私が得るエラーです:
??? Operands to the || and && operators must be convertible to
logical scalar values.
Error in ==> fzero at 323
elseif ~isfinite(fx) || ~isreal(fx)
Error in ==> HW4 at 52
rt(cnt) = fzero(f,ii); % Pass each guess to FZERO.
\thetaの最初の解を取得したいと思います。ありがとう。