数学实验报告.docVIP

  1. 1、本文档共15页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 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)

ww90055 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档