数学实验 期末作业

发布 2022-07-02 02:42:28 阅读 3144

第一题。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];

endsubplot(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,进行下次循环。

endend

第二题。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];

endxn(end迭代结果与求根结果比较。

p=sym2poly(y);

roots(p)ans =

ans =-0.5771 + 1.9998i

-0.5771 - 1.9998i

第三题。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))]

endplot(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 =

ans =2.5190 - 0.0000i

-8.3507 - 0.0000i

-2.9158 + 0.0000i

第四题。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)

第五题。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));

endend

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);

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));

endend

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

enddisp(e)

clear,clc,clf

**eage_time

**eage_timebuy

**eage_wait_line

**eage_wait_time

for j=1:1000

buyers

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

things

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

things=randsrc(1,200,x);

total=sum(things);

buything

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

buything_all=su=0; %

for i=1:200

su=0;for n=1:things(i)

su=su+normrnd(1.2,0.3);

endbuything_all=[buything_all,su];

endbuyersend=0

buyerswait=0;buyersmax=0; %

time=0

timebuy=0

数学实验期末作业

第一题 可口可乐 雪碧 健力宝等销量极大的饮料罐 易拉罐 顶盖的直径和从顶盖到底部的高之比为多少?为什么?它们的形状为什么是这样的?制造一只易拉罐不仅要求在保证满足各种物理条件,还要易于制造 美观大方易于使用 从而能吸引顾客 便于运输等,而且要求在这些前提下成本要最低。全世界每年要制造几十亿只易拉罐...

2019工程数学实验C期末大作业

成绩。轻工学院。2014工程数学实验c报告。学院。姓名。班级。学号。email 展示图形之美篇。要求 要求用中文宋体五号字输入文字,mathematica程序默认格式即可。用word自带公式编辑器输入所有数学公式,在运行结果后输入结果分析,可以是题目的结论,也可以是对题目的理解与认识。数学实验一 题...

数学实验作业

学院 理学院。班级 统计11 1 姓名 吴。学号 201111051026 实验1一 问题的提出。已知方程组,其中,定义为。通过迭代法求解方程组。1 选取不同的初始向量,和不同的方程组右端向量,给定迭代误差要求,用雅克比和高斯赛德尔迭代法计算,观察得到的迭代向量序列是否收敛?2 取定右端向量和初始向...