MATLAB中有壹個內置函數是humps,該函數是下列函數的特例:
當q=0.3,r=0.9和s=6的時候就是humps函數了。humps函數到底有什麽 用呢?
根據MATLAB的幫助,如下:
>> help humps
HUMPS ?A function used by QUADDEMO, ZERODEMO and FPLOTDEMO.?Y = HUMPS(X) is a function with strong maxima near x = .3?
?and x = .9. ?
?[X,Y] = HUMPS(X) also returns X. ?With no input arguments,
?HUMPS uses X = 0:.05:1.
?Example:
plot(humps)
即該函數通常用於幾個函數的測試工作,這個函數有個顯著的特點就是在0.3和0.9處有很強的值變化,在相對較短的區間內展現出平滑和陡峭的曲線。
繪制humps函數可以使用如下簡單的語句:
>> x=0:0.001:1;
>> y=humps(x);
>> plot(x,y);
>> grid on
圖形如下:
首先,這個程序是可積的,而且更美妙的是可以有解析公式,這樣得到精確結果就很簡單了,然而,真正需要 使用這個程序可以測試數值積分的穩定性和計算精度,尤其是像quad和quadl函數的測試工作,如果妳覺得自己寫出什麽比較高明的數值積分公式或者方法、函數,不妨用它來測試壹把。~