数学建模matlab作业.docVIP

数统 应数 20121323003 王妍 数统 应数 20121323022 胡可旺 1 生成5阶矩阵，使其元素满足均值为1，方差为4的正态分布； 代码：y=1+sqrt(4)*randn(5) 结果: 2，生成一个20行5列矩阵A，其元素分别以概率0.7、0.1、0.2取自0、1、4这三个整数， 然后将“删除A的全零行”之后的结果赋给矩阵B； for i = 1:20 for j = 1:5 p = rand(); if p=0.7 A(i,j) = 0; elseif p0.7p=0.8 A(i,j) = 1; else A(i,j) = 4; end end end B = []; for i = 1:20 flag = true; for j = 1:5 if A(i,j)~=0 flag=false; break; end end if ~flag B = [B;A(i,:)]; end end B 3, （函数文件）对于输入的正整数，判断其是否为素数 function [ output_args ] = u3( n ) %U3 Summary of this function goes here % Detailed explanation goes here y=1; for i=2:1:sqrt(n) if rem(n,i)==0 y=0; break end end if(y==1) fprintf(是素数) else fprintf(不是素数) end 4，（函数文件）对于输入的正整数，求出其中0的个数； function [ i ] = u4( n ) %U4 Summary of this function goes here % Detailed explanation goes here i=0; while n/10~=0 if rem(n,10)==0 i=i+1; end n=fix(n/10); end 5，（函数文件）有n个人围成一圈，按序列编号。从第一个人开始报数，数到m时该人退出， 并且下一个从1开始重新报数，那么最后剩下的那个人的原始编号为多少？（nm，例如n=20，m=7 function [ x ] = u5( n,m ) %U5 Summary of this function goes here % Detailed explanation goes here x=1:n; y=zeros(1,n-m+1); t=0; p=size(x,2); s=fix(p/m); while s0 j=1; while j=s y(t+j)=x(m*j); j=j+1; end x([m:m:m*s])=[]; if m*sp x=x([((m-1)*s+1):(p-s),1:((m-1)*s)]); end t=t+s; p=size(x,2); s=fix(p/m); end 6, 对曲线 | x^2+y^2+z^2=4 | x+y+z=1 作图，并修饰 ezplot(x^2+y^2+(1-x-y)^2-4,[-2,2,-2,2]) 曲面作图：z=x^2+y^2+sin(xy)，定义域自拟； x=0:100;y=0:100; [x,y]=meshgrid(x,y); z=x.^2+y.^2+sin(x.*y); plot3(x,y,z) 8，函数f(x)定义为如下的行列式： | x x^2 x^3 | | 1 2x 3x^2 | | 0 2 6x | 求f(x)在x=1处的数值导数，并画图比较数值导数与真实导函数之间的差别 syms x; A=[ x x^2 x^3 ; 1 2*x 3*x^2; 0 2 6*x ]; B=det(A); y=diff(B,x) subs(y,1) lims=[-5,5]; fplot(y,lims); hold on; fpl