matlab中求多项式的积分

把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)

通过矩形法估计面积：>> sum(yy(xx>400 & xx<600))*.1ans =  9.8796e

symsabcdx;p=a*x^2+c*x+c+d*x;t=coeffs(p,x);t(2)%%%%输出结果=c+d即为所得.coeffs(p,x)的结果是按照变量的幂来排列的.如上t(1)为常系数c

symsx;m=1;n=463;p1=[];while(m

使用quad有问题>>quad('0.0732*exp(-(x+0.02287).^2/64.90)',0,inf)Warning:InfiniteorNot-a-Numberfunctionvalu

for(i=1:n)x=rand;y=rand*(exp(1)-1)+1;if(exp(x^2)>=y)N=N+1;endends=(N/n)*(exp(1)-1)+1s=1.4627

最好给出具体的函数表达式.一般的做法差不多可以写成：fminunc(@(yz)quadl(@(x)f(x,yz(1),yz(2)),a,b),yz0)其中,yz0是y和z的初值,a和b是积分上下限.&

fori=1:3q=quadv(@(x)f(x,3),0,t(i))这样就能得出结果了,积分上下线应该都是标量值；你可以再仔细看看quadv的帮助文档

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)

clearall;closeall;clcsigma1=1;mu1=0;mu2=0;sigma2=1;p=0.5;y=-4:0.1:4;symsxf=1/(2*pi*sigma1*sigma2*sqr

symskcv>>k=2;>>c=8.5;>>int(k/c*(v/c)^(k-1)*exp(-(v/c)^k),'v',8,22)ans=1/exp(256/289)-1/exp(1936/289)

这条语句太长,分成几个表达式写,是括号错误.我调了一下,有两个地方：1.mu=0.4l;%滑动摩擦系数,即μ应为0.41；你是0.4L.2.倒数第二句：f=dW*ds,这是symbol对象,改为inl

程序改成如下：A=[13631];s=poly2str(A,'x');x=unifrnd(0,5,150,1)%产生150个[05]上的均匀分布随机数y=polyval(A,x)%5阶多项式拟合p=p

可以用复合的形式,比如说:symsaxy;f=x^2+y^2;I=int(int(f,x,(y-a),y),y,a,3*a)%先对x积分,积分限为y-a到y,再对y积分,积分限为a到3a.执行结果：I

套公式,代数!再问：呵呵呵，无语再答：那你要给个具体数据吧！

设被积函数为f=ax+b,积分区间为（c,d)MATLAB代码如下：symsabcdx;f=a*x+b;answer=int(f,x,c,d).

symsxa>>b=int(a*exp(-a*x),x,0,10)b=1-1/exp(10*a)>>aa=solve(b-1)%题有问题,若积分值为1,则a无穷大,所以无解Warning:Explic

polyder的对应积分函数为polyint

按springehome的做法,是得不出结果的,这是MATLAB给的结果:int=int(-((73*(C/T)^(1/2))/100+(503*(C/T)^(1/2))/(4*T)-53/2000)

poly2sym(poly([x1,x2,x3,...]))x1,x2,x3...是根.