24小时热门版块排行榜    

查看: 1207  |  回复: 0

丫头丫头2014

银虫 (小有名气)

[求助] matlab数值积分出现divide by zero 无法求出结果该如何解决

matlab数值求积分求面积,程序如下,报错,divide by zero,结果为 s1 = NaN +    NaNi,请高手帮忙看看

clear;
clc;
z1=6;
z2=z1-1;
R=48.78;
r=8.13;
Q=R-r;
e=7.05;
syms theta
k=12.6;
f=r/e;
a=sin(z1*theta);
b=f+cos(z1*theta);
m=atan(a./b);
p=m-theta;       %theta
c=sin(theta+p);
n=asin(f*c);
o=n-p;          %phi
g=(z1/z2)*o;    %psi

x2t=Q*sin(theta)+e*sin(z2*theta)-k*sin(p);       %原始齿形
y2t=Q*cos(theta)-e*cos(z2*theta)+k*cos(p);
x3t=x2t.*cos(o-g)-y2t.*sin(o-g)-e*sin(g);        %共轭齿形
y3t=x2t.*sin(o-g)+y2t.*cos(o-g)-e*cos(g);

d1=diff(x3t);
ds1=inline(y3t.*d1);
s1=quadl(ds1,0,2*pi/5)

下面是要求面积的图形

matlab数值积分出现divide by zero 无法求出结果该如何解决
转子.png
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 丫头丫头2014 的主题更新
信息提示
请填处理意见