東北大學(xué)matlab實(shí)驗(yàn)參考題答案-資料下載頁(yè)
2025-06-23 19:02本頁(yè)面
【正文】 0 1 0 0。 0 1 0 0 1 0。 0 0 1 0 0 1]。B=[30。 18。 10。 18。 30]。 intlist=[1。1。1。1。1]。ctype=[0。0。0。1。1]。 xm=zeros(5,1)。 xM=inf*ones(5,1)。[res,b]=ipslv_mex(f,A,B,intlist,xM,xm,ctype)。 resres =08221008 Aeq=[1 1 1 0 0 0。 0 0 0 1 1 1。 1 0 0 1 0 0]。 Beq=[30。 18。 10]。A=[0 1 0 0 1 0。 0 0 1 0 0 1]。 B=[18。 30]。intlist=ones(6,1)。 xm=zeros(6,1)。 xM=20000*ones(6,1)。 x0=xm。[errmsg,f,x]=bnb20(39。exc6fun339。,x0,intlist,xm,xM,A,B,Aeq,Beq)。if length(errmsg)==0, x=round(x), endx =08221008 試求出微分方程的解析解通解，并求出滿足邊界條件的解析解。 syms xy=dsolve(39。D2y(21/x)*Dy+(11/x)*y=x^2*exp(5*x)39。,39。x39。) y = exp(x)*C2+exp(x)*log(x)*C1+1/216*Ei(1,6*x)*exp(x)+11/1296*exp(5*x)+5/216*exp(5*x)*x+1/36*x^2*exp(5*x) syms xy=dsolve(39。D2y(21/x)*Dy+(11/x)*y=x^2*exp(5*x)39。,...39。y(1)=sym(pi)39。,39。y(sym(pi))=139。,39。x39。) y = 1/1296*exp(x)*(1296*sym(pi)*exp(5)6*exp(6)*Ei(1,6)77)/exp(1)/exp(5)1/1296*exp(x)*log(x)*(1296*exp(1)*exp(5)+1296*exp(sym(pi))*sym(pi)*exp(5)6*exp(sym(pi))*exp(6)*Ei(1,6)77*exp(sym(pi))+6*exp(5*sym(pi))*exp(6*sym(pi))*Ei(1,6*sym(pi))*exp(1)*exp(5)+11*exp(5*sym(pi))*exp(1)*exp(5)+30*exp(5*sym(pi))*sym(pi)*exp(1)*exp(5)+36*exp(5*sym(pi))*sym(pi)^2*exp(1)*exp(5))/exp(sym(pi))/log(sym(pi))/exp(1)/exp(5)+1/1296*(6*exp(6*x)*Ei(1,6*x)+11+30*x+36*x^2)*exp(5*x) vpa(y,10) ans = .1912617421e5*exp(x)*(*sym(pi))*exp(x)*log(x)*(+*exp(sym(pi))*sym(pi)*exp(sym(pi))+*exp(5.*sym(pi))*exp(6.*sym(pi))*Ei(1.,6.*sym(pi))+*exp(5.*sym(pi))+*exp(5.*sym(pi))*sym(pi)+*exp(5.*sym(pi))*sym(pi)^2)/exp(sym(pi))/log(sym(pi))+.7716049383e3*(6.*exp(6.*x)*Ei(1.,6.*x)+11.+30.*x+36.*x^2)*exp(5.*x) 試求出下面微分方程的通解。（1） ；（2）。(1) syms t。x=dsolve(39。D2x+2*t*Dx+t^2*x=t+139。) x = exp(t1/2*t^2)*C2+exp(t1/2*t^2)*C11/2*i*pi^(1/2)*2^(1/2)*erf(1/2*i*2^(1/2)*(1+t))*exp(1/2+t1/2*t^2)(2) syms xy=dsolve(39。Dy+2*x*y=x*exp(x^2)39。,39。x39。) y = 1/2*(x^2+2*C1)*exp(x^2)1 考慮著名的化學(xué)反應(yīng)方程組，選定，且，繪制仿真結(jié)果的三維相軌跡，并得出其在xy平面上的投影。在實(shí)際求解中建議將作為附加參數(shù)，同樣的方程若設(shè)，時(shí)，繪制出狀態(tài)變量的二維圖和三維圖。 f=inline(39。[x(2)x(3)。 x(1)+a*x(2)。 b+(x(1)c)*x(3)]39。,...39。t39。,39。x39。,39。flag39。,39。a39。,39。b39。,39。c39。)。[t,x]=ode45(f,[0,100],[0。0。0],[],)。plot3(x(:,1),x(:,2),x(:,3))。 grid [t,x]=ode45(f,[0,100],[0。0。0],[],10)。plot3(x(:,1),x(:,2),x(:,3))。 grid1 試選擇狀態(tài)變量，將下面的非線性微分方程組轉(zhuǎn)換成一階顯式微分方程組，并用 MATLAB對(duì)其求解，繪制出解的相平面或相空間曲線。 f=inline([39。[x(2)。 x(1)x(3)(3*x(2))^2+(x(4))^3+6*x(5)+2*t。 39。,...39。x(4)。 x(5)。 x(5)x(2)exp(x(1))t]39。],39。t39。,39。x39。)。[t1,x1]=ode45(f,[1,0],[2, 4, 2, 7, 6]39。)。[t2,x2]=ode45(f,[1,2],[2, 4, 2, 7, 6]39。)。t=[t1(end:1:1)。 t2]。 x=[x1(end:1:1,:)。 x2]。plot(t,x)figure。 plot(x(:,1),x(:,3))1 考慮簡(jiǎn)單的線性微分方程，且方程的初值為，，試用Simulink搭建起系統(tǒng)的仿真模型，并繪制出仿真結(jié)果曲線。 [t,x,y]=sim(39。exc7m539。,[0,10])。 plot(t,x)figure。 plot(t,y)1 用生成一組較稀疏的數(shù)據(jù)，并用一維數(shù)據(jù)插值的方法對(duì)給出的數(shù)據(jù)進(jìn)行曲線擬合，并將結(jié)果與理論曲線相比較。 t=0::2。y=t.^2.*exp(5*t).*sin(t)。 plot(t,y,39。o39。) ezplot(39。t.^2.*exp(5*t).*sin(t)39。,[0,2])。 hold onx1=0::2。 y1=interp1(t,y,x1,39。spline39。)。plot(x1,y1)1. 若不給自己設(shè)限，則人生中就沒有限制你發(fā)揮的藩籬。2. 若不是心寬似海，哪有人生風(fēng)平浪靜。在紛雜的塵世里，為自己留下一片純靜的心靈空間，不管是潮起潮落，也不管是陰晴圓缺，你都可以免去浮躁，義無反顧，勇往直前，輕松自如地走好人生路上的每一步3. 花一些時(shí)間，總會(huì)看清一些事。用一些事情，總會(huì)看清一些人。有時(shí)候覺得自己像個(gè)神經(jīng)病。既糾結(jié)了自己，又打擾了別人。努力過后，才知道許多事情，堅(jiān)持堅(jiān)持，就過來了。4. 歲月是無情的，假如你丟給它的是一片空白，它還給你的也是一片空白。歲月是有情的，假如你奉獻(xiàn)給她的是一些色彩，它奉獻(xiàn)給你的也是一些色彩。你必須努力，當(dāng)有一天驀然回首時(shí)，你的回憶里才會(huì)多一些色彩斑斕，少一些蒼白無力。只有你自己才能把歲月描畫成一幅難以忘懷的人生畫卷。
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1