如何用matlab求解二阶微分方程,以和程序实例.docVIP

如何用matlab求解二阶微分方程，以及程序实例 (2010-03-25 12:56:36)  HYPERLINK javascript:; 转载 标签：  HYPERLINK /c.php?t=blogk=matlabts=bpoststype=tag \t _blank matlab  HYPERLINK /c.php?t=blogk=%B3%CC%D0%F2%CA%B5%C0%FDts=bpoststype=tag \t _blank 程序实例  HYPERLINK /c.php?t=blogk=%B6%FE%BD%D7%CE%A2%B7%D6%B7%BD%B3%CCts=bpoststype=tag \t _blank 二阶微分方程  HYPERLINK /c.php?t=blogk=dsolvets=bpoststype=tag \t _blank dsolve  HYPERLINK /c.php?t=blogk=%CE%C4%BB%AFts=bpoststype=tag \t _blank 文化微分方程的求解主要用到一个dsolve函数，如下面的“Uc1=dsolve(D2y+1000*Dy+10^6*y=0,y(0)=10,Dy(0)=0,t);”，可以看出，函数的第一部分是所要求解的微分方程，其次是初始条件，最后是对自变量的说明。 下面给出的程序实例是用于分析一个最简单零输入的二阶电路。其中C=1uf,L=1H。R是不确定的，他的值的选取将会直接影响到方程解的形式以及最后画出的曲线形状，在此我取R的值分别为1000,2000,3000欧姆。 R=1000;while (R=3000)if R2000Uc1=dsolve(D2y+1000*Dy+10^6*y=0,y(0)=10,Dy(0)=0,t);It1=-1*diff(Uc1)*(1e-6);Ul1=diff(It1);elseif R==2000Uc2=dsolve(D2y+2000*Dy+10^6*y=0,y(0)=10,Dy(0)=0,t);It2=-1*diff(Uc2)*(1e-6);Ul2=diff(It2);elseUc3=dsolve(D2y+3000*Dy+10^6*y=0,y(0)=10,Dy(0)=0,t);It3=-1*diff(Uc3)*(1e-6);Ul3=diff(It3);endR=R+1000;end%while i=3figure (1)xlabel(t)ylabel(Uc(t))hold on;p1=ezplot(Uc1,[0,0.02]);p2=ezplot(Uc2,[0,0.02]);p3=ezplot(Uc3,[0,0.02]);title(电容电压Uc(t)变化曲线)set(p1,LineStyle,:);set(p2,LineStyle,--);legend(uc1,uc2,uc3);set(p1,Color,r);set(p2,Color,k);legend(uc1,uc2,uc3);hold off;axis auto;figure (2)xlabel(t)ylabel(I(t))hold on;p4=ezplot(It1,[0,0.02]);p5=ezplot(It2,[0,0.02]);p6=ezplot(It3,[0,0.02]);title(电流I（t）变化曲线)set(p4,LineStyle,:);set(p5,LineStyle,--);legend(uc1,uc2,uc3);set(p4,Color,r);set(p5,Color,k);legend(It1,It2,It3);hold off;axis auto;figure (3)xlabel(t)ylabel(Ul(t))hold on;p7=ezplot(Ul1,[0,0.02]);p8=ezplot(Ul2,[0,0.02]);p9=ezplot(Ul3,[0,0.02]);title(电感电压Ul(t)变化曲线)set(p7,LineStyle,:);set(p8,LineStyle,--);legend(uc1,uc2,uc3);set(p7,Color,r);set(p8,Color,k);legend(Ul1,Ul2,Ul3);hold off;axis auto;