《Matlab上机实验报告(共4页).docx》由会员分享,可在线阅读,更多相关《Matlab上机实验报告(共4页).docx(4页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精选优质文档-倾情为你奉上Matlab上机实验报告姓名:王梅 学号: 上机时间:2013-5-141用符号工具箱计算:01(x4*(1-x)4)/(1+x2)dx这个积分能用数值积分计算吗?如果能,应该如何做?方法:符号计算和数值计算。 程序流程图:新建一个指令窗口编辑指令将指令复制到操作窗口得出结果程序清单及运行结果: syms xa=int(x4*(1-x)4)/(1+x2),0,1) a = 22/7 - pi f=inline(x.4.*(1-x).4./(1+x.*x),x)f = Inline function: f(x) = x.4.*(1-x).4./(1+x.*x) quad
2、(f,0,1)ans =0.00132.误差函数erf(x)的定义是: Erf(x)=2/(pi0.5)0xe-x2dx用指令quad计算当x=0.1,0.21.0时相应的erf函数值。经结果与MATLAB的内部函数erf(x)的结果进行比较。方法:数值计算。程序流程图:新建一个指令窗口编辑指令将指令复制到操作窗口得出结果程序清单及运行结果: x=0.1:0.1:1.0f=inline(2./pi.0.5.*exp(-x.*x),x)erf(x)for x=0.1:0.1:1.0quad(f,0,x)endx = 0.1000 0.2000 0.3000 0.4000 0.5000 0.600
3、0 0.7000 0.8000 0.9000 1.0000f = Inline function: f(x) = 2./pi.0.5.*exp(-x.*x)ans = 0.1125 0.2227 0.3286 0.4284 0.5205 0.6039 0.6778 0.7421 0.7969 0.8427ans = 0.1125ans = 0.2227ans = 0.3286ans = 0.4284ans = 0.5205ans = 0.6039ans = 0.6778ans = 0.7421ans = 0.7969ans =0.84273.求微分方程xy+y-ex=0在初值条件y1=2e下的特接,并画出解函数的图形(解析解与数值解)程序流程图:新建一个指令窗口编辑指令将指令复制到操作窗口得出结果程序按及运行结果:y=dsolve(x*Dy+y-exp(x)=0,y(1)=2*exp(1),x)y =(exp(1) + exp(x)/xezplot(y)专心-专注-专业