Matlab不繪製，無限循環

Question

我試圖用Matlab中的二分法找到方程的根，因爲q以10步間隔從2000-3000變化。然而，我的代碼並沒有打印出一個圖表，即使我有一個plot聲明，我認爲它創建了一個無限循環，因爲當我運行它時，matlab說繁忙，除非強制關閉，否則我無法關閉該程序。我看不到任何會導致此問題的代碼，有人可以幫我解決嗎？

function myFunction 

a = 20; 
b = 40; 
tol = 1e-4; 
q = 2000:10:3000; 
t = zeros(101,1); 

for i=(1:length(q)) 
    f = @(x) (((1800).*log((160000)./(160000 - (x.*q(i)))) - (9.812).*x)./750) - 1; 
    t(i) = bisect(f,a,b,tol); 
end 

figure(1) 
plot(q,t) 

    function c=bisect(f,a,b,tol) 
     k=0; 
     while b-a > tol 
      c = (a-b)/2; 
      if sign(f(c)) == sign(f(b)) 
       b=c; 
      else 
       a=c; 
      end 
      k=k+1; 
     end 
    end 
end 

還應當指出的是，我已經使用這個bisect方法之前和它的工作，所以我不認爲這個問題是與該功能。

Answer 1

你的錯誤是在這裏：

c = (a-b)/2; 

您初始化a=20和b=40。 c最初設置爲-10。但是，你真的想c是半路a和b之間，這意味着你要：

c = (a+b)/2; 

Answer 2

此外，添加一個drawnow後，您繪製語句強制MATLAB繪製圖形。