搜索
Matlab上机实验报告姓名：王梅学号：11116010224上机时间：2013-5-141用符号工具箱计算：∫这个积分能用数值积分计算吗？如果能，应该如何做？方法：符号计算和数值计算。程序流程图：新建一个指令窗口编辑指令将指令复制到操作窗口得出结果程序清单及运行结果：symsxa=int((x^4*(1-x)^4)/(1+x^2),0,1)a=22/7-pif=inline('x.^4.*(1-x).^4./(1+x.*x)','x')f=Inlinefunction:f(x)=x.^4.*(1-x).^4./(1+x.*x)quad(f,0,1)ans=0.00132.误差函数erf(x)的定义是：Erf(x)=2/(pi^0.5)∫()用指令quad计算当x=0.1,0.2^1.0时相应的erf函数值。经结果与MATLAB的内部函数erf(x)的结果进行比较。方法：数值计算。程序流程图：新建一个指令窗口编辑指令将指令复制到操作窗口得出结果程序清单及运行结果：x=[0.1:0.1:1.0]f=inline('2./pi.^0.5.*exp(-x.*x)','x')erf(x)forx=0.1:0.1:1.0quad(f,0,x)endx=0.10000.20000.30000.40000.50000.60000.70000.80000.90001.0000f=Inlinefunction:f(x)=2./pi.^0.5.*exp(-x.*x)ans=0.11250.22270.32860.42840.52050.60390.67780.74210.79690.8427ans=0.1125ans=0.2227ans=0.3286ans=0.4284ans=0.5205ans=0.6039ans=0.6778ans=0.7421ans=0.7969ans=0.84273.求微分方程xy’+y-e^x=0在初值条件y1=2e下的特接，并画出解函数的图形（解析解与数值解）程序流程图：新建一个指令窗口编辑指令将指令复制到操作窗口得出结果程序按及运行结果：y=dsolve('x*Dy+y-exp(x)=0','y(1)=2*exp(1)','x')y=(exp(1)+exp(x))/xezplot(y)-6-4-20246-30-20-1001020304050x(exp(1)+exp(x))/x