未定義的函數或方法「INT」類型的輸入參數「雙重」

Question

我已經寫了代碼在蠕動泵送來評估流率的臨界值

dp = -(((theta - 1)+(1 + phi.*cos(2.*pi.*x))).*M.^3.*(cosh((1 + phi.*cos(2.*pi.*x)).*M)+ M.*beta.*sinh((1 + phi.*cos(2.*pi.*x)).*M)))./((1 + phi.*cos(2.*pi.*x)).*M.*cosh((1 + phi.*cos(2.*pi.*x)).*M)+(-1 +(1 + phi.*cos(2.*pi.*x)).*M.^2.*beta).*sinh((1 + phi.*cos(2.*pi.*x)).*M))-(alpha.*(M.^2.*((theta - 1)+(1 + phi.*cos(2.*pi.*x))).*(-1+2.*(1 + phi.*cos(2.*pi.*x)).^2.*M.^2+ cosh(2.*(1 + phi.*cos(2.*pi.*x)).*M)-2.*(1 + phi.*cos(2.*pi.*x)).*M.*sinh(2.*(1 + phi.*cos(2.*pi.*x)).*M)))./(8.*((1 + phi.*cos(2.*pi.*x)).*M.*cosh((1 + phi.*cos(2.*pi.*x)).*M)+(-1+(1 + phi.*cos(2.*pi.*x)).*M.^2.*beta).*sinh((1 + phi.*cos(2.*pi.*x)).*M)).^2)); 

和p = int(y,x,0,1 + phi.*cos(2.*pi.*x));

什麼是關鍵在θ= 0的臨界值時theta的值; 我寫的代碼：

M = 2; 
alpha = 0.2; 
beta = 0.03; 
phi = 0.6; 
theta = 0; 
x = 0.5; 
y = -(((theta - 1)+(1 + phi.*cos(2.*pi.*x))).*M.^3.*(cosh((1 + phi.*cos(2.*pi.*x)).*M)+ M.*beta.*sinh((1 + phi.*cos(2.*pi.*x)).*M)))./((1 + phi.*cos(2.*pi.*x)).*M.*cosh((1 + phi.*cos(2.*pi.*x)).*M)+(-1 +(1 + phi.*cos(2.*pi.*x)).*M.^2.*beta).*sinh((1 + phi.*cos(2.*pi.*x)).*M))-(alpha.*(M.^2.*((theta - 1)+(1 + phi.*cos(2.*pi.*x))).*(-1+2.*(1 + phi.*cos(2.*pi.*x)).^2.*M.^2+ cosh(2.*(1 + phi.*cos(2.*pi.*x)).*M)-2.*(1 + phi.*cos(2.*pi.*x)).*M.*sinh(2.*(1 + phi.*cos(2.*pi.*x)).*M)))./(8.*((1 + phi.*cos(2.*pi.*x)).*M.*cosh((1 + phi.*cos(2.*pi.*x)).*M)+(-1+(1 + phi.*cos(2.*pi.*x)).*M.^2.*beta).*sinh((1 + phi.*cos(2.*pi.*x)).*M)).^2)); 
p = int(y,x,0,1 + phi.*cos(2.*pi.*x)); 
sym theta 
solve('p','theta'); 

檢查一次，但我正在逐漸

Undefined function or method 'int' for input arguments of type 'double'. 

你能幫忙嗎？謝謝。

Answer 1

如果您要做的是符號集成，則需要符號工具箱即INT函數所屬的。變量應該是符號的（用sym或syms定義）。用文檔中的示例檢查函數語法。