for i=1:1:100
y(i)=sin(pi*i/10);
end
t=1:1:100;
plot(t,y,'r.-');for i=0:0.01:1I=40I1=dsolve('DI1=(1/7200)*(40-I1)','I1(0)=0','i')%I1ÐèÒª¸ø³õÖµ0Q=dsolve('DQ=-40','Q(0)=300','i')%QÐèÒª¸ø³õÖµ300DOC=1-Q*(1-0.11*(I/51.5)^1.75)/352.8DOD=1-DOC这程序还有一半发不过去了 能不能把你qq给我 我给你发过去帮我看看 这个分给你了谢谢啦