有如下程序:t=0:0.001:0.2sym Tgb=870000,Vm=0.00249,Cp=20.8,f=0.3,Vt=0.0000325,At=0.0005,Ac=0.0028,R=8.3145,P1=100000,P0= ,Th=300,Tl=60,L=0.06,
有如下程序:t=0:0.001:0.2sym Tgb=870000,Vm=0.00249,Cp=20.8,f=0.3,Vt=0.0000325,At=0.0005,Ac=0.0028,R=8.3145,P1=100000,P0= ,Th=300,Tl=60,L=0.06,F=20,Tr=300,Lt=0.065M=(Th-Tl)/Lw=2*pi*3.14Pt=P0+P1.*sin(w*t)j=(w.*Pt*Vt./(6*R*Ac*Tg)).*sin(w*t)+(w*P1*At*Lt./(1.6*R*Ac*Tg)).*cos(w*t)Tg=dsolve('DTg+(b*Vm/((1-f)*Cp))*Tg+j*Vm*M/(1-f)-(P1*w*Vm/Cp)*cos(w*t)-b*Vm*Tr/((1-f)*Cp)=0','Tg(0)=Th','t')simplify(Tg)目的是求解Tg和t的关系,解出来的Tg=exp(b*Vm/(-1+f)/Cp*t)*(-Tr*b^3*Vm^2+Vm^2*w*P1*b^2*f+j*M*Cp*b^2*Vm^2+2*Tr*b*w^2*Cp^2*f+j*M*Cp^3*w^2*f^2-Tr*b*w^2*Cp^2*f^2-Tr*b*w^2*Cp^2-2*j*M*Cp^3*w^2*f+j*M*Cp^3*w^2-Vm^2*w*P1*b^2+Th*b^3*Vm^2+Th*b*w^2*Cp^2-2*Th*b*w^2*Cp^2*f+Th*b*w^2*Cp^2*f^2)/(b^2*Vm^2+w^2*Cp^2-2*w^2*Cp^2*f+w^2*Cp^2*f^2)/b-(j*M*Cp^3*w^2*f^2-Tr*b*w^2*Cp^2*f^2-Vm*w^2*P1*Cp*sin(w*t)*b*f^2+2*Tr*b*w^2*Cp^2*f-2*j*M*Cp^3*w^2*f+2*Vm*w^2*P1*Cp*sin(w*t)*b*f-Tr*b^3*Vm^2-Tr*b*w^2*Cp^2+Vm^2*w*P1*b^2*cos(w*t)*f+j*M*Cp*b^2*Vm^2-Vm^2*w*P1*b^2*cos(w*t)+j*M*Cp^3*w^2-Vm*w^2*P1*Cp*sin(w*t)*b)/(b^2*Vm^2+w^2*Cp^2-2*w^2*Cp^2*f+w^2*Cp^2*f^2)/b之前我已经对某些参数赋值了,为什么matlab不计算出来最后结果呢?一大堆东西很烦人
展开
全部