Matlab上機(jī)練習(xí)二答案

1、1、 求的共軛轉(zhuǎn)置。>> x=4+8i 3+5i 2-7i 1+4i 7-5i;3+2i 7-6i 9+4i 3-9i 4+4i; >> xans = 4.0000 - 8.0000i 3.0000 - 2.0000i 3.0000 - 5.0000i 7.0000 + 6.0000i 2.0000 + 7.0000i 9.0000 - 4.0000i 1.0000 - 4.0000i 3.0000 + 9.0000i 7.0000 + 5.0000i 4.0000 - 4.0000i2、計(jì)算與的數(shù)組乘積。>> a=6 9 3;2 7 5;>>

2、b=2 4 1;4 6 8;>> a.*bans = 12 36 3 8 42 403、 對(duì)于，如果，求解X。>> A=4 9 2;7 6 4;3 5 7;>> B=37 26 28;>> X=ABX = -0.5118 4.0427 1.33184、 ，觀察a與b之間的六種關(guān)系運(yùn)算的結(jié)果。>> a=1 2 3;4 5 6;>> b=8 7 4;3 6 2;>> a>bans = 0 1 0 1 0 1>> a>=bans = 0 1 0 1 0 1>> a<bans

3、= 1 0 1 0 1 0>> a<=bans = 1 0 1 0 1 0>> a=bans = 0 0 0 0 0 0>> a=bans = 1 1 1 1 1 15、，在進(jìn)行邏輯運(yùn)算時(shí)，a相當(dāng)于什么樣的邏輯量。相當(dāng)于a=1 1 0 1 1。6、 角度，求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>> ta

4、n(x1)ans = 0.5774 1.0000 1.7321>> cot(x1)ans = 1.7321 1.0000 0.57747、 用四舍五入的方法將數(shù)組2.4568 6.3982 3.9375 8.5042取整。>> b=2.4568 6.3982 3.9375 8.5042;>> round(b)ans = 2 6 4 98、設(shè)，求C1=A*B;C2=A*B;C3=A.*B,并求上述所有方陣的逆陣。>> A=1 4 8 13;-3 6 -5 -9;2 -7 -12 -8;>> B=5 4 3 -2;6 -2 3 -8;-1

5、 3 -9 7;>> C1=A*B'C1 = 19 -82 30 12 27 3 -38 54 29>> C2=A'*BC2 = -15 16 -24 36 63 -17 93 -105 22 6 117 -60 19 46 84 -10>> C3=A.*BC3 = 5 16 24 -26 -18 -12 -15 72 -2 -21 108 -56>> inv(C1)ans = 0.0062 0.0400 -0.0106 -0.0046 0.0169 0.0030 0.0168 0.0209 0.0150>> inv(

6、C2)Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 8.997019e-019.ans = 1.0e+015 * -0.9553 -0.2391 -0.1997 0.2700 0.9667 0.2420 0.2021 -0.2732 -0.4473 -0.1120 -0.0935 0.1264 -1.1259 -0.2818 -0.2353 0.3182>> inv(C3)? Error using => invMatrix must be sq

7、uare.9、設(shè)x=rcost+3t,y=rsint+3,分別令r=2,3,4，畫出參數(shù)t=010區(qū)間生成的xy曲線。>> t=linspace(0,10);>> r1=2;>> x1=(r1*cos(t)+3*t);>> y1=r1*sin(t)+3;>> r2=3;>> x2=(r2*cos(t)+3*t);>> y2=r2*sin(t)+3;>> r3=4;>> x3=(r3*cos(t)+3*t);>> y3=r3*sin(t)+3;>> plot(x1,

8、y1,'r',x2,y2,'b',x3,y3,'m')10、設(shè)f(x)=x5- 4x4 +3x2- 2x+ 6(1) 在x=-2,8之間取100個(gè)點(diǎn)，畫出曲線，看它有幾個(gè)過零點(diǎn)。（提示：用polyval 函數(shù)）>> x=linspace(2,8,100);>> y=polyval(1 0 -4 3 -2 6,x);>> plot(x,y,'b',x,0,'y')(2) 用roots函數(shù)求此多項(xiàng)式的根。t=1 0 -4 3 -2 6p=roots(t)12、 求解多項(xiàng)式x3-7x2

9、+2x+40的根。 >> r=1 -7 2 40; >> p=roots(r); -0.2151 0.4459 0.7949 0.270713、符號(hào)函數(shù)繪圖法繪制函數(shù)x=sin(3t)cos(t)，y=sin(3t)sin(t)的圖形，t的變化范圍為0,2p。 >> syms t >> ezplot(sin(3*t)*cos(t),sin(3*t)*sin(t),0,pi)14、設(shè),求x=sym('x'); y=(sin(x)4+(cos(x)4;diff(y,10)15、x=sym('x');int(sqrt(9

10、-x2)3)/x6),x)16、x=sym('x');int(x+x3,x,1,2)17、求級(jí)數(shù)的和: sym(n);symsum(n+2/2n,n,1,inf)18、利用函數(shù)int計(jì)算二重不定積分x=sym('x'); y=sym('y');z=(x+y)*exp(-x*y);a=int(z,x);int(a,y)19、試求出如下極限。（1）； （2）；（1）>> syms x;f=(3x+9x)(1/x);l=limit(f,x,inf) ans= 9（2）>> syms x y;f=x*y/(sqrt(x*y+1)-

11、1);limit(limit(f,x,0),y,0) ans = 2（3）>> syms x y;f=(1-cos(x2+y2)*exp(x2+y2)/(x2+y2);limit(limit(f,x,0),y,0) ans = 020、已知參數(shù)方程，試求出和>> syms t;x=log(cos(t);y=cos(t)-t*sin(t);diff(y,t)/diff(x,t) ans = -(-2*sin(t)-t*cos(t)/sin(t)*cos(t) >> f=diff(y,t,2)/diff(x,t,2);subs(f,t,sym(pi)/3) an

12、s = 3/8-1/24*pi*3(1/2)21、假設(shè)，試求>> syms x y t>> s=int(exp(-t2),t,0,x*y);>> x/y*diff(f,x,2)-2*diff(diff(f,x),y)+diff(f,y,2) ans =2*x2*y2*exp(-x2*y2)-2*exp(-x2*y2)-2*x3*y*exp(-x2*y2)22試求出下面的極限。（1）；（2）>> syms k n;symsum(1/(2*k)2-1),k,1,inf) ans = 1/2>> limit(symsum(1/(2*k)2-

13、1),k,1,n),n,inf) ans = 1/2（2）>> limit(n*symsum(1/(n2+k*pi),k,1,n),n,inf) ans = 123、假設(shè)一曲線數(shù)據(jù)點(diǎn)為x=0:2:4*pi y=sin(x).*exp(-x/5)試將x的間距調(diào)成0.1，并用下列方法進(jìn)行內(nèi)插：（1） 線性內(nèi)插法（method=linear）（2） 樣條內(nèi)插法（method=spline）（3） 三次多項(xiàng)式內(nèi)插法（method=cubic）（4） 多項(xiàng)式擬合法：直接利用6次多項(xiàng)式去通過7個(gè)數(shù)據(jù)點(diǎn)請(qǐng)將這些內(nèi)插法的結(jié)果及原先的數(shù)據(jù)點(diǎn)畫在同一個(gè)圖上25、閱讀以下程序并在指定位置上寫注釋（6分）

14、t=(0:0.02:2)*pi;x=sin(t);y=cos(t);z=cos(2*t);plot3(x,y,z,'r-',x,y,z,'bd') %_(1)_繪制一條紅色連續(xù)曲線和藍(lán)色菱形離散點(diǎn)view(-82,58) %_(2)_以方位角-82度，俯角58度觀察視圖box on %顯示坐標(biāo)軸的矩形框legend('圖形','視角') %_ (3)_標(biāo)注圖例_hold on %_(4)_保持圖形_x=-4:4;y=x;X,Y=meshgrid(x,y); %_(5)_生成平面網(wǎng)格坐標(biāo)矩陣_Z=X.2+Y.2;surf(X,Y,Z

15、); %_ (6)_繪制三維曲面_26、編寫一個(gè)M函數(shù)文件fun_es(x)，計(jì)算如下函數(shù)：，其中參數(shù)可以為標(biāo)量，也可以為向量。在Matlab指令窗口鍵入如下命令檢驗(yàn)該函數(shù)： fun_es(3)function fun_es(x)y=0.5.*exp(x./3)-x.2.*sin(x)27、編寫m函數(shù)文件求半徑為r的圓的面積和周長(zhǎng)。function f=circle(r)S=pi*r2L=2*pi*rx,yMatlab上機(jī)實(shí)驗(yàn)四28、編制m文件，等待鍵盤輸入，輸入密碼123，密碼正確，顯示輸入密碼正確，程序結(jié)束；否則提示，重新輸入。function f=zy2()t=input('In

16、put key(key is 123) ')if t=123 'good'else 'You are fool, input again' zy2; /從新輸入；end4.3已知y=。當(dāng)n=100時(shí)，求y的值步驟：建立函數(shù)，命名為Cbb.m求平方倒數(shù)的和。其程序如下function F=Cbb(n)F=0;k=1;while k<=n F=F+(1/k)/k; k=k+1;endend在工作空間中求y值 y= Cbb(100)y=1.63504.5編制程序，計(jì)算1+2+n<2000時(shí)的最大n值步驟：建立函數(shù)文件Dbb.m,求累加和小于n的最大

17、數(shù)，使n=2000程序如下function k=Dbb(n)k=1;f=0;while f<n f=f+k; k=k+1;endk=k-2;end將2000作為輸入值，求n.在工作空間中輸入如下： Dbb(2000)ans =624.6分別用for和while循環(huán)結(jié)構(gòu)編寫程序，求出K=分別用for和while循環(huán)結(jié)構(gòu)編寫命令文件Ebb,Fbb.并在工作空間中求KEbb程序如下f=1;J(1)=1;for i=1:63J(i+1)=J(i)*2;f=f+J(i+1);End求K程序：Ebb.m ff = 1.8447e+019Fbb程序如下I=1;H=1;i=1;while i<=6

18、3 H=H*2; I=I+H; i=i+1;end在工作空間中求KFbb.mII = 1.8447e+0191.倘若要是生成一個(gè)4行5列的矩陣a=rand(4,5)2.倘若要是對(duì)x進(jìn)行賦值，從5到25，間隔是0.1x=5:0.1:25;倘若要是對(duì)x進(jìn)行賦值，從2到20，中間間隔100個(gè)點(diǎn)x=linspace(2,20,100);4. 用3種方法將圖畫在同一個(gè)窗口中x=5:0.1:25;y1=sin(x).*cos(2*x);y2=sin(2*x).*cos(x);plot(x,y1,x,y2)有一點(diǎn)必須注意：在式子與式子之間，必須要用點(diǎn)乘要想對(duì)所畫的圖形進(jìn)行修飾，可以用以下的語句，要點(diǎn)是必須是

19、先畫出圖形，再寫有關(guān)對(duì)圖形修飾的語句plot(x,y1,x,y2);plot(x,y1,'r-',x,y2,'b-.');xlabel('橫坐標(biāo)');title('標(biāo)題');legend('sinx*cosx','sin2x*cosx');("加示例")怎樣改變線寬？以及怎樣給一個(gè)中間空著的正方形的邊緣及內(nèi)部填充顏色plot(x,y,'markeredge','g');plot(x,y,'markerfacecolor','

20、b');5關(guān)于三維圖的畫法 餅圖 pie 火柴桿圖 stem(x,y) 階梯圖6. 關(guān)于排序問題a=rand(4,5); 隨機(jī)生成4*5的矩陣sort(a); 按列排序，且每列由上到下是由小到大sort(a,1); 按列排序，且每列由上到下是由小到大sort(a,2); 按行排序，且每行由左到右是由小到大倘若想要改變由大到小還是想改變由小到大，需要用到增長(zhǎng) ascend減小 descend例如：按行進(jìn)行排序，且每一行自左至右是遞減的,即：按照行從高到低排序a=1 2 3;4 5 6;7 8 9sort(a,2,'descend')a = 1 2 3 4 5 6 7 8

21、9ans = 3 2 1 6 5 4 9 8 77. 關(guān)于矩陣問題求矩陣中每一列的平均值mean(a)mean(a,1)求矩陣中每一行的平均值mean(a,2)從矩陣a中找出最大值max(max(a) 內(nèi)層的max相當(dāng)于找出每一列的最大值，外層的max相當(dāng)于找出每一列最大值的最大值 或者 a=rand(3,4)max(a(:) 即a的所有行，所有列倘若對(duì)矩陣a前兩行的 2，3 ,4 列求最大值max(a(1:2,2:4)求矩陣a的行列式值|a|b=det(a) det語句是求n*n的行列式的值求矩陣a的逆矩陣b=inv(a) inv是矩陣求逆的意思求矩陣a的轉(zhuǎn)置b=a'對(duì)于矩陣a:求出矩陣a中每一列中的最大值，最小值a=1 2 3;4 5 6;7 8 9b=max(a)c=min(a)求出矩陣a中每一行中的最大值，最小值a=1 2 3;4 5 6;7 8 9b=max(a,2)c=min(a,2)8. 關(guān)于取整問題a=5/2b=ceil(a) 向大取c=floor(a) 向小取9. 關(guān)于多形式的乘法和除法 例如：求x2+2x+3與4x2+5x+6的乘積 a=1 2 3b=4 5 6c=conv(a,b)求x2+2x+