2012-04-08 103 views
1

我已经写了代码在蠕动泵送来评估流率的临界值未定义的函数或方法“INT”类型的输入参数“双重”

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'. 

你能帮忙吗?谢谢。

回答

0

如果您要做的是符号集成,则需要符号工具箱INT函数所属的。变量应该是符号的(用symsyms定义)。用文档中的示例检查函数语法。

相关问题