MATLAB實(shí)習(xí)一答案.doc

實(shí)習(xí)一 數(shù)值計(jì)算和符號(hào)計(jì)算1、“左除”與“右除”有什么區(qū)別？ 在通常情況下，左除x=ab是a*x=b的解，右除x=b/a是x*a=b的解，一般情況下，abb/a。 2、，在進(jìn)行邏輯運(yùn)算時(shí)，a相當(dāng)于什么樣的邏輯量。 =相當(dāng)于3、在sin(x)運(yùn)算中，x是角度還是弧度？ 在sin(x)運(yùn)算中，x是弧度，MATLAB規(guī)定所有的三角函數(shù)運(yùn)算都是按弧度進(jìn)行運(yùn)算。 4、角度，求x的正弦、余弦、正切和余切 x=30 45 60; x1=x/180*pi; sin(x1) ans = 0.5000 0.7071 0.8660 cos(x1) ans = 0.8660 0.7071 0.5000 tan(x1) ans = 0.5774 1.0000 1.7321 cot(x1) ans = 1.7321 1.0000 0.5774 5、用四舍五入的方法將數(shù)組2.4568 6.3982 3.9375 8.5042取整。 b=2.4568 6.3982 3.9375 8.5042; round(b) ans = 2 6 4 9 6、將矩陣、和組合成兩個(gè)新矩陣： =（1）組合成一個(gè)4x3的矩陣，第一列為按列順序排列的a矩陣元素，第二列為按列順序排列的b矩陣元素，第三列為按列順序排列的c矩陣元素，即 （2）按照a、b、c的列順序組合成一個(gè)行矢量，即 a=4 2;5 7; b=7 1;8 3; c=5 9;6 2; % (1) d=a(:) b(:) c(:) d = 4 7 5 5 8 6 2 1 9 7 3 2 % (2) e=a(:);b(:);c(:) e = 4 5 2 7 7 8 1 3 5 6 9 2 或利用（1）中產(chǎn)生的d e=reshape(d,1,12) ans = 4 5 2 7 7 8 1 3 5 6 9 2 7 將(x-6)(x-3)(x-8)展開為系數(shù)多項(xiàng)式的形式。 a=6 3 8; pa=poly(a); ppa=poly2sym(pa) ppa = x3-17*x2+90*x-144 8求解多項(xiàng)式x3-7x2+2x+40的根。 r=1 -7 2 40; p=roots(r); 5.0000 4.0000 -2.0000 9求解在x=8時(shí)多項(xiàng)式(x-1)(x-2) (x-3)(x-4)的值。 p=poly(1 2 3 4); polyvalm(p,8) ans = 840 10計(jì)算多項(xiàng)式乘法(x2+2x+2)(x2+5x+4)。 c=conv(1 2 2,1 5 4) c = 1 7 16 18 8 11計(jì)算多項(xiàng)式除法(3x3+13x2+6x+8)/(x+4)。 d=deconv(3 13 6 8,1 4) d = 3 1 2 12對(duì)下式進(jìn)行部分分式展開： a=1 3 4 2 7 2; b=3 2 5 4 6; r,s,k=residue(b,a) r = 1.1274 + 1.1513i 1.1274 - 1.1513i -0.0232 - 0.0722i -0.0232 + 0.0722i 0.7916 s = -1.7680 + 1.2673i -1.7680 - 1.2673i 0.4176 + 1.1130i 0.4176 - 1.1130i -0.2991 k = 13計(jì)算多項(xiàng)式的微分和積分。 p=4 12 14 5; pder=polyder(p); pders=poly2sym(pder) pint=polyint(p); pints=poly2sym(pint) pders = 12*x2-24*x-14 pints = x4-4*x3-7*x2+5*x 14解方程組。 a=2 9 0;3 4 11;2 2 6; b=13 6 6; x=ab x = 7.4000 -0.2000 -1.4000 15矩陣，計(jì)算a的行列式和逆矩陣。 a=4 2 -6;7 5 4 ;3 4 9; ad=det(a) ai=inv(a) ad = -64 ai = -0.4531 0.6562 -0.5937 0.7969 -0.8437 0.9062 -0.2031 0.1562 -0.0937 16 y=sin(x)，x從0到2p，Dx=0.02p，求y的最大值、最小值、均值和標(biāo)準(zhǔn)差。 x=0:0.02*pi:2*pi; y=sin(x); ymax=max(y) ymin=min(y) ymean=mean(y) ystd=std(y) ymax = 1 ymin = -1 ymean = 2.2995e-017 ystd = 0.7071 17在不加專門指定的情況下，以下符號(hào)表達(dá)式中的哪一個(gè)變量是獨(dú)立自由變量？sym(sin(w*t) wsym(a*exp(-X) asym(z*exp(j*theta) z18求以下符號(hào)矩陣的行列式值和逆，所得結(jié)果應(yīng)用“子表達(dá)式替換”簡潔化。A=sym(a11,a12,a13;a21,a22,a23;a31,a32,a33);det(A)ans =a11*a22*a33-a11*a23*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a22inv(A)ans =(a22*a33-a23*a32)/(a11*a22*a33-a11*a23*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a22), -(a12*a33-a13*a32)/(a11*a22*a33-a11*a23*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a22), (a12*a23-a13*a22)/(a11*a22*a33-a11*a23*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a22) -(a21*a33-a23*a31)/(a11*a22*a33-a11*a23*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a22), (a11*a33-a13*a31)/(a11*a22*a33-a11*a23*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a22), -(a11*a23-a13*a21)/(a11*a22*a33-a11*a23*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a22)(a21*a32-a22*a31)/(a11*a22*a33-a11*a23*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a22), -(a11*a32-a12*a31)/(a11*a22*a33-a11*a23*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a22), (a11*a22-a12*a21)/(a11*a22*a33-a11*a23*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a22)subexpr(ans,w)w = a11*a22*a33-a11*a23*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a22ans = (a22*a33-a23*a32)/w, -(a12*a33-a13*a32)/w, (a12*a23-a13*a22)/w -(a21*a33-a23*a31)/w, (a11*a33-a13*a31)/w, -(a11*a23-a13*a21)/w (a21*a32-a22*a31)/w, -(a11*a32-a12*a31)/w, (a11*a22-a12*a21)/w19(1)通過符號(hào)計(jì)算求的導(dǎo)數(shù)；(2)然后根據(jù)此結(jié)果，求和。（1）syms t y = abs(sin(t);d=diff(y,t)d =abs(1,sin(t)*cos(t)(2)subs(d,t,0)ans = 1subs(d,t,pi/2)ans = 020在區(qū)間，畫出曲線，并計(jì)算。syms t x zy = sin(t)/t;z = int(y,t,0,x);ezplot(z,x,0,2*pi)subs(z,x,4.5)ans = 1.6541404143792421求方程的解。x,y = solve(x2+y2=1,x*y=2,x,y)x = -1/2*(1/2*5(1/2)+1/2*i*3(1/2)3+1/4*5(1/2)+1/4*i*3(1/2) -1/2*(1/2*5(1/2)-1/2*i*3(1/2)3+1/4*5(1/2)-1/4*i*3(1/2) -1/2*(-1/2*5(1/2)+1/2*i*3(1/2)3-1/4*5(1/2)+1/4*i*3(1/2) -1/2*(-1/2*5(1/2)-1/2*i*3(1/2)3-1/4*5(1/2)-1/4*i*3(1/2)y = 1/2*5(1/2)+1/2*i*3(1/2) 1/2*5(1/2)-1/2*i*3(1/2) -1/2*5(1/2)+1/2*i*3(1/2) -1/2*5(1/2)-1/2*i*3(1/2)22求一階微分方程，的解。s = dsolve(Dx=a*t2+b*t,x(0)=2)s =1/3*a*t3+1/2*b*t2+223求邊值問題，的解。f,g = dsolve(Df=3*f+4*g,Dg=-4*f+3*g,f(0)=0,g