数学建模matlab作业

《数学建模matlab作业》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、数统应数20121323003王妍数统应数20121323022胡可旺1生成5阶矩阵，使其元素满足均值为1，方差为4的正态分布;代码：y=l+sqrt(4)*randn(5)结果：y=2.34303.06942.77683.87680.79551.41502.4538-1.29411.65040.51712.43450.3931-1.1377-0.50991.63844.26051.5877-0.61903.74061.62571.9778-0.5746-4.8886-2.4230-0.72982,生成一个20行5

2、列矩阵A,其元素分别以概率0.7、0.1、0.2取自0、1、4这三个整数,然后将“删除A的全零行”之后的结果赋给矩阵B;fori=1:20forj=1:5p=rand();ifp<=0.7A(ij)=0;elseifp>0.7&&p<=0.8A(i,j)=l；elseA(i,j)=4;endendendB=[]；fori=1:20flag=true;forj=1:5ifA(i，j)〜=0flag=false;break;endendif〜flagendendBB=40000114000144000404001000

3、04000410004440001000040104401400411000000104440000010404404fx»I3，(函数文件)对于输入的正整数，判断其是否为素数function[output_args]=u3(n)%U3Summaryofthisfunctiongoeshere%Detailedexplanationgoesherey=l;fori=2:1:sqrt(n)ifrem(n,i)==0y=o；breakendendif(y==l)fprintf('足素数1〉elsefprintf（’不是

4、素数’）endCommandWindow»u3⑸是羞教>>»u3(10)fx不是羞数>>4,(函数文件)对于输入的正整数，求出其中0的个数;function[i]=u4(n)%U4Summaryofthisfunctiongoeshere%Detailedexplanationgoesherei=0;whilen/10〜=◦ifrem(n,10)==0i=i+l;endn=fix(n/10);endCommandWindow»u4(10)arts=>>u4(66)arts=05,(函数文件)有n个人围成一圈，按序列

5、编号。从第一个人开始报数，数到m时该人退出，并且下一个从1开始重新报数，那么最后剩下的那个人的原始编号为多少？(n>m,例如n=20,m=7function[x]=u5(n,m)%U5Summaryofthisfunctiongoeshere%Detailedexplanationgoesherex=l:n;y=zeros(1,n-m+1);t=0;p=size(x,2);s=fix(p/m);whiles〉0j=l;whilej<=sy(t+j)=x(m*j);endx(])=[];ifm*s

6、-1)*s+l):(p-s),l:((m-1)*s)]);endt=t+s;p=size(x,2);s=fix(p/m);endCommandWindow»u5(5,2)arts=3»u5(20,9)ans=13151612356A»4,对曲线丨xA2+yA2+zA2=4

7、x+y+z=1作閔，并修饰ezplot(，xA2+y八2+(1-x-y)八2-4，，[一2,2,-2,2])CZ3lure1FileEditViewInsertToolsDesktopWindowHelpQQB给I，丨昆I口j^+Z+O-x-y)

8、^=021iiiiiii1.5-_1--0.5--0--•0.5---1---1.5--21IIIII曹I-2-1.5-1-0.500.511.52曲面作图：z=xA2+yA2+sin(xy)，定义域自拟;x=0:100;y=0:100;[x,y]=meshgrid(x,y);z=x.A2+y.A2+sin(x.*y);plot3(x,y,z)8,函数f(x)定义为如下的行列式：

9、XX八2X八3

10、

11、12x3xA2

12、

13、026x

14、求f(x)在x=l处的数值导数，并画图比较数值导数与真实导函数之间的差别symsx;A=[

15、xxA2xA3;12*x3*xA2;026*x];B=det(A);y=diff(B,xsubs(y,1))lims=[—5,5];fplot(，，lims);holdon;fplot(*6*xA2*,lims)CommandWindowy=6*xA2arts二6A»8,当x取值pi/6,2pi/6,".,5pi/6时，求(sin(x))A2+(cos(2x