| 查看: 236 | 回复: 1 | ||
[求助]
错误出在哪里了？
|
|
>> a=5.0;L=10.0;B=10.0; rho=1000.0;gra=9.81;W=100000.0; ratio=L/a; theta=linspace(0.0,pi/2.0); VoverVw=cos(theta)+0.5*ratio*cos(theta).*(theta); options=optimset('display','off'); figure(1) subplot(1,2,1) plot(theta*180/pi,VoverVw,'k') axis([0.0,90.0,0.0,1.5]) ylabel('V/V_w') xlabel('/theta(degrees)') grid on hold on MaxTheta=inline('1-(cos(theta)+0.5*ratio*cos(theta).*sin(theta))','theta','ratio'); ThetaMaxDeg=fzero(MaxTheta,[0.01,pi/2.0],options,ratio)*180/pi; plot(ThetaMaxDeg,1.0,'sk') text(19,0.95,['\theta_{max}=',num2str(ThetaMaxDeg,4)'\'circ']) ThetaMax Vol=fiminbnd(MaxTheta,0.0,ThetaMaxDeg*pi/180,options,ratio); MaxVol=1-MaxTheta(ThetaMaxVol,ratio); plot(ThetaMaxVol*180.0/pi,MaxVol,'ks') text(10,MaxVol+0.1,['max vol='num2str(MaxVol,4)]) text(10,MaxVol+0.05,['at\theta='num2str(ThetaMaxVol*180/pi,4)'\circ']) subplot(1,2,2) theta=linspace(0.01,ThetamaxDeg*pi/180); h=zeros(1,length(theta)); for i=1:length(theta) r=roots([tan(theta(i))*0.5,a,-a*L]); h(i)=r(2); end plot(theta*180/pi,h) ylabel('h(meters)') xlabel('\theta(degrees)') grid on figure(2) Frod=((B*rho*grav*h.^3)./(6*cos(theta).^2)+0.5*W*L*sin(theta))./(L*sin(theta+acos(cos(theta)/sqrt(2)))); plot(theta*180.0/pi,Frod*le-6 ylabel('F_{rod}(MN)') xlabel('\theta(degrees)') grid on Matrix dimensions must agree.（运行后所显示的）怎么来解决啊？？？ |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有9人回复
-
救命帖
已经有4人回复
-
限项规定
已经有5人回复
-
为什么nbs上溴 没有产物点出现呢
已经有9人回复
-
招博士
已经有3人回复
-
存款400万可以在学校里躺平吗
已经有35人回复
-
最失望的一年
已经有18人回复
-
求推荐博导
已经有4人回复
-
求助一下有机合成大神
已经有4人回复
-
求推荐英文EI期刊
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求大神帮我改一下下面的程序！matlab迭代法的 不知道错在哪里！
已经有7人回复
- gromacs水中的溶菌酶案例,哪出错了?
已经有7人回复
- 《用正确的方法解决问题》【转载】
已经有217人回复
- arcgis10不能导入excel数据，出现“找不到可安装的ISAM”错误
已经有5人回复
- 求高人帮忙看一下我这个Fortran程序的错误出在哪里？
已经有12人回复
- 烧陶瓷的三相加热电炉有一相没有示数，是哪里出问题了？
已经有9人回复
- 在家上不去小木虫，哪出问题了？
已经有13人回复
- 层错能与什么有关？
已经有5人回复
- 高手帮助分析一下流式细胞仪双染测细胞凋亡的问题出在哪？
已经有8人回复
- 网页打开错误，是什么情况？
已经有20人回复
- 【求助】锂电池循环伏安是一条水平直线，问题出在哪里
已经有17人回复
- 【求助】高斯安装后多用户使用、多核运行及服务器配置问题
已经有15人回复
- 【求助】水热法制备磷酸铁锂，大家帮我看看 问题出在哪
已经有35人回复
- 【求助】运行错误，大家看看问题在哪？
已经有7人回复
【答案】应助回帖
★ ★ ★ ★
感谢参与，应助指数 +1
hpuwyb: 金币+2 2013-05-10 10:41:04
ben_ladeng: 金币+2, 谢谢交流 2013-05-11 17:21:11
感谢参与，应助指数 +1
hpuwyb: 金币+2 2013-05-10 10:41:04
ben_ladeng: 金币+2, 谢谢交流 2013-05-11 17:21:11
|
楼主是在出题考大家吗？ 第18行应为： text(19,0.95,['\theta_{max}=', num2str(ThetaMaxDeg,4),'\circ']) 第19行： ThetaMax Vol是神马东东，中间加个空格表示什么？ fiminbnd是什么，自定义的函数？是fminbnd吧？ 第23行： text(10,MaxVol+0.05,[num2str(ThetaMaxVol*180/pi,4),'\circ']) 第36行： grav是神马东东？ [ Last edited by xiegangmai on 2013-5-9 at 22:50 ] |
2楼2013-05-09 22:42:44