数学建模上机练习2011.docVIP

PAGE PAGE 7 练习1 基础练习 一、矩阵及数组操作： 1．利用基本矩阵产生3×3和15×8的单位矩阵、全1矩阵、全0矩阵、均匀分布随机矩阵（[-1，1]之间）、正态分布矩阵（均值为1，方差为4）。- 2．利用fix及rand函数生成[0，10]上的均匀分布的10×10的整数随机矩阵a，然后统计a中大于等于5的元素个数。 3．在给定的矩阵中删除含有 整行内容全为0的行，删除整列内容全为0的列。 二、绘图：` 4．在同一图形窗口画出下列两条曲线图像： ?????????????? y1=2x+5； y2=x^2-3x+1， 并且用legend标注。 5．画出下列函数的曲面及等高线： ????z=x^2+y^2+sin(by)． 三、程序设计： 6．编写程序计算（x在[-3，3]，间隔0.01） 7．有一列分数序列： 求前15项的和。 8．用至少三种方法编写函数实现求任意整数n的阶乘。 9．将任意大于6的偶数m写成两个素数p1、p2的和（试着写出所有的m=p1+p2的可能形式）。 10．是否任意3的倍数m可以写成两个素数p1、p2、p3的和（试着写出所有的m=p1+p2+p3 的可能形式）？ 四、数据处理与拟合初步： 11．通过测量得到一组数据： t 1 2 3 4 5 6 7 8 9 10 y 4.842 4.362 3.754 3.368 3.169 3.038 3.034 3.016 3.012 3.005 分别采用y=c1+c2e^(-t)和y=d1+d2te^(-t)进行拟合，并画出拟合曲线进行对比。 t = 1:10; y = [4.842 4.362 3.754 3.368 3.169 3.038 3.034 3.016 3.012 3.005]; p = polyfit(exp(-t),y,1) p = 5.2165 3.1564 xt=1:0.5:10; yt=polyval(p,exp(-xt)) yt = Columns 1 through 7 5.0755 4.3204 3.8624 3.5846 3.4161 3.3139 3.2520 Columns 8 through 14 3.2144 3.1916 3.1777 3.1694 3.1643 3.1612 3.1593 Columns 15 through 19 3.1582 3.1575 3.1571 3.1568 3.1567 subplot(1,2,1) plot(t,y,r*,xt,yt,go’) p1=polyfit(t.*exp(-t),y,1) p1 = 5.0273 2.9973 t1=1:0.5;10; y2=polyval(p1,t1.*exp(-t1)) y2 = Empty matrix: 1-by-0 subplot(1,2,2) plot(t,y,r*,t1,y2,go) 12．计算下列定积分： （1）t1=0:0.01:2; z1=exp(-2*t); trapz(t1,z1) ans = 0.4909 （2）t2=0:0.01:2; z2=exp(2*t); trapz(t2,z2) ans = 26.8000 （3）x=-1:0.01:1; z3=x.^2-3.*x+0.5; trapz(x,z3) ans = 1.6667 13．微分方程组 当t=0时，x1(0)=1，x2(0)=-0.5，求微分方程t在[0，25]上的解，并画出相空间轨道图像。 t=0:25; [x1,x2]=dsolve(Dx1=0.5-x1,Dx2=x1-4*x2,x1(0)=1,x2(0)=-0.5,t) x1 =1/2+1/2*exp(-t) x2 =1/8+1/6*exp(-t)-19/24*exp(-4*t) 14．设通过测量得到时间t与变量y的数据： ??????? ?t=[0 0.3 0.8 1.1 1.6 2.3]; ???????????y=[0.5 0.82 1.14 1.25 1.35 1.41]; 分别采用多项式： ? y=a0+a1t+a2t2 和指数函数 ? y=b0+b1e^t+b2te^t 进行拟合，并计算均方误差、画出拟合效果图进行比较。 t=[0 0.3 0.8 1.1 1.6 2.3]; y=[0.5 0.82 1.14 1.25 1.35 1.41]; p1=polyfit(t,y,2) p1 = -0.2346 0.9134 0.5326 t1=1:0.01:2.3; y1=polyval(p1,t1) y1