作业帮 > 综合 > 作业

matlab int积分太慢,积不出来,怎么办?

来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/05/17 08:10:29
matlab int积分太慢,积不出来,怎么办?
clear
clc
syms x
format long
H=[0:0.01:1.2];
v=pi*4.1/180;a=0.4;b=2.05;n=1.2;m=1.78;
for i=1:15
h=H(i);
S2=(pi/2+(h+(a-x)*tan(v)-n/2)/n*2.*sqrt(1-((h+(a-x)*tan(v)-n/2)/n*2).^2)+asin((h+(a-x)*tan(v)-n/2)/n*2))*m*n/4;
v2(i)=int(S2,x,0,a+h*cot(v));
end
for i=16:118
h=H(i);
S2=(pi/2+((a-x)*tan(v)-n/2+h)/n*2.*sqrt(1-(((a-x)*tan(v)-n/2+h)/n*2).^2)+asin(((a-x)*tan(v)-n/2+h)/n*2))*m*n/4;
v2(i)=int(S2,x,0,a+b);
end
for i=119:121
h=H(i);
S2=(pi/2+(x*tan(v)-n/2)/n*2.*sqrt(1-((x*tan(v)-n/2)/n*2).^2)+asin((x*tan(v)-n/2)/n*2))*m*n/4;
vv2(i)=int(S2,x,0,(n-h)*cot(v)+b);
v2(i)=pi/4*m*n*(a+b)-vv2(i);
end
V2=eval(v2);
V1=V2';
这是计算CUMCM2010 A (1)题的一个程序,但是程序写对了,积分却积得很慢.2个小时都没积出来.
matlab 一直busy
电脑是win7 64位 家庭高级版,
用的是R2009A(32位) matlab
用quad(数值积分)来求.(int是符号积分)
编写函数M文件比如函数名f
那么quad(@f,1,2)就能求出1到2的积分