matlab复合梯形法求积分及结果
来源:学生作业帮助网 编辑:作业帮 时间:2024/09/20 09:13:23
把f=a^2*(1-e^2)*(cos(x))/((1-((e^2)*(sin(x))^2))^2);改为f=a^2*(1-e^2)*(cos(x))./((1-((e^2)*(sin(x)).^2)
1.数值积分(梯形法)>>x=0:0.01:1;>>y=1./(1+x.^4);>>sy=trapz(x,y)sy=0.86702.符号积分,对比>>symsx>>y=1/(1+x^4);>>yi=e
s1='(x+1)^10*exp(-x)'s1=(x+1)^10*exp(-x)>>int(s1,0,1)ans=-26813184*exp(-1)+9864101>>numeric(int(s1,0
使用quad有问题>>quad('0.0732*exp(-(x+0.02287).^2/64.90)',0,inf)Warning:InfiniteorNot-a-Numberfunctionvalu
没有定义积分表达式f,可以用inline()来定义.如f=inline('x.*sin(x)','x')
楼上的回答对了一半.1、在近似估算中,梯形法比矩形法精确.估算的难度大一些.2、在无限分割的极限情况下,两种方法得到的结果是一样的.都是100%准确的.3、无论什么函数,包括sinx,都可用两种估算或
sysmx;int(积分表达式,要对其积分的变量,积分下限,积分上限)例如要计算f(x)=x^2在[0,1]区间的积分sysmx;result=int(x^2,x,0,1);result则为结果
详细的源程序,程序设计原理及流程图.
设梯形上边长、下边长和高分别为a、b、h坐标系:取下底边建x轴,y轴沿高向,则有:先由中性面静矩为0求中性面位置:0=∫((h-y)/h*(b-a)+a)*(y-y0)dy再求中性面的惯性矩:Iyy0
我给您解了一下f=Inf3.30695.986114.968323.875764.4201109.3522231.1783449.1438533.5097602.3386509.2590381.845
最好给出具体的函数表达式.一般的做法差不多可以写成:fminunc(@(yz)quadl(@(x)f(x,yz(1),yz(2)),a,b),yz0)其中,yz0是y和z的初值,a和b是积分上下限.&
symsaf=((2.72.^a-a-1).^(-1/2)-1);g=int(f,a,a,3)结果:g=int(1/((68/25)^a-a-1)^(1/2)-1,a=a..3)
这条语句太长,分成几个表达式写,是括号错误.我调了一下,有两个地方:1.mu=0.4l;%滑动摩擦系数,即μ应为0.41;你是0.4L.2.倒数第二句:f=dW*ds,这是symbol对象,改为inl
%bydynamic%2009.2.%积分Matlab中提供了数值和解析解法两种1.数值方法>>fun=@(x)738.5*x.^2.*exp(-0.00105*x.^2);>>quadl(fun,0
d(2x+3)=2dx所以dx=1/2d(2x+3)所以原式=1/2*∫(2x+3)^100d(2x+3)=1/2*(2x+3)^101/101(0,1)=1/2*(2+3)^101/101-1/2*
我的建议是把e,p的值代入int()中,可以计算,这样是有值的
n=0.01;x=0:n:1;y=4./(1+x.^2);z=trapz(x,y)n取得越小,积分结果越接近精确值pi.再问:老师的意思是把n的取值看做一个函数直接出来不同n的不同结果不用自己一个一个
积分是微分的逆运算,只要把需要微分的参量按原样返还使积分后的式子求导后得原式就可以了!关于积分的计算相信数学书上是有公式的.再利用复合函数积分法很容易就可以得到答案了!打不出式子,我截个图给你~
想给你做,可是宿舍电脑上不了网了、有空可以商讨以下啊
f=inline('sin(x)')>>g=@(x)x.*f(x)>>quad(g,1,2)