2012-09-05
60=7.18*0.001*dI/dt+1/36.8*30/pi*df/dt*93+1.4*I
8.74/93*(d2f/dt2)*1.15=259*0.001*I
初始条件:
t=0时
f=15/180*pi;
df/dt=11.4028;
I=0
评论1
楼主 2012/9/5 19:43:07
clc;clear
[I,f]=dsolve('60=7.18*0.001*DI+1/36.8*30/pi*Df*93+1.4*I','8.74/93*D2f*1.15=259*0.001*I','f(0)=15/180*pi','Df(0)=11.4028','I(0)=0')
t=0:0.001:0.1;
I=subs(I);f=subs(f);
[AX,H1,H2]= plotyy(t,I,t,f)
xlabel('t')
set(get(AX(1),'Ylabel'),'String','I')
set(get(AX(2),'Ylabel'),'String','\phi')
figure(2)
subplot(2,1,1)
plot(t,I),xlabel('t'),ylabel('I')
subplot(2,1,2)
plot(t,f),xlabel('t'),ylabel('\phi')
图形参见：
