matlab求解定积分上下限的问题？
初学matlab，遇到了一点求定积分的问题。 等式左边的那个τ是个常数，等于0.65，要求定积分的上限θ。用符号积分方法，也就是int方法，可能是积分的式子过于复杂，表达不出具体的式子，算出来的结果还是带有y的。syms x t yy=int(t/(exp(t)-1),0,x)solve('(1+(4/x)*((1/x)*y-1))','x')Warning: Explicit integral could not befound.
y = int(t/(exp(t) - 1), t == 0..x)ans =
(2*y)/((1 - y)^(1/2) + 1) -(2*y)/((1 - y)^(1/2) - 1)所以考虑用数值分析法，看了很多关于这方面的解答，但是关于“@”的用法不是很懂。还有quadl和quad在用法上有什么不同？话说这两个函数不是只能用常数作为上下限吗？如何在不确定上下限的时候用quad或者quald解出这个方程呢？如果可以，最好能给出解法，谢谢知乎的各位大神。
按投票排序
matlab说quad函数改成了integral。定义一个theta是自变量的函数func1，返回tau-0.65。"@"就是一个函数的句柄，跟c里边的地址差不多。fsolve可以解数值解建立一个func1.m文件到当前目录function y=func1(theta)func2=@(t) t./(exp(t)-1);y=1+4/theta*(1/theta*integral(func2,0,theta)-1)-0.65;然后输入fsolve(@func1,10)，就可以解出9.4376
我的妈呀。知乎上连问这个问题的都有。不知道楼主是不是在算Frank Copula函数时Tau与Theta的关系时遇到的。想问最后这个问题怎么解决的？小白我对二楼的回答看的不是很懂~~~扫扫二维码，随身浏览文档
手机或平板扫扫即可继续访问
Matlab定积分的近似计算
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由：
将文档分享至：
分享完整地址
文档地址：
粘贴到BBS或博客
flash地址：
支持嵌入FLASH地址的网站使用
html代码：
&embed src='/DocinViewer-4.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布，请您等待！
3秒自动关闭窗口Matlab课程设计_交互界面求定积分_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
评价文档：
Matlab课程设计_交互界面求定积分
m​a​t​l​a​b​课​程​设​计​,​交​互​界​面​求​定​积​分​,​电​气​学​院
阅读已结束，如果下载本文需要使用
想免费下载本文？
你可能喜欢如何用matlab计算定积分_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
评价文档：
如何用matlab计算定积分
使​用​m​a​t​l​a​b​计​算​定​积​分
阅读已结束，如果下载本文需要使用
想免费下载本文？
你可能喜欢59如何用matlab计算定积分
上亿文档资料，等你来发现
59如何用matlab计算定积分
用matlab计算积分；4．1积分的有关理论；定积分：积分是微分的无限和，函数f(x)在区间[；I=∫f(x)dx=；max(?xi)→0；lim；∑f(ξ)?x；i=1；其中；a=x0&x1&?&xn=b,；于[a,b]上非负函数f(x)，记分值I是曲线y；微积分基本定理（Newton-Leibniz公式；F'(x)=f(x),x∈[a,
用matlab计算积分4．1积分的有关理论定积分：积分是微分的无限和，函数f(x)在区间[a,b]上的积分定义为I=∫f(x)dx=abmax(?xi)→0lim∑f(ξ)?xii=1ni 其中a=x0&x1&?&xn=b,?xi=xi?xi?1,ξi∈(xi?1,xi),i=1,2,?,n.从几何意义上说，对于[a,b]上非负函数f(x)，记分值I是曲线y=f(x)与直线x=a,x=b及x轴所围的曲边梯形的面积。有界连续（或几何处处连续）函数的积分总是存在的。微积分基本定理（Newton-Leibniz公式）：f(x)在[a,b]上连续，且F'(x)=f(x),x∈[a,b]，则有∫baf(x)dx=F(b)?F(a) 这个公式表明导数与积分是一对互逆运算，它也提供了求积分的解析方法：为了求f(x)的定积分，需要找到一个函数F(x)，使F(x)的导数正好是f(x)，我们称F(x)是f(x)的原函数或不定积分。不定积分的求法有学多数学技巧，常用的有换元积分和分部积分法。从理论上讲，可积函数的原函数总是存在的，但很多被积函数的原函数不能用初等函数表示，也就是说这些积分不能用解析方法求解，需用数值积分法解决。在应用问题中，常常是利用微分进行分析，而问题最终归结为微分的和（即积分）。一些更复杂的问题是含微分的方程，不能直接积分求解。多元函数的积分称为多重积分。二重积分的定义为∫∫f(x,y)dxdy=Gmax(?xi2+?yi2)→0lim∑∑f(ξ,ηiijj)?xi?yj 当f(x,y)非负时，积分值表示曲顶柱体的体积。二重积分的计算主要是转换为两次单积分来解决，无论是解析方法还是数值方法，如何实现这种转换，是解决问题的关键。4．2积分的数值方法梯形法：将[a,b]划分为若干小区间ba=x0&x1&?&xn=b,.则nxiI=∫f(x)dx=∑∫ai=1xi?1f(x)dx 在每一小区间[xi?1,xi]上f(x)近似为一直线，用弦代替，有∫xixi?1f(x)dx≈xi?xi?1(f(xi?1)+f(xi))2从而I≈∑i=1nxi?xi?1(f(xi?1)+f(xi))2h=b?a,xi=a+ihn，称为梯形公式。通常将区间[a,b]n等分，f(b)+f(a)n?1+∑f(xi?1))I≈Tn=h(2i=1可以证明，当n→∞时由上述公式给出的梯形法是收敛的。重积分：重积分的数值计算可通过若干次单积分的组合实现，如对于二重积分I=∫∫f(x,y)dxdyG 先化为二次计分I=∫dx∫abd(x)c(x)f(x,y)dy 利用梯形法，先将[a,b]区间m等分，公式可得hx=b?a,xi=a+ihx,i=0,1,?,m.m利用梯形积分m?1d(xi)1I≈hx((G(a)+G(b))+∑G(xi)),G(xi)=∫f(xi,y)dy.c(xi)2i=1再将[c(xi),d(xi)]区间n等分，hy(i)=d(xi)?c(xi),yij=a+jhy(i),j=0,1,?,n.n利用梯形积分公式可得n?11G(xi)≈hy(i)((f(xi,c(xi))+f(xi,d(xi)))+∑f(xi,yij)).2j=1 4．3积分的MATLAB命令MATLAB中主要用int进行符号积分，用trapz,dblquad,quad,quad8等进行数值积分。
R=int(s,v)
%对符号表达式s中指定的符号变量v计算不定积分.表达式R只是表达式函数s的一个原函数,后面没有带任意常数C.R=int(s)
%对符号表达式s中确定的符号变量计算计算不定积分. R=int(s,a,b)
%符号表达式s的定积分，a,b分别为积分的上、下限R=int(s,x,a,b) %符号表达式s关于变量x的定积分，a,b分别为积分的上、下限 trapz(x,y) 梯形积分法，x时表示积分区间的离散化向量，y是与x同维数的向量，表示被积函数，z返回积分值。fblquad(‘fun’,a,b,c,d)
矩形区域二重数值积分，fun表示被积函数的M函数名，a,b分别为x的上、下限，c,d分别为y的上、下限.可以用help int, help trapz, help quad等查阅有关这些命令的详细信息 例1
用符号积分命令int计算积分MATLAB代码为：&& &&int(x^2*sin(x))2x∫sinxdx.结果为ans =-x^2*cos(x)+2*cos(x)+2*x*sin(x)如果用微分命令diff验证积分正确性，MATLAB代码为： &&&&diff(-x^2*cos(x)+2*cos(x)+2*x*sin(x))结果为ans =x^2*sin(x)例2
计算数值积分∫2?2x4dx.2先用梯形积分法命令trapz计算积分∫?2x4dx，MATLAB代码为：&& x=-2:0.1:2; y=x.^4;
%积分步长为0.1 &&trapz(x,y)结果为ans = 12.853364=12.8xdx∫5?2实际上，积分的精确值为。如果取积分步长为0.01, MATLAB代码为：24&& x=-2:0.01:2; y=x.^4;
%积分步长为0.01&&trapz(x,y)结果为ans =12.8005可用不同的步长进行计算，考虑步长和精度之间的关系。一般说来,trapz是最基本的数值积分方法,精度低,适用于数值函数和光滑性不好的函数.如果用符号积分法命令int计算积分输入MATLAB代码为：&& &&int(x^4,x,-2,2) 结果为ans =64/52∫2?2x4dx，例3
计算数值积分x+y2≤1∫∫(1+x+y)dxdy,可将此二重积分转化为累次积分x2+y2≤1∫∫(1+x+y)dxdy=∫∫11?x2?1??x2(1+x+y)dy 输入MATLAB代码为： &&&&iy=int(1+x+y,y,-sqrt(1-x^2),sqrt(1-x^2)); &&int(iy,x,-1,1) 结果为ans =pix2I=∫exp(sinx?dx?∞50。
例4（广义积分） 计算广义积分+∞输入MATLAB代码为：&&&&y=int(exp(sin(x)-x^2/50),-inf,inf); &&vpa(y,10) 结果为15.。包含各类专业文献、专业论文、生活休闲娱乐、中学教育、应用写作文书、高等教育、59如何用matlab计算定积分等内容。
　用递推公式计算定积分(matlab版)_理学_高等教育_教育专区。完整的实验报告哦～～用递推公式计算定积分实验目的: 1.充分理解不稳定的计算方法会造成误差的积累,在... 　求定积分问题就分解为求和问题。 2.数值积分的实现方法 数值积分的实现方法 基于变步长辛普生法, 函数来求定积分。 基于变步长辛普生法,MATLAB 给出了 quad ... 　定积分的近似计算 MATLAB R2012b 利用牛顿―莱布尼兹公式虽然可以精确地计算定积分的值, 但它仅适用于被 积函数的原函数能用初等函数表达出来的情形.如果这点办不... 　matlab中求定积分的几种方法_互联网_IT/计算机_专业资料。求函数'exp(-x*x)的定积分,积分下限为 0,积分上限为 1。 程序: M-文件:建立内联函数 function y... 　实验目的: 1 加深理解积分理论中分割、近似、求和、取极限的思想方法; 2 了解定积分近似计算的矩形法、梯形法与抛物线法; 3 会用 MATLAB 语言编写求定积分近似... 　MATLAB实验三 定积分的近似计算_理学_高等教育_教育专区。MATLAB实验三 定积分的近似计算实验三 实验三 定积分的近似计算一、问题背景与实验目的利用牛顿―莱布尼兹公... 　积分的计算及其简单应用摘要: 本文简要的概述了 MATLAB 在高等数学中积分的计算及应用; 利用 MATLAB 中符号积分和数值积分的命令,计算定积分和不定积分。同时,也... 　matlab实验报告--定积分的近似计算_理学_高等教育_教育专区。定积分的近似计算-...如何用matlab计算定积分 4页 1下载券 matlab实验报告之定积分... 13页 4下载...