Matlab的集成

我寫在MATLAB代碼來計算的整體使用高斯 - 切比雪夫正交，但它不工作：Matlab的集成

function int = chebquad('-1i*exp(x+3)',1e-8,-1,1); 
f=inline('-1i*exp(x+3)','x') 
old_int = inf; 
for n=1:1000  
    x = cos(((2*(1:n) - 1)/(2*n))*pi);  
    w = pi/n;  
    fx = f(x);  
    int = sum(w.*fx);  
    if abs(int_old-int) < tol 
     break  
    end  
    old_int = int; 
end

有什麼建議？

謝謝！

來源

2012-05-12 Monique

定義「不起作用」。你是否要求它返回一個「int」值？浮點數不是一個更好的主意嗎？你將函數與exp或餘弦整合了嗎？請解釋你正在嘗試做什麼。 – duffymo

我正在試圖計算從-1到1的定積分-i次int e ^（x + 3）/ \ sqrt（1-x^2），也就是說我正在試圖積分指數函數，餘弦是重量函數，它不工作，因爲matlab給這個錯誤錯誤：文件：int.m行：1列：27 意外的MATLAB表達式。 – Monique

對不起，是「我」等於負數的平方根還是循環計數器？ – duffymo

以供將來參考，這將幫助我們，你沒有使用指定的變量作爲常量顯示你的函數：

這樣顯示：

function hv= someName(firstVar, secondVar)

，而不是：

function hv= someName(1, 'some string')

我不知道你爲什麼有-1和1，但我假設tol = 1e-8。話雖這麼說，你在你的代碼犯了個小錯誤使用int_old代替old_int.

編輯：所以起初我以爲，你只顯示功能「我們的利益」，現在我想你甚至沒有定義Matlab函數正確。請致電read this，瞭解基本的Matlab編碼。切比雪夫高斯積分定義了-1 to 1，因此並不需要在功能代碼，下面是修改後的代碼：

intV = chebquad('-1i*exp(x+3)', 1e-8)

：

function intV = chebquad(funC,tol) 

f=inline(funC,'x'); 
old_int = inf; 
for n=1:1000  
    x = cos(((2*(1:n) - 1)/(2*n))*pi);  
    w = pi/n;  
    fx = f(x);  
    intV = sum(w.*fx);  
    if abs(old_int - intV) < tol 
     break  
    end  
    old_int = intV ; 
end

在命令提示符下，你把這個

來源

2012-05-12 17:22:01 Rasman

+1 - 尼斯答案。做得好。 – duffymo

謝謝拉斯曼;） – Monique