matlab70x課后習(xí)題答案

1、1、利用基本矩陣產(chǎn)生3*3和15*8的單位矩陣、全1矩陣、全0矩陣、均勻分布隨機(jī)陣（-1,1之間）、正態(tài)分布隨機(jī)陣（均值為1，方差為4）。 解：A1=eye(3); A2=ones(3); A3=zeros(3); A4=2*rand(3)-1; A5=2*randn(3)+1; B1=eye(15,8); B2=ones(15,8); B3=zeros(15,8); B4=2*rand(15,8)-1; B5=2*randn(15,8)+1; 結(jié)果：由于數(shù)據(jù)是隨機(jī)產(chǎn)生的，所以在沒(méi)有給出運(yùn)行結(jié)果。2、利用diag等函數(shù)產(chǎn)生下列矩陣：a=0 0 8;0 -7 5;2 3 0 b=2 0 4;0

2、5 0;7 0 8然后利用reshape函數(shù)將它們變換成行向量。解：產(chǎn)生a的程序：b=diag(8 -7 2);c=b+diag(5 3,-1);a=fliplr(c)產(chǎn)生b的程序：s=2 2 8;t=4 3 7;v=diag(s);p=diag(t)+fliplr(v);b=fliplr(p)運(yùn)行結(jié)果：a = 0 0 8 0 -7 5 2 3 0b = 2 0 4 0 5 0 7 0 8利用reshape函數(shù)將它們變換成行向量：reshape(a,1,9)ans = 0 0 2 0 -7 3 8 5 03、產(chǎn)生一均勻分布在（-5,5）之間的隨機(jī)陣（50*2），要求精確到小數(shù)點(diǎn)后一位。解： A

3、=5-round(100*rand(50,2)/10部分?jǐn)?shù)據(jù)結(jié)果：A = 2.4000 4.2000 -0.1000 2.7000 -4.6000 -3.3000 -0.5000 -0.40003.5000 4.20004、編程實(shí)現(xiàn)當(dāng)t-,間隔為1時(shí)求解正弦和余弦值。解：t=(-1*pi:1/180:pi);y1=sin(t)y2=cos(t)部分?jǐn)?shù)據(jù)結(jié)果： Columns 10 through 18（y1）-0.0500 -0.0555 -0.0611 -0.0666 -0.0722 -0.0777 -0.0832 -0.0888 -0.09435、利用rand函數(shù)產(chǎn)生（0,1）間的均勻分布

4、的10*10隨機(jī)矩陣A，然后統(tǒng)計(jì)A中大于等于0.6的元素的個(gè)數(shù)。解：A=rand(10);B=A = 0.6;C=sum(B);count=sum(C)運(yùn)行結(jié)果（每次運(yùn)行結(jié)果是不同的，僅作參考）：count=326、利用randn函數(shù)產(chǎn)生均值為0，方差為1的10*10隨機(jī)矩陣A，然后統(tǒng)計(jì)A中大于-0.5且小于0.5的元素的個(gè)數(shù)。解：A=randn(10);B=(A-0.5);C=sum(sum(B)運(yùn)行結(jié)果（每次運(yùn)行結(jié)果是不同的，僅作參考）：C=48解：if and(a1,b=0.5)語(yǔ)句1;elseif and(a0.5)語(yǔ)句2;elseif and(a=1,b=0.5)語(yǔ)句3;else語(yǔ)句

5、4;有一矩陣A，找出矩陣中值等于1的元素，并將它們重新排列成列向量B。解：A=2*rand(4);k=find(A100);A(i,:)= %刪去存在奇異值的行運(yùn)行結(jié)果（每次運(yùn)行結(jié)果是不同的，僅作參考）：A = 49.5355 -23.7550 -73.0269 -118.4354 39.3214 -88.4472 91.1482 -28.5962 -209.9855 21.2336 -74.0239 -9.5871 -15.8184 72.1322 22.0444 71.4429 11.0770 34.8948在給定的100*100矩陣中，刪去整行為0的行，刪去整列為0的列。解：A=diag

6、(1 2 3 4,1)B=any(A)i,j=find(B=0)A(:,i)= %刪除全為0的列B=any(A)i,j=find(B=0)A(j,:)= %刪除全為0的行運(yùn)行結(jié)果：初始值：A = 0 1 0 0 0 0 0 2 0 0 0 0 0 3 0 0 0 0 0 4 0 0 0 0 0操作后：A = 1 0 0 0 0 2 0 0 0 0 3 0 0 0 0 4將窗口分割成四格，分別繪制正弦、余弦、正切和余切函數(shù)曲線，并加上適當(dāng)?shù)臉?biāo)注。程序?yàn)椋簒=0:pi/50:2*pi;k=1 26 51 76 101;x(k)=;%刪除正切和余切的奇異點(diǎn)figure(1)subplot(2,2,1

7、)plot(x,sin(x),k-),grid onlegend(ity=sin(x)title(y=sin(x)xlabel(x), ylabel(y) subplot(2,2,2)plot(x,cos(x),r-),grid onlegend(ity=cos(x)title(y=con(x)xlabel(x), ylabel(y) subplot(2,2,3)plot(x,tan(x),k),grid onlegend(ity=tan(x)title(y=tan(x)xlabel(x), ylabel(y) subplot(2,2,4)plot(x,cot(x),b-),grid onle

8、gend(ity=cot(x)title(y=cot(x)xlabel(x), ylabel(y)運(yùn)行如下：2、繪制多峰函數(shù)peaks和三角函數(shù)多條曲線。 多峰函數(shù)peaks：x,y=meshgrid(-3:0.15:3);z=peaks(x,y);x1=x(1,:);figure(1)plot(x1,z),grid ontitle(二維多峰函數(shù))圖形為：x,y=meshgrid(-3:0.15:3);z=peaks(x,y);figure(1)plot3(x,y,z),grid ontitle(三維多峰函數(shù))三角函數(shù)多條曲線： 程序?yàn)椋簍=-pi:pi/20:pi;y1=sinh(t); %

9、雙曲正弦y2=cosh(t); %雙曲余弦figure(1)subplot(2,1,1)plot(t,y1,r-,t,y2,k-),grid onlegend(ity1=sinh(t),ity2=cosh(t)title(三角函數(shù)1)xlabel(t), ylabel(y) subplot(2,1,2)plot(t,sin(t),k-),grid onhold on %保持原有圖像函數(shù)plot(t,cos(t),r-)legend(ity2=cos(t),ity1=sin(t)title(三角函數(shù)2)xlabel(t), ylabel(y) 運(yùn)行圖形為：3、將圖形窗口分成兩個(gè)，分別繪制以下函數(shù)

10、在-3,3區(qū)間上的曲線，并利用axis調(diào)整軸刻度，使他們具有相同縮放尺度。y1=2x+5； y2=x2-3x+1。程序?yàn)椋簒=-3:0.1:3;y1=2*x+5; y2=x.2-3*x+1;figure(1)subplot(2,2,1)plot(x,y1,r-),grid onlegend(ity1=2*x+5)title(y1=2x+5)xlabel(x), ylabel(y1)subplot(2,2,2)plot(x,y2,k-),grid onlegend(ity2=x.2-3*x+1)title(y2=x2-3x+1)xlabel(x), ylabel(y2) subplot(2,2,

11、3)plot(x,y1,r-),grid onlegend(ity1=2*x+5)title(調(diào)整后的y1=2x+5)axis(-3 3 -10 10)xlabel(x), ylabel(y1) subplot(2,2,4)plot(x,y2,k-),grid onlegend(ity2=x.2-3*x+1)title(調(diào)整后的y2=x2-3x+1)axis(-3 3 -10 10) %調(diào)整坐標(biāo)軸xlabel(x), ylabel(y2) 運(yùn)行后的圖形：繪制餅圖。程序?yàn)椋簒=190 33 45 42 45;explode=0 1 0 0 0;figure(1)subplot(2,1,1)col

12、ormap hsvpie(x,explode)gtext(生活費(fèi))gtext(資料費(fèi))gtext(電話費(fèi))gtext(衣服)gtext(其它)title(二維餅圖) subplot(2,1,2)colormap hsvpie3(x,explode)title(三維餅圖)圖形為：5、畫(huà)出函數(shù)z=（x-2）2+（y-1.2）2+sin(xy)的三維曲線和網(wǎng)格曲線。程序?yàn)椋簒,y=meshgrid(0:0.5:10); %為三維繪圖產(chǎn)生x，y數(shù)據(jù)矩陣z=(x-2).2+(y-1.2).2;figure(1)subplot(2,1,1) mesh(x,y,z),grid on %繪制網(wǎng)格曲線title

13、(網(wǎng)格曲線)subplot(2,1,2)plot3(x,y,z),grid ontitle(三維曲線)運(yùn)行后的圖形：6、畫(huà)出下列函數(shù)的曲面及等高線圖z=x2+y2+sin(xy)。程序?yàn)椋簒,y=meshgrid(0:pi/10:2*pi);z=x.2+y.2+sin(x*y);figure(1)subplot(2,1,1)surfc(x,y,z), grid ontitle(曲面和等高線) subplot(2,1,2)c,h=contour(x,y,z);set(h,showtext,on,textstep,get(h,levelstep)*2);title(等高線)運(yùn)行后的圖形：1、將圖形

14、窗口分成兩個(gè)，分別繪制正割和余割曲線，并加上標(biāo)注。程序?yàn)椋簒1=0:pi10:2*pi;figure(1)subplot(2,1,1)plot(x,sec(x),k-),grid onlegend(ity=sec(x)title(y=sec(x)xlabel(x), ylabel(y) subplot(2,1,2)plot(x,csc(x),k-),grid onlegend(ity=csc(x)title(y=csc(x)xlabel(x), ylabel(y)運(yùn)行后圖形為：2、畫(huà)出對(duì)數(shù)和指數(shù)曲線并加上標(biāo)注。x=0.01:0.1:10;y1=log10(x);y2=exp(x);figure

15、(1)subplot(2,1,1)plot(x,y1,k-),grid onlegend(ity1=log-10(x)title(y1=log-10(x)xlabel(x), ylabel(y1) subplot(2,1,2)plot(x,y2,k-),grid onlegend(ity2=exp(x)title(y2=exp(x)xlabel(x), ylabel(y2) 運(yùn)行后圖形為：3、設(shè)有函數(shù)y=exp(x+5)+x.3，在半對(duì)數(shù)坐標(biāo)系中繪制曲線。程序?yàn)椋簒=1:0.01:10;y=exp(x+5)+x.3;figure(1)subplot(2,1,1)plot(x,y,r-),gri

16、d onlegend(ity=exp(x+5)+x.3)title(平面坐標(biāo))xlabel(x), ylabel(y) subplot(2,1,2)semilogx(x,y,k-),grid on %半對(duì)數(shù)坐標(biāo)軸legend(ity=exp(x+5)+x.3)title(半對(duì)數(shù)坐標(biāo))xlabel(x), ylabel(y) 運(yùn)行后圖形為：4、畫(huà)出各種大小和形狀的球和柱體。繪制柱體的程序?yàn)椋簍=0:pi/10:2*pi;figure(1)subplot(2,1,1)x,y,z=cylinder(2+cos(t);surf(x,y,z),axis squaretitle(復(fù)雜柱面體) subplo

17、t(2,1,2)cylinder, axis squaretitle(簡(jiǎn)單柱體)繪制球的程序?yàn)椋篺igure(1)subplot(2,1,1)sphereaxis equaltitle(半徑為1的球) subplot(2,1,2)x,y,z=sphere;x=2*x;y=2*y;z=2*z;surf(x,y,z),axis squaretitle(半徑為2的球)運(yùn)行后的圖形： 5、繪制三維條形圖：程序?yàn)椋篩=cool(7);figure(1)subplot(2,2,1),bar3(Y,detached),title(Detached)subplot(2,2,2),bar3(Y,0.25,det

18、ached),title(Width=0.25)subplot(2,2,3),bar3(Y,grouped),title(Grouped)subplot(2,2,4),bar3(Y,stacked),title(Stacked)運(yùn)行后的圖形為：6、繪制二維條形圖程序?yàn)椋篩=round(rand(5,3)*10);figure(1)subplot(2,2,1),bar(Y,group),title(Group)subplot(2,2,2),bar(Y,stack),title(Stack)subplot(2,2,3),barh(Y,stack),title(Stack)subplot(2,2,4

19、),bar(Y,1.5),title(Width=1.5)運(yùn)行后的圖形：編寫(xiě)M函數(shù)實(shí)現(xiàn)：求一個(gè)數(shù)是否為素?cái)?shù)，在編寫(xiě)一主程序，要求通過(guò)鍵盤(pán)輸入一個(gè)整數(shù)，然后完成判斷其是否為素?cái)?shù)。解：function prime(x)n=fix(sqrt(x);for i=2:n if rem(x,i)=0 a=fasle return else a=true end end運(yùn)行結(jié)果： x=56; prime(x)a =fasle編寫(xiě)程序完成從表示字符的響亮中刪去空格，并求出字符個(gè)數(shù)。解：function nstr,n=del(str)nstr=;k=find(str= );nstr=str(k);n=lengt

20、h(nstr);end運(yùn)行后為： str=dr hy fgh gtesd hgfds; nstr,n=del(str)nstr =drhyfghgtesdhgfdsn = 173、編寫(xiě)M函數(shù)統(tǒng)計(jì)十進(jìn)制數(shù)值中0的個(gè)數(shù)，然后編寫(xiě)腳本文件，實(shí)現(xiàn)統(tǒng)計(jì)所有自然數(shù)12006中0的個(gè)數(shù)。解：M函數(shù)為：function y=geshu(x)s=num2str(x);n=length(s);m=0; if s(1)=0 disp(x is error); return end for i=2:n if s(i)=0 m=m+1; end end y=m;腳本文件為 jiu4：sum=0;for x=1:2006

21、 y=geshu(x); sum=sum+y;enddisp(sum)運(yùn)行結(jié)果為： jiu4 5044、利用menu函數(shù)輸入選擇參數(shù)ch。當(dāng)ch=1時(shí)，產(chǎn)生-10,10之間均勻分布的隨機(jī)數(shù)；當(dāng)ch=2時(shí)，產(chǎn)生-5,5之間均勻分布的隨機(jī)數(shù)；當(dāng)ch=3時(shí)，產(chǎn)生-1,1之間均勻分布的隨機(jī)數(shù)；當(dāng)ch=4時(shí)，產(chǎn)生均值為0，方差為1的正態(tài)分布隨機(jī)數(shù)。要求使用switch函數(shù)。解：s=menu(ch,1,2,3,4);n=;switch s case 1,n=20*rand(3)-10 case 2,n=10*rand(3)-5 case 3,n=2*rand(3)-1 case 4,n=randn(3)

22、otherwise disp(error)end運(yùn)行后：按下2后：n = 4.2274 0.4366 3.3897 3.0037 4.8478 -0.6674 -2.1405 2.1568 -0.2938求陣列x的平均值和標(biāo)準(zhǔn)差解：function mean1,stdev=stat2(x)m,n=size(x);if m=1m=n;ends1=sum(x);s2=sum(x.2);mean1=s1/m;stdev=sqrt(s2/m-mean1.2);運(yùn)行后： x=rand(4,4)+2; mean1,stdev=stat2(x)mean1 = 2.5207 2.3922 2.6498 2.2

23、539stdev = 0.1713 0.1892 0.1725 0.2027測(cè)試程序執(zhí)行時(shí)間% tech1.mtici=0;for t=0:.01:100 i=i+1; y(i)=sin(t);endtoc % tech2.mtict=0:.01:100;y=sin(t);Toc運(yùn)行后：Elapsed time is 0.015217 seconds.Elapsed time is 0.000508 seconds.1、產(chǎn)生menu選擇輸出顏色解：s=menu(color selection,red,green,blue,yellow,black)switch s case 1,scolor=

24、red; case 2,scolor=green; case 3,scolor=blue; case 4,scolor=yellow; case 5,scolor=black; otherwise disp(error)end Scolor運(yùn)行后：按下red后：s = 1scolor =red2、企業(yè)發(fā)放的獎(jiǎng)金按個(gè)人完成的利潤(rùn)（I）提成。分段提成比例為=即如王某完成25萬(wàn)元利潤(rùn)時(shí)，個(gè)人可得y=10 x 10% + 10 x 5% + 5 x 2% (萬(wàn)元)據(jù)此編寫(xiě)程序，求企業(yè)職工的獎(jiǎng)金。解：function bonus=bon(I)n=fix(I/100000)if(n4) n=4;endbo

25、n1=100000*0.1;bon2=0.05*(200000-100000);bon3=0.02*(400000-200000);switch n case 0,bonus=I*100000; case 1 bonus=bon1+0.05*(I-100000); case 2,3 bonus=bon1+bon2+0.02*(I-200000); case 4,bonus=bon1+bon2+bon3+0.01*(I-400000);end運(yùn)行后： I=1700000; bonus=bon(I)n = 17bonus = 32000有一分?jǐn)?shù)序列2/1，3/2，5/3/，8/5求前15項(xiàng)和。解：

26、s=1;t=2;sum=0;x=t/s;sum=sum+x;for i=1:15z=t;t=s+t;s=z;x=t/s;sum=sum+x;endsum運(yùn)行后： qiuhesum = 26.1881約瑟夫環(huán)解：n=input(please input n:);m=input(please input m:);b=1:n;i=1;c=0;s=0;while s yueseplease input n:12please input m:3a = Columns 1 through 8 3 6 9 12 4 8 1 7 Columns 9 through 16 2 11 5 10 3 16 5 20

27、 Columns 17 through 23 11 9 2 10 19 15 1編寫(xiě)程序計(jì)算x在（-3,3）上，并畫(huà)出曲線。解：function y=func2(x)n=length(x);for i=1:n; if (x(i)=-3)&(x(i)=-1)&(x(i)=1)&(x(i) v,d=eig(a,b)v = -0.4330 -0.2543 -0.1744 -0.5657 0.9660 -0.6091 -0.7018 0.0472 0.7736d = 13.5482 0 0 0 4.8303 0 0 0 3.6216 a=9 1 2;5 6 3;8 2 7; u,s,v=svd(a)u

28、 = -0.5601 0.5320 -0.6350 -0.4762 -0.8340 -0.2788 -0.6779 0.1462 0.7204s = 15.5234 0 0 0 4.5648 0 0 0 3.3446v = -0.8275 0.3917 -0.4023 -0.3075 -0.9156 -0.2592 -0.4699 -0.0907 0.8781 l,u=lu(a)l = 1.0000 0 0 0.5556 1.0000 0 0.8889 0.2041 1.0000u = 9.0000 1.0000 2.0000 0 5.4444 1.8889 0 0 4.8367 q,r=qr

29、(a)q = -0.6903 0.3969 -0.6050 -0.3835 -0.9097 -0.1592 -0.6136 0.1221 0.7801r = -13.0384 -4.2183 -6.8260 0 -4.8172 -1.0807 0 0 3.7733求解微分方程。解：function dy=funf(t,y)dy=5*y(1)-5*y(2)-6*y(3);3*y(1)-2*y(2)+5*y(3);2*y(1)-y(2)-4*y(3)；腳本文件：x0=1,-4,5;tspan=30,100;t,x=ode45(funf,tspan,x0);plot3(x(:,1),x(:,2),x

30、(:,3),grid ontitle(微分方程曲線)運(yùn)行后：微分方程組x=10(-x+y);y=28x-y-xz;z=xy-8z/3,x0=12,2,9,求微分方程在0,30上的解，并畫(huà)出系統(tǒng)軌跡。解：腳本文件： 二維圖：三維圖：2、分別用多項(xiàng)式和指數(shù)函數(shù)進(jìn)行擬合。y1(t)=c1+c2t+c3t2,y2(t)=d1+d2exp(t)解：t=0 0.2 0.4 0.6 0.8 1.0 2.0 5.0;y=1.0 1.51 1.88 2.13 2.29 2.40 2.60 -4.00;B1=ones(size(t) t t.*t;B2=ones(size(t) exp(t);A=B1y;C=B2

31、y;T=0:.1:6;Y1=ones(size(T) T T.*T*A;Y2=ones(size(T) exp(T)*C;plot(T,Y1,-,T,Y2,-,t,y,o)legend(itY1,itY2)運(yùn)行后：將(x-6)(x-3)(x-8)展開(kāi)為系數(shù)多項(xiàng)式的形式。解： a=6 3 8; pa=poly(a); ppa=poly2sym(pa) ppa = x3-17*x2+90*x-144求解多項(xiàng)式x3-7x2+2x+40的根。解： r=1 -7 2 40; p=roots(r); -0.2151 0.4459 0.7949 0.2707求解在x=8時(shí)多項(xiàng)式(x-1)(x-2) (x-3

32、)(x-4)的值。解： p=poly(1 2 3 4); polyvalm(p,8) ans = 840計(jì)算多項(xiàng)式乘法(x2+2x+2)(x2+5x+4)。解： c=conv(1 2 2,1 5 4) c = 1 7 16 18 8計(jì)算多項(xiàng)式除法(3x3+13x2+6x+8)/(x+4)。解： d=deconv(3 13 6 8,1 4) d = 3 1 2求雙重積分。9、微分方程組當(dāng)t=0，=1；=-0.5，求微分方程組t【0,25】上的解，并畫(huà)出x1-x2的系統(tǒng)軌跡。解：function dy=fund(t,y)dy=0.5-y(1);y(1)-4*y(2);腳本文件：x0=1,-0.5;

33、tspan=0,20;T,Y=ode23(fund,tspan,x0);figure(1)plot(T,Y(:,1),r-,T,Y(:,2)legend(itx1,itx2)運(yùn)行后：利用下標(biāo)建立多維陣列。產(chǎn)生一個(gè)332的多維矩陣AA=5 7 2; 0 1 2; 3 4 2; %產(chǎn)生一個(gè)3*3矩陣A(:, :, 2)=2 7 3; 4 2 8; 2 0 3A(:,:,1) = 5 7 2 0 1 2 3 4 2A(:,:,2) = 2 7 3 4 2 8 2 0 3利用MATLAB函數(shù)產(chǎn)生多維陣列。 利用MATLAB的函數(shù)(如rand、randn、ones、zeros等)都可直接產(chǎn)生多維陣列，在

34、函數(shù)調(diào)用時(shí)可指定每一維的尺寸。例如，為產(chǎn)生10032維的正態(tài)分布隨機(jī)數(shù)R，可輸入R=randn(100, 3, 2);A=5*ones(3, 4, 2); %產(chǎn)生元素相同的多維陣列B=repmat(5, 3 4 2); %產(chǎn)生元素相同的多維陣列利用cat函數(shù)建立多維陣列A=2 8; 0 5; B=1 8; 2 4; C=cat(3,A,B); D=cat(4,A,B); size(C)ans = 2 2 2 size(D)ans = 2 2 1 2這說(shuō)明得到的C為222維，而D為2212維。利用賦值語(yǔ)句建立結(jié)構(gòu)陣列 =John Doe;patient.billing=1

35、27.00;patient.test=79 75 73;180 178 177.5;220 210 205; 這時(shí)就建立起了一個(gè)具有三個(gè)域的結(jié)構(gòu)patient，當(dāng)輸入 patientpatient = name: John Doe billing: 127 test: 3x3 double size(patient)ans =1 1 5、利用struct函數(shù)建立結(jié)構(gòu)陣列利用單元陣列還可以一次輸入多個(gè)結(jié)構(gòu)元素，例如輸入 n=John Doe Ann Lane Alan Johnson這是一個(gè)單元陣列n = John Doe Ann Lane Alan Johnson b=127 28.5 95.

36、8; t1=79 75 73; 180 178 177.5; 220 210 205; t2=68 70 68; 118 118 119; 172 170 169; t3=37 38 36; 119 121 120; 165 166 159; patient2=struct(name,n,billing,b,test,t1 t2 t3); patient2patient2 = 13 struct array with fields: name billing test 結(jié)構(gòu)嵌套A=struct(data,3 4 7;8 0 1,nest, struct(testnum,Test 1,xdata

37、,4 2 8,ydata,7 1 6);A(2).data=9 3 2;7 6 5;A(2).nest.testnum=Test 2;A(2).nest.xdata=3 4 2;A(2).nest.ydata=5 0 9 ; A(1).dataans = 3 4 7 8 0 1 A(2).nestans = testnum: Test 2 xdata: 3 4 2 ydata: 5 0 9 A(2).nest.testnumans =Test 2 A(2).nest.ydata(1,3)ans = 97利用賦值語(yǔ)句建立單元陣列利用賦值語(yǔ)句可對(duì)一個(gè)單元陣列進(jìn)行賦值，與一般陣列賦值時(shí)使用方括號(hào)不同

38、，它應(yīng)使用花括號(hào)，花括號(hào)可出現(xiàn)在賦值號(hào)右邊，也可以出現(xiàn)在左邊。例如，建立22單元陣列A，可輸入A(1,1)=1 4 3; 0 5 8; 7 2 9;A(1,2)=Anne Smith;A(2,1)=3+7i;A(2,2)= -pi:pi/10:pi; 為得到同樣的單元陣列A，也可以輸入A1,1=1 4 3; 0 5 8; 7 2 9;A1,2=Anne Smith;A2,1=3+7i;A2,2=-pi:pi/10:pi; 8利用cell函數(shù)定義單元陣列cell函數(shù)可預(yù)分配單元陣列，但其內(nèi)容為空陣列。例如，建立23的單元陣列，可輸入B=cell(2,3);然后利用賦值語(yǔ)句可以給各個(gè)單元賦值，如輸

39、入 B(1,3)=1:4; B(2,2)=rand(3,3); BB = 1x4 double 3x3 double 7、 綜合設(shè)計(jì)示例設(shè)某高等學(xué)校某系兩個(gè)年級(jí)共有八個(gè)班(29912994，29812984)，某年度由五位教師共開(kāi)設(shè)了九門(mén)課程，設(shè)置情況如表1所示。 表1 課程設(shè)置情況學(xué)分29912992299329942981298229832984高等數(shù)學(xué)5普通物理4英語(yǔ)(1)3英語(yǔ)(2)3計(jì)算機(jī)應(yīng)用基礎(chǔ)3電路分析基礎(chǔ)4低頻電子線路4高頻電子線路2數(shù)字信號(hào)處理3 表2 教師教學(xué)情況29912992299329942981298229832984高等數(shù)學(xué)3112普通物理2332英語(yǔ)(1)545

40、4英語(yǔ)(2)4545計(jì)算機(jī)應(yīng)用基礎(chǔ)1234電路分析基礎(chǔ)3215低頻電子線路1143高頻電子線路5134數(shù)字信號(hào)處理2243表3 考試平均成績(jī)統(tǒng)計(jì)表 29912992299329942981298229832984高等數(shù)學(xué)76.2 68.9 70.1 71.3普通物理60.7 59.2 55.4 61.2英語(yǔ)(1)60.2 65.4 58.1 59.2英語(yǔ)(2)62.3 68.3 71.2 63.9計(jì)算機(jī)應(yīng)用基礎(chǔ)75.4 80.1 73.2 74.2電路分析基礎(chǔ)66.1 65.4 68.2 70.0低頻電子線路56.9 59.1 52.3 57.2高頻電子線路66.1 60.9 61.2 64.

41、9數(shù)字信號(hào)處理77.1 78.9 69.9 72.1要求設(shè)計(jì)MATLAB程序保存這些結(jié)構(gòu)，計(jì)算各個(gè)班的學(xué)分加權(quán)平均成績(jī)，計(jì)算公式為 各教師的教學(xué)工作量，計(jì)算公式為 教學(xué)工作量=學(xué)時(shí)(1+w) 其中 x =人數(shù)-40 MATLAB程序如下： clear課程基本信息C=高等數(shù)學(xué); Course.score=5; Course.time=90;Course(2).name=普通物理; Course(2).score=4;Course(2).time=80; Course(3).name=英語(yǔ)(1); Course(3).score=3;Course(3).time=64;Cour

42、se(4).name=英語(yǔ)(2); Course(4).score=3;Course(4).time=64;Course(5).name=計(jì)算機(jī)應(yīng)用基礎(chǔ); Course(5).score=3;Course(5).time=64;Course(6).name=電路分析基礎(chǔ); Course(6).score=4;Course(6).time=80;Course(7).name=低頻電子線路; Course(7).score=4;Course(7).time=80;Course(8).name=高頻電子線路; Course(8).score=2;Course(8).time=46;Course(9)

43、.name=數(shù)字信號(hào)處理; Course(9).score=3; Course(9).time=64;教師基本信息T=王立偉; Teacher.posit=講師; Teacher.age=30;Teacher(2).name=趙杰; Teacher(2).posit=副教授;Teacher(2).age=36;Teacher(3).name=李志范; Teacher(3).posit=教授; Teacher(3).age=45;Teacher(4).name=陳根生; Teacher(4).posit=教授; Teacher(4).age=50;Teacher(5).na

44、me=孫興國(guó); Teacher(5).posit=講師; Teacher(5).age=26;學(xué)生班基本信息 Student.class=2991; Student.tol=130;Student(2).class=2992; Student(2).tol=56;Student(3).class=2993; Student(3).tol=45;Student(4).class=2994; Student(4).tol=118;Student(5).class=2981; Student(5).tol=125;Student(6).class=2982; Student(6).tol=53;St

45、udent(7).class=2983; Student(7).tol=40;Student(8).class=2984; Student(8).tol=106;教學(xué)基本信息(用單元陣顯示)TS1=高等數(shù)學(xué),2991 2992 2993 2994,3 1 1 2,76.2 68.9 70.1 71.3;TS2=普通物理,2991 2992 2993 2994,2 3 3 2,60.7 59.2 55.4 61.2; TS3=英語(yǔ)(1),2991 2992 2993 2994,5 4 5 4,60.2 65.4 58.1 59.2;TS4=英語(yǔ)(2),2981 2982 2983 2984,4

46、5 4 5,62.3 68.3 71.2 63.9;TS5=計(jì)算機(jī)應(yīng)用基礎(chǔ),2991 2992 2993 2994,1 2 3 4,75.4 80.1 73.2 74.2;TS6=電路分析基礎(chǔ),2981 2982 2983 2984,3 2 1 5,66.1 65.4 68.2 70.0;TS7=低頻電子線路,2991 2992 2993 2994,1 1 4 3,56.9 59.1 52.3 57.2;TS8=高頻電子線路,2981 2982 2983 2984,5 1 3 4,66.1 60.9 61.2 64.9;TS9=數(shù)字信號(hào)處理,2981 2982 2983 2984,2 2 4

47、3,77.1 78.9 69.9 72.1;TS=cat(3,TS1,TS2,TS3,TS4,TS5,TS6,TS7,TS8,TS9);計(jì)算8個(gè)班的學(xué)分加權(quán)平均disp1=班的學(xué)分加權(quán)平均為; for i=1:8 c=Student(i).class; s=0;v=0; for j=1:9 k=find(TS1,2,j=c); if (isempty(k) v1=TS1,4,j(k); s1=Course(j).score; s=s+s1; v=v+v1*s1; end end av=v/s; disp(num2str(c) disp1 num2str(av)end計(jì)算教師工作量disp2=的

48、教學(xué)工作量為;disp3=小時(shí);for i=1:5 c=Teacher(i).name; c1=Teacher(i).posit; g=0; for j=1:9 k=find(TS1,3,j=i); if (isempty(k) for l=1:length(k) cn=TS1,2,j(l) -2990; if cnx(j) temp=x(i); x(i)=x(j); x(j)=temp; end end end y=x;運(yùn)行結(jié)果： x=12 34 654 2 5 76 23; y=bubblesort(x)y = 2 5 12 23 34 76 654以上為按照升序排列的，若要降序，則 if x(i) x=12 21 2 4 5 19 45 30; y=bubblesort(x)y =45 30 21 19 12 5 4 2傅里葉變換應(yīng)用付立葉