jcst.net
当前位置:首页 >> mAtlAB梯形法求F(x)在[A,B]上的定积分 >>

mAtlAB梯形法求F(x)在[A,B]上的定积分

给你一个例子。 x=-1:0.05:1的积分区间划成步长为0.05的41个等分 y=exp(-x.^2) %被积函数f(x) trapz(x,y) %梯形积分法函数 运行结果: ans=1.4933

1.数值积分(梯形法) >> x=0:0.01:1; >> y=1./(1+x.^4); >> sy=trapz(x,y) sy = 0.8670 2.符号积分,对比 >> syms x >> y=1/(1+x^4); >> yi=eval(int(y,0,1)) yi = 0.8670 >>

我能帮你写,今晚就可以完成,请私信我看具体要求。 function y = jifen(f, a, b, N) % 矩形法算积分 f = sym(f); f = matlabFunction(f); D = b - a; dx = D./N; yn = zeros(1,N); for i = 1:N yn(i) = f(a+dx*i); end y = sum(yn .* dx); 将以...

for N=10:10:100 x=0:2*pi/(N-1):2*pi; y=sin(x); z(N/10)=0.5*sum((y(1:end-1)+y(2:end)).*2*pi/(N-1)); end N=10:10:100; plot(N,z,'-o',N,N*0,'*')

可以利用matlab的trapz函数命令 x=0:0.00001:1;%x用来储存积分点 y=(x+1).*sin(x);%y用来求解积分点x处的函数值 I=trapz(x,y) I = 0.7608663730793 验证该问题的解析解 syms x y=(x+1)*sin(x);%被积函数表达式 II=int(y,0,1) II = sin(1) - 2*co...

1、使用int函数,函数由integrate缩写而来,int 函数表达式,变量,积分上限,积分下限。 2、比如求一个Fx = a*x^2,在区间(m,n)对x进行积分, 首先要将 m,x,a,b 这四个变量定义为符号变量 syms m x a b; Fx = a*x^2; int(Fx,x,m,n) 3、通过上...

k=[1 2 3 4 5 6 7]; xk=[0.3 0.5 0.7 0.9 1.1 1.3 1.5]; fxk=[0.3895 0.6598 0.9147 1.1611 1.3971 1.6212 1.8325]; s1=sum(fxk(1:end-1).*(diff(xk))); %矩形 s2=trapz(xk,fxk); %梯形 矩形公式就是把将微分的面积元看成是矩形 ΔS(i)=fxk(i)*Δx...

>> clear >> syms a b c d e; >> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e') y = (2*a*e + b*d + d*(b^2 - 4*a*c)^(1/2))/(2*exp((t*(b - (b^2 - 4*a*c)^(1/2)))/(2*a))*(b^2 - 4*a*c)^(1/2)) - (2*a*e + b*d - d*(b^2 - 4*a*c)^(1/2))/(...

比如你要处理的信号叫x(n),是一个N点的序列。理想低通滤波器最简单,先对x(n)做FFT,得到频域特性X(k),然后把高频部分的数据改为0,最后在用IFFT变换转换到时域,就是理想滤波了,因为高频完全被抹掉了,而低频信息丝毫不变。其他滤波器就不能...

%Compound trapezoid formula复化梯形法 function y = traint(a,b,n,f) h = (b - a) / n; x = linspace(a,b,n+1); y1 = h * feval(f,x); y1(1) = y1(1) / 2; y1(n+1) = y1(n+1) / 2; y = sum(y1); %compound Simpson formula复化辛普森公式 func...

网站首页 | 网站地图
All rights reserved Powered by www.jcst.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com