matlab求变步长积分的例子

同样步长的放在一个向量里,然后把全部向量合在一起例如：a=1:2:5;b=6:1:9;c=[a,b];

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

大致的意思是：变量zeors的数据类型有问题,你定义zeors用的什么类型?再问：zeros不是系统函数么？上面是我左右的思想·我是初级菜鸟·再答：你的t变量如何定义的？还是没有定义？function

积分区间[a,b]等分为n段,积分步长h＝(b-a)/n

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

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

cleara=0:0.5:2;n=5;fori=1:n;a(i)*quadl(@(x)exp(-x.^2),0,1)end再问：若是表达式变成exp(-a(i)*x*x),求它的积分，同样是a=0:0

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

我有excel的程序,是用复化simpson求积的工具.要启用宏.前几天刚考完计算方法,复习完了就做了个小程序.你要的话回我一个邮件,我发给你.

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

http://numericalmethods.eng.usf.edu/mws/gen/07int/mws_gen_int_txt_gaussquadrature.pdf

%假设X1的时间步长比X2大%下面的x1x2y1y2模拟的x1=0:1.2:10;y1=sin(x1);x2=0:0.1:10;y2=sin(x2)+0.1;y1_linear=interp1(x1,

%bydynamic%2009.2.%积分Matlab中提供了数值和解析解法两种1.数值方法>>fun=@(x)738.5*x.^2.*exp(-0.00105*x.^2);>>quadl(fun,0

其实你可以考虑尝试用数组的算法,毕竟MATLAB并不擅长循环运算.如果不是特别复杂的运算,数组的算法可能并不难想出来.如果你非要用变步长的循环,我想到一个折中的办法,因为你提前知道了m应该取那些数值,

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

楼主是不是编译的时候显示sin没定义?编译的时候给最后面加个-lm就行了

set(gca,'XLim',[01.5]);%X轴的数据显示范围set(gca,'XTick',[0:0.1:1.5]);%设置要显示坐标刻度set(gca,'XTickLabel',[0:0.1:

symsxtf=abs(sin(cos(t)));%需要积分的公式I=int(f,x,0,2*pi)%int(f,x,a,b)表示定积分,f是方程式,x表示参数,a表示下限,b表示上限sim=simp

有20个v,那应该有20个结果吧?而且你的那个函数的定积分是没有解析解的.所以只能求数值解.建议你利用trapz函数,具体的看help.其实也可以用quad但是它需要函数的定义,所以有点麻烦,我不推荐