-
-
手机阅读
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
-
-
-
-
lxw19881017 | 当前状态：离线
总积分：0 2024年可用积分：0
注册时间： 0001-01-01
最后登录时间： 0001-01-01
-
-
lxw19881017 发表于 2012/9/5 19:43:07
2楼 回复本楼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')
图形参见：
引用 lxw19881017 2012/9/5 19:43:07 发表于2楼的内容
-
