原创力文档- 1、本文档共12页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数学建模上机练习20110815.
练习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 =
Colu
文档评论（0）