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]的定积分

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

于无法求得exp(x^2)的原函数,我们只能用数值算法来求解,可以用复化梯形公式、Romberg公式、Gauss公式等,有好多种。我用Matlab编了一个用Gauss公式求解积分的函数。 function S=GaussIntegrate() %运用Gauss求积公式计算数值积分 %f为被积函...

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

咱班群里有快去看!

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

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