建模练习题9166.doc

1．利用基本矩阵产生3×3和15×8的单位矩阵、全1矩阵、全0矩阵、均匀分布随机矩阵（[-1，1]之间）、正态分布矩阵（均值为1，方差为4）。A=eye(3) b=eye(15,8) c=ones(3) d=zero(15,8) e=1+2*rand(3) 2．利用fix及rand函数生成[0，10]上的均匀分布的10×10的整数随机矩阵a，然后统计a中大于等于5的元素个数。A=fix(rand(10)*10) find(A=5) size(find(A=5)) 3．在给定的矩阵中删除含有整行内容全为0的行，删除整列内容全为0的列。 function a=question3(x) if size(x,1)==1 a=x(any(x),:); elseif size(x,2)==1 a=x(:,any(x)); else a=x(any(x),:); a=a(:,any(a)); end 要想生成（A,B）上的6阶随机均匀分布矩阵，就用A+B*RAND(6) 4．在同一图形窗口画出下列两条曲线图像： ?????????????? y1=2x+5； y2=x^2-3x+1， 并且用legend标注。x=0:pi/100:2*pi; y1=2*x+5; y2=x.^2-3*x+1; plot(x,y1,x,y2) legend(y1,y2) 5．画出下列函数的曲面及等高线： ????z=x^2+y^2+sin(xy)． [x,y]=meshgrid(-5:0.05:5); z=x.^2+y.^2+sin(x.*y); contour3(x,y,z,20,g)编写程序计算（x在[-3，3]，间隔0.01） x=input(请输入x=); if -3=xx-1 y=(-x^2-4*x-3)/2 elseif -1=xx1 y=-x^2+1 else 1=xx=3 y=(-x^2+4*x-3)/2 end 7．有一列分数序列： 求前15项的和。 for i=1:15 if i=2 f(1)=1;f(2)=2;f(3)=3; else f(i+1)=f(i)+f(i-1); end a(i)=f(i+1)/f(i); end s=sum(a) 8．用至少三种方法编写函数实现求任意整数n的阶乘 n=input(n=); s=1; for i=1:n s=s*i; end s 直接用FACTORIAL(N)或者 PROD(1:N) 或者GAMMA(N+1)(此处是求N的阶乘) function y=jiecheng(x) if( x==1)||(x==0) y=1; elseif x1 y=x*jiecheng(x-1); end 9．将任意偶数m写成两个素数p1、p2的和（试着写出所有的m=p1+p2的可能形式）。 m=input(m=); p=1:m;p(1)=0; for i=2:sqrt(m) for j=2*i:i:m p(j)=0; end end n=find(p~=0); p(n) k=length(p); for q=1:k for w=1:k if p(q)+p(w)==m if p(w)=p(q) disp([p(q),p(w)]) end end end end 10．是否任意3的倍数m可以写成两个素数p1、p2、p3的和（试着写出所有的m=p1+p2+p3m=input(m=); p=1:m;p(1)=0; for i=2:sqrt(m) for j=2*i:i:m p(j)=0; end end n=find(p~=0); p(n) p=p(p0) k=length(p); for q=1:k for w=1:k for t=1:(k/3+1) if p(q)+p(w)+p(t)==m if p(t)=p(w)p(w)=p(q) disp([p(q),p(w),p(t)]) end end end end end [4.842,4.362,3.754,3.368,3.