匿名用户
一级
2016-09-18答复
试试这个
clear all；
clc；
L1=1； L2=1； L=1； C1=1； C2=1； C=1； m=0.5； w=50；
a=sqrt(2)/2； b=sqrt(2)/2； d=sqrt(2)/2； c=-sqrt(2)/2；
k=1.381*10^(-23； h=6.63*10^(-34；
T=273；
b=1/(k*t )；
m1=1/(L-M )； m2=1/(Lm )；
R1=1/sqrt(C* ) L-M )； R2=1/sqrt(c* ) Lm )；
a1=0； b1=0；
t=0:0.1:10；
a2=-sqrt(2) ) cos(wR2 ) t )/(2) ) wR2 ) t )/(2) ) w-R2 ) ) w-R2 ) )1/(2) ) ) wR2 ) ) )
B2=--(sqrt(2) *M2/r2 ) ) sin ) ) wR2 ) )-sin ) ) w-R2 ) ) t ) )2) ) ) ) ) w ) w ) w ) w ) w ) w ) w )
a1=a1*cos(R1*T ) B1*R1*sin ) R1*T )/M1；
B1=B1*cos(R1*T )-A1*M1*sin ) R1*T )/r1；
A2=A2.*cos(R2*T ) B2*R2.*sin ) R2*T )/M2；
B2=B2.*cos(R2*T )-A2*M2.*sin ) R2*T )/r2；
q1=a*B1 b*B2； q2=c*B1 d*B2；
I1=(m* ) a*a2-b*a1 ) L2* ) c*a2-d*a1 ) )/)/(m^2-L1*L2 )；
I2=(m*(d*a1-c*a2 )-L1* ) a*a2-b*a1 ) )/)/(m^2-L1*L2 )；
k1=(a^2) h*m1/)2*R1*tanh ) h*R1*B/2 ) ) b^2) h*M2/)2*R2*tanh )2) a
k2=(c^2) h*m1/)2*R1*tanh ) h*R1*B/2 ) ) d^2) h*M2/(2*R2*tanh ) h*R2*B/2 ) )2) c
S1=K1-(Q1 ).^2；
S2=K2-(Q2 ).^2；
subplot (2，1，1 )； plot(t，I1 )； 保持打开； 网格开； title ()温度为273K时电流I1的热力学平均值随时间t变化的曲线)； xlabel(t/小时)； ylabel(I1 )；
subplot (2，1，2 )； plot(t，I2 )； 保持打开； 网格开； title ()温度为273K时电流I2的热力学平均值随时间t变化的曲线)； xlabel(t/小时)； ylabel(I2 )；