wanjiadi
2010-05-12, 16:36
2487
各位大虾 实在做不出来 在此请教
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的应该是虚数 =。=
各位大虾 实在做不出来 在此请教
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的应该是虚数 =。=