matlab习题——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:); %分支结构 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。%产生测量阵 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. %对数曲线 x=0.1:0.01:100; y=log2(x); semilogx(x,y) title(对数曲线) xlabel(log(x) ylabel(y) gtext(leftarrow对数曲线) %指数曲线 x=0:0.1:10; y=exp(x); figure(2) plot(x,y) title(指数曲线) xlabel(x) ylabel(y) gtext(指数曲线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. %对数坐标 x=0.01:0.01:100; y=log10(x); figure(1) subplot(2,1,1) plot(x,y,k-) grid on title(迪卡尔坐标系中y=log_10(x)曲线),ylabel(y) subplot(2,1,2) semilogx(x,y,k-),grid on title(半对数坐标系中y=log_10(x)曲线) xlabel(x),ylabel(y) %极坐标 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(迪卡尔坐标系中特征值的表示) subplot(2,1,2) polar(c2,c1,rx) gtext(极坐标中的特征值的表示)4. %多峰函数多条曲线 x,y=meshgrid(-3:0.15:3); z=peaks(x,y); x=x(1,:); plot(x,z),grid on %三角函数多条曲线 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(王某月生活费调查) 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(请输入一个整数 x:); if myprime(x) disp(您输入的整数x是一个素数。) else disp(您输入的数x不是一个素数。) end %myprime函数 function y=myprime(x) y=1; for i=2:fix(sqrt(x) if mod(x,i)=0 y=0; end第四章2.编写程序完成从表示字符的向量中删去空格，并求出字符的个数a=wer rty inx=a= b=a(inx)num=length(b);num=sum(inx);3.编写M函数统计十进制数值中“0”的个数，然后编写脚本文件，统计1-2006中“0”的个数。s=0; for i=1:1999 s=s+my_total_O(i); end s %my_total_0函数 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函数输入选择参数ch，当ch=1时，产生【-10,10】之间均匀分布的随机数；当ch=2时，产生【-5,5】之间均匀分布的随机数；当ch=3时，产生【-1,1】之间均匀分布的随机数；当ch=4时，产生均值为0，方差为1的正太分布随机数。要求使用switch函数。N=3; %n为产生矩阵的大小 ch=menu(选择矩阵形式,-10,10均匀分布,-5,5均匀分布,. -1,1均匀分布,正态分布); 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个人围城一圈，按序列编号。从第一个人起数到m时该人退出，并且从1开始报数，求出圈人顺序。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. % 设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. % 将下列多项式进行因式分解，也即计算出多项式的根。 % （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.% 对函数y=10*exp(-|x|)取x-5,-4,-3,，3，4，5点的值作为粗值， % 分别采用最邻近内插，线性内插，三次样条内插和三次曲线内插方法， % 对-5,5内的点进行内插，比较其结果。 % 程序如下： %产生低分辨率峰值函数 x,y=meshgrid(-5:1:5); z=peaks(x,y); %定义内差点 xi,yi=meshgrid(-5:0.25:5); %利用四种方法进行内插 %最邻近内插 zi1=interp2(x,y,z,xi,yi,nearest); %线性内插 zi2=interp2(x,y,z,xi,yi,linear); %三次样条内插 zi3=interp2(x,y,z,xi,yi,spline); %三次曲线内插 zi4=interp2(x,y,z,xi,yi,cubic); %绘出曲面图进行比较 figure(1) subplot(2,3,1) surf(x,y,z),title(原始数据曲线) subplot(2,3,2) surf(xi,yi,zi1),title(最邻近内插) subplot(2,3,3) surf(xi,yi,zi2),title(线性内插) subplot(2,3,4) surf(xi,yi,zi3),title(三次样条内插) subplot(2,3,5) surf(xi,yi,zi4),title(三次曲线内插) %绘制轮廓图进行比较 figure(2) subplot(2,3,1) contour(x,y,z),title(原始数据曲线) subplot(2,3,2) contour(xi,yi,zi1),title(最邻近内插) subplot(2,3,3) contour(xi,yi,zi2),title(线性内插) subplot(2,3,4) contour(xi,yi,zi3),title(三次样条内插) subplot(2,3,5) contour(xi,yi,zi4),title(三次曲线内插)16.% 产生三个信号： % x1sin(kt)+randn(size(t) % x2cos(kt)+randn(size(t) % x3sin(kt)+randn(size(t) % 试计算x1与x2，x1与x3之间的相关系数，从中可得出什么结论? % 如果信号不含正余弦信号分量，结论有如何？ %程序如下： x0=2*pi*20; t=0:1/119:1; x1=sin(x0*t)+randn(size(t);