24小时热门版块排行榜    

查看: 395  |  回复: 1

wzxchars

铁杆木虫 (著名写手)

[求助] matlab积分计算问题求助已有1人参与

很早以前写的一段计算小程序,在matlab8的时候能够运行,前几天需要再计算一下,发现在2014里面无法运行,
提示错误为:
未定义与 'double' 类型的输入参数相对应的函数 'int'。

出错         A2=int(A1,0,4);

求助高人修改下。

程序如下:


% A=int(int f(x,y)*df(x,y)/dx)dxdy;

clear;
format long g;

a=4;
b=4;



f1='(4^3+2*x^3-3*4*x^2)*(4^3+2*y^3-3*4*y^2)';
f2='4*x*(x-4)^2*(4^3+2*y^3-3*4*y^2)';
f3='(4^3+2*x^3-3*4*x^2)*4*y*(y-4)^2';
f4='4*4*x*y*(x-4)^2*(y-4)^2';
f5='(3*4*x^2-2*x^3)*(4^3+2*y^3-3*4*y^2)';
f6='4*(x^3-4*x^2)*(4^3+2*y^3-3*4*y^2)';
f7='4*(3*4*x^2-2*x^3)*y*(y-4)^2';
f8='4*4*x*y*(x^2-4*x)*(y-4)^2';
f9='(3*4*x^2-2*x^3)*(3*4*y^2-2*y^3)';
f10='4*(x^3-4*x^2)*(3*4*y^2-2*y^3)';
f11='4*(3*4*x^2-2*x^3)*(y^3-4*y^2)';
f12='4*4*x*y*(x^2-4*x)*(y^2-4*y)';
f13='(4^3+2*x^3-3*4*x^2)*(3*4*y^2-2*y^3)';
f14='4*x*(x-4)^2*(3*4*y^2-2*y^3)';
f15='4*(4^3+2*x^3-3*4*x^2)*(y^3-4*y^2)';
f16='4*4*x*y*(x-4)^2*(y^2-4*y)';


df1=diff(f1);
df2=diff(f2);
df3=diff(f3);
df4=diff(f4);
df5=diff(f5);
df6=diff(f6);
df7=diff(f7);
df8=diff(f8);
df9=diff(f9);
df10=diff(f10);
df11=diff(f11);
df12=diff(f12);
df13=diff(f13);
df14=diff(f14);
df15=diff(f15);
df16=diff(f16);


f=str2mat(f1,f2,f3,f4,f5,f6,f7,f8,f9,f10,f11,f12,f13,f14,f15,f16);
df=[df1,df2,df3,df4,df5,df6,df7,df8,df9,df10,df11,df12,df13,df14,df15,df16];

Aero=zeros(16,16);

for i=1:16;
    for j=1:16
        A1=f(i,*df(j);
        A2=int(A1,'0','4');
        A3=int(A2,'0','4');
        Aero(i,j)=A3/(4^11);
    end
end
回复此楼
静坐常思自己过,闲谈莫论他人非
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +1
https://www.mathworks.com/help/symbolic/integration.html
CODE:
syms y
f = y^(-1);

int(f)
ans =
log(y)

2楼2016-10-13 22:16:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wzxchars 的主题更新
信息提示
请填处理意见