matlab积分画图问题,急!!

这个问题的困难不在于积分限里有变量，而在于积分限里有无穷大。如果用数值积分，quad系列函数都不能处理积分区域为无界的情况，但2012a引入了integral系列新的数值积分函数可以处理。不知道你的函数F(x,y)表达式，这里举个例子来说明。假设F(x,y)是二元平方指数函数F(x,y)=exp(-(x^2+y^2))（...

书上的思路应该是先对y积分，然后再对x积分。不过他实际操作int命令错了，对比下面的程序就能看出，而且图中标注的积分范围的选取也确实是错的。先对y积分，再对x积分，应是下面的结果。代码如下：syms x y;f = x;B = fzero('2-x-0.8*x',8)F1 = int(f,y,0.8*x,3*x) %要指明对...

plot（）使用错误 你应该按下列步骤来绘图。x=-6:0.1:6;y=(-2/5)*x+4;z=1000*9.8*pi*x.*y.^2;plot(x,z)xlabel('x'),ylabel('y');grid on

clear;clc;f=@(x)(1./x-1).^(1/3);F=@(x)quadl(f,x,1);fplot(@(x)F(x)-x,[eps 1])

用Matlab求积分，可以直接使用int（）函数求解。（1）不定积分 syms x int(x*sin(x)) %计算结果，sin(x) - x*cos(x)int(x*exp(-x)) %计算结果，-exp(-x)*(x + 1)（2）定积分 syms x int(x*exp(-x),0,1) %计算结果，1 - 2*exp(-1)int(x*cos(x),0,pi/2) ...

那就存在一个问题：y是被积分的变量，积分之后得到的结果与y无关，只能是一个常数而已，怎么画图？参考代码（各常数的值随便取的）：M = 1;G = 2;H = 3;a = 4;b = 5;F = 1/M * ( G*quadl(@(y)1./(H/(M*a)*y.^b+1),0,1) - 1)得到的结果：F = 0.8184 ...

1、确定变上限的范围，如x0=[1,2,3,4,5,6,7,8,9,10,50,100,500,1000,5000,10000,50000,100000,500000,1000000];2、用for循环语句，求出每个x值对应的y值 3、再用plot函数，绘出其积分图像 4、求解代码 x0=[1,2,3,4,5,6,7,8,9,10,50,100,500,1000,5000,10000,50000,100000,...

问题主要出现在f0这个函数。在s比较小而x比较大的时候，指数函数的值在数值意义上为0（小于realmin=2.2251e-308），而Bessel函数的值为无穷大（大于realmax=1.7977e+308），导致出现0*Inf的情况，结果为非数NaN。被积函数的值一旦出现NaN，数值积分就会失败，返回的结果也是NaN，最后得到f是一串NaN，...

这里再对两个问题做进一步探讨：一是对出现NaN的原因做更深入分析，二是把积分下限换成0.1的误差分析。1、结果中出现NaN的原因 之前分析过（参见问题2051722037141864067），之所以不能画图，归根到底是由于f0在某些条件下计算结果出现NaN引起的，而NaN又是由于指数项为0、Bessel函数为无穷大导致的。■ ...

plot作图基本的两种输入方式 （1）plot(y),这个y是一个矩阵，n行m列的话，一共会画出m条，每条n个点 （2）plot(x,y),x与y是同维矩阵，比如 同是1行n列，或者同是n行1列，都能作图，只要能找到x与y一一对应的点就行，如果x与y的元素个数不能一一对应（比如不同维），那就做不出来图...