Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
![]() |
#1 |
初级会员
注册日期: 2010-05-12
帖子: 1
声望力: 0 ![]() |
![]()
QQ截图未命名.jpg
各位大虾 实在做不出来 在此请教 C0=1.2 A=60 Lp=5 B=C=125 AP=12.5 M=225.25 c=4.2 T0=25 m=40 Pabsorb=62 实现: 画出关于t Mt的图像 我的代码: syms t JD1=(10^4)*exp(AP-0.1351*(M^(2/3))+0.003*M-(10454*c*m)/(c*m*T0+Pabsorb*t)) %D(t)表达式 %JD=cumtrapz(t,JD1) JD=int(JD1,t) %对D(t)求积分 syms n G=symsum('exp((-(2*n+1)^2*pi^2*JD)/5^2)',n,0,inf) %G为累加求和的那一部分 %G=double(G) %vpa(G,7) JD2=JD1*G JD3=int(JD2,t) %JD3=cumtrapz(t,JD2) f=((4*A*C0)/Lp)*JD3-B*exp(t/C) %plot(t,f) %这里是重点 不知道怎么画 t的应该是虚数 =。= |
![]() |
![]() |