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

dx=0.1;x=0:dx:4;y=x.*sin(x);s=cumtrapz(y)*dx; %梯形法求累计积分 plot(x,y,x,s) 计算xsin(x)在[0,4]的定积分

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,'*')

我能帮你写,今晚就可以完成,请私信我看具体要求。 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); 将以...

可以利用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、通过上...

程序24行 漏了个字母n fprintf('随机生成分割梯形的面积和')

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))/(...

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

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