关于matlab积分对于这样一个积分,采用matlab编程syms t r a;R1=@(t,r,a)(((35-9*t)^2+(r*cos(a))^2+(r*sin(a)-14)^2)^0.5);R2=@(t,r,a)(((35-9*t)^2+(r*cos(a))^2+(14-r*sin(a))^2)^0.5);f=(35-9*t)./(2880*pi)*((r*cos(a)-14)/(R1)^3+(1.88*(r*cos(a)-
来源：学生作业帮助网 编辑：六六作业网 时间：2024/05/05 09:35:07
关于matlab积分对于这样一个积分,采用matlab编程syms t r a;R1=@(t,r,a)(((35-9*t)^2+(r*cos(a))^2+(r*sin(a)-14)^2)^0.5);R2=@(t,r,a)(((35-9*t)^2+(r*cos(a))^2+(14-r*sin(a))^2)^0.5);f=(35-9*t)./(2880*pi)*((r*cos(a)-14)/(R1)^3+(1.88*(r*cos(a)-
关于matlab积分
对于这样一个积分,采用matlab编程
syms t r a;
R1=@(t,r,a)(((35-9*t)^2+(r*cos(a))^2+(r*sin(a)-14)^2)^0.5);
R2=@(t,r,a)(((35-9*t)^2+(r*cos(a))^2+(14-r*sin(a))^2)^0.5);
f=(35-9*t)./(2880*pi)*((r*cos(a)-14)/(R1)^3+(1.88*(r*cos(a)-14)/(R1)^3)+1.2672/(R2*(R2+14-r*sin(a))));
S=int((int(f,r,0,3)),a,0,2*pi)
出错,Undefined function 'mpower' for input arguments of type 'function_handle'
clear all； clc；
syms t r a;
f=(35-9*t)./(2880*pi)*((r*cos(a)-14)/(((35-9*t)^2+(r*cos(a))^2+(r*sin(a)-14)^2)^0.5)^3+(1.88*(r*cos(a)-14)/(((35-9*t)^2+(r*cos(a))^2+(r*sin(a)-14)^2)^0.5)^3)+1.2672/((((35-9*t)^2+(r*cos(a))^2+(14-r*sin(a))^2)^0.5)*((((35-9*t)^2+(r*cos(a))^2+(14-r*sin(a))^2)^0.5)+14-r*sin(a))));
S=int((int(f,r,0,3)),a,0,2*pi)
警告：Warning:Explicit integral could not be found.
这个积分到底能不能求解,
关于matlab积分对于这样一个积分,采用matlab编程syms t r a;R1=@(t,r,a)(((35-9*t)^2+(r*cos(a))^2+(r*sin(a)-14)^2)^0.5);R2=@(t,r,a)(((35-9*t)^2+(r*cos(a))^2+(14-r*sin(a))^2)^0.5);f=(35-9*t)./(2880*pi)*((r*cos(a)-14)/(R1)^3+(1.88*(r*cos(a)-
syms类型就没必要定义函数了.积分有三个参数,所以最后的S应该是三重积分才对.
f中有个除号是./ 把点去掉,另外sin抄成了cos
clear;
clc;
syms t r a;
R1=(((35-9*t)^2+(r*cos(a))^2+(r*sin(a)-14)^2)^0.5);
R2=(((35-9*t)^2+(r*cos(a))^2+(14-r*sin(a))^2)^0.5);
f=(35-9*t)/(2880*pi)*((r*sin(a)-14)/(R1)^3+(1.88*(r*sin(a)-14)/(R1)^3)+1.2672/(R2*(R2+14-r*sin(a))));
更正后就和你下面一样了.可能是不可积分类型的,因此没有解析式.
但是使用简单的解析式实验可以求解,用subs()可以解除数值解.简单解析式用上面方法可以通过.
clear;
clc;
syms t r a;
f=sin(t)*cos(r)*cos(a);
g=f*r;
temp=int(int(int(int(g,t,0.01,1)),r,0.01,1),a,0.01,pi);
subs(temp);
就上面不太复杂的都需要计算一会.
综上,1.可能是积分本身无解,2.可能计算量超大,matlab计算一会就不算了.