数学实验-期末作业

第一题 clc,clear

p0=0.5;pn=[p0];n=15;x=1:n;k=[1.5,2.5,3.2,3.8];

for j=1:4; %建立循环 for a=k(j)

for i=1:n-1;

p=k(j)*pn(i)*(1-pn(i)); %分别计算pn pn=[pn;p]; end

subplot(2,2,j) %制图排版 scatter(x,pn,10,'filled');title(['k=' num2str(k(j))]) %分别定义标题 axis([0,15,0,1]) %定义图像区间 set(gca,'ytick',0.1:0.1:1) %建立合适的网格线 grid on;

pn=[p0]; %重新定义pn,进行下次循环 end end

510

15

0.10.20.30.40.50.60.70.80.9

1k=1.5

510

15

0.10.20.30.40.50.60.70.80.9

1k=2.5

5

10

15

0.10.20.30.40.50.60.70.80.9

1k=3.2

5

10

15

0.10.20.30.40.50.60.70.80.9

1k=3.8

第二题

clc,clear

syms x

y=x^3+3*x-5;

y1=diff(y); %求导函数

xn=[1];n=100;

for i=1:n;

s=xn(i)-subs(y,xn(i))/subs(y1,xn(i)); %利用牛顿切线法迭代xn=[xn;s];

end

xn(end) %迭代结果与求根结果比较

P=sym2poly(y);

roots(P)

ans =

1.1542

ans =

-0.5771 + 1.9998i

-0.5771 - 1.9998i

1.1542

第三题

clc,clear

x=0:0.001:18.2;alpha=pi/9;theta=[];

for i=1:length(x);%计算每个x对应的theta值

a=((2.74+x(i)*cos(alpha))^2+(9.45-x(i)*sin(alpha))^2)^(1/2);

b=((2.74+x(i)*cos(alpha))^2+(x(i)*sin(alpha)-1.83)^2)^(1/2);

theta=[theta;acos((a^2+b^2-7.62^2)/(2*a*b))];

end

plot(x,theta),grid on

set(gca,'xtick',0:0.91:18.2)

axis([0,20,0.35,1])

text(2.2,0.86,'第四排','fontsize',12,'color','r') %标记最大值与最小值的座位排数text(17.8,0.4,'第二十一排','fontsize',12,'color','r')

[x1 x2]=max(theta); %求theta取最大值时的x值

x(x2)

syms x

a=((2.74+x*cos(alpha))^2+(9.45-x*sin(alpha))^2)^(1/2);

b=((2.74+x*cos(alpha))^2+(x*sin(alpha)-1.83)^2)^(1/2);

y=acos((a^2+b^2-7.62^2)/(2*a*b));

y1=diff(y);

a=solve(y1); %利用求根命令找零点

double(a)

ans =

2.5190

ans =

2.5190 - 0.0000i

-8.3507 - 0.0000i

-2.9158 + 0.0000i

00.91 1.82 2.73 3.64 4.55 5.46 6.377.288.199.110.0110.9211.8312.7413.6514.5615.4716.3817.2918.2

第四题

function dx=zhuitu(t,x) %定义微分方程 u=3;v=5;

s=sqrt((70-x(1))^2+(75-u*t-x(2))^2); %以猎狗为原点,兔子在右列出追击方程 dx=[v*(70-x(1))/s;v*(75-u*t-x(2))/s];

clc,clear

ts=0:0.01:20;x0=[0 0];

[t,x]=ode45('zhuitu',ts,x0); %用obe45求解微分方程 plot(x(:,1),x(:,2)),grid on,hold on scatter(70,75,'b','filled'); scatter(70,15,'r','filled'); legend('y(x)',2)

01020304050607080

1020304050607080

第五题

clc,clear

n=200;

a=normrnd(2,0.5); %计算第一个人来的时间

b=a;

for x1=randperm(10); %随机产生顾客买东西的个数并计算总花费时间

x=x1(1);

if x>5

s=normrnd(1.2,0.3);

elseif x>3&&x<6

s=sum(normrnd(1.2,0.3,1,2));

elseif x>1&&x<4

s=sum(normrnd(1.2,0.3,1,3));

else

s=sum(normrnd(1.2,0.3,1,4));

end

end

e=b+s; %第一个顾客服务完的时间

for i=1:n-1;%产生循环，从第二个人到第两百个人

a=normrnd(2,0.5);

b=b+a; %第i+1个顾客到达的时间

c=max(b,e); %从第i+1个顾客到达的时间和服务完第i个顾客的时间取较大值for x1=randperm(10);

x=x1(1);

ifx>5

s=normrnd(1.2,0.3);

elseif x>3&&x<6

s=sum(normrnd(1.2,0.3,1,2));

elseif x>1&&x<4

s=sum(normrnd(1.2,0.3,1,3));

else

s=sum(normrnd(1.2,0.3,1,4));

end

end

e=c+s; %第i+1个顾客服务完的时间

end

disp(e)

clear,clc,clf

aveage_time=[]; %??′??￡?a×ü·t??ê±?????ù?μ

aveage_timebuy=[]; %??′??￡?aè?2?1??í1o??ê±??μ????ù?μ

aveage_wait_line=[]; %??′??￡?a×?′ó?ó3¤μ????ù?μ

aveage_wait_time=[]; %??′??￡?aμè′yê±??μ????ù?μ

for j=1:1000

buyers =[]; %1??íμ?′?ê±??????

buyers=normrnd(2,0.5,1,200);

things=[]; %1??í1o?òéì?·êyá?

x=[1,2,3,4:0.5,0.2,0.2,0.1];

things=randsrc(1,200,x);

total=sum(things);

buything=[]; %1??í1o?ò??ò??t???·?ù·?ê±??

buything=normrnd(1.2,0.3,1,total);

buything_all=[];su=0; %??ò???1??í1o??×üê±??

for i=1:200

su=0;

for n=1:things(i)

su=su+normrnd(1.2,0.3);

end

buything_all=[buything_all,su];

end

buyersend=0; %1o??íêμ?1??íêyá?

buyerswait=0;buyersmax=0; %???óμ?1??íêyá?

time=0; %×üê±??

timebuy=0; %μ±?°1o??1??í?ùê￡ê±??

time_next=0;

buyers_waittime=[];

timebuy=buything_all(1);

time_next=buyers(2);

buyerswaittime(1)=0;

while (buyersend<200)

if timebuy

aveage_wait_time=[average_wait_time,mean( buyers_waittime-buything_al l)];

end

end

end

disp(mean(average_time))

数学软件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)

重庆大学数学模型数学实验作业四讲解

开课学院、实验室：数统学院实验时间：2015年10月28日 课程名称数学实验实验项目 名称 种群数量的状态转移—— 微分方程 实验项目类型 验证演示综合设计其他 指导 教师 肖剑成绩 实验目的 [1] 归纳和学习求解常微分方程(组)的基本原理和方法； [2] 掌握解析、数值解法，并学会用图形观察解的形态和进行解的定性分析； [3] 熟悉MATLAB软件关于微分方程求解的各种命令； [4] 通过范例学习建立微分方程方面的数学模型以及求解全过程； 通过该实验的学习，使学生掌握微分方程(组)求解方法（解析法、欧拉法、梯度法、改进欧拉法等），对常微分方程的数值解法有一个初步了解，同时学会使用MATLAB软件求解微分方程的基本命令，学会建 立微分方程方面的数学模型。这对于学生深入理解微分、积分的数学概念，掌握数学的分析思维方法，熟 悉处理大量的工程计算问题的方法是十分必要的。 实验内容 1．微分方程及方程组的解析求解法； 2．微分方程及方程组的数值求解法——欧拉、欧拉改进算法； 3．直接使用MATLAB命令对微分方程(组)进行求解(包括解析解、数值解)； 4．利用图形对解的特征作定性分析； 5．建立微分方程方面的数学模型，并了解建立数学模型的全过程。 基础实验 一、问题重述 1．求微分方程的解析解, 并画出它们的图形， y’= y + 2x, y(0) = 1, 0

大学数学实验

大学数学实验 项目一 矩阵运算与方程组求解 实验1 行列式与矩阵 实验目的 掌握矩阵的输入方法. 掌握利用Mathematica (4.0以上版本) 对矩阵进行转置、加、减、数乘、相乘、乘方等运算, 并能求矩阵的逆矩阵和计算方阵的行列式. 基本命令 在Mathematica 中, 向量和矩阵是以表的形式给出的. 1. 表在形式上是用花括号括起来的若干表达式, 表达式之间用逗号隔开. 如输入 {2,4,8,16} {x,x+1,y,Sqrt[2]} 则输入了两个向量. 2. 表的生成函数 (1) 最简单的数值表生成函数Range, 其命令格式如下: Range[正整数n]—生成表{1,2,3,4,…,n }; Range[m, n]—生成表{m ,…,n }; Range[m, n, dx]—生成表{m ,…,n }, 步长为d x . (2) 通用表的生成函数Table. 例如,输入命令 Table[n^3,{n,1,20,2}] 则输出 {1,27,125,343,729,1331,2197,3375,4913,6859} 输入 Table[x*y,{x,3},{y,3}] 则输出 {{1,2,3},{2,4,6},{3,6,9}} 3. 表作为向量和矩阵 一层表在线性代数中表示向量, 二层表表示矩阵. 例如,矩阵 ??? ? ??5432 可以用数表{{2,3},{4,5}}表示. 输入 A={{2,3},{4,5}} 则输出 {{2,3},{4,5}} 命令MatrixForm[A]把矩阵A 显示成通常的矩阵形式. 例如, 输入命令: MatrixForm[A] 则输出 ??? ? ??5432 但要注意, 一般地, MatrixForm[A]代表的矩阵A 不能参与运算. 输入 B={1,3,5,7} 输出为 {1,3,5,7} 输入 MatrixForm[B] 输出为

数学实验期末考试

数学实验 期末上机考核 学号201519030102 姓名曹欣辉年级专业2015级水产养殖学 学号201519030103 姓名陈妙珊年级专业2015级水产养殖学 学号201519030104 姓名杜日臻年级专业2015级水产养殖学 学号姓名年级专业 学号姓名年级专业

注意事项： 1、考核方式：以组（3~4人）为单位，请于指定时间内开卷完成所布置的任务，地点为实验 室机房或课室。 2、发题时间为6月25号早上8:30，请到如下邮箱提取题目，账号：nongkeshuxue@https://www.wendangku.net/doc/674085226.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])

清华大学数学实验报告4

清华大学数学实验报告4

————————————————————————————————作者: ————————————————————————————————日期： ?

电13 苗键强2011010645

一、实验目的 1.掌握用 MATLAB 软件求解非线性方程和方程组的基本用法, 并对结果作初步分析； 2.练习用非线性方程和方程组建立实际问题的模型并进行求解。 二、实验内容 题目1 【问题描述】 (Q1）小张夫妇以按揭方式贷款买了1套价值２0万元的房子，首付了5万元，每月还款1000元，15年还清。问贷款利率是多少？ (Q2)某人欲贷款50 万元购房,他咨询了两家银行,第一家银行 开出的条件是每月还45０0元,15 年还清；第二家银行开出的条件是每年还４5０00 元,2０年还清。从利率方面看,哪家银行较优惠(简单假设:年利率=月利率×12）? 【分析与解】 假设初始贷款金额为x０,贷款利率为ｐ,每月还款金额为x,第i 个月还完当月贷款后所欠银行的金额为x i,(i=1,2，3,.．....，n）。由题意可知： x1=x0(1+p)?x x2=x0(1+p)2?x(1+p)?x x3=x0(1+p)3?x(1+p)2?x(1+p)?x ……

x n=x0(1+p)n?x(1+p)n?1???x(1+p)?x =x0(1+p)n?x (1+p)n?1 p =0 因而有： x0(1+p)n=x (1+p)n?1 p (1) 则可以根据上述方程描述的函数关系求解相应的变量。 （Ｑ1) 根据公式(１)，可以得到以下方程： 150p(1+p)180?(1+p)180+1=0 设 f(p)=150p(1+p)180?(1+p)180+1，通过计算机程序绘制ｆ（p）的图像以判断解ｐ的大致区间,在Ｍatlab中编程如下： ｆorｉ = 1:25 t = ０．0001＊i; p(i) = t； f(i) =１50*ｔ＊(１+t).^180-(１+t).^180+1; ｅｎd; ｐｌｏｔ(ｐ,ｆ),hｏld ｏn,grid on; 运行以上代码得到如下图像：

清华大学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]

下学期数学实验作业

实验一 图形的画法 1. 做出下列函数的图像： （1）)2sin()(2 2--=x x x x y ，22≤≤-x （分别用plot 、fplot ） （2）2 2 /9/251x y +=（用参数方程） (3) 在同一图形窗口中，画出四幅不同图形（用subplot 命令）： 1cos()y x =，2sin(/2)y x pi =-，23cos()y x x pi =-，sin()4x y e =（]2,0[π∈x ） 2 作出极坐标方程为)cos 1(2t r -=的曲线的图形. 3 作出极坐标方程为10/t e r =的对数螺线的图形. 4 绘制螺旋线?? ? ??===t z t y t x ,sin 4,cos 4在区间［0，π4］上的图形.在上实验中，显示坐标轴名称. 5 作出函数2 2 y x xye z ---=的图5形. 6 作出椭球面11 942 22=++z y x 的图形. (该曲面的参数方程为 ,cos ,sin sin 3,cos sin 2u z v u y v u x === (π π20,0≤≤≤≤v u ).) 7 作双叶双曲面13 .14.15.122 2222-=-+z y x 的图形. (曲面的参数方程是 ,csc 3.1,sin cot 4.1,cos cot 5.1u z v u y v u x === 其中参数πππ<<-≤

北理工数学实验作业

一． 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的图形特点。

大学数学实验心得体会

大学数学实验心得体会 [模版仅供参考，切勿通篇使用] 大学数学实验心得体会（一） 数学，在整个人类生命进程中至关重要，从小学到中学，再到大学，乃至更高层次的科学研究都离不开数学，随着时代的发展，人们越来越重视数学知识的应用，对数学课程提出了更高层次的要求，于是便诞生了数学实验。 学期最初，大学数学实验对于我们来说既熟悉又陌生，在我们的记忆中，我们做过物理实验、化学实验、生物实验，故然我们以为数学实验与它们一样，当我们在网上搜索有关数学实验的信息时，我们才知道，大学数学实验作为一门新兴的数学课程在近十年来取得了迅速的发展。数学实验以计算机技术和数学软件为载体，将数学建模的思想和方法融入其中，现在已经成为一种潮流。 当我们怀着好奇的心情走进屈静国老师的数学实验课堂时，我们才渐渐懂得，数学实验是一门有关计算机软件的课程，就像c语言一样，需要编辑运行程序，从而进行数学运算，它不需要自己来运算，就像计算器一样，只要我们自己记下重要程序语句，输入运行程序，便可得到运行结果，大大降低了我们的运算量，

给我们生活带来许多便捷，在大一时，我学过c语言，由于这样的基础，让我能够更快的学会并应用此软件。 时间飞逝，转眼间，我们就要结课了，这学期我们学习了mathematics的基础，微积分实验，线性代数实验，概率论与数理统计实验，数值计算方法及实验。通过这学期的学习，我也积累了些自己的学习方法和心得。首先，我们要在平时上课牢记那些mathematics语言和公式，那些东西就想单词和公式一样，只需要背诵;然后，我们要看几遍书，并多看一下例题;最后，我们要多应用mathematics软件去练习。正所谓熟能生巧，我坚信，只要我们能够做到这三步，我们就能很好的掌握这门课程。 通过学习使用数学软件，数学实验建模，使我们能够从实际问题出发，认真分析研究，建立简单数学模型，然后借助先进的计算机技术，最终找出解决实际问题的一种或多种方案，从而提高了我们的数学思维能力，为我们参加数学竞赛和数学建模打下了坚实的基础，同时也为我们进一步深造和参加工作打下一定的实践基础! 大学数学实验心得体会（二） 在此期间我充分利用研修活动时间学习，感到既有辛苦，又有收获。既有付出，又有新所得。这次远程研修让我有幸与专家和各地的数学精英们交流，面对每次探讨的主题，大家畅所欲言，

数学实验期末考试上机考试

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分）

山东建筑大学数学实验期末作业matlab

数学实验 期 末 作 业 学号： 班级： 姓名：

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”。

数学实验作业题目(赛车跑道)

数学实验报告实验题目：赛车车道路况分析问题 小组成员： 填写日期 2012 年 4 月 20 日

一．问题概述 赛车道路况分析问题 现要举行一场山地自行车赛，为了了解环行赛道的路况，现对一选手比赛情况进行监测，该选手从A地出发向东到B，再经C、D回到A地（如下图）。现从选手出发开始计时，每隔15min观测其位置，所得相应各点坐标如下表（假设其体力是均衡分配的）： 由D→C→B各点的位置坐标（单位：km） 假设：1. 车道几乎是在平原上，但有三种路况（根据平均速度(km/h)大致区分）： 平整沙土路（v>30）、坑洼碎石路（10

2．估计车道的长度和所围区域的面积； 3．分析车道上相关路段的路面状况（用不同颜色或不同线型标记出来）； 4．对参加比赛选手提出合理建议. 二．问题分析 1.模拟比赛车道的曲线：因为赛道散点分布不规则，我们需要用光滑曲线来近似 模拟赛道。由于数据点较多，为了避免龙格现象，应采用三次样条插值法来对曲线进行模拟（spline命令）。全程曲线为环路，我们需要对上下两部分分别 模拟，设模拟出的曲线为P：。 2.把A到B点的曲线分成若干小段： 赛道的路程L：取dL=，对模拟出的整条曲线求线积分，即 所围区域的面积：用上下部分曲线的差值对求定积分，即 3.用样条插值法模拟出比赛车道曲线后，根据曲线分别计算出原数据中每两点 （）间的路程，即求线积分 由于每两点间时间间隔相同且已知（15min），故可求出每段路程的平均速度 易知即为的积分中值 将此速度近似作为两点间中点时刻的速度，然后再次采用样条插值法，模拟出全过程的图像。而根据求出的与之间的关系，再次采用样条插值法，即可模拟出全过程的图像 4. 由赛道曲线可求出赛道上任一点到点的路程

数学实验 作业10

实验十三回归分析 电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)

数学实验8月13日作业

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 π θπ ?>??=??+

东南大学高等数学数学实验报告上

高等数学数学实验报告实验人员：院（系） ___________学号_________姓名____________ 实验地点：计算机中心机房 实验一 一、实验题目： 根据上面的题目，通过作图，观察重要极限：lim（1+1/n）n=e 二、实验目的和意义 方法的理论意义和实用价值。 利用数形结合的方法观察数列的极限，可以从点图上看出数列的收敛性，以及近似地观察出数列的收敛值；通过编程可以输出数列的任意多项值，以此来得到数列的收敛性。通过此实验对数列极限概念的理解形象化、具体化。 三、计算公式（1+1/n）n 四、程序设计 五、程序运行结果 六、结果的讨论和分析 当n足够大时，所画出的点逐渐接近于直线，即点数越大，精确度越高。对于不同解题方法最后均能获得相同结果，因此需要择优，从众多方法中尽可能选择简单的一种。程序编写需要有扎实的理论基础，因此在上机调试前要仔细审查细节，对程序进行尽可能的简化、改进与完善。 实验二 一、实验题目 制作函数y=sin cx的图形动画，并观察参数c对函数图形的影响。 二、实验目的和意义 本实验的目的是让同学熟悉数学软件Mathematica所具有的良好的作图功能，并通过函数图形来认识函数，运用函数的图形来观察和分析函数的有关性态，建立数形结合的思想。 三、计算公式：y=sin cx 四、程序设计 五、程序运行结果

六、结果的讨论和分析 c 的不同导致函数的区间大小不同。 实验三 一、实验题目 观察函数f(x)=cos x 的各阶泰勒展开式的图形。 二、实验目的和意义 利用Mathematica 计算函数)(x f 的各阶泰勒多项式，并通过绘制曲线图形，来进一步掌握泰勒展开与函数逼近的思想。 三、计算公式 四、程序设计 五、程序运行结果 六、结果的讨论和分析 函数的泰勒多项式对于函数的近似程度随着阶数的提高而提高，但是对于任一确定次数的多项式，它只在展开点附近的一个局部范围内才有较好的近似精确度。 实验四 一、实验题目 计算定积分的黎曼和 二、实验目的和意义 在现实生活中许多实际问题遇到的定积分，被积函数往往不能用算是给出，而通过图像或表格给出；或虽然给出，但是要计算他的原函数却很困难，甚至原函数非初等函数。本实验目的，就是为了解决这些问题，进行定积分近似计算。 三、计算公式 四、程序设计 五、程序运行结果 六、结果的讨论和分析 本实验求的近似值由给出的n 的值的不同而不同。给出的n 值越大，得到的结果越接近准确的

大学数学实验—期末考试试题6

数学实验试题 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米/平方秒。）

数学实验作业 韩明版

练习6.7 1.有两个煤厂A,B，每月进煤不少于60t,100t，它们担负供应三个居 民区的用煤任务，这三个居民区每月用煤量分别为45t,75t和45t.A 厂离这三个居民区的距离分别为10km,5km,6km,B厂离这三个居民区的距离分别为4km,8km,15km.问这两个煤厂如何分配供煤量能使总运输量（t.km）最小。 解：设甲对三个居民区的供煤量分别为：x1,x2,x3,乙对三个居民区的供煤量分别为x4,x5,x6.由已知有： y=10x1+5x2+6x3+4x4+8x5+15x6 -x1-x2-x3<=-60, -x4-x5-x6<=-100, x1+x4=45,x2+x5=75,x3+x6=40, X1>=0,x2>=0,x3>=0,x4>=0,x5>=0,x6>=0. 输入命令： > c=[10 5 6 4 8 15];A=[-1 -1 -1 0 0 0;0 0 0 -1 -1 -1;0 0 0 0 0 0;0 0 0 0 0 0;0 0 0 0 0 0;0 0 0 0 0 0]; >> b=[-60;-100;0;0;0;0];Aeq=[1 0 0 1 0 0;0 1 0 0 1 0;0 0 1 0 0 1;0 0 0 0 0 0;0 0 0 0 0 0;0 0 0 0 0 0]; >> beq=[45 75 40 0 0 0]; >> lb=ones(6,1); >> [x,fval]=linprog(c,A,b,Aeq,beq,lb) Optimization terminated.

结果为： x = 1.0000 20.0000 39.0000 44.0000 55.0000 1.0000 fval =975.0000 这说明甲乙两个煤厂分别对三个居民区输送1t 20t 39t,44t 55t 1t的煤才能使总运输量最小，且总运输量为975t.km 2．某银行经理计划用一笔资金进行有价证券的投资，可供购进的证券及其信用等级、到期年限、税前收益如下表所示。按照规定，市政证券的收益可以免税，其他证券的收益需按40%的税率纳税。此外还有以下限制： （1）政府及待办机构的证券总共至少购进400万元； （2）所构证券的平均信用等级不超过1.4（信用等级数字越小，信用程度越高）； （3）所构证券的平均到期年限不超过5年。

电子科技大学《数学实验》2008-2009学年期末试题(含答案)

电子科技大学二零零八到二零零九学年第二学期期末考试《数学实验》课程考试题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(yx.^2);的功能是( ) (A) 统计2000个随机点中落入特殊区域的点的索引值； (B) 统计1000个随机点落入特殊区域的点的索引值； (C) 模拟2000个随机点落入特殊区域的过程； (D) 模拟1000个随机点落入特殊区域的过程。 5、MATLAB计算二项分布随机变量分布律的方法是( ) (A) binocdf(x,n,p); (B) normpdf(x,mu,s); (C)binopdf(x,n,p); (D) binornd(x,n,p)。 6、MATLAB命令syms e2;f=sqrt(1-e2*cos(t)^2);S=int(f,t,0,pi/2)功能是() (A) 计算f(x)在[0,pi/2]上的积分；(B) 计算f(t)不定积分符号结果； (C) 计算f(x)积分的数值结果；(D) 计算f(t)定积分的符号结果。 7、y=dsolve(‘Dy=1/(1+x^2)-2*y^2’,’y(0)=0’,’x’);ezplot(y)的功能是( ) (A) 求微分方程特解并绘图；(B) 解代数方程(C) 求定积分；(D)求微分方程通解。 8、X=10000 ;0.5*asin(9.8*X/(515^2))的功能是计算关于抛射体问题的() (A) 十公里发射角；(B) 十公里飞行时间；(C)最大飞行时间；(D)最大射程。 9、theta=linspace(0,2*pi,100) ;r=cos(4*theta) ;polar(theta,r,’k’)功能是() (A) 绘四叶玫瑰线；(B)绘三叶玫瑰线；(C)绘心脏线；(D) 绘八叶玫瑰线。 10、北京和纽约的经度分别是：东经118和西经76，根据经度差计算时差用() (A) fai1=118;fai2=-76;Dfai=(fai1+fai2)/24; (B) fai1=118;fai2=-76;Dfai=(fai1+fai2)/15; (C) fai1=118;fai2=-76;Dfai=(fai1-fai2)/24; (D) fai1=118;fai2=-76;Dfai=(fai1-fai2)/15。 二、程序阅读题(40分) 1、直方图功能是将数据分为n个类，统计各个类的数据量并绘图。借用现有的直方图命令hist，编写新直方图程序如下。 function m=myhist(data,n) if nargin==1,n=7;end Xmin=min(data);Xmax=max(data);h=(Xmax-Xmin)/n; m=hist(data,n)/length(data)/h;