東北大學(xué)matlab實(shí)驗(yàn)參考題答案(參考版)
2025-06-26 19:02本頁(yè)面
【正文】 。你必須努力，當(dāng)有一天驀然回首時(shí)，你的回憶里才會(huì)多一些色彩斑斕，少一些蒼白無(wú)力。4. 歲月是無(wú)情的，假如你丟給它的是一片空白，它還給你的也是一片空白。既糾結(jié)了自己，又打擾了別人。用一些事情，總會(huì)看清一些人。2. 若不是心寬似海，哪有人生風(fēng)平浪靜。)。 y1=interp1(t,y,x1,39。,[0,2])。) ezplot(39。 plot(t,y,39。 t=0::2。 plot(t,x)figure。exc7m539。 plot(x(:,1),x(:,3))1 考慮簡(jiǎn)單的線性微分方程，且方程的初值為，試用Simulink搭建起系統(tǒng)的仿真模型，并繪制出仿真結(jié)果曲線。 x2]。 t2]。)。)。)。,39。],39。 x(5)。,...39。 x(1)x(3)(3*x(2))^2+(x(4))^3+6*x(5)+2*t。 f=inline([39。plot3(x(:,1),x(:,2),x(:,3))。0。plot3(x(:,1),x(:,2),x(:,3))。0。)。,39。,39。,39。,39。,39。,...39。 x(1)+a*x(2)。 f=inline(39。) y = 1/2*(x^2+2*C1)*exp(x^2)1 考慮著名的化學(xué)反應(yīng)方程組，選定，且，繪制仿真結(jié)果的三維相軌跡，并得出其在xy平面上的投影。,39。) 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。x=dsolve(39。（1） ；（2）。x39。y(sym(pi))=139。y(1)=sym(pi)39。D2y(21/x)*Dy+(11/x)*y=x^2*exp(5*x)39。x39。D2y(21/x)*Dy+(11/x)*y=x^2*exp(5*x)39。if length(errmsg)==0, x=round(x), endx =08221008 試求出微分方程的解析解通解，并求出滿足邊界條件的解析解。exc6fun339。 x0=xm。 xm=zeros(6,1)。 30]。 0 0 1 0 0 1]。 10]。 Beq=[30。 0 0 0 1 1 1。[res,b]=ipslv_mex(f,A,B,intlist,xM,xm,ctype)。 xm=zeros(5,1)。1。0。1]。1。 intlist=[1。 18。 18。 0 0 1 0 0 1]。 1 0 0 1 0 0。A=[1 1 1 0 0 0。 function y=exc6fun2(x)y=(592*x(1)+381*x(2)+273*x(3)+55*x(4)+48*x(5)+37*x(6)+23*x(7))。,ff)Maximum number of function evaluations exceeded。,x0,A,B,Aeq,Beq,xm,xM,39。x=fmincon(39。 =1e10。x0=(xm+xM)/2。 xM=[10。 xm=[10。 Aeq=[]。 A=[]。 x(1)*x(2)x(1)x(2)+。 function [c,ce]=exc6fun6a(x)ce=[]。exc6ff39。 y=int((exp(x)c*x)^2,x,0,1) y = 1/22*c+1/2*exp(2)+1/3*c^2function y=exc6ff(c)y=1/2*exp(1)^2+1/3*c^21/22*c。) 試求出使得取得極小值的值。),(2) ezsurf(39。(1) ezplot(39。 F=ztrans(f) F = 1/a*(a*T*z/(z1)^2z/(z1)+z/exp(a*T)/(z/exp(a*T)1)) iztrans(F) ans = ((1/exp(a*T))^n1+a*T*n)/a 用數(shù)值求解函數(shù)求解下述一元和二元方程的根，并對(duì)得出的結(jié)果進(jìn)行檢驗(yàn)。 F=ztrans(f) F = T^2*z*exp(a*T)*(z+exp(a*T))/(zexp(a*T))^3 f1=iztrans(F) f1 = T^2*(1/exp(a*T))^n*n^2(3) syms a k T。 F=ztrans(f) F = (zcos(a*T))*z/(z^22*z*cos(a*T)+1) f1=iztrans(F) f1 = cos(a*T*n)(2) syms k T a。(1) syms k a T。 F=fourier(f) F = 2*pi*(4*i*pi*dirac(3,w)4*pi^2*dirac(2,w)+dirac(4,w)) ifourier(F) ans = x^2*(2*pi+x)^2 請(qǐng)將下述時(shí)域序列函數(shù)進(jìn)行Z變換，并對(duì)結(jié)果進(jìn)行反變換檢驗(yàn)。 F=fourier(f) F = 6*(4+pi^2*dirac(2,w)*w^4)/w^4 ifourier(F) ans = x^2*(4*x*heaviside(x)+3*pi+2*x)(2) syms t。(1) syms x。 ilaplace(F) ans = 1/t*(exp(b*t)exp(a*t)) 試求出下面函數(shù)的Fourier變換，對(duì)得出的結(jié)果再進(jìn)行Fourier反變換，觀察是否能得出原來(lái)函數(shù)。 ilaplace(F) ans = 1/2/t^(3/2)/pi^(1/2)*(exp(b*t)exp(a*t))(3) syms a b s。 ilaplace(F) ans = 1/2/b^2/a^3/(a^2b^2)*(2*t*a*b^3+2*(1b*texp(b*t))*a^3+(2*a+exp(a*t)*(ab)+(a+b)*exp(a*t))*b^2)(2) syms s a b。(1) syms s a b。 laplace(f) ans = 20160/(si*a)^9+20160/(s+i*a)^9 對(duì)下面的式進(jìn)行Laplace反變換。 laplace(f) ans = 60*i*(1/(si*a)^6+1/(s+i*a)^6)(3) syms t a。 laplace(f) ans = atan(a/s)(2) syms t a。(1) syms a t。練習(xí)題： 對(duì)下列的函數(shù)進(jìn)行Laplace變換。syms x t。 0,1,1,3]。 ,4,。syms x t。 0,1,1,3]。 ,4,。 syms t。 0,1,1,3]。 ,4,。1 假設(shè)已知矩陣如下，試求出。 2,1,2,2]。 0,。 A=vander([a b c d e])A =[ a^4, a^3, a^2, a, 1][ b^4, b^3, b^2, b, 1][ c^4, c^3, c^2, c, 1][ d^4, d^3, d^2, d, 1][ e^4, e^3, e^2, e, 1]det(A), simple(ans)ans =(
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1