下面积分运算表达式错在什么地方？

ssg

>> x=0.26;
y=13.3333*((0.2/x).^(3)-(0.2/x).^(9));
>> int(y,0.26,inf)
??? Function 'int' is not defined for values of class 'double'.
本想对图示函数积分进行运算，尝试对其编辑，但是老是出错，请朋友们指教。

lin2009

“int(y,0.26,inf)”
x＝0.26已赋值，y已有确切的值，int(y,0.26,inf)是对常数积分的，而不是你要的y关于x的表达式。

ssg

x是用来描述h0,该如何定义它呢

hill711

你可以直接积分，然后再计算的

时间有限，难免有误，但是思路是没有问题的

（可怜俺的手机，就这效果了，凑乎看吧）

hill711

还能看清楚
貌似第四行少个中括号
看来着急弄是不行的

ssg

十分感谢！好久没用高数积分，把这么重要的处理方法都忘在一边了。

hill711

十分感谢！好久没用高数积分，把这么重要的处理方法都忘在一边了。
ssg 发表于 2010-7-20 20:03

要温故而知新哦~
问题这不就迎刃而解了！

lin2009

直接的方法：
clear all;
clc;
syms x y
y = 40/3*((0.2/x)^3-(0.2/x)^9);
yfcn = inline(char(y));
x1 = 0.26:0.1:2.0;
hc = zeros(size(x1));
k = 1;
for x1 = 0.26:0.1:2.0
    hc(k) = double(x1 + int(y,x1,inf)/yfcn(x1));
    k = k + 1;
end
hc
% int  不支持矢量运算

间接的方法（先推导不定积分）
clear all;
clc;
syms x y
y = 40/3*((0.2/x)^3-(0.2/x)^9);
yfcn = vectorize(inline(char(y)));
yint = vectorize(inline(char(int(y,x))));
yint2 = vectorize(inline(char(yint(inf)-yint(x))));
x = 0.26:0.1:2.0;
hc = double(x + yint2(x)./yfcn(x))

ssg

感谢大家帮忙！

luoyecolin

:) 谢谢！！！！！！！！！！！！！！！！！！！！