第一题。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 取定右端向量和初始向...