matlab習(xí)題——2.doc

第二章1。%3X3矩陣 a1=eye(3) a2=ones(3) a3=zeros(3) a4=1-2*rand(3) a5=2*randn(3)+1 %15X8矩陣 b1=eye(15,8) b2=ones(15,8) b3=zeros(15,8) b4=1-2*rand(15,8) b5=2*randn(15,8)+1 2。%生成a矩陣 x=diag(2,-7,8); xx=diag(3,5,-1); a=rot90(x+xx) %生成b矩陣 x=diag(2 5 8); xx=rot90(diag(7 0 4); b=x+xx %變成行向量 ah=reshape(a,1,9) bh=reshape(b,1,9)3。a=5-10*rand(50,2); b=round(10*a)/10; c=num2str(b,%+10.1f)4。alpha=-pi:pi/180:pi; x=sin(alpha); y=cos(alpha); %畫曲線 plot(alpha,x,b-) hold on lot(alpha,y,r-）5。a=rand(10); b=a=0.6; total=sum(sum(b)6。a=randn(10); b=and(a-.5); total=sum(sum(b)7。%輸入a和b a=input(please input a:); b=input(please input b:); %分支結(jié)構(gòu) if and(a1,b=.5) disp(語句1) elseif and(a.5) disp(語句2) elseif and(a1,b1); b=a(i)8b。a=rand(10)+0.2; x=a1; =a.*x; i,j,b=find(c); b9。%產(chǎn)生測量陣 a=30*rand(100,3)+5; a(97)=103; a(156)=280; %查找刪除 i,j=find(a100); a(i,:)=10。%生成10X10矩陣，包含全零行，全零列 a=round(10*rand(10); a(4,:)=0; a(:,7)=0; %查找全零行，全零列 hang=any(a); lie=any(a); hang_hao=find(hang); lie_hao=find(lie); %刪除！ a(hang_hao,:)= a(:,lie_hao)=第三章1. %對數(shù)曲線 x=0.1:0.01:100; y=log2(x); semilogx(x,y) title(對數(shù)曲線) xlabel(log(x) ylabel(y) gtext(leftarrow對數(shù)曲線) %指數(shù)曲線 x=0:0.1:10; y=exp(x); figure(2) plot(x,y) title(指數(shù)曲線) xlabel(x) ylabel(y) gtext(指數(shù)曲線rightarrow)2. x=0:pi/50:4*pi; k=1:25:201; x(k)=; subplot(1,2,1) plot(x,sec(x),r-) grid on title(正割曲線) subplot(1,2,2) plot(x,csc(x),g-) grid on title(余割曲線) 3. %對數(shù)坐標(biāo) x=0.01:0.01:100; y=log10(x); figure(1) subplot(2,1,1) plot(x,y,k-) grid on title(迪卡爾坐標(biāo)系中y=log_10(x)曲線),ylabel(y) subplot(2,1,2) semilogx(x,y,k-),grid on title(半對數(shù)坐標(biāo)系中y=log_10(x)曲線) xlabel(x),ylabel(y) %極坐標(biāo) a=randn(2,2);b=eig(a) c1=abs(b),c2=angle(b) figure(1) subplot(2,1,1) plot(b,rx),grid on title(迪卡爾坐標(biāo)系中特征值的表示) subplot(2,1,2) polar(c2,c1,rx) gtext(極坐標(biāo)中的特征值的表示)4. %多峰函數(shù)多條曲線 x,y=meshgrid(-3:0.15:3); z=peaks(x,y); x=x(1,:); plot(x,z),grid on %三角函數(shù)多條曲線 t=-pi:pi/20:pi; y1=sin(t);y2=2*cos(t); plotyy(t,y1,t,y2),grid on title(雙y軸正余弦曲線) text(0,0,leftarrow sin(t) text(pi/2,0,leftarrow 2cos(t)5. x=-3:0.1:3; y1=2*x+5; y2=x.2-3*x+1; subplot(1,2,1) plot(x,y1,r-) grid on title(y=2x+1) axis(-3,3,-5,20) subplot(1,2,2) plot(x,y2,g-) grid on title(y=x2-3x+1)6. camp=colormap; L=length(camp); x=1:L; y=x*ones(size(x); figure(1) subplot(6,1,3) bar(x(1:2),y(1:2,:) title(hsv顏色板) colormap(hsv)7. x=190 33 45 42 45; ex=0 1 0 0 0; pie(x,ex) text(0.3,-.8,資料費) text(0.5,.2,服裝費) text(0.6,-.1,電話費) text(0.1,.8,其它費) text(-1.1,-.8,生活費) title(王某月生活費調(diào)查) colormap(autumn)8. x,y=meshgrid(-10:.1:14,-10:.1:12.4); z=(x-2).2+(y-1.2).2; mesh(x,y,z),grid on9. x,y=meshgrid(-4:0.1:4); z=x.2+y.2+sin(x.*y); surfc(x,y,z) colormap(spring).%腳本文件 x=input(請輸入一個整數(shù) x:); if myprime(x) disp(您輸入的整數(shù)x是一個素數(shù)。) else disp(您輸入的數(shù)x不是一個素數(shù)。) end %myprime函數(shù) function y=myprime(x) y=1; for i=2:fix(sqrt(x) if mod(x,i)=0 y=0; end第四章2.編寫程序完成從表示字符的向量中刪去空格，并求出字符的個數(shù)a=wer rty inx=a= b=a(inx)num=length(b);num=sum(inx);3.編寫M函數(shù)統(tǒng)計十進(jìn)制數(shù)值中“0”的個數(shù)，然后編寫腳本文件，統(tǒng)計1-2006中“0”的個數(shù)。s=0; for i=1:1999 s=s+my_total_O(i); end s %my_total_0函數(shù) function y=my_total_0(x) y=0; n=fix(log10(x)+1; for i=1:n y=y+mod(x,10); x=fix(x/10); end5. 編寫程序計算（步長0.01）并且畫出曲線。 y=0;step=0.01; x1=-3:step:-1-step; x2=-1:step:1-step; x3=1:step:3; y1=(-x1.2-4*x1-3)/2; y2=-x2.2+1; y3=(-x3.2+4*x3-3)/2; x=x1,x2,x3; yy=y1,y2,y3; plot(x,yy,linestyle)6. 利用menu函數(shù)輸入選擇參數(shù)ch，當(dāng)ch=1時，產(chǎn)生【-10,10】之間均勻分布的隨機(jī)數(shù)；當(dāng)ch=2時，產(chǎn)生【-5,5】之間均勻分布的隨機(jī)數(shù)；當(dāng)ch=3時，產(chǎn)生【-1,1】之間均勻分布的隨機(jī)數(shù)；當(dāng)ch=4時，產(chǎn)生均值為0，方差為1的正太分布隨機(jī)數(shù)。要求使用switch函數(shù)。N=3; %n為產(chǎn)生矩陣的大小 ch=menu(選擇矩陣形式,-10,10均勻分布,-5,5均勻分布,. -1,1均勻分布,正態(tài)分布); switch ch case 1,a=10-20*rand(n); case 2,a=5-10*rand(n); case 3,a=1-2*rand(n); case 4,a=randn(n); end 9. 有一分布序列，求其前15項和。n=15 x(1)=1; x(2)=2; %求出費波納契序列 for i=3:n+1 x(i)=x(i-2)+x(i-1); end %求出各項 for i=1:n a(i)=x(i+1)/x(i); end %求出前n項和 y=sum(a);10.有n個人圍城一圈，按序列編號。從第一個人起數(shù)到m時該人退出，并且從1開始報數(shù)，求出圈人順序。function y=chap4_10(m,n) x=1:n; for i=1:n N=n+1-i; t=mod(m,N); if t=0 t=N;end y(i)=x(t); x=x(t+1:N),x(1:t-1); end第五章4. % 設(shè)A=11.9 12.81 15.66;15.2 4.18 13.61;10.59 7.59 9.22 , % B=16.00 4.41 -10.37 -21.61;0.88 -20.04 12.86 8.56;-1.43 10.71 18.81 - 5.99;-12.48 24.35 -23.9 10.34 % 分別求出這兩個矩陣的LU和QR的分解。 %程序如下： % A的LU分解 A=11.9 12.81 15.66;15.2 4.18 13.61;10.59 7.59 9.22 L,U=lu(A) % B的LU分解 B=16.00 4.41 -10.37 -21.61;0.88 -20.04 12.86 8.56;-1.43 10.71 18.81 -5.99; -12.48 24.35 -23.9 10.34 L,U=lu(B) % A的QR分解 A=11.9 12.81 15.66;15.2 4.18 13.61;10.59 7.59 9.22 Q,R=qr(A) % B的QR分解 B=16.00 4.41 -10.37 -21.61;0.88 -20.04 12.86 8.56;-1.43 10.71 18.81 -5.99; -12.48 24.35 -23.9 10.34 Q,R=qr(B)8. % 將下列多項式進(jìn)行因式分解，也即計算出多項式的根。 % （1）：p1(x)=x4-2*x3-3*x2+4*x+2 % （2）：p2(x)=x4-7*x3+5*x2+31*x-30 % （3）：p3(x)=x3-x2-25*x+25 % （4）：p4(x)=-2*x5+3*x4+x3+5*x2+8*x %程序如下： p1=1 -2 -3 4 2; r1=roots(p1) p2=1 -7 5 31 -30; r2=roots(p2) p3=1 -1 -25 25; r3=roots(p3) p4=-2 3 1 5 8 0; r4=roots(p4)12.% 對函數(shù)y=10*exp(-|x|)取x-5,-4,-3,，3，4，5點的值作為粗值， % 分別采用最鄰近內(nèi)插，線性內(nèi)插，三次樣條內(nèi)插和三次曲線內(nèi)插方法， % 對-5,5內(nèi)的點進(jìn)行內(nèi)插，比較其結(jié)果。 % 程序如下： %產(chǎn)生低分辨率峰值函數(shù) x,y=meshgrid(-5:1:5); z=peaks(x,y); %定義內(nèi)差點 xi,yi=meshgrid(-5:0.25:5); %利用四種方法進(jìn)行內(nèi)插 %最鄰近內(nèi)插 zi1=interp2(x,y,z,xi,yi,nearest); %線性內(nèi)插 zi2=interp2(x,y,z,xi,yi,linear); %三次樣條內(nèi)插 zi3=interp2(x,y,z,xi,yi,spline); %三次曲線內(nèi)插 zi4=interp2(x,y,z,xi,yi,cubic); %繪出曲面圖進(jìn)行比較 figure(1) subplot(2,3,1) surf(x,y,z),title(原始數(shù)據(jù)曲線) subplot(2,3,2) surf(xi,yi,zi1),title(最鄰近內(nèi)插) subplot(2,3,3) surf(xi,yi,zi2),title(線性內(nèi)插) subplot(2,3,4) surf(xi,yi,zi3),title(三次樣條內(nèi)插) subplot(2,3,5) surf(xi,yi,zi4),title(三次曲線內(nèi)插) %繪制輪廓圖進(jìn)行比較 figure(2) subplot(2,3,1) contour(x,y,z),title(原始數(shù)據(jù)曲線) subplot(2,3,2) contour(xi,yi,zi1),title(最鄰近內(nèi)插) subplot(2,3,3) contour(xi,yi,zi2),title(線性內(nèi)插) subplot(2,3,4) contour(xi,yi,zi3),title(三次樣條內(nèi)插) subplot(2,3,5) contour(xi,yi,zi4),title(三次曲線內(nèi)插)16.% 產(chǎn)生三個信號： % x1sin(kt)+randn(size(t) % x2cos(kt)+randn(size(t) % x3sin(kt)+randn(size(t) % 試計算x1與x2，x1與x3之間的相關(guān)系數(shù)，從中可得出什么結(jié)論? % 如果信號不含正余弦信號分量，結(jié)論有如何？ %程序如下： x0=2*pi*20; t=0:1/119:1; x1=sin(x0*t)+randn(size(t);