数学实验作业
练习2﹒1
画出下列常见曲线的图形（其中a=1，b=2,c=3）。 1.
立方抛物线y =
解： x=-4:0.1:4; y=x.^(1/3);
plot(x,y)
-4
-3-2-101234
0.20.40.60.811.21.4
1.6
2.高斯曲线2
x
y e -=
解：
fplot('exp(-x^2)',[-4,4])
-4
-3
-2
-1
1
2
3
4
00.10.20.30.40.50.60.70.80.9
1
3、笛卡儿曲线23
3
2
2
33,(3)11at at x y x y axy t
t
=
=
+=++
解：ezplot('x^3+y^3-3*x*y',[-4,4])
-4
-3-2-1
01234
-4-3-2-10123
4x
y
x 3+y 3-3 x y = 0
或：t=-4:0.1:4; x=3*t./(1+t.^2); y=3*t.^2./(1+t.^2); plot(x,y)
-1.5
-1-0.500.51 1.5
00.5
1
1.5
2
2.5
3
4、蔓叶线233
2
2
2
,()11at
at
x
x y y t
t
a x
=
=
=
++-
解：t=-4:0.1:4; x=t.^2./(1+t.^2); y=t.^3,/(1+t.^2); y=t.^3./(1+t.^2); plot(x,y)
00.10.20.30.40.50.60.70.80.91
-4
-3-2-10123
4
或：
ezplot('y .^2-x.^3/(1-x)',[-4,4])
-4
-3-2-1
01234
-4-3-2-10123
4x
y
y.2-x.3/(1-x) = 0
5、摆线(sin ),(1cos )x a t t y b t =-=-
解：t=-4:0.1:4; x=t-sin(t); y=2-2*cos(t); plot(x,y)
-5
-4-3-2-1012345
00.511.522.533.5
4
6、星形线2
2
2
3
3
333(cos )
,(sin )
()
x a t y a t x y a ==+=
解：t=0:0.1:2*pi; x=cos(t).^3; y=sin(t).^3; plot(x,y)
-1
-0.8
-0.6
-0.4
-0.2
0.2
0.4
0.6
0.8
1
-1-0.8-0.6-0.4-0.200.20.40.60.8
1
或: ezplot('x.^(2/3)+y .^(2/3)-1',[-4,4])
-4
-3-2-1
01234
-4-3-2-10123
4x
y
x.2/3+y.2/3-1 = 0
7、螺旋线cos ,sin ,x a t y b t z ct ===
解：t=0:0.1:2*pi; x=cos(t); y=2*sin(t); z=3*t; plot3(x,y ,z)
-1
1
8、阿基米德螺线r a θ=
解：θ=0:0.1:2*pi; r=;θ polar(θ，r)
90
270
9、对数螺线a
r eθ
=
θ=0:0.1:2*pi;
r=exp(θ);
polar(θ,r)
90
270
1800
10、双纽线2
2
222222
cos 2(()())
r a x y a x y θ
=+=-
解：θ=0:0.1:2*pi;
r=sqrt(cos(2*θ));
90
270
或：ezplot('(x.^2+y.^2).^2-(x.^2-y.^2)',[-4,4])
hold on;
grid
on
-4
-3-2-1
01234
-4-3-2-10123
4x
y
(x.2+y.2).2-(x.2-y.2) = 0
11、双纽线2
2
2
2
2
2
sin 2(()2)r a x y a xy θ
=+=
解：t=0:0.1:2*pi; r=sqrt(sin(2*t)); polar(t,r)
90
270
1800
或：ezplot('(x.^2+y^2).^2-2*x*y',[-4,4])
-4
-3-2-1
01234
-4
-3-2-101234x
y
(x.2+y 2).2-2 x y = 0
12、心形线(1cos )r a θ=+
解：t=0:0.1:2*pi; >> r=1+cos(t); >> polar(t,r)
90270
180
练习2.2
1、求出下列极限值。
(1)
lim
n →∞
解：syms n
limit((n^3+3^n)^(1/n),n,inf) ans =3
（2）lim n →∞
解：syms x
>> limit(sqrt(x+2)-2*(sqrt(x+1))+sqrt(x),x,inf)
ans =0
（3）0
lim cot 2n x x → 解：syms x
limit(x*cot(2*x),x,0) ans =1/2
（4）lim (cos
)
x
n m x
→∞
解：syms x m
limit((cos(m/x))^x,x,inf) ans =1
（5）1
11lim (
)1
x
n x
e →-
-
解：syms x
limit((1/x)-1/(exp(x)-1),x,1) ans =(exp(1)-2)/(exp(1)-1)
（6）lim
)n x →∞
解：syms x
limit(sqrt(x^2+x)-x,x,inf) ans =1/2
2、有个客户看中某套面积为1802m ，每平方米7500元。他计划首付30%，其余70%用20年按揭贷款（贷款年利率5.04%），按揭贷款中还有10万元为公积金贷款（贷款年利率4.05%），请问他的房屋总价、首付款额和月付还款额分别为多少？
解：（1）房屋总价：18075001350000
S
=?= (元)
（2）首付款额：13500000.3405000N =?= （元）
（3）房屋未付钱：945000M S N =-=
(元)
设揭贷款的年利率为x ，则20
(1)
240
a x y ?+=
其中a 为本金，y 为
每月所付的钱。
解：当a =945000-100000=845000， 5.04%x =时， syms x y
y=845000*(1+x)^20/240; x=0.0504; eval(y)
ans = 9.4133e+003
当a =100000， 4.05%x =时； syms x y
y=100000*(1+x)^20/240; x=0.0405; eval(y) ans =921.7867 即每月付还款额为9413.3921.786710335.0867
Z
=+= (元)
3、作出下列函数及其导函数的图形，观察极值点、最值点的位置点的位置并求出，求出所有驻点以及对应的二阶导函数，求出函数的单调区间。 （1）2
2
()sin(2),[2,2];f x x x x =
---
解：函数图象：
fplot('x.^2*sin(x.^2-x-2)',[-2,2])
-2
-1.5-1-0.500.51 1.52
-4-3
-2
-1
1
2
3
原函数在-1附近的极小值：
[x,f]=fminsearch('x.^2*sin(x.^2-x-2)',-1) x = -0.7315 f =-0.3582
原函数在1.5附近的极小值： [x,f]=fminsearch('x.^2*sin(x.^2-x-2)',1.5) x =1.5951 f =-2.2080
原函数在-1.5附近的极大值：
[x,f]=fminsearch('-x.^2*sin(x.^2-x-2)',-1.5)
x =-1.5326
f =2.2364
原函数在0附近的极大值：
[x,f]=fminsearch('-x.^2*sin(x.^2-x-2)',0) x =0
f =0
原函数在[-2,2]上的最小值：
x=-2:0.1:2;
y=x.^2.*sin(x.^2-x-2);
[m,k]=min(y)
m =-3.0272
k =1
原函数在[-2,2]上的最大值：
x=-2:0.1:2;
y=x.^2.*sin(x.^2-x-2);
[m,k]=max(y)
m =2.2140
k =6
原函数的导函数图像：
syms x
y=x^2*sin(x^2-x-2);
diff(y,x)
ans =2*x*sin(x^2-x-2)+x^2*cos(x^2-x-2)*(2*x-1)
导函数在-1.5附近的极小值：
[x,f]=fminsearch('2*x*sin(x^2-x-2)+x^2*cos(x^2-x-2)*(2*x-1)',-1.5) x =-1.2650
f =-5.5890
导函数在1.5附近的极小值：
[x,f]=fminsearch('2*x*sin(x^2-x-2)+x^2*cos(x^2-x-2)*(2*x-1)',1.5)
x =1.2404
f =-2.7572
导函数在-2附近的极大值：
[x,f]=fminsearch('-(2*x*sin(x^2-x-2)+x^2*cos(x^2-x-2)*(2*x-1))',-2) x =-1.9240
f =17.6746
导函数在-0.5附近的极大值：
[x,f]=fminsearch('-(2*x*sin(x^2-x-2)+x^2*cos(x^2-x-2)*(2*x-1))',-0.5) x =-0.4742
f =0.7973
导函数在[-2,2]上的最大值：
x=-2:0.1:2;
y=2*x.*sin(x.^2-x-2)+x.^2.*cos(x.^2-x-2).*(2*x-1);
[m,k]=max(y)
数学软件MATLAB实验作业
数学软件与数学实验作业 一.《数学软件》练习题（任选12题，其中19－24题至少选2题）： 3.对下列各式进行因式分解. (1). syms x y >> factor(x^5-x^3) (2). syms x y >> factor(x^4-y^4) (3). syms x >> factor(16-x^4) (4). syms x >> factor(x^3-6*x^2+11*x-6) (5). syms x y >> factor((x+y)^2-10*(x+y)+25) (6). syms x y >> factor(x^2/4+x*y+y^2) (7). syms x y a b >> factor(3*a*x+4*b*y+4*a*y+3*b*x) (8). syms x >> factor(x^4+4*x^3-19*x^2-46*x+120) 5.解下列方程或方程组. (1).solve('(y-3)^2-(y+3)^3=9*y*(1-2*y)') (2). solve('3*x^2+5*(2*x+1)') (3). solve('a*b*x^2+(a^4+b^4)*x+a^3*b^3','x') (4). solve('x^2-(2*m+1)*x+m^2+m','x') (5). [x,y]=solve('4*x^2-9*y^2=15','2*x-3*y=15') 6.计算极限. (1). syms x f=(exp(x)-exp(-x))/sin(x); limit(f,x,0) (2) syms x >> f=(x/(x-1)-1/log(x)); >> limit(f,x,1) (3). syms x >> f=(1-cos(x))/x^2; >> limit(f,x,0)
数学实验上机汇总未完成
数学实验上机作业整理∈hyd 实验一 1. 计算球体体积（半径r=5） r=5;v=(4/3)*pi*r^3 v =523.5988 2.设矩阵1234567891023416A ?? ? = ? ??? （1）提取A 的第2列赋值给B; A=[1 2 3 4 5;6 7 8 9 10;2 3 4 1 6];B=A(:,2) B = 2 7 3 （2）提取A 的第2行前3个数给C ； A=[1 2 3 4 5;6 7 8 9 10;2 3 4 1 6];C=A(2,[1,2,3]) C = 6 7 8 （3）提取A 第1,3行和2, 4列相交位置元素构成子矩阵D ； A=[1 2 3 4 5;6 7 8 9 10;2 3 4 1 6];D=A([1,3],[2,4]) D = 2 4 3 1 （4）构造矩阵E 使得E 的结构为：132213C E D C ???? ?= ? ?? A=[1 2 3 4 5;6 7 8 9 10;2 3 4 1 6];E=[D [C;C]] E = 2 4 6 7 8 3 1 6 7 8
（5）把A 中间的8换为0； A(2,3)=0;A A = 1 2 3 4 5 6 7 0 9 10 2 3 4 1 6 （6）去掉A 的第2行； A=[1 2 3 4 5;6 7 8 9 10;2 3 4 1 6]; A(2,:)=[] A = 1 2 3 4 5 2 3 4 1 6 3.写出完成下列操作的命令 （1） 建立10阶单位矩阵A; A=eye(10) （2）建立5×6的随机矩阵A ，其元素为[100，200]范围内的随机数； A=rand(5,6)*100+100 （3）将A 对角线元素加30 A+eye(5,6)*30 4.（选做题）设有分块矩阵333223E R A O S ????? =? ??? ，其中E,R,O,S 分别为单位矩阵、随机矩阵、零矩阵和对角矩阵，试通过数值计算验证2 2 E R RS A O S +?? =? ??? 。 S=[1 1;1 1]; E=eye(3);R=rand(3,2); O=zeros(2,3); [E R;O S]^2 [E R+R*S;O S^2] 实验二 1.设矩阵1215346562A -?? ? = ? ?-?? （1）求A 的秩、A 的每个元素3次方; A=[1 2 -1;5 34 6;-5 6 2];
数学实验期末考试
数学实验 期末上机考核 学号201519030102 姓名曹欣辉年级专业2015级水产养殖学 学号201519030103 姓名陈妙珊年级专业2015级水产养殖学 学号201519030104 姓名杜日臻年级专业2015级水产养殖学 学号姓名年级专业 学号姓名年级专业
注意事项： 1、考核方式：以组（3~4人）为单位，请于指定时间内开卷完成所布置的任务，地点为实验 室机房或课室。 2、发题时间为6月25号早上8:30，请到如下邮箱提取题目，账号：nongkeshuxue@https://www.wendangku.net/doc/075488855.html,， 密码：shuxueshiyan。 3、关于试卷提交时间： （1）电子版提交时间于6月26晚上12:00前，发送本班任课老师给定的email地址，任课老师以收到信件的时间为准，提交文件的同学可通过收到任课老师回复的邮件接收函作为提交信息。 （2）纸质版提交时间于6月27日早上11:30前，由学委收齐后交与任课老师。 4、每小组同学可以使用无生命的数据或材料：如计算机、软件、参考文献、网络、图书等。 5、除小组成员内相互讨论，队伍成员不可以向老师及其他人员寻求帮助。任何从非小组成员 内得到的帮助都是被严格禁止的，这包括通过邮件，电话交谈，聊天，网络聊天等其他交流工具得到的他人的帮助。 6、每位同学需在承诺书上签字，如无签字，可视为放弃该科目考试，并且一经发现抄袭作弊 等行为，将取消该组所有同学的答卷分数。 7、每组同学完成答题后，请组内同学根据所作贡献协商讨论后进行评价打分，每组同学贡献 值总分为100。 8、请在下列表中有学生姓名的地方填上相应的名字。 组内同学互评后贡献值表：注：贡献值≤100 学生姓名曹欣辉陈妙珊杜日臻张照明 贡献值 老师评分表： 题号 1 2 3 4 5 7 总分得分 签名 学生成绩： 学生姓名 成绩 注：表中每位学生成绩得分计算公式如下： 该学生贡献分 卷面总分 该组最高贡献分
数学实验作业
练习2﹒1 画出下列常见曲线的图形（其中a=1，b=2,c=3）。 1. 立方抛物线y = 解： x=-4:0.1:4; y=x.^(1/3); plot(x,y) -4 -3-2-101234 0.20.40.60.811.21.4 1.6 2.高斯曲线2 x y e -= 解： fplot('exp(-x^2)',[-4,4])
-4 -3 -2 -1 1 2 3 4 00.10.20.30.40.50.60.70.80.9 1 3、笛卡儿曲线23 3 2 2 33,(3)11at at x y x y axy t t = = +=++ 解：ezplot('x^3+y^3-3*x*y',[-4,4])
-4 -3-2-1 01234 -4-3-2-10123 4x y x 3+y 3-3 x y = 0 或：t=-4:0.1:4; x=3*t./(1+t.^2); y=3*t.^2./(1+t.^2); plot(x,y)
-1.5 -1-0.500.51 1.5 00.5 1 1.5 2 2.5 3 4、蔓叶线233 2 2 2 ,()11at at x x y y t t a x = = = ++- 解：t=-4:0.1:4; x=t.^2./(1+t.^2); y=t.^3,/(1+t.^2); y=t.^3./(1+t.^2); plot(x,y)
00.10.20.30.40.50.60.70.80.91 -4 -3-2-10123 4 或： ezplot('y .^2-x.^3/(1-x)',[-4,4])
清华大学2002至2003学年第二学期数学实验期末考试试题A
清华大学2002至2003学年第二学期数学实验期末考试试题A 数学实验试题 2003.6.22 上午 (A卷；90分钟) 一. 某两个地区上半年6个月的降雨量数据如下（单位：mm）： 月份123456 地区A259946337054 地区B105030204530 在90%的置信水平下，给出A地区的月降雨量的置信区 间： 在90%的置信水平下，A地区的月降雨量是否不小于70（mm）？ 在90%的置信水平下，A、B地区的月降雨量是否相同？ A地区某条河流上半年6个月对应的径流量数据如下（单位：m3）：110，184，145，122，165，143。该河流的径流量y与当地的降雨量x的线性回归方程为；若当地降雨量为55mm，该河流的径流量的预测区间为（置信水平取90%）。 答案：（程序略） (1) [32.35，76.65] (2) 是 (3) 否 (4) y=91.12+0.9857x (5) [130.9,159.7] 二．（10分） （1）（每空1分）给定矩阵，如果在可行域上考虑线性函数，其中，那么的最小值是，最小点为；最大值是，最大点为。 （每空2分）给定矩阵，，考虑二次规划问题，其最优解为，（2） 最优值为，在最优点处起作用约束 为 。 答案：（1）最小值为11/5，最大值为7/2，最小点为（0，2/5，9/5），最大点为（1/2，0，3/2）。 （2）最优解为（2.5556，1.4444），最优值为–1.0778e+001，其作用约束为。 三．（10分）对线性方程组：，其中Ａ＝，b= （3分）当时，用高斯—赛德尔迭代法求解。取初值为，写出迭代第4步的结果=____________________。 （4分）当时，用Jacobi 迭代法求解是否收敛？__________ , 理由是_________________________________________________ 。 （3分）求最大的c, 使得对任意的，用高斯—赛德尔迭代法求解一定收敛，则c应为__________。 答案：（1）x = [ -1.0566 1.0771 2.9897]
数学实验作业题目(赛车跑道)
数学实验报告实验题目：赛车车道路况分析问题 小组成员： 填写日期2012 年 4 月20 日
一．问题概述 赛车道路况分析问题 现要举行一场山地自行车赛，为了了解环行赛道的路况，现对一选手比赛情况进行监测，该选手从A地出发向东到B，再经C、D回到A地（如下图）。现从选手出发开始计时，每隔15min观测其位置，所得相应各点坐标如下表（假设其体力是均衡分配的）： 由D→C→B各点的位置坐标（单位：km） 假设：1. 车道几乎是在平原上，但有三种路况（根据平均速度v(km/h)大致区分）： 平整沙土路（v>30）、坑洼碎石路（10 2．估计车道的长度和所围区域的面积； 3．分析车道上相关路段的路面状况（用不同颜色或不同线型标记出来）； 4．对参加比赛选手提出合理建议. 二．问题分析 1.模拟比赛车道的曲线：因为赛道散点分布不规则，我们需要用光滑曲线来近 似模拟赛道。由于数据点较多，为了避免龙格现象，应采用三次样条插值法来对曲线进行模拟（spline命令）。全程曲线为环路，我们需要对上下两部分分别 模拟，设模拟出的曲线为P：。 2.把A到B点的曲线分成若干小段： 赛道的路程L：取dL=，对模拟出的整条曲线求线积分，即 所围区域的面积：用上下部分曲线的差值对求定积分，即 3.用样条插值法模拟出比赛车道曲线后，根据曲线分别计算出原数据中每两点 （）间的路程，即求线积分 由于每两点间时间间隔相同且已知（15min），故可求出每段路程的平均速度 易知即为的积分中值 将此速度近似作为两点间中点时刻的速度，然后再次采用样条插值法，模拟出全过程的图像。而根据求出的与之间的关系，再次采用样条插值法，即可模拟出全过程的图像 4. 由赛道曲线可求出赛道上任一点到点的路程 同时图像也可以求出赛道上任一点到点的路程 数学建模作业——实验1 学院：软件学院 姓名： 学号： 班级：软件工程2015级 GCT班 邮箱： 电话： 日期：2016年5月10日 基本实验 1.椅子放平问题 依照1.2.1节中的“椅子问题”的方法，将假设中的“四腿长相同并且四脚连线呈正方形”，改为“四腿长相同并且四脚连线呈长方形”，其余假设不变，问椅子还能放平吗？如果能，请证明；如果不能，请举出相应的例子。 答：能放平，证明如下： 如上图，以椅子的中心点建立坐标，O为原点，A、B、C、D为椅子四脚的初始位置，通过旋转椅子到A’、B’、C’、D’，旋转的角度为α，记A、B两脚，C、D两脚距离地面的距离为f(α)和g(α)，由于椅子的四脚在任何位置至少有3脚着地，且f(α)、g(α)是α的连续函数，则f(α)和g(α)至少有一个的值为0，即f(α)g(α)=0，f(α)≥0，g(α)≥0，若f(0)＞0，g(0)=0， 则一定存在α’∈（0，π），使得 f(α’)=g(α’)=0 令α=π（即椅子旋转180°,AB 边与CD 边互换），则 f(π)=0，g(π)＞0 定义h(α)= f(α)-g(α)，得到 h(0)=f(0)-g(0)＞0 h(π)=f(π)-g(π) ＜0 根据连续函数的零点定理，则存在α’∈（ 0，π），使得 h(α’)= f(α’)-g(α’)=0 结合条件f(α’)g(α’)=0,从而得到 f(α’)=g(α’)=0，即四脚着地，椅子放平。 2. 过河问题 依照1.2.2节中的“商人安全过河”的方法，完成下面的智力游戏：人带着猫、鸡、米过河，船除需要人划之外，至多能载猫、鸡、米之一，而当人不在场时，猫要吃鸡、鸡要吃米，试设计一个安全过河的方案，并使渡河的次数尽量的少。 答： 用i =1,2,3,4分别代表人，猫，鸡，米。1=i x 在此岸，0 =i x 在对岸， ()4321,,,x x x x s =此岸状态，()43211,1,1,1x x x x D ----=对岸状态。安全状态集合为 : 一． 1. 1/e 2. 3 3.1 4.e3 5. ∞ 6. 0 7.∞ 8.0 9.1/2 10.0 11.e2c12.不存在13. 1/12 Matlab实验过程： 1.1/exp(1) syms n; f=(1-1/n)^n; limit(f,n,inf) ans = 1/exp(1) 2.3 syms n; f=(n^3+3^n)^(1/n); limit(f,n,inf) ans = 3 3. 1 syms n; f=(1+sin(2*n))/(1-cos(4*n)); limit(f,n,pi/4) ans = 1 4.e^3 syms x; f=(1+cos(x))^(3*sec(x)); limit(f,x,pi/2) ans = exp(3) 5.inf syms x; f=(x^2)*exp(1/(x^2)); limit(f,x,0) ans = Inf 6.0 syms x; f=(x^2-2*x+1)/(x^3-x); limit(f,x,1) ans = 7.inf syms x; f=((2/pi)*atan(x))^x; limit(f,x,+inf) ans = Inf 8.0 syms x y; f=(1-cos(x^2+y^2))/((x^2+y^2)*exp(x^2+y^2)); limit(limit(f,x,0),y,0) ans = 9.1/2 syms x; f=(1-cos(x))/(x*sin(x)); limit(f,x,0) ans = 1/2 10.0 syms x; f=atan(x)/(2*x); limit(f,x,inf) ans = 11.exp(2*c) syms c; f=sym('((x+c)/(x-c))^x'); limit(f,'x',inf) ans = exp(2*c) 12.极限不存在 syms x; f=cos(1/x); limit(f,x,0) ans = limit(cos(1/x), x = 0) 13.1/12 syms x; f=1/(x*log(x)^2)-1/(x-1)^2; limit(f,x,1) ans = 1/12 二．观察函数logbx,当b=1/2,1/3,1/4和b=2,3,4时函数的变化特点，总结logbx的图形特点。 2014-2015学年第一学期数学实验上机试卷 一、上机操作题 1. 画出以下函数图形（要求写出程序和结果）： ⑴ 3411()2 1 x x x f x x x ?++≥=? +>求 ⑵22(sin )(1cos )x a t t d y y a t dx =-??=-? 求 ⑶ 2cos (sin )'x y x y =求 ⑷22 2''02 1 (,),x y x y x y u u x y e u u x +==?=?求及 结果：⑴ ⑵ ⑶ ⑷ 4. 计算下列积分（要求写出程序和计算结果）： ⑴ 211ln 11x dx x x +--? ⑵2220sec 2tan x dx x π+? (3) 2,02}x x ≤≤≤?? 2其中D={(x,y):y (4) 2221 L dl x y z ++?其中L 为空间螺旋线cos ,sin ,(02,0)x a t y a t z bt t b π===≤≤> . (5) 222()S x y z dS ++?? 其中S 是球2222x y z az ++=. (6) 22S x dydz y dzdx +?? 其中S 为球面2222()()()x a y b z c R -+-+-=的外侧. 结果：⑴ ⑵ ⑶ ⑷ ⑸ ⑹ 5. 判断以下级数的敛散性： ⑴ 1()21n n n n +∞ =+∑ ⑵ 12n n n x n +∞ =∑ 结果：⑴ ⑵ 6. 用两种以上的方法求解下列方程组： 1234234 1242342344331733 x x x x x x x x x x x x x -+-=??-+=-?? ++=??-++=-? 结果： 二、写出解题的思想，计算过程和程序，结果及分析等内容. ⑴在某化学反应里，由实验得到生成物的浓度y 与时间t 有如下关系，求浓度与时间的关系的拟合函数．（30分） ⑵某公司刊登广告：“现有一栋住宅楼，每套只需自备七万元，其余由公司贷款，贷款可分期偿还，每月只需800元，十年还清。”现在的问题是如果一次性付清该付多少（即该房屋的实际价格是多少）？如果贷款，买房人实际借了多少钱？（假设月利率为1%）（40分） 数学实验 期 末 作 业 学号： 班级： 姓名： 1. 求函数x x y 2sin 3=的5阶导数。 2. 使用sparse 命令描述? ? ???? ? ? ??30001 020******* 01020 10003。 3. 求解边值问题 1)0(,0)0(,34,43==+-=+=g f g f dx dg g f dx df 。 4. 建立函数1 2sin )(3-=x x f x 的M-文件，并计算)2(f 和)10(f 。 5. 计算二重积分dy dx x y ??211 0][。 6. 已知数列满足2,11 01=+= +a ka a k k ，求5a ，并要求最后结果分别以小数点后两位和有理数这两种数据显示格式输出。 7. 大约在1500年前，《孙子算经》中就记载了这个有趣的问题：“今有鸡兔同笼，上有三十五头，下有九十四足，问鸡兔各几何？”请根据你的思路编程求解。 8. 绘制以下方程所表示的图形。 （1）x x y -=23 2 （2）y z cos =绕z 轴的旋转曲面 （3）)40(,) 2sin(sin )]2cos(4[cos )]2cos(4[π< 10.根据中华人民共和国个人所得税法规定：公民的个人工资、薪金应依法缴纳个人所得税。所得税计算办法为：在每个人的月收入中超过2000元以上的部分应该纳税，这部分收入称为应纳税所得额。应纳税所得额实行分段累计税率，按下列税率表计算： 个人所得税税率表： 等级全月应纳税所得额税率（%） 1 不超过500元的部分 5 2 超过500元，不到2000元的部分10 3 超过2000元，不到5000元的部分15 4 超过5000元，不到20000元的部分20 5 超过20000元，不到40000元的部分25 6 超过40000元，不到60000元的部分30 7 超过60000元，不到80000元的部分35 8 超过80000元，不到100000元的部分40 9 超过100000元的部分45 若某人的工资是x元，试建立税款y与收入x之间的M-文件，并要求程序运行时可以告知操作者“please input the number of your wage”。 实验十三回归分析 电61 张俊翔2016010891 13.5 （1）首先对于所给数据，分别画出y关于三个因素x1、x2、x3的散点图如下：犯罪率y关于年收入低于5000美元家庭的百分比x1： 犯罪率y关于失业率x2： 犯罪率y关于人口总数x3： 由上图可以看出，y关于x1、x2应该有线性关系，而与x3无明显的相关性。 由此选取y关于x1、x2、x3的线性模型进行拟合。即 Y=β0+β1*x1+β2*x2+β3*x3 首先选取x1、x2作拟合，程序如下： n=20; X=[ones(n,1),x1',x2']; [b,bint,r,rint,s]=regress(y',X); b,bint,s 三者比较可知，最好的模型是只选择x1、x2的情况，此时决定系数最大，剩余方差最小，而且不存在系数的置信区间包含零的情况。 β3的置信区间包含零点，说明x3对y几乎没有什么影响，因此包含3个自变量的模型并没有比只含x1、x2的模型好。 因此选择最终模型是只含x1、x2的模型。 表达式为y=-34.0725+1.2239*x1+4.3989*x2 （3）对最终模型用rcoplot命令观察残差，可得下面的图形： 可见剩余方差和决定系数都有了明显的改进。此时的残差图如下： 这时不再有异常数据点，表达式为：y=-35.7095+1.6023*x1+3.3926*x2 13.10 首先假设风险偏好度对人寿保险额没有二次效应，两个自变量对人寿保险额也没有交互效应，来看已经确定的影响因素的系数： 由于已知经理的年均收入和人寿保险额之间存在着二次关系，而风险偏好度对人寿保险额有线性效应，因此模型为： Y=β0+β1*x1+β2*x2+β3*x1^2 程序如下（数据输入略）： n=18; xx1=x1.^2; xx2=x2.^2; xx=x1.*x2; X=[ones(n,1),x1',x2',xx1']; [b,bint,r,rint,s]=regress(y',X); b,bint,s rcoplot(r,rint) 数学实验作业汇总 文稿归稿存档编号：[KKUY-KKIO69-OTM243-OLUI129-G00I-FDQS58- （1）产生一个5阶魔方矩阵M：M=magic(5) （2）将矩阵M的第3行4列元素赋值给变量t：t=M(3,4) （3）将由矩阵M第2，3，4行第2,5列构成的子矩阵赋给变N：N=M(2:4,2:3:5) （4）将由矩阵M的前3行赋给变量N：N=M(1:3,:) （5）将由矩阵M的后3列赋给变量N：N=M(:,end:-1:end-2) （6）提取M的主对角线元素，并以这些对角线元素构成对角矩阵N：N=diag(diag(M))或N=tril(triu(M)) (7)随机产生1000个100以内的整数赋值给变量t：t=round(rand(1,1000)*100) （8）随机产生100*5个100以内的实数赋值给变量M：M=rand(100,5)*100 （1）删除矩阵M的第7个元素M(7)=[] （2）将含有12个元素的向量t转换成3*4的矩阵：reshape(t,3,4) （3）产生和M同样大小的单位矩阵：eye(size(M)) （4）寻找向量t中非零元素的下标：find(t) （5）逆序显示向量t中的元素：t(end:-1:1) （6）显示向量t偶数位置上的元素：t(2:2:end) （7）利用find函数，将向量t中小于10的整数置为0：t(find(t<10&rem(t,1)==0))=0（8）不用find函数，将向量t中小于10的整数置为0：t(t<10&rem(t,1)==0)=0 （9）将向量t中的0元素用机器0（realmin）来代替：t(find(t=0))=realmin （10）将矩阵M中小于10的整数置为0：M(find(M<10)&rem(M,1)==0)=0 数学实验试题 2003.6.22 上午 班级姓名学号得分 说明： （1）第一、二、三题的答案直接填在试题纸上； （2）第四题将数学模型、简要解题过程和结果写在试题纸上；卷面空间不够时，可写在背面； （3）考试时间为90分钟。 一．（10分，每空2分）（计算结果小数点后保留4位有效数字） 地区的月降雨量的置信区间： （2）在90%的置信水平下，A地区的月降雨量是否不小于70（mm）？ （3）在90%的置信水平下，A、B地区的月降雨量是否相同？ （4）A地区某条河流上半年6个月对应的径流量数据如下（单位：m3）：110，184，145，122，165，143。该河流的径流量y与当地的降雨量x的线性回归方程为；若当地降雨量为55mm，该河流的径流量的预测区间为（置信水平取90%）。 二．（10分） （1）（每空1分）给定矩阵，如果在可行域上考虑线性函数，其中，那么的最小值是，最小点为；最大值是，最大点为。 （2）（每空2分）给定矩阵，，考虑二次规划问题，其最优解 为，最优值为，在最优点处起作用约束为。 三．（10分）对线性方程组：，其中Ａ＝，b= (1)（3分）当时，用高斯—赛德尔迭代法求解。取初值为， 写出迭代第4步的结果=____________________。 (2)（4分）当时，用Jacobi 迭代法求解是否收敛？__________ , 理由是_________________________________________________ 。 (3)（3分）求最大的c, 使得对任意的，用高斯—赛德尔迭代法求解一 定收敛，则c应为__________。 四．（20分）一个二级火箭的总重量为2800公斤。第一级火箭的重量为1000公斤，其中燃料为800公斤。第一级火箭燃料燃烧完毕后自动脱落，第二级火箭立即继续燃烧。第二级火箭中的燃料为600公斤。假设火箭垂直向上发射，两级火箭中的燃料同质，燃烧率为15公斤/秒，产生的推力为30000牛顿。火箭上升时空气阻力正比于速度的平方，比例系数为0.4公斤/米。 （1）建立第一级火箭燃烧时火箭运行的数学模型，并求第一级火箭脱落时的高度、速度和加速度； （2）建立第二级火箭燃烧时火箭运行的数学模型，并求火箭所有燃料燃烧完毕瞬间的高度、速度、和加速度。 （提示：牛顿第二定律f=ma，其中f为力，m为质量，a为加速度。重力加速度9.8米/平方秒。） 数学实验作业一 对以下问题,编写M文件: (1)用起泡法对10个数由小到大排序. 即将相邻两个数比较,将小的调到前头. 解： 代码如下： zuoye1 clear all;clc; a=[7 2 1 0 9 4 5 -3 8 6]; n=length(a); for ii=1:n-1 if a(ii+1)>=a(ii) t1=a(ii); a(ii)=a(ii+1); a(ii+1)=t1; end for jj=1:n-1 if a(jj+1)>=a(jj) t2=a(jj); a(jj)=a(jj+1); a(jj+1)=t2; end end end a 运行结果显示如下： a = 9 8 7 6 5 4 2 1 0 -3 (2)有一个 矩阵,编程求出其最大值及其所处的位置. 解： 代码如下：zuoye2.m clear; clc; a=[1 2 3 4 5 3 4 5 6 9 6 7 8 8 0 1 2 4 5 6] max=-1; flage1=0; flage2=0 for i=1:4 for j=1:5 if (a(i,j)>max) t=max; max=a(i ,j); a(i,j)=t; flage1=i; flage2=j ； end end end max flage1 flage2 运行结果显示如下： a = 1 2 3 4 5 3 4 5 6 9 6 7 8 8 0 1 2 4 5 6 flage2 = max = 45′ 9 flage1 = 2 flage2 = 5 结果： （3）编程求∑=20 1 !n n 。 解： 代码如下：zuoye3.m clear; clc; sum=0; for i=2:11 sum=sum+gamma(i); end sum 电子科技大学二零零八到二零零九学年第二学期期末考试《数学实验》课程考试题A卷(120分钟) 考试形式：闭卷考试日期：2009年7月8日 一、单项选择题(20分) 1、三阶幻方又称为九宫图，提取三阶幻方矩阵对角元并构造对角阵用( ) (A) diag(magic(3)); (B) diag(magic); (C) diag(diag(magic(3))); (D) diag(diag(magic))。 2、MATLAB命令P=pascal(3)将创建三阶帕斯卡矩阵，max(P)的计算结果是( ) (A) 1 2 3 (B) 1 2 1 (C) 3 6 10 (D) 1 3 6 3、命令J=*1;1;1+**1,2,3+;A=j+j’-1将创建矩阵( ) (A) 123 234 345 ?? ?? ?? ?? ?? ; (B) 234 345 456 ?? ?? ?? ?? ?? (C) 123 123 123 ?? ?? ?? ?? ?? (D) 111 222 333 ?? ?? ?? ?? ?? 4、data=rand(1000,2);x=data(:,1);y=data(:,2);II=find(y 1.取不同的初值计算下列平方和形式的非线性规划，尽可能求出所有局部极小点，进 而找出全局极小点，并对不同算法（搜索方向、搜索步长、数值梯度与分析梯度等）的结 果进行分析、比较。 (2). ( )( ) 2 2 2 22 121212min 12114949812324681x x x x x x +-++++-, (4).()()212222 23 12123min10010,1x x x x x x θ??????-++-+?????????????? ，其中 ()()()21112211 1 arc ,02,11arc ,0 22tg x x x x x tg x x x π θπ ?>??=??+ 《数学实验》报告 实验名称 matlab拟合与插值学院机械工程学院 专业班级 姓名 学号 2011年 10月 一、【实验目的】 掌握Matlab关于采用最小二乘法拟合曲线的方法。学会使用matlab求实际中得到数据的插值曲线。 二、【实验任务】 P130第8、10、12题 三、【实验程序】 P130第8题： x=[0.10,0.30,0.40,0.55,0.70,0.80,0.95]; y=[15,18,19,21,22.6,23.8,26]; p1=polyfit(x,y,1); p3=polyfit(x,y,3); p5=polyfit(x,y,5); disp('一阶拟合函数'),f1=poly2str(p1,'x') disp('三阶拟合函数'),f3=poly2str(p3,'x') disp('五阶拟合函数'),f5=poly2str(p5,'x') x1=0.1:0.0017:0.95; y1=polyval(p1,x1); y3=polyval(p3,x1); y5=polyval(p5,x1); plot(x,y,'rp',x1,y1,'--',x1,y3,'k-.',x1,y5); legend('拟合点','一次拟合','三次拟合','七次拟合') P130第10题 x=[10,15,20,25,30]; y=[25.2,29.8,31.2,31.7,29.4]; xi=10:.5:30; yi1=interp1(x,y,xi,'*nearest'); yi2=interp1(x,y,xi,'*linear'); yi3=interp1(x,y,xi,'*spline'); yi4=interp1(x,y,xi,'*cubic'); plot(x,y,'ro',xi,yi1,'--',xi,yi2,'-',xi,yi3,'k.-',xi,yi4,'m:') ,grid on 《数学实验》期末样题(90分钟) 学号：0906010036 姓名：马德祥 一、符号计算（40'） 1.符号表达式的因式分解：factor 题目：x 2 – 1 2.符号表达式的展开：expand 题目：(x-2)(x+3) 3.求复合函数：compose 题目：已知：f(x)= 2 11 x +,g(y)=e y 求：f(g(y)) 4.求反函数：finverse 题目：求y=e x 的反函数 5.求极限：limit 题目：求x x x sin lim 0→ 的值 6.求导数：diff 题目：求函数y=e x 2 的一阶导数y ’ 7.求积分：int 题目：计算：S=? 5 xdx 的值 8：计算级数的和：symsum 题目：求级数： ∑∞ =1 21 n n 的和 9：泰勒级数展开：taylor 题目：求函数f(x)=lnx 在x=1处的taylor 展开式的前5项 syms x y y=log(x) taylor(y,x,5,1) y = log(x) ans = x-1-1/2*(x-1)^2+1/3*(x-1)^3-1/4*(x-1)^4 10.解方程：solve 题目：求方程：ax 2+bx+c=0的根，a 不为0. x=solve('a*x^2 + b*x + c') x = 1/2/a*(-b+(b^2-4*a*c)^(1/2)) 1/2/a*(-b-(b^2-4*a*c)^(1/2)) 二、基本数据输入与计算(20') 1.复数 题目：保存复数3+4i ，输出其实部、虚部，并计算其模 z=3+4i real(z) imag(z) abs(z) z = 3.0000 + 4.0000i ans = 3 ans = 4 ans = 5 2.矩阵输入 题目：把矩阵???? ? ??087654321保存给变量A A=[1 2 3;4 5 6;7 8 0] A = 1 2 3 4 5 6 7 8 0 3.矩阵的基本运算（转置'、乘*、除/\、逆inv 、行列式det 、秩rank ） 题目：把矩阵的转置赋值给变量B ； 4. 问题： 某公司将3种不同含硫量的液体原料（分别记为甲、乙、丙）混合生产两种产品（分别记为A,B ）。按照生产工艺的要求，原料甲、乙必须首先导入混合池中混合，混合后的液体再分别与原料丙混合生产A,B 。一直原料甲、乙、丙的含硫量分别是3%，1%，2%，进货价格分别为6千元/t ，16千元/t ，10千元/t ；产品A,B 的含硫量分别不能超过2.5%，1.5%，售价分别为9千元/t ，15千元/t 。根据市场信息，原料甲、乙、丙的供应量都不能超过500t ；产品A,B 的最大市场需求量分别为100t ，200t 。 (1) 应如何安排生产？ (2) 如果产品A 的最大市场需求量增长为600t ，应如何安排生产？ (3) 如果乙的进货价格下降为13千元/t ，应如何安排生产？分别对(1)、(2)两种情况进 行讨论。 模型： （只考虑问题1，问题2,3只需改变一些约束条件） 设生产时使用原料甲、乙分别为12,x x t ，分别取混合后的液体34,x x t 再加入原料丙 56,x x t 生产产品A,B 。 有质量守恒，可得 1234x x x x +=+ 甲乙混合后的液体的含硫量可表示为 12 12 3%x x x x ++，根据含硫量的要求，可得 12 353512 124646 12 3%*2%* 2.5%*()3%*2%* 1.5%*() x x x x x x x x x x x x x x x x +?+≤+?+?? +?+≤+?+? 根据市场的限制，易得 12563546500 500500100200 x x x x x x x x ≤?? ≤?? +≤??+≤??+≤? 当然还有非负约束 123456,,,,,0x x x x x x ≥ 公司的净利润为（单位：千元）： 35461256123456 9()15()61610()6169155z x x x x x x x x x x x x x x =+++---+=--++-+ 电子科技大学成都学院 2013—2014学年第一学期期末试卷 数学实验 课程考试题B （120 、单项选择题（共 40分，每题4分） 1、提取5阶幻方矩阵对角元并求对角元之和 （） (C) sum(diag(diag(magic(5)))); (D) diag(diag sum(diag(magic(5)))) 2、data=rand(20000,2);x=data(:,1);y=data(:,2);n= find(y数学建模作业——实验1
北理工数学实验作业
数学实验期末考试上机考试
山东建筑大学数学实验期末作业matlab
数学实验 作业10
数学实验作业汇总终审稿)
大学数学实验—期末考试试题6
数学实验作业一
电子科技大学《数学实验》2008-2009学年期末试题(含答案)
数学实验8月13日作业
《数学实验》报告matlab-第五次作业
数学实验期末考试试题2
数学实验第七次作业
2013-2014期末数学实验考试题
- 最新实验数据与处理大作业题目及答案
- 数学实验大作业---数学曲面..
- 大作业任务实验报告
- 数学实验大作业
- c++大作业学生实验报告
- 哈工大 数学实验 大作业
- MATLAB实验练习题计算机南邮MATLAB数学实验大作业答案
- 《实验设计与数据处理》大作业xx
- c++大作业学生实验报告
- 数学实验大作业
- matlab与数学实验大作业
- MB实验练习题计算机南邮MAAB数学实验大作业答案
- MATLAB大作业实验设计报告
- 中科大模式识别大作业实验
- MATLAB实验练习题(计算机)-南邮-MATLAB-数学实验大作业答案
- MATLAB实验练习题(计算机)-南邮-MATLAB-数学实验大作业答案
- MTLB实验练习题计算机南邮MATLAB数学实验大作业答案
- 北京科技大学数学大作业物理组
- 数学实验大作业(贝壳)
- Matlab实验大作业.doc