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

用matlab绘制这类复杂的积分函数图形,应采用循环语句求出,当T为某值时的C值.实现方法:x=20:1:70; for i=1:length(x) T=x(i); C(i)=fun_ex(T); end T=x; plot(x,C,'*-'),grid on; title('T--C(T)曲线图'); xlabel('T');ylabel('C(T)'); 运行结果

举个例子如下:clc,clear g=inline('cos(x^2)'); i=0; for x=0:0.05:5; i=i+1; f=@(y)exp(-y.^4+x).*sin(x.^2); Q(i)=quad(f,0,g(x)); end plot(0:0.05:5,Q)

下面是求sin(t)/t的一个函数:变上限积分,另外提醒你:这个函数没有对应的初等函数形式的原函数,也就是说不定积分没有解析解!! clear,clc t=linspace(0,10,100); f=@(t) sin(t)./(t+eps); F=t; n=length(t); for i=1:n F(i)=quad(f,0,t(i)); end plot(t,F)

提示你从符号变量变到双精度数值出现错误.还有你注意有些地方需要用到点乘或点除.给你一个可以运行的程序:clear all; syms x; M=2; fx=exp(-x*x/2)/sqrt(2*pi); sn=-20:1:100; t=length(sn); Q0=zeros(length(sn)); Q=Q0(1:length(sn)); for k=1:t Q(k)=eval(int(fx,x,sn(k),inf)); end P=2*(M-1)*Q/M; plot(sn,P);

matlab说他也不会做.如果你只是想求解这个积分的值,可按如下方法做1. 首先看看被积函数是什么形态,有没有奇点 ezplot('(x^2*cos(x)^2 + 1)^(1/2)'); 你这个函数是连续的,所以可积,然后2. 用简单的数值积分 step = 1e-2; x = 0:step:1; y = (x.^2.*cos(x).^2 + 1).^(1/2); step*sum(y) % 积分值

plot作图基本的两种输入方式 (1)plot(y),这个y是一个矩阵,n行m列的话,一共会画出m条,每条n个点 (2)plot(x,y),x与y是同维矩阵,比如 同是1行n列,或者同是n行1列,都能作图,只要能找到x与y一一对应的点就行,如果x与y的元素个数不能一一对应(比如不同维),那就做不出来图的 你现在workplace里看一看fp与y的储存形式,起码应该是矩阵.要是函数绘图 ,有自己单独的绘图指令:fplot,基本用法是 fplot('fun',[xmin,xmax])

函数Q不是数组,需要赋值不采用sym定义,Conversion to double from “sym” is not possible.T=300:1:600;U=293./(273+T)*2100/15; Q=-(x-U).^2./(U.^2)/0.0032;p=int(1/(sqrt(2*pi)*8.4)*exp(Q),x,-inf,sqrt(2/3)*110);plot(T,Q);

(1)首先绘制积分函数曲线 (2)利用fill命令填充积分函数曲线在积分区间内和x轴包围的区域.如下命令绘制了积分函数exp(-(x.^2)) 在[0 1]区间内的定积分,即蓝色区域的面.

我在你前面的提问里回答过了,但答案无法正常显示,说是“提交中”,唉,恶心的度娘. 分析 1、你要求的那个方程无法求出解析解; 2、求解微分方程需要初始.