24小时热门版块排行榜    

查看: 250  |  回复: 2

yujunhui

银虫 (著名写手)

[求助] 积分程序报错,求指点?

global phi1 alfa1 Rsw Rsr L a b P1 n P gamma alfaL delta r deltar;
global zeta zetaa zetab phi2 phi3 phi4 phi11 phi10 phi9 ;
phi1=0*pi/180;%计算时需要修改的参数
alfa1=-32.915*pi/180;
Rsw=110.5*0.0001;Rsr=102.5*0.0001;%若是不同的螺杆,则需要修改
L=49*0.0001;a=164*0.0001;b=28*0.0001;%若是不同的螺杆,则需要修改
P1=101325;n=2970;P=11/6;gamma=2*pi/11;%固定几何参数
alfaL=acos((a-Rsr)/Rsw);delta=asin(b/2/Rsw);r=Rsw-L;deltar=asin(b/2/r);%固定几何参数
zeta=0.00001;zetaa=0.000005;zetab=zeta-zetaa;
phi2=gamma+phi1;phi3=gamma+phi2;phi4=gamma+phi3;% 求其余齿星轮转角
phi11=phi1-gamma;phi10=phi1-2*gamma;phi9=phi1-3*gamma;% 求其余齿星轮转角
%**************************介质参数*****************************************
global rho mu;
rho=850;mu=0.0256;
%************************压力函数******************************************
syms eta alfa;
f=sym('P*(Rsr^2-((a-cos(alfa)*sqrt(Rsw*Rsw-eta*eta)+eta*sin(alfa))^2))','eta','alfa');
eta1=-b/2;
eta2=b/2;
alfa3=alfa1+delta;
alfa4=alfaL-delta;
V2=dblquad(f,'eta','alfa',eta1,eta2,alfa3,alfa4);


系统提示  
??? Error using ==> sym
Too many input arguments
求指点?
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 积分.jpg
  • 2012-12-11 21:19:03, 9.79 K

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

luomia

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
xiegangmai: 金币+2, 谢谢参与 2012-12-13 23:42:01
yujunhui: 金币+30, ★★★★★最佳答案, 谢谢指点! 2012-12-15 11:20:29
方程表达式有问题,包含参数的运算符应该是数组运算
global phi1 alfa1 Rsw Rsr L a b P1 n P gamma alfaL delta r deltar;
global zeta zetaa zetab phi2 phi3 phi4 phi11 phi10 phi9 ;
phi1=0*pi/180;%计算时需要修改的参数
alfa1=-32.915*pi/180;
Rsw=110.5*0.0001;Rsr=102.5*0.0001;%若是不同的螺杆,则需要修改
L=49*0.0001;a=164*0.0001;b=28*0.0001;%若是不同的螺杆,则需要修改
P1=101325;n=2970;P=11/6;gamma=2*pi/11;%固定几何参数
alfaL=acos((a-Rsr)/Rsw);delta=asin(b/2/Rsw);r=Rsw-L;deltar=asin(b/2/r);%固定几何参数
zeta=0.00001;zetaa=0.000005;zetab=zeta-zetaa;
phi2=gamma+phi1;phi3=gamma+phi2;phi4=gamma+phi3;% 求其余齿星轮转角
phi11=phi1-gamma;phi10=phi1-2*gamma;phi9=phi1-3*gamma;% 求其余齿星轮转角
%**************************介质参数*****************************************
global rho mu;
rho=850;mu=0.0256;
%************************压力函数******************************************

f=@(alfa,eta) P*(Rsr^2-((a-cos(alfa).*sqrt(Rsw*Rsw-eta.*eta)+eta.*sin(alfa)).^2));
eta1=-b/2;
eta2=b/2;
alfa3=alfa1+delta;
alfa4=alfaL-delta;
V2=dblquad(f,eta1,eta2,alfa3,alfa4);
一次次萍水相逢
2楼2012-12-13 16:26:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yujunhui

银虫 (著名写手)

送鲜花一朵
谢谢指点!
安静
3楼2012-12-15 11:19:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yujunhui 的主题更新
信息提示
请填处理意见