24小时热门版块排行榜    

查看: 670  |  回复: 1

xiaozhixiao

铁虫 (初入文坛)

[求助] 利用quad2dggen求解二重积分的问题

我编写了一个函数,可以有三个变量,求解斜圆锥(底面被斜切)体积,
??? Error: File: D:\MATLAB7\work\三维分析\jifen3.m Line: 12 Column: 35
Missing MATLAB operator.
错误提示的意思是我 strf  里面写f0 表达式是弄错了,求那个高手给看下,万分感激!被积分函数  fo=tan(arf)+k-1/tan(fai)*sqrt(2^2+y^2)

function  y=jifen3(afj,lin,fai)     %求解斜圆锥的体积
  
   strfl=['-(-(1-tand(',num2str(afj),')^2*tand(',num2str(fai),')^2)*(y^2)+2*',num2str(lin),'*tand(',num2str(afj),')*tand(',num2str(fai),')^2*y+',num2str(lin),'^2*tand(',num2str(fai),')^2)^0.5'];
   fl=inline(strfl,'y');
   
   strfh=['(-(1-tand(',num2str(afj),')^2*tand(',num2str(fai),')^2)*(y^2)+2*',num2str(lin),'*tand(',num2str(afj),')*tand(',num2str(fai),')^2*y+',num2str(lin),'^2*tand(',num2str(fai),')^2)^0.5'];
   fh=inline(strfh,'y');
   
   yn=-(-2*lin*tand(afj)*tand(fai)^2)/2/(1-tand(afj)^2*tand(fai)^2)-((-2*lin*tand(afj)*tand(fai)^2)^2-4*(1-tand(afj)^2*tand(fai)^2)*(-lin^2*tand(fai)^2))^0.5/2/(1-tand(afj)^2*tand(fai)^2);
   ym=-(-2*lin*tand(afj)*tand(fai)^2)/2/(1-tand(afj)^2*tand(fai)^2)+((-2*lin*tand(afj)*tand(fai)^2)^2-4*(1-tand(afj)^2*tand(fai)^2)*(-lin^2*tand(fai)^2))^0.5/2/(1-tand(afj)^2*tand(fai)^2);
  
   strf=['tand(',num2str(fai),')*'y'+',num2str(lin),''-1/tand(',num2str(fai),')*sqrt(x^2+y^2)'];
   f0=inline(strf,'x','y');
   
   y=quad2dggen(f0,fl,fh,yn,ym,eps);
end
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : NIT.rar
  • 2012-04-17 15:30:35, 21.22 K
  • 附件 2 : jifen3.m
  • 2012-04-17 15:30:44, 1.05 K

» 猜你喜欢

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

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

xiaozhixiao

铁虫 (初入文坛)

好心人给看看啊,解决了我一定给以重谢啊!
我来了
2楼2012-04-17 15:34:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaozhixiao 的主题更新
信息提示
请填处理意见