- 1、本文档共15页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数学实验报告
1 阻滞增长logistic模型
代码
x=0:0.1:1
n=length(x);
xm=1;
r=1;
x(1)=0;
for k=1:n
y(k)=r*x(k)*(1-x(k)/xm)
end
plot(x,y,r)
2.实验结果
clear;
t=0:0.01:3
n=length(t);
xm=1;
x0=0.1;
r=2;
for k=1:n
x(k)=xm/(1+(xm/x0-1)*exp(-r*t(k)))
end
plot(t,x,r)
实验结果
2 拉格朗日
1)M1:
function y=lagr(x,x0,y0) (保存时 自动命名 为 lagr 不能更改)
m=length(x0);
n=length(y0);
if m~=n
return
else
s=0
for l=1:n
t=1
for k=1:n
if l~=k
t=t*(x-x0(k))/(x0(l)-x0(k));
end
end
s=s+y0(l)*t
end
y=s
end
M2:
x0=-5:5
y0=1./(1+x0.^2);
x=-5:0.1:5;
n=length(x);
y=[ ];
for k=1:n;
y1=lagr(x(k),x0,y0);
y=[y y1]
end
y
3 最小二乘法拟合
t=[20.5 32.5 51 73 95.7 ];
r=[765 826 873 942 1032];
aa=polyfit(t,r,1);
a=aa(1)
b=aa(2)
y=polyval(aa,t);
plot(t,r,’k+’,t,y,’r’)
xlabel(‘t’),ylabel(‘R’)
结果
4 机床加工
x0=[0 3 5 7 9 11 12 13 14 15]
y0=[0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6]
x=0:0.1:15;
y1=lagr(x0,y0,x);
y2=interp1(x0,y0,x);
y3=interp1(x0,y0,x,spline);
[x,y1,y2,y3];
subplot(3,1,1)
plot(x0,y0,k+,x,y1,r)
grid,title(lagrange)
subplot(3,1,2)
plot(x0,y0,k+,x,y2,r)
grid,title(piecewise linear)
subplot(3,1,3)
plot(x0,y0,k+,x,y3,r)
grid,title(spline)
5 随机模拟方法 求圆锥的 V
n=10000;
r=1;
h=1;
a=r*rand(2,n)
b=h*rand(1,n);
sum(sum(a.^2)=(1-r/n*b).^2)/n*4
ans =
3.1304
6 梯形公式 辛普森公式
h=pi/200;
x=0:h:pi/2
y=sin(x)
z1=sum(y(1:100))*h
z2=sum(y(2:101))*h
z3=trapz(y)*h
z4=quad(sin(x),0,pi/2)
z1 =
0.9921
z2 =
1.0078
z3 =
1.0000
z4 =
1.0000
18
t=0:pi/200:pi/2
y1=weixing(t)
format long
l1=4*trapz(t,y1)
l2=4*quad(weixing,0,pi/2)
7 卫星轨道
t=0:pi/10:pi/2
y1=weixing(t)
format long
l1=4*trapz(t,y1)
l2=4*quad(weixing,0,pi/2)
l1 =4.870744099902405e+004
l2 =4.870744099903280e+004
8 人口模型
t=[1900 1910 1920 1930 1940 1950 1960 1970 1980 1990 2000]
x=[76 92 106.5 123.2 131.7 150.7 179.3 179.3 204 226.5 251.4 281.4]
r(1)=(-3*x(1)+4*x(2)-x(3))/20/x(1)
for k=2:10
r(k)=(x(k+1)-x(k-1)/20/x(k));
end
r(11)=(x(9)-4*x(10)+3*x(11))/20
文档评论(0)