版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、数学实验报告实验人:范晓玉(201004058)院 系: 管理与经济学部 专 业: 经济系 班 级: 管经1009班 实验目的1、结合实际问题展现matlab在生活与学习的广泛应用;2、学会利用matlab编程并求解实际问题;3、培养自学及用学过的只是独立解决问题能力。实验地点大连理工大学大黑楼c2副楼实验语言matlab c语言实验问题一:用梯形法和抛物线法求定积分用梯形法和抛物线法求(误差不超过)实验过程:(1)梯形法梯形法的基本思路n=4,x=0,0.25,0.5,0.75; x=0.25i=0.5*(1/4)*f1(0)+2*f1(0.25)+2*f1(0.5)+2*f1(0.75)+
2、f1(1);n=8;x=0,1/8,2/8,3/8,4/8,5/8,6/8,7/8,8/8; x=1/8i=0.5*(1/8)*f1(0)+2*f1(1/8)+2*f1(2/8)+2*f1(3/8)+2*f1(4/8)+2*f1(5/8)+2*f1(6/8)+2*f1(7/8)+f1(1);再通过c语言编程进行循环计算得到近似计算结果。源程序:#include #include #include #include #define n 100#define e 2.7 double fun(double x) return pow(e,-x*x); void main() double low,
3、up,s,h; int i; low= 0.0; up= 1.0; s=(fun(low)+fun(up)/2.0; h=(up-low)/n; for(i=1;in;i+) s+=fun(low+i*h); printf(the result of integral is: %f.n,s*h); getch();运行结果:the result of integral is 0.748099(2)抛物线法实验语言:matlab抛物线法基本思路:n=4;x=0, 1/4, 1/2, 3/4, 1 ; x=1/4i=(1/12)*f(0)+4f(1/4)+2f(1/2)+4f(3/4)+f(1);
4、 0.8671137 n=8;x=0,1/8,1/4,3/8,1/2,5/8,3/4,7/8,1;x=1/8i=(1/24)*f(0)+4f(1/8)+2f(1/4)+4f(3/8)+2f(1/2)+4f(5/8)+2f(3/4)+4f(7/8)+f(1);源程序:format longring n=100;a=0;b=1;inum=0;e=2.7;syms x fxfx=e(-x*x);for i=1:nxj=a+(i-1)*(b-a)/n; %左点xi=a+i*(b-a)/n; %右点xk=(xi+xj)/2; %中点fxj=subs(fx,x,xj); fxi=subs(fx,x,xi)
5、; fxk=subs(fx,x,xk); inum=inum+(fxj+4*fxk+fxi)*(b-a)/(6*n);endinumintegrate=int(fx,0,1)integrate=double(integrate)fprintf(the relative error between inum and real-value is about: %dnn,.abs(inum-integrate)/integrate)运行结果:inum =0.74810502168631实验问题二:求微分方程数值解并作图在【1,4】上求数值解,并作图。实验过程:y=dsolve(dy=x*y*sin(
6、x)-1/x,y(1)=1)function numericclear;clc; x,y=ode45(fun,1 4,1)plot(x,y,*,x,exp(-x.2)%*为数值解,线段为解析解function yhat=fun(x,y)yhat=x*y*sin(x)-1/x实验问题三:如图a bd(m)c d 10m在cd上找一点p,使p的光照强度最小的点。(已知单光源对某点光照强度与光源强度成正比与光源到该点距离平方成反比。)实验过程:#include #include #include float fx(float x) return (exp(-x*x);float tixing(flo
7、at(*fun)(float x),float a,float b,int n) float s,h,y; int i; s=(fun(a)+fun(b)/2; h=(b-a)/n; for(i=1;in;i+) s=s+fun(a+i*h); y=s*h; return y;float juxing(float(*fun)(float x),float a,float b,int n) double h; double fx; double x; int i; double sum=0; h=fabs(a-b)/100; x=a; for(i=0;i=eps)p=0.0; for (k=0;
8、k=n-1;k+) x=a+(k+0.5)*h; p=p+fun(x); t2=(t1+h*p)/2.0; s2=(4.0*t2-t1)/3.0; ep=fabs(s2-s1); t1=t2; s1=s2; n=n+n; h=h/2.0; return(s2);void menu() system(cls); printf(1.梯形法n2.矩形法n3.辛普森法n4.退出n);int main() int c=0; menu(); while(1) scanf(%d,&c); switch(c) case 1:printf(梯形法:%fn,tixing(fx,0,1,1000); break; c
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2024届山东省淄博市临淄区第一中学中考英语适应性模拟试题含答案
- 《柿柿如意》幼儿园小学少儿美术教育绘画课件创意教程教案
- 2024年机器人及具有独立功能专用机械行业营销策略方案
- 2024年泌尿肛肠外科手术器械相关行业营销方案
- 一级造价工程师《造价管理》秭归县2024年高分冲刺试卷含解析
- 中国箱式运输半挂车行业市场现状分析及竞争格局与投资发展研究报告2024-2029版
- 中国书报架行业市场现状分析及竞争格局与投资发展研究报告2024-2029版
- 中国轻型飞机行业调研分析及发展趋势预测研究报告2024-2029版
- 一年级上册美术《看谁涂得更好看》说课稿和教案
- 中国唇彩行业供需趋势及投资风险研究报告2024-2029版
- 《无粘结钢绞线体外预应力束》(JT∕T 853-2013)
- 医院常见不安全因素及防范课件
- 研发部门薪酬福利制度样本
- 高压电工复审培训课件
- 办公软件常用操作培训课件
- 2024年《生物化学》教案
- 中药药剂学与制药工艺
- 标准机房施工工艺要求
- 中医中药防治冠心病培训课件
- 关于欧亚牛奶的分析报告
- 现在进行时(全国一等奖)
评论
0/150
提交评论