xcxd.net
当前位置:首页 >> 积分上限的mAtlAB求解 >>

积分上限的mAtlAB求解

clear allclc;syms t;f = sin(t)/t;x =0:10/1000:10;for i =1:length(x); y(i) =int(sin(t)/t,t,0,x(i));endplot(x,y);y1 = int(sin(t)/t,t,0,4.5)

x=0.6时,y=68.581359956511627574834151414947 close allclear,clcx=0:0.05:1;y=zeros(1,length(x));syms y Yex1=int(int('exp(-(t-50)^2/200)',-inf,y),0,Y);ex2=int(int('exp(-(t-50)^2/200)',-inf,0.6*y),0,0.6*Y);ex1=ex1/(75*sqrt(2)*pi^(1/2));ex2=ex2/(75*

直接使用CFTOOL工具箱命令行输入cftool即可,然后选择拟合类型x=[6.69:0.02:7.53];y=[4.2,3.7,3.3,2.95,2.63,2.33,2.11,1.87,1.65,1.47,1.32,1.17,1.04,0.925,0.82,0.735,0.653,0.582,0.52,0.462,0.412,0.366,0.325,0.289,0.258,0.23,0.205,0.182,0.

fun=@(x)x.^2;% a1=[1 3 4 5]; for i=1:length(a1) result(i)=quadl(fun,0,a1(i)); end

syms x Q solve(int(1/100*(x/100)^4*exp(-x/100)/24,x,0,Q)-4/5,Q) ans =672.09787874865558455510741713836

int(f,x,a,b),f是被积函数,x是变量,如果你函数里只有x这一个变量,其他都是数字的话那可以省略这一项,a、b是积分限,a是下限,b是上限,无穷大是inf

>> syms t r m x>> y=1/(t+2*(r-sqrt(r^2-(m/2-x)^2)));>> f=int(y,0,m);f=simple(f);f1/2*(-i*t*atanh((4*r^2+(-t*(t+4*r))^(1/2)*m)/((t+2*r)^2)^(1/2)/(-m^2+4*r^2)^(1/2))+2*t*atan(m/(t^2+4*t*r)^(1/2))-i*t*atanh((-4*r^2+(-t*(t+4*r))^(1/2)*m)/((t+2*r)^2)^(1/2)/(-m^2+4*r^2

数值试算法:编一个程序即可,不小心我试对了,上限是0.3. m=0.3 f=quadl('1./x.^ 0.6131./(1-x).^ 1.0633',0.0001,m)% -1.7266结果:f = 1.7266

可以求的.使用命令前,先执行这几个命令,clear all;close all;clc

按你的方程积分不出来吧 syms x z int(pi^2/(sqrt(1-x^2)*sqrt(1-(x-z)^2)),x,-1,z)Warning: Explicit integral could not be found.

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