matlab中求多项式的积分
来源:学生作业帮助网 编辑:作业帮 时间:2024/09/20 17:33:48
把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...是根.