资源描述:
《数学建模上机练习2011.doc》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、练习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.画
2、出下列函数的曲面及等高线: 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.通过测量得到一组数据:t12345678910y4.8424.3
3、623.7543.3683.1693.0383.0343.0163.0123.005分别采用y=c1+c2e^(-t)和y=d1+d2te^(-t)进行拟合,并画出拟合曲线进行对比。t=1:10;y=[4.8424.3623.7543.3683.1693.0383.0343.0163.0123.005];p=polyfit(exp(-t),y,1)p=5.21653.1564xt=1:0.5:10;yt=polyval(p,exp(-xt))yt=Columns1through75.07554.32043.862
4、43.58463.41613.31393.2520Columns8through143.21443.19163.17773.16943.16433.16123.1593Columns15through193.15823.15753.15713.15683.1567subplot(1,2,1)plot(t,y,'r*',xt,yt,'go’)p1=polyfit(t.*exp(-t),y,1)p1=5.02732.9973t1=1:0.5;10;y2=polyval(p1,t1.*exp(-t1))y2=Empty
5、matrix:1-by-0subplot(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.666713.微分方程组当t=0时,x1(0)=1,x2(0)=-0.5,求微分方程t在[0,
6、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=[00.30.81.11.62.3]; y=[0.50.821.141.251.351.41];分别采用多项式: y=a0+a1t+a2t2和指
7、数函数 y=b0+b1e^t+b2te^t进行拟合,并计算均方误差、画出拟合效果图进行比较。t=[00.30.81.11.62.3];y=[0.50.821.141.251.351.41];>>p1=polyfit(t,y,2)p1=-0.23460.91340.5326t1=1:0.01:2.3;y1=polyval(p1,t1)y1=Columns1through71.21131.21571.22011.22441.22871.23291.2371Columns8through141.24131.24531.
8、24941.25341.25731.26121.2651Columns15through211.26891.27271.27641.28001.28371.28721.2908Columns22through281.29421.29771.30101.30441.30771.31091.3141Columns29through351.31731.32041.32341.3