在Matlab中的兩個函數來近似積分 - 沒有足夠的輸入參數？

Question

我想寫一個近似梯形法則積分的函數。

我首先在一個文件中定義的函數：

function[y] = integrand(x) 
y = x*exp(-x^2); %This will be integrand I want to approximate 
end 

然後我寫函數近似定積分與下界和上界B（也在另一文件）：

function [result] = trapez(integrand,a,b,k) 

sum = 0; 
h = (b-a)/k; %split up the interval in equidistant spaces 

for j = 1:k 
    x_j = a + j*h; %this are the points in the interval 
    sum = sum + ((x_j - x_(j-1))/2) * (integrand(x_(j-1)) + integrand(x_j)); 

end 
result = sum 
end 

但是當我想從命令窗口調用這個函數時，例如使用result = trapez(integrand,0,1,10)，我總是會得到一個錯誤'沒有足夠的輸入參數'。我不知道我做錯了什麼？

Answer 1

有你的代碼許多問題：

• x_(j-1)沒有定義，而不是真的是一個有效的Matlab語法（假設你希望它是一個變量）。
• 通過調用trapez(integrand,0,1,10)你實際上調用integrand函數沒有輸入參數。如果您想要傳遞句柄，請改用@integrand。但在這種情況下，根本不需要傳遞它。
• 您應該避免與Matlab函數一致的變量名稱，例如sum。如果您也嘗試使用sum作爲函數，這很容易導致難以調試的問題。

這裏有一個工作版本（注意也更好的代碼風格）：

function res = trapez(a, b, k) 
    res = 0; 
    h = (b-a)/k; % split up the interval in equidistant spaces 

    for j = 1:k 
     x_j1 = a + (j-1)*h; 
     x_j = a + j*h; % this are the points in the interval 
     res = res+ ((x_j - x_j1)/2) * (integrand(x_j1) + integrand(x_j)); 
    end 
end 

function y = integrand(x) 
    y = x*exp(-x^2); % This will be integrand I want to approximate 
end 

，並呼籲事情是這樣的：result = trapez(0, 1, 10);

Answer 2

你integrand函數需要輸入參數x，你是不是在你的命令行函數調用提供